diff --git a/solver/compressibleInterFlow/compressibleInterFlow.C b/solver/compressibleInterFlow/compressibleInterFlow.C
index 09700486..5252c7f3 100644
--- a/solver/compressibleInterFlow/compressibleInterFlow.C
+++ b/solver/compressibleInterFlow/compressibleInterFlow.C
@@ -133,13 +133,13 @@ int main(int argc, char *argv[])
                 {
                     MRF.update();
 
-                    advector->surf().mapAlphaField();
-                    alpha2 = 1.0 - alpha1;
-                    alpha2.correctBoundaryConditions();
-                    rho == alpha1*rho1 + alpha2*rho2;
-                    rho.correctBoundaryConditions();
-                    rho.oldTime() = rho;
-                    alpha2.oldTime() = alpha2;
+                    // advector->surf().mapAlphaField();
+                    // alpha2 = 1.0 - alpha1;
+                    // alpha2.correctBoundaryConditions();
+                    // rho == alpha1*rho1 + alpha2*rho2;
+                    // rho.correctBoundaryConditions();
+                    // rho.oldTime() = rho;
+                    // alpha2.oldTime() = alpha2;
 
                     Info<< "Execution time for mesh.update() = "
                         << runTime.elapsedCpuTime() - timeBeforeMeshUpdate
diff --git a/solver/interFlow/interFlow.C b/solver/interFlow/interFlow.C
index 33a65d3c..9009d6ad 100644
--- a/solver/interFlow/interFlow.C
+++ b/solver/interFlow/interFlow.C
@@ -86,13 +86,13 @@ int main(int argc, char *argv[])
                     // gets recompute by surfaces forces
                     // gh = (g & mesh.C()) - ghRef;
                     // ghf = (g & mesh.Cf()) - ghRef;
-                    advector->surf().mapAlphaField();
-                    alpha2 = 1.0 - alpha1;
-                    alpha2.correctBoundaryConditions();
-                    rho == alpha1*rho1 + alpha2*rho2;
-                    rho.correctBoundaryConditions();
-                    rho.oldTime() = rho;
-                    alpha2.oldTime() = alpha2;
+                    // advector->surf().mapAlphaField();
+                    // alpha2 = 1.0 - alpha1;
+                    // alpha2.correctBoundaryConditions();
+                    // rho == alpha1*rho1 + alpha2*rho2;
+                    // rho.correctBoundaryConditions();
+                    // rho.oldTime() = rho;
+                    // alpha2.oldTime() = alpha2;
 
                     MRF.update();
 
@@ -159,4 +159,4 @@ int main(int argc, char *argv[])
 }
 
 
-// ************************************************************************* //
\ No newline at end of file
+// ************************************************************************* //
diff --git a/solver/multiRegionPhaseChangeFlow/fluid/solveFluid.H b/solver/multiRegionPhaseChangeFlow/fluid/solveFluid.H
index bc8644c6..a47f018c 100644
--- a/solver/multiRegionPhaseChangeFlow/fluid/solveFluid.H
+++ b/solver/multiRegionPhaseChangeFlow/fluid/solveFluid.H
@@ -32,13 +32,13 @@ if (finalIter)
         {
             MRF.update();
 
-            advector.surf().mapAlphaField();
-            alpha2 = 1.0 - alpha1;
-            alpha2.correctBoundaryConditions();
-            rho == alpha1*rho1 + alpha2*rho2;
-            rho.correctBoundaryConditions();
-            rho.oldTime() = rho;
-            alpha2.oldTime() = alpha2;
+            // advector.surf().mapAlphaField();
+            // alpha2 = 1.0 - alpha1;
+            // alpha2.correctBoundaryConditions();
+            // rho == alpha1*rho1 + alpha2*rho2;
+            // rho.correctBoundaryConditions();
+            // rho.oldTime() = rho;
+            // alpha2.oldTime() = alpha2;
 
             Info<< "Execution time for mesh.update() = "
                 << runTime.elapsedCpuTime() - timeBeforeMeshUpdate