diff --git a/applications/solvers/compressible/steadyUniversalMRFFoam/pEqn.H b/applications/solvers/compressible/steadyUniversalMRFFoam/pEqn.H
index 8eec3318e..1d4b2dc94 100644
--- a/applications/solvers/compressible/steadyUniversalMRFFoam/pEqn.H
+++ b/applications/solvers/compressible/steadyUniversalMRFFoam/pEqn.H
@@ -4,9 +4,10 @@
     surfaceScalarField psisf = fvc::interpolate(psis);
     surfaceScalarField rhof = fvc::interpolate(rho);
 
-    // Needs to be outside of loop since p is changing, but psi and rho are not
+    // Needs to be outside of loop since p is changing, but psi and rho are not.
     surfaceScalarField rhoReff = rhof - psisf*fvc::interpolate(p);
 
+    // --- PISO loop
     while (pimple.correct())
     {
         U = rUA*UEqn.H();
@@ -18,6 +19,11 @@
 
         surfaceScalarField phid("phid", psisf/rhof*phi);
 
+        // Make fluxes relative within the MRF zone
+        mrfZones.relativeFlux(rhoReff, phi);
+        mrfZones.relativeFlux(psisf, phid);
+        mrfZones.relativeFlux(rhoReff, phid2);
+
         p.storePrevIter();
 
         while (pimple.correctNonOrthogonal())