diff --git a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.H b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.H
index c2ed55c2c082ad7103350dbcea1a091048120dde..f90c86b9181f4b2514e5225b49c48dc2ad81bd8e 100644
--- a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.H
+++ b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcel.H
@@ -341,6 +341,9 @@ public:
             //  routine
             inline label faceInterpolation() const;
 
+            //- Cell owner mass
+            inline scalar massCell(const label cellI) const;
+
             //- Particle mass
             inline scalar mass() const;
 
diff --git a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcelI.H b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcelI.H
index edf4b0771e4cd81aa86c04fe727efe1a58d0631a..bf4545032af3ffa0b9ee2708bc762ab1e31bdf47 100644
--- a/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcelI.H
+++ b/src/lagrangian/intermediate/parcels/Templates/KinematicParcel/KinematicParcelI.H
@@ -290,6 +290,16 @@ inline Foam::label Foam::KinematicParcel<ParcelType>::faceInterpolation() const
 }
 
 
+template <class ParcelType>
+inline Foam::scalar Foam::KinematicParcel<ParcelType>::massCell
+(
+    const label cellI
+) const
+{
+    return rhoc_*this->cloud().pMesh().cellVolumes()[cellI];
+}
+
+
 template <class ParcelType>
 inline Foam::scalar Foam::KinematicParcel<ParcelType>::mass() const
 {