interDyMFoam: Update surface fields following mesh change

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
     (

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)