ENH: rhoPimple[DyM]Foam - updated treatment for dpdt

applications/solvers/compressible/rhoPimpleFoam/createFields.H

@@ -35,8 +35,7 @@ volVectorField U
         IOobject::MUST_READ,
         IOobject::AUTO_WRITE
     ),
-    mesh,
-    true
+    mesh
 );
 
 #include "compressibleCreatePhi.H"
@@ -57,20 +56,7 @@ autoPtr<compressible::turbulenceModel> turbulence
     )
 );
 
-Info<< "Creating field dpdt\n" << endl;
-volScalarField dpdt
-(
-    IOobject
-    (
-        "dpdt",
-        runTime.timeName(),
-        mesh,
-        IOobject::READ_IF_PRESENT,
-        IOobject::AUTO_WRITE
-    ),
-    mesh,
-    dimensionedScalar("dpdt", p.dimensions()/dimTime, 0)
-);
+#include "createDpdt.H"
 
 #include "createK.H"
 

src/finiteVolume/cfdTools/compressible/createDpdt.H

+Info<< "Creating field dpdt\n" << endl;
+volScalarField dpdt
+(
+   IOobject
+    (
+        "dpdt",
+        runTime.timeName(),
+        mesh,
+        IOobject::NO_READ,
+        IOobject::NO_WRITE
+    ),
+    fvc::ddt(p)
+);
+
+if (!thermo.dpdt())
+{
+    dpdt == dimensionedScalar("0", dpdt.dimensions(), 0);
+}