diff --git a/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcel.C b/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcel.C
index 70d04a910e..7692c2011f 100644
--- a/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcel.C
+++ b/src/lagrangian/intermediate/parcels/Templates/ReactingMultiphaseParcel/ReactingMultiphaseParcel.C
@@ -453,7 +453,7 @@ Foam::scalar Foam::ReactingMultiphaseParcel<ParcelType>::calcDevolatilisation
         dMassDV[i] = YVolatile[i]*dMassTot;
     }
 
-    td.cloud().addToMassDevolatilisation(dMassTot);
+    td.cloud().addToMassDevolatilisation(this->nParticle_*dMassTot);
 
     return -dMassTot*td.constProps().LDevol();
 }
@@ -515,7 +515,8 @@ Foam::scalar Foam::ReactingMultiphaseParcel<ParcelType>::calcSurfaceReactions
 
     td.cloud().addToMassSurfaceReaction
     (
-        sum(dMassSRGas) + sum(dMassSRLiquid) + sum(dMassSRSolid)
+        this->nParticle_
+       *(sum(dMassSRGas) + sum(dMassSRLiquid) + sum(dMassSRSolid))
     );
 
     // Add enthalpy of consumed components to the carrier phase enthalpy
diff --git a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
index e7a2a71152..237b0dda39 100644
--- a/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
+++ b/src/lagrangian/intermediate/parcels/Templates/ReactingParcel/ReactingParcel.C
@@ -274,7 +274,7 @@ Foam::scalar Foam::ReactingParcel<ParcelType>::calcPhaseChange
     scalar dMassTot = sum(dMassPC);
 
     // Add to cumulative phase change mass
-    td.cloud().addToMassPhaseChange(dMassTot);
+    td.cloud().addToMassPhaseChange(this->nParticle_*dMassTot);
 
     // Effective latent heat of vaporisation
     scalar LEff = td.cloud().composition().L(idPhase, YComponents, pc_, T);