Added missing flux correction code

src/finiteVolume/cfdTools/general/MRF/MRFZones.C

@@ -50,12 +50,42 @@ Foam::MRFZones::MRFZones(const fvMesh& mesh)
             IOobject::NO_WRITE
         ),
         MRFZone::iNew(mesh)
-    )
+    ),
+    mesh_(mesh)
 {}
 
 
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 
+Foam::tmp<Foam::surfaceScalarField> Foam::MRFZones::fluxCorrection() const
+{
+    tmp<surfaceScalarField> tMRFZonesPhiCorr
+    (
+        new surfaceScalarField
+        (
+            IOobject
+            (
+                "MRFZonesPhiCorr",
+                mesh_.time().timeName(),
+                mesh_,
+                IOobject::NO_READ,
+                IOobject::NO_WRITE
+            ),
+            mesh_,
+            dimensionedScalar("zero", dimVelocity*dimArea, 0)
+        )
+    );
+    surfaceScalarField& MRFZonesPhiCorr = tMRFZonesPhiCorr();
+
+    forAll(*this, i)
+    {
+        operator[](i).relativeFlux(MRFZonesPhiCorr);
+    }
+
+    return tMRFZonesPhiCorr;
+}
+
+
 void Foam::MRFZones::addCoriolis(fvVectorMatrix& UEqn) const
 {
     forAll(*this, i)

src/finiteVolume/cfdTools/general/MRF/MRFZones.H

@@ -53,6 +53,11 @@ class MRFZones
 :
     public IOPtrList<MRFZone>
 {
+    // Private data
+
+        //- Mesh reference
+        const fvMesh& mesh_;
+
 
     // Private Member Functions