diff --git a/applications/solvers/multiphase/interFoam/alphaEqnSubCycle.H b/applications/solvers/multiphase/interFoam/alphaEqnSubCycle.H
index 428876fd2297f7e39ee228d688bc51307dc49c00..ac40a9bc9e1d163440129b15fc4eb53d49e85cf6 100644
--- a/applications/solvers/multiphase/interFoam/alphaEqnSubCycle.H
+++ b/applications/solvers/multiphase/interFoam/alphaEqnSubCycle.H
@@ -3,7 +3,17 @@
 if (nAlphaSubCycles > 1)
 {
     dimensionedScalar totalDeltaT = runTime.deltaT();
-    surfaceScalarField rhoPhiSum(0.0*rhoPhi);
+    surfaceScalarField rhoPhiSum
+    (
+        IOobject
+        (
+            "rhoPhiSum",
+            runTime.timeName(),
+            mesh
+        ),
+        mesh,
+        dimensionedScalar("0", dimDensity, 0)
+    );
 
     for
     (
diff --git a/applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C b/applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C
index 4c986ef2d8ab48263f7f520814aabded1a11df2d..032f2db7303cc5004bc0e01355c78bceb346dfa6 100644
--- a/applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C
+++ b/applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C
@@ -100,6 +100,8 @@ int main(int argc, char *argv[])
 
             // Make the flux relative to the mesh motion
             fvc::makeRelative(phi, U);
+
+            interface.correct();
         }
 
         if (mesh.changing() && checkMeshCourantNo)