diff --git a/src/finiteVolume/fvMatrices/solvers/isoAdvection/isoAdvection/isoAdvection.C b/src/finiteVolume/fvMatrices/solvers/isoAdvection/isoAdvection/isoAdvection.C
index 1bb244e84aa94686654054c7f3bc657cae58f56d..003226f682810274d21085571deda97f7f7f40c7 100644
--- a/src/finiteVolume/fvMatrices/solvers/isoAdvection/isoAdvection/isoAdvection.C
+++ b/src/finiteVolume/fvMatrices/solvers/isoAdvection/isoAdvection/isoAdvection.C
@@ -479,7 +479,7 @@ void Foam::isoAdvection::limitFluxes()
     // Get time step size
     const scalar dt = mesh_.time().deltaT().value();
 
-    volScalarField alphaNew = alpha1_ - fvc::surfaceIntegrate(dVf_);
+    volScalarField alphaNew(alpha1_ - fvc::surfaceIntegrate(dVf_));
     const scalar aTol = 1.0e-12;          // Note: tolerances
     scalar maxAlphaMinus1 = gMax(alphaNew) - 1;      // max(alphaNew - 1);
     scalar minAlpha = gMin(alphaNew);           // min(alphaNew);