ENH: modified laplacianFoam to use dimScalar or volScalar for DT (#1388)

applications/solvers/basic/laplacianFoam/createFields.H

@@ -14,8 +14,24 @@ volScalarField T
 );
 
 
-Info<< "Reading transportProperties\n" << endl;
+Info<< "Reading diffusivity DT\n" << endl;
 
+volScalarField DT
+(
+    IOobject
+    (
+        "DT",
+        runTime.timeName(),
+        mesh,
+        IOobject::READ_IF_PRESENT,
+        IOobject::AUTO_WRITE
+    ),
+    mesh,
+    dimensionedScalar(dimViscosity, Zero)
+);
+
+if (!DT.headerOk())
+{
     IOdictionary transportProperties
     (
         IOobject
@@ -27,10 +43,7 @@ IOdictionary transportProperties
             IOobject::NO_WRITE
         )
     );
-
+    DT = dimensionedScalar("DT", dimViscosity, transportProperties);
-
+}
-Info<< "Reading diffusivity DT\n" << endl;
-
-dimensionedScalar DT("DT", dimViscosity, transportProperties);
 
 #include "createFvOptions.H"

applications/solvers/basic/laplacianFoam/write.H

@@ -41,6 +41,18 @@
             gradT.component(vector::Z)
         );
 
+        volVectorField DTgradT
+        (
+            IOobject
+            (
+                "flux",
+                runTime.timeName(),
+                mesh,
+                IOobject::NO_READ,
+                IOobject::AUTO_WRITE
+            ),
+            DT*gradT
+        );
 
         runTime.write();
     }