diff --git a/src/engine/engineMesh/fvMotionSolverEngineMesh/fvMotionSolverEngineMesh.C b/src/engine/engineMesh/fvMotionSolverEngineMesh/fvMotionSolverEngineMesh.C
index 4f43e0d822..8c7db2230b 100644
--- a/src/engine/engineMesh/fvMotionSolverEngineMesh/fvMotionSolverEngineMesh.C
+++ b/src/engine/engineMesh/fvMotionSolverEngineMesh/fvMotionSolverEngineMesh.C
@@ -68,9 +68,11 @@ void Foam::fvMotionSolverEngineMesh::move()
     motionSolver_.cellMotionU().boundaryField()[pistonIndex_] == deltaZ;
 
     {
-        scalarField linerPoints =
+        scalarField linerPoints
+        (
             motionSolver_.cellMotionU()
-           .boundaryField()[linerIndex_].patch().Cf().component(vector::Z);
+           .boundaryField()[linerIndex_].patch().Cf().component(vector::Z)
+        );
 
         motionSolver_.cellMotionU().boundaryField()[linerIndex_] ==
             deltaZ*pos(deckHeight_.value() - linerPoints)
diff --git a/src/engine/engineMesh/layeredEngineMesh/layeredEngineMesh.C b/src/engine/engineMesh/layeredEngineMesh/layeredEngineMesh.C
index 41285ea67a..abe366a5c5 100644
--- a/src/engine/engineMesh/layeredEngineMesh/layeredEngineMesh.C
+++ b/src/engine/engineMesh/layeredEngineMesh/layeredEngineMesh.C
@@ -64,7 +64,7 @@ void Foam::layeredEngineMesh::move()
     // Position of the top of the static mesh layers above the piston
     scalar pistonPlusLayers = pistonPosition_.value() + pistonLayers_.value();
 
-    pointField newPoints = points();
+    pointField newPoints(points());
 
     forAll(newPoints, pointi)
     {
diff --git a/src/engine/include/StCorr.H b/src/engine/include/StCorr.H
index 80df7f0af4..c3bdc5771f 100644
--- a/src/engine/include/StCorr.H
+++ b/src/engine/include/StCorr.H
@@ -80,8 +80,10 @@
 
             // Calculate kernel area from b field consistent with the
             // discretisation of the b equation.
-            volScalarField mgb =
-                fvc::div(nf, b, "div(phiSt,b)") - b*fvc::div(nf) + dMgb;
+            const volScalarField mgb
+            (
+                fvc::div(nf, b, "div(phiSt,b)") - b*fvc::div(nf) + dMgb
+            );
             dimensionedScalar AkEst = gSum(mgb*mesh.V().field());
 
             StCorr.value() = max(min((Ak/AkEst).value(), 10.0), 1.0);
diff --git a/src/errorEstimation/errorDrivenRefinement/errorDrivenRefinement.C b/src/errorEstimation/errorDrivenRefinement/errorDrivenRefinement.C
index e7c08e25af..feb0219954 100644
--- a/src/errorEstimation/errorDrivenRefinement/errorDrivenRefinement.C
+++ b/src/errorEstimation/errorDrivenRefinement/errorDrivenRefinement.C
@@ -130,7 +130,7 @@ void Foam::errorDrivenRefinement::setRefinement(polyTopoChange& ref) const
     const volVectorField& resError =
         runTime.lookupObject<volVectorField>(errorField_);
 
-    const volScalarField magResError = Foam::mag(resError);
+    const volScalarField magResError(Foam::mag(resError));
 
     scalar min = Foam::min(magResError).value();
     scalar max = Foam::max(magResError).value();