diff --git a/applications/solvers/basic/laplacianFoam/createFields.H b/applications/solvers/basic/laplacianFoam/createFields.H
index dbddccb904936402d36f7330d82e00cfae71992d..6f7ac045291086079d3d6e0df2769f6cbab5123a 100644
--- a/applications/solvers/basic/laplacianFoam/createFields.H
+++ b/applications/solvers/basic/laplacianFoam/createFields.H
@@ -14,23 +14,36 @@ volScalarField T
 );
 
 
-Info<< "Reading transportProperties\n" << endl;
+Info<< "Reading diffusivity DT\n" << endl;
 
-IOdictionary transportProperties
+volScalarField DT
 (
     IOobject
     (
-        "transportProperties",
-        runTime.constant(),
+        "DT",
+        runTime.timeName(),
         mesh,
-        IOobject::MUST_READ_IF_MODIFIED,
-        IOobject::NO_WRITE
-    )
+        IOobject::READ_IF_PRESENT,
+        IOobject::AUTO_WRITE
+    ),
+    mesh,
+    dimensionedScalar(dimViscosity, Zero)
 );
 
-
-Info<< "Reading diffusivity DT\n" << endl;
-
-dimensionedScalar DT("DT", dimViscosity, transportProperties);
+if (!DT.headerOk())
+{
+    IOdictionary transportProperties
+    (
+        IOobject
+        (
+            "transportProperties",
+            runTime.constant(),
+            mesh,
+            IOobject::MUST_READ_IF_MODIFIED,
+            IOobject::NO_WRITE
+        )
+    );
+    DT = dimensionedScalar("DT", dimViscosity, transportProperties);
+}
 
 #include "createFvOptions.H"
diff --git a/applications/solvers/basic/laplacianFoam/write.H b/applications/solvers/basic/laplacianFoam/write.H
index 6b14ba49955810d81208f587331c05ac991b0f70..c4d405d449eaa9d619010916a824b42ba7bbf370 100644
--- a/applications/solvers/basic/laplacianFoam/write.H
+++ b/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();
     }