Test on second order restart using backward time scheme

applications/solvers/compressible/rhoPimpleFoam/createFields.H

@@ -35,7 +35,8 @@ volVectorField U
         IOobject::MUST_READ,
         IOobject::AUTO_WRITE
     ),
-    mesh
+    mesh,
+    true
 );
 
 #include "compressibleCreatePhi.H"
@@ -63,7 +64,9 @@ volScalarField dpdt
     (
         "dpdt",
         runTime.timeName(),
-        mesh
+        mesh,
+        IOobject::READ_IF_PRESENT,
+        IOobject::AUTO_WRITE
     ),
     mesh,
     dimensionedScalar("dpdt", p.dimensions()/dimTime, 0)

src/OpenFOAM/fields/GeometricFields/GeometricField/GeometricField.C

@@ -770,10 +770,11 @@ void Foam::GeometricField<Type, PatchField, GeoMesh>::storeOldTimes() const
     )
     {
         storeOldTime();
+        timeIndex_ = this->time().timeIndex();
     }
 
     // Correct time index
-    timeIndex_ = this->time().timeIndex();
+    //timeIndex_ = this->time().timeIndex();
 }
 
 

src/finiteVolume/finiteVolume/ddtSchemes/backwardDdtScheme/backwardDdtScheme.H

@@ -92,27 +92,13 @@ public:
         backwardDdtScheme(const fvMesh& mesh)
         :
             ddtScheme<Type>(mesh)
-        {
-            // Ensure the old-old-time cell volumes are available
-            // for moving meshes
-            if (mesh.moving())
-            {
-                mesh.V00();
-            }
-        }
+        {}
 
         //- Construct from mesh and Istream
         backwardDdtScheme(const fvMesh& mesh, Istream& is)
         :
             ddtScheme<Type>(mesh, is)
-        {
-            // Ensure the old-old-time cell volumes are available
-            // for moving meshes
-            if (mesh.moving())
-            {
-                mesh.V00();
-            }
-        }
+        {}
 
 
     // Member Functions

src/finiteVolume/fvMesh/fvMesh.C

@@ -197,7 +197,7 @@ void Foam::fvMesh::storeOldVol(const scalarField& V)
                     "V0",
                     time().timeName(),
                     *this,
-                    IOobject::READ_IF_PRESENT,//NO_READ,
+                    IOobject::NO_READ,
                     IOobject::AUTO_WRITE,
                     true
                 ),
@@ -281,8 +281,8 @@ Foam::fvMesh::fvMesh(const IOobject& io)
                 time().timeName(),
                 *this,
                 IOobject::MUST_READ,
-                IOobject::NO_WRITE,
-                false
+                IOobject::AUTO_WRITE,
+                true
             ),
             *this
         );

src/finiteVolume/fvMesh/fvMeshGeometry.C

@@ -258,13 +258,14 @@ const Foam::volScalarField::Internal& Foam::fvMesh::V00() const
                 "V00",
                 time().timeName(),
                 *this,
-                IOobject::READ_IF_PRESENT,//NO_READ,
-                IOobject::AUTO_WRITE,
+                IOobject::NO_READ,
+                IOobject::NO_WRITE,
                 true
             ),
             V0()
         );
 
+
         // If V00 is used then V0 should be stored for restart
         V0Ptr_->writeOpt() = IOobject::AUTO_WRITE;
     }