diff --git a/applications/solvers/basic/potentialFoam/potentialFoam.C b/applications/solvers/basic/potentialFoam/potentialFoam.C
index 33931f051e..fbe29e652f 100644
--- a/applications/solvers/basic/potentialFoam/potentialFoam.C
+++ b/applications/solvers/basic/potentialFoam/potentialFoam.C
@@ -59,9 +59,10 @@ int main(int argc, char *argv[])
     // function objects so do it ourselves
     runTime.functionObjects().start();
 
+    fvOptions.relativeFlux(phi);
+
     adjustPhi(phi, U, p);
 
-    fvOptions.relativeFlux(phi);
 
     for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
     {
diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/fluid/pEqn.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/fluid/pEqn.H
index 3e566a21cf..161dff5447 100644
--- a/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/fluid/pEqn.H
+++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/chtMultiRegionSimpleFoam/fluid/pEqn.H
@@ -19,12 +19,12 @@
         fvc::interpolate(rho)*(fvc::interpolate(HbyA) & mesh.Sf())
     );
 
+    fvOptions.relativeFlux(fvc::interpolate(rho), phiHbyA);
+
     bool closedVolume = adjustPhi(phiHbyA, U, p_rgh);
 
     phiHbyA += phig;
 
-    fvOptions.relativeFlux(fvc::interpolate(rho), phiHbyA);
-
     dimensionedScalar compressibility = fvc::domainIntegrate(psi);
     bool compressible = (compressibility.value() > SMALL);
 
diff --git a/applications/solvers/incompressible/pimpleFoam/pEqn.H b/applications/solvers/incompressible/pimpleFoam/pEqn.H
index 1d133108d0..f085e548ce 100644
--- a/applications/solvers/incompressible/pimpleFoam/pEqn.H
+++ b/applications/solvers/incompressible/pimpleFoam/pEqn.H
@@ -13,10 +13,10 @@ surfaceScalarField phiHbyA
   + fvc::ddtPhiCorr(rAU, U, phi)
 );
 
-adjustPhi(phiHbyA, U, p);
-
 fvOptions.relativeFlux(phiHbyA);
 
+adjustPhi(phiHbyA, U, p);
+
 // Non-orthogonal pressure corrector loop
 while (pimple.correctNonOrthogonal())
 {
diff --git a/applications/solvers/incompressible/simpleFoam/pEqn.H b/applications/solvers/incompressible/simpleFoam/pEqn.H
index 79f4f56076..9bd829d3f1 100644
--- a/applications/solvers/incompressible/simpleFoam/pEqn.H
+++ b/applications/solvers/incompressible/simpleFoam/pEqn.H
@@ -5,10 +5,11 @@
     UEqn.clear();
 
     surfaceScalarField phiHbyA("phiHbyA", fvc::interpolate(HbyA) & mesh.Sf());
-    adjustPhi(phiHbyA, U, p);
 
     fvOptions.relativeFlux(phiHbyA);
 
+    adjustPhi(phiHbyA, U, p);
+
     // Non-orthogonal pressure corrector loop
     while (simple.correctNonOrthogonal())
     {
diff --git a/applications/solvers/incompressible/simpleFoam/porousSimpleFoam/pEqn.H b/applications/solvers/incompressible/simpleFoam/porousSimpleFoam/pEqn.H
index b05a9fe270..bb906cc8c6 100644
--- a/applications/solvers/incompressible/simpleFoam/porousSimpleFoam/pEqn.H
+++ b/applications/solvers/incompressible/simpleFoam/porousSimpleFoam/pEqn.H
@@ -10,10 +10,11 @@ else
 
 UEqn.clear();
 surfaceScalarField phiHbyA("phiHbyA", fvc::interpolate(HbyA) & mesh.Sf());
-adjustPhi(phiHbyA, U, p);
 
 mrfZones.relativeFlux(phiHbyA);
 
+adjustPhi(phiHbyA, U, p);
+
 while (simple.correctNonOrthogonal())
 {
     tmp<fvScalarMatrix> tpEqn;