From f8673223237f39b2c144d47afe3a12c656c1bb80 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Thu, 20 Jun 2013 16:30:07 +0100
Subject: [PATCH] interDyMFoam: phi.oldTime() is not needed but U.oldTime()
 needs to be created

---
 .../solvers/multiphase/interFoam/interDyMFoam/correctPhi.H     | 1 -
 .../solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C   | 3 +++
 2 files changed, 3 insertions(+), 1 deletion(-)

diff --git a/applications/solvers/multiphase/interFoam/interDyMFoam/correctPhi.H b/applications/solvers/multiphase/interFoam/interDyMFoam/correctPhi.H
index 8d3e60da40..75560b6337 100644
--- a/applications/solvers/multiphase/interFoam/interDyMFoam/correctPhi.H
+++ b/applications/solvers/multiphase/interFoam/interDyMFoam/correctPhi.H
@@ -52,7 +52,6 @@
             phiAbs = phi;
             phiAbs.oldTime() = phi;
             fvc::makeRelative(phi, U);
-            phi.oldTime() = phi;
         }
     }
 }
diff --git a/applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C b/applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C
index 3ccfa73fee..0c22dd01b7 100644
--- a/applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C
+++ b/applications/solvers/multiphase/interFoam/interDyMFoam/interDyMFoam.C
@@ -81,6 +81,9 @@ int main(int argc, char *argv[])
         scalar timeBeforeMeshUpdate = runTime.elapsedCpuTime();
 
         {
+            // Ensure old-time U exists for mapping
+            U.oldTime();
+
             // Calculate the relative velocity used to map the relative flux phi
             volVectorField Urel("Urel", U);