diff --git a/applications/utilities/errorEstimation/estimateScalarError/estimateScalarError.C b/applications/utilities/errorEstimation/estimateScalarError/estimateScalarError.C
index 5c89107d79fcfadf89b6a0756d71cabdacffe412..4daa1c116d370d3a3b21be2acca6ab9cf419f89e 100644
--- a/applications/utilities/errorEstimation/estimateScalarError/estimateScalarError.C
+++ b/applications/utilities/errorEstimation/estimateScalarError/estimateScalarError.C
@@ -112,7 +112,7 @@ int main(int argc, char *argv[])
             );
 
             ee.residual()().write();
-            volScalarField e = ee.error();
+            volScalarField e(ee.error());
             e.write();
             mag(e)().write();
         }
diff --git a/applications/utilities/errorEstimation/icoErrorEstimate/icoErrorEstimate.C b/applications/utilities/errorEstimation/icoErrorEstimate/icoErrorEstimate.C
index edeb80a4b25a2527aed9026f69b150ac5ab2c910..ec6f5ddcd0d2e4bb3f6c0679bdc2365e4617980f 100644
--- a/applications/utilities/errorEstimation/icoErrorEstimate/icoErrorEstimate.C
+++ b/applications/utilities/errorEstimation/icoErrorEstimate/icoErrorEstimate.C
@@ -108,7 +108,7 @@ int main(int argc, char *argv[])
                -fvc::grad(p)
             );
 
-            volVectorField e = ee.error();
+            volVectorField e(ee.error());
             e.write();
             mag(e)().write();
         }
diff --git a/applications/utilities/errorEstimation/icoMomentError/icoMomentError.C b/applications/utilities/errorEstimation/icoMomentError/icoMomentError.C
index b3e9f927ce087022260a5207f90da5bacd1d88fe..fbe6561f94735d663f24c3b4c26597952c1a17d4 100644
--- a/applications/utilities/errorEstimation/icoMomentError/icoMomentError.C
+++ b/applications/utilities/errorEstimation/icoMomentError/icoMomentError.C
@@ -101,8 +101,8 @@ int main(int argc, char *argv[])
 
 #           include "createPhi.H"
 
-            volScalarField ek = 0.5*magSqr(U);
-            volTensorField gradU = fvc::grad(U);
+            volScalarField ek(0.5*magSqr(U));
+            volTensorField gradU(fvc::grad(U));
 
             // Divergence of the error in U squared
 
diff --git a/applications/utilities/errorEstimation/momentScalarError/momentScalarError.C b/applications/utilities/errorEstimation/momentScalarError/momentScalarError.C
index ce38d1dfb4ac2ae077bd7cfe3a2022d5b79a5842..dd6a0cb25a820feb10404072204244bfda33cc04 100644
--- a/applications/utilities/errorEstimation/momentScalarError/momentScalarError.C
+++ b/applications/utilities/errorEstimation/momentScalarError/momentScalarError.C
@@ -106,9 +106,9 @@ int main(int argc, char *argv[])
 
 #           include "createPhi.H"
 
-            volVectorField gradT = fvc::grad(T);
+            volVectorField gradT(fvc::grad(T));
 
-            volScalarField TE = 0.5*sqr(T);
+            volScalarField TE(0.5*sqr(T));
 
             volScalarField L
             (
diff --git a/src/errorEstimation/errorEstimate/resErrorDiv.C b/src/errorEstimation/errorEstimate/resErrorDiv.C
index d374b180ea6ba93afbf342905bffac8a53272d21..b1a0d908818947b6ee115c85bae2b5ef22c8030e 100644
--- a/src/errorEstimation/errorEstimate/resErrorDiv.C
+++ b/src/errorEstimation/errorEstimate/resErrorDiv.C
@@ -57,13 +57,14 @@ div
     scalarField aNorm(vols.size(), 0.0);
 
     // Get sign of flux
-    const surfaceScalarField signF = pos(flux);
+    const surfaceScalarField signF(pos(flux));
 
     // Calculate gradient of the solution
     GeometricField
     <
         typename outerProduct<vector, Type>::type, fvPatchField, volMesh
-    > gradVf = fvc::grad(vf);
+    >
+    gradVf(fvc::grad(vf));
 
     // Internal faces
     forAll(owner, faceI)
diff --git a/src/errorEstimation/errorEstimate/resErrorLaplacian.C b/src/errorEstimation/errorEstimate/resErrorLaplacian.C
index c310a583b304434be3ef9d727ad954e35904d6b0..3d895e62fbf3a74de5a062160f99d7edd1551d9c 100644
--- a/src/errorEstimation/errorEstimate/resErrorLaplacian.C
+++ b/src/errorEstimation/errorEstimate/resErrorLaplacian.C
@@ -122,7 +122,7 @@ laplacian
 
     const scalarField& vols = mesh.V();
     const surfaceVectorField& Sf = mesh.Sf();
-    const surfaceScalarField magSf = mesh.magSf();
+    const surfaceScalarField magSf(mesh.magSf());
     const fvPatchList& patches = mesh.boundary();
     const labelUList& owner = mesh.owner();
     const labelUList& neighbour = mesh.neighbour();
@@ -137,7 +137,8 @@ laplacian
     GeometricField
     <
         typename outerProduct<vector, Type>::type, fvPatchField, volMesh
-    > gradVf = fvc::grad(vf);
+    >
+    gradVf(fvc::grad(vf));
 
     // Internal faces
     forAll(owner, faceI)