applications/solvers/multiphase: for MRF make phiHbyA relative BEFORE adjustPhi

Resolves bug-report http://openfoam.org/mantisbt/view.php?id=1665

applications/solvers/multiphase/driftFluxFoam/pEqn.H

@@ -11,8 +11,8 @@
         (fvc::interpolate(HbyA) & mesh.Sf())
       + fvc::interpolate(rho*rAU)*fvc::ddtCorr(U, phi)
     );
-    adjustPhi(phiHbyA, U, p_rgh);
     fvOptions.makeRelative(phiHbyA);
+    adjustPhi(phiHbyA, U, p_rgh);
 
     surfaceScalarField phig
     (

applications/solvers/multiphase/interFoam/pEqn.H

@@ -11,8 +11,8 @@
         (fvc::interpolate(HbyA) & mesh.Sf())
       + fvc::interpolate(rho*rAU)*fvc::ddtCorr(U, phi)
     );
-    adjustPhi(phiHbyA, U, p_rgh);
     fvOptions.makeRelative(phiHbyA);
+    adjustPhi(phiHbyA, U, p_rgh);
 
     surfaceScalarField phig
     (

applications/solvers/multiphase/potentialFreeSurfaceFoam/pEqn.H

@@ -16,9 +16,8 @@ surfaceScalarField phiHbyA
   + rAUf*fvc::ddtCorr(U, phi)
 );
 
-adjustPhi(phiHbyA, U, p_gh);
-
 fvOptions.makeRelative(phiHbyA);
+adjustPhi(phiHbyA, U, p_gh);
 
 // Update the fixedFluxPressure BCs to ensure flux consistency
 setSnGrad<fixedFluxPressureFvPatchScalarField>