diff --git a/src/finiteVolume/cfdTools/general/MRF/MRFZoneTemplates.C b/src/finiteVolume/cfdTools/general/MRF/MRFZoneTemplates.C
index 55821dbfde6ab56fa5ff4370bb808722bd0b7597..23cae40603c5b710c7fbe7db9c12bb48e3d06974 100644
--- a/src/finiteVolume/cfdTools/general/MRF/MRFZoneTemplates.C
+++ b/src/finiteVolume/cfdTools/general/MRF/MRFZoneTemplates.C
@@ -44,11 +44,15 @@ void Foam::MRFZone::relativeRhoFlux
     const vector& origin = origin_.value();
     const vector& Omega = Omega_.value();
 
+    const vectorField& Cfi = Cf.internalField();
+    const vectorField& Sfi = Sf.internalField();
+    scalarField& phii = phi.internalField();
+
     // Internal faces
     forAll(internalFaces_, i)
     {
         label facei = internalFaces_[i];
-        phi[facei] -= rho[facei]*(Omega ^ (Cf[facei] - origin)) & Sf[facei];
+        phii[facei] -= rho[facei]*(Omega ^ (Cfi[facei] - origin)) & Sfi[facei];
     }
 
     // Included patches
@@ -91,11 +95,15 @@ void Foam::MRFZone::absoluteRhoFlux
     const vector& origin = origin_.value();
     const vector& Omega = Omega_.value();
 
+    const vectorField& Cfi = Cf.internalField();
+    const vectorField& Sfi = Sf.internalField();
+    scalarField& phii = phi.internalField();
+
     // Internal faces
     forAll(internalFaces_, i)
     {
         label facei = internalFaces_[i];
-        phi[facei] += rho[facei]*(Omega ^ (Cf[facei] - origin)) & Sf[facei];
+        phii[facei] += rho[facei]*(Omega ^ (Cfi[facei] - origin)) & Sfi[facei];
     }
 
     // Included patches