diff --git a/src/finiteVolume/finiteVolume/fvc/fvcSurfaceIntegrate.C b/src/finiteVolume/finiteVolume/fvc/fvcSurfaceIntegrate.C
index faaf05c..425074d 100644
--- a/src/finiteVolume/finiteVolume/fvc/fvcSurfaceIntegrate.C
+++ b/src/finiteVolume/finiteVolume/fvc/fvcSurfaceIntegrate.C
@@ -72,7 +72,7 @@ void surfaceIntegrate
         }
     }
 
-    ivf /= mesh.V();
+    ivf /= mesh.Vsc();
 }
 
 
diff --git a/src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C b/src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C
index a7a5282..f0ee13c 100644
--- a/src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C
+++ b/src/finiteVolume/fvMatrices/solvers/MULES/MULESTemplates.C
@@ -60,7 +60,7 @@ void Foam::MULES::explicitSolve
             mesh.Vsc0()().field()*rho.oldTime().field()
            *psi0*rDeltaT/mesh.Vsc()().field()
           + Su.field()
-          - psiIf
+          - psiIf*(mesh.V()/mesh.Vsc()().field())
         )/(rho.field()*rDeltaT - Sp.field());
     }
     else