diff --git a/src/dynamicFvMesh/include/meshCourantNo.H b/src/dynamicFvMesh/include/meshCourantNo.H
index 501611db18d2c635adf4454a5a6383eb9c42c7d6..ad8f679657382f98a3af8d7639439342f0c96140 100644
--- a/src/dynamicFvMesh/include/meshCourantNo.H
+++ b/src/dynamicFvMesh/include/meshCourantNo.H
@@ -34,8 +34,10 @@ scalar meanMeshCoNum = 0.0;
if (mesh.nInternalFaces())
{
- scalarField sumPhi =
- fvc::surfaceSum(mag(mesh.phi()))().internalField();
+ scalarField sumPhi
+ (
+ fvc::surfaceSum(mag(mesh.phi()))().internalField()
+ );
meshCoNum = 0.5*gMax(sumPhi/mesh.V().field())*runTime.deltaTValue();
diff --git a/src/finiteVolume/cfdTools/compressible/compressibleCourantNo.H b/src/finiteVolume/cfdTools/compressible/compressibleCourantNo.H
index c96ddeb22ffa9f33765ef6e71e32d434c5096728..cd5e92ed3dfae1325f6c8bde75917a38a8de66cc 100644
--- a/src/finiteVolume/cfdTools/compressible/compressibleCourantNo.H
+++ b/src/finiteVolume/cfdTools/compressible/compressibleCourantNo.H
@@ -34,9 +34,11 @@ scalar meanCoNum = 0.0;
if (mesh.nInternalFaces())
{
- scalarField sumPhi =
+ scalarField sumPhi
+ (
fvc::surfaceSum(mag(phi))().internalField()
- /rho.internalField();
+ / rho.internalField()
+ );
CoNum = 0.5*gMax(sumPhi/mesh.V().field())*runTime.deltaTValue();
diff --git a/src/finiteVolume/cfdTools/general/SRF/derivedFvPatchFields/SRFVelocityFvPatchVectorField/SRFVelocityFvPatchVectorField.C b/src/finiteVolume/cfdTools/general/SRF/derivedFvPatchFields/SRFVelocityFvPatchVectorField/SRFVelocityFvPatchVectorField.C
index 177a6ce049072f686e059ffd61d8cb127ea387a4..2f00ad181eea512529a8619e125093539c7bf921 100644
--- a/src/finiteVolume/cfdTools/general/SRF/derivedFvPatchFields/SRFVelocityFvPatchVectorField/SRFVelocityFvPatchVectorField.C
+++ b/src/finiteVolume/cfdTools/general/SRF/derivedFvPatchFields/SRFVelocityFvPatchVectorField/SRFVelocityFvPatchVectorField.C
@@ -137,7 +137,7 @@ void Foam::SRFVelocityFvPatchVectorField::updateCoeffs()
db().lookupObject<SRF::SRFModel>("SRFProperties");
// Determine patch velocity due to SRF
- const vectorField SRFVelocity = srf.velocity(patch().Cf());
+ const vectorField SRFVelocity(srf.velocity(patch().Cf()));
operator==(-SRFVelocity + inletValue_);
}
diff --git a/src/finiteVolume/cfdTools/incompressible/CourantNo.H b/src/finiteVolume/cfdTools/incompressible/CourantNo.H
index 8b917e1f89fb2dc15e207c42d1a1c0b6d0320d52..d5198141a327d0fb8aacca67bad4d3934265ccb6 100644
--- a/src/finiteVolume/cfdTools/incompressible/CourantNo.H
+++ b/src/finiteVolume/cfdTools/incompressible/CourantNo.H
@@ -34,8 +34,10 @@ scalar meanCoNum = 0.0;
if (mesh.nInternalFaces())
{
- scalarField sumPhi =
- fvc::surfaceSum(mag(phi))().internalField();
+ scalarField sumPhi
+ (
+ fvc::surfaceSum(mag(phi))().internalField()
+ );
CoNum = 0.5*gMax(sumPhi/mesh.V().field())*runTime.deltaTValue();
diff --git a/src/finiteVolume/cfdTools/incompressible/continuityErrs.H b/src/finiteVolume/cfdTools/incompressible/continuityErrs.H
index 181cd251eb72ef1d945bf4e738a0292fe8e03a56..06287750bd09f8ae7034589c13871d0c02ae38d9 100644
--- a/src/finiteVolume/cfdTools/incompressible/continuityErrs.H
+++ b/src/finiteVolume/cfdTools/incompressible/continuityErrs.H
@@ -30,7 +30,7 @@ Description
\*---------------------------------------------------------------------------*/
{
- volScalarField contErr = fvc::div(phi);
+ volScalarField contErr(fvc::div(phi));
scalar sumLocalContErr = runTime.deltaTValue()*
mag(contErr)().weightedAverage(mesh.V()).value();
diff --git a/src/finiteVolume/fields/fvPatchFields/derived/movingWallVelocity/movingWallVelocityFvPatchVectorField.C b/src/finiteVolume/fields/fvPatchFields/derived/movingWallVelocity/movingWallVelocityFvPatchVectorField.C
index 89e9cf9dc3dde447aeabe5854c280d10f1d92b08..aefa37b8f5f6f2e77985aa1283f9c788943f2aab 100644
--- a/src/finiteVolume/fields/fvPatchFields/derived/movingWallVelocity/movingWallVelocityFvPatchVectorField.C
+++ b/src/finiteVolume/fields/fvPatchFields/derived/movingWallVelocity/movingWallVelocityFvPatchVectorField.C
@@ -120,8 +120,10 @@ void Foam::movingWallVelocityFvPatchVectorField::updateCoeffs()
const vectorField Up((pp.faceCentres() - oldFc)/mesh.time().deltaTValue());
const volVectorField& U = db().lookupObject<volVectorField>(UName_);
- scalarField phip =
- p.patchField<surfaceScalarField, scalar>(fvc::meshPhi(U));
+ scalarField phip
+ (
+ p.patchField<surfaceScalarField, scalar>(fvc::meshPhi(U))
+ );
const vectorField n(p.nf());
const scalarField& magSf = p.magSf();
diff --git a/src/finiteVolume/finiteVolume/ddtSchemes/CoEulerDdtScheme/CoEulerDdtScheme.C b/src/finiteVolume/finiteVolume/ddtSchemes/CoEulerDdtScheme/CoEulerDdtScheme.C
index 691522364c304d2b7583dbf5a5eec228d2967523..2c0e5f1e0b11339f41e86f7b068b3ca361b7baef 100644
--- a/src/finiteVolume/finiteVolume/ddtSchemes/CoEulerDdtScheme/CoEulerDdtScheme.C
+++ b/src/finiteVolume/finiteVolume/ddtSchemes/CoEulerDdtScheme/CoEulerDdtScheme.C
@@ -382,7 +382,7 @@ CoEulerDdtScheme<Type>::fvmDdt
fvMatrix<Type>& fvm = tfvm();
- scalarField rDeltaT = CorDeltaT()().internalField();
+ scalarField rDeltaT(CorDeltaT()().internalField());
fvm.diag() = rDeltaT*mesh().V();
@@ -417,7 +417,7 @@ CoEulerDdtScheme<Type>::fvmDdt
);
fvMatrix<Type>& fvm = tfvm();
- scalarField rDeltaT = CorDeltaT()().internalField();
+ scalarField rDeltaT(CorDeltaT()().internalField());
fvm.diag() = rDeltaT*rho.value()*mesh().V();
@@ -454,7 +454,7 @@ CoEulerDdtScheme<Type>::fvmDdt
);
fvMatrix<Type>& fvm = tfvm();
- scalarField rDeltaT = CorDeltaT()().internalField();
+ scalarField rDeltaT(CorDeltaT()().internalField());
fvm.diag() = rDeltaT*rho.internalField()*mesh().V();
diff --git a/src/finiteVolume/finiteVolume/ddtSchemes/SLTSDdtScheme/SLTSDdtScheme.C b/src/finiteVolume/finiteVolume/ddtSchemes/SLTSDdtScheme/SLTSDdtScheme.C
index 567df96452cf979d97be73cf9c49e25216977fab..af50dc93cead4e2318b76b243d8f594817ade56c 100644
--- a/src/finiteVolume/finiteVolume/ddtSchemes/SLTSDdtScheme/SLTSDdtScheme.C
+++ b/src/finiteVolume/finiteVolume/ddtSchemes/SLTSDdtScheme/SLTSDdtScheme.C
@@ -385,7 +385,7 @@ SLTSDdtScheme<Type>::fvmDdt
fvMatrix<Type>& fvm = tfvm();
- scalarField rDeltaT = SLrDeltaT()().internalField();
+ scalarField rDeltaT(SLrDeltaT()().internalField());
Info<< "max/min rDeltaT " << max(rDeltaT) << " " << min(rDeltaT) << endl;
@@ -422,7 +422,7 @@ SLTSDdtScheme<Type>::fvmDdt
);
fvMatrix<Type>& fvm = tfvm();
- scalarField rDeltaT = SLrDeltaT()().internalField();
+ scalarField rDeltaT(SLrDeltaT()().internalField());
fvm.diag() = rDeltaT*rho.value()*mesh().V();
@@ -459,7 +459,7 @@ SLTSDdtScheme<Type>::fvmDdt
);
fvMatrix<Type>& fvm = tfvm();
- scalarField rDeltaT = SLrDeltaT()().internalField();
+ scalarField rDeltaT(SLrDeltaT()().internalField());
fvm.diag() = rDeltaT*rho.internalField()*mesh().V();
diff --git a/src/finiteVolume/fvMatrices/fvMatrix/fvMatrixSolve.C b/src/finiteVolume/fvMatrices/fvMatrix/fvMatrixSolve.C
index 0b6f3f7d69614348d412ffb98f7cf74a73594f1b..bc320829cdd4ff89d5829e5b24fcf31aa33d6046 100644
--- a/src/finiteVolume/fvMatrices/fvMatrix/fvMatrixSolve.C
+++ b/src/finiteVolume/fvMatrices/fvMatrix/fvMatrixSolve.C
@@ -71,9 +71,9 @@ Foam::lduMatrix::solverPerformance Foam::fvMatrix<Type>::solve
psi.name()
);
- scalarField saveDiag = diag();
+ scalarField saveDiag(diag());
- Field<Type> source = source_;
+ Field<Type> source(source_);
// At this point include the boundary source from the coupled boundaries.
// This is corrected for the implict part by updateMatrixInterfaces within
@@ -227,7 +227,7 @@ Foam::tmp<Foam::Field<Type> > Foam::fvMatrix<Type>::residual() const
// Loop over field components
for (direction cmpt=0; cmpt<Type::nComponents; cmpt++)
{
- scalarField psiCmpt = psi_.internalField().component(cmpt);
+ scalarField psiCmpt(psi_.internalField().component(cmpt));
scalarField boundaryDiagCmpt(psi_.size(), 0.0);
addBoundaryDiag(boundaryDiagCmpt, cmpt);
diff --git a/src/finiteVolume/fvMatrices/fvScalarMatrix/fvScalarMatrix.C b/src/finiteVolume/fvMatrices/fvScalarMatrix/fvScalarMatrix.C
index e877b51af547d8902e257d39eca11ae55eb286b4..f0c6404abacd414a9c1df306d8076f955a1e694e 100644
--- a/src/finiteVolume/fvMatrices/fvScalarMatrix/fvScalarMatrix.C
+++ b/src/finiteVolume/fvMatrices/fvScalarMatrix/fvScalarMatrix.C
@@ -66,7 +66,7 @@ Foam::fvMatrix<Foam::scalar>::solver
<< endl;
}
- scalarField saveDiag = diag();
+ scalarField saveDiag(diag());
addBoundaryDiag(diag(), 0);
autoPtr<fvMatrix<scalar>::fvSolver> solverPtr
@@ -102,10 +102,10 @@ Foam::lduMatrix::solverPerformance Foam::fvMatrix<Foam::scalar>::fvSolver::solve
const_cast<GeometricField<scalar, fvPatchField, volMesh>&>
(fvMat_.psi());
- scalarField saveDiag = fvMat_.diag();
+ scalarField saveDiag(fvMat_.diag());
fvMat_.addBoundaryDiag(fvMat_.diag(), 0);
- scalarField totalSource = fvMat_.source();
+ scalarField totalSource(fvMat_.source());
fvMat_.addBoundarySource(totalSource, false);
// assign new solver controls
@@ -145,10 +145,10 @@ Foam::lduMatrix::solverPerformance Foam::fvMatrix<Foam::scalar>::solve
GeometricField<scalar, fvPatchField, volMesh>& psi =
const_cast<GeometricField<scalar, fvPatchField, volMesh>&>(psi_);
- scalarField saveDiag = diag();
+ scalarField saveDiag(diag());
addBoundaryDiag(diag(), 0);
- scalarField totalSource = source_;
+ scalarField totalSource(source_);
addBoundarySource(totalSource, false);
// Solver call
diff --git a/src/thermophysicalModels/radiation/derivedFvPatchFields/wideBandDiffusiveRadiation/wideBandDiffusiveRadiationMixedFvPatchScalarField.C b/src/thermophysicalModels/radiation/derivedFvPatchFields/wideBandDiffusiveRadiation/wideBandDiffusiveRadiationMixedFvPatchScalarField.C
index e18684d9f40cf8d1bf24f2e64f79fef7fa4de313..8de1c8a2b21faf79fd2b2ff74a0d3e65a7394a9c 100644
--- a/src/thermophysicalModels/radiation/derivedFvPatchFields/wideBandDiffusiveRadiation/wideBandDiffusiveRadiationMixedFvPatchScalarField.C
+++ b/src/thermophysicalModels/radiation/derivedFvPatchFields/wideBandDiffusiveRadiation/wideBandDiffusiveRadiationMixedFvPatchScalarField.C
@@ -169,8 +169,10 @@ updateCoeffs()
ray.Qr().boundaryField()[patchI] += Iw*(n & ray.dAve());
- const scalarField Eb =
- dom.blackBody().bLambda(lambdaId).boundaryField()[patchI];
+ const scalarField Eb
+ (
+ dom.blackBody().bLambda(lambdaId).boundaryField()[patchI]
+ );
forAll(Iw, faceI)
{
diff --git a/src/turbulenceModels/incompressible/LES/Smagorinsky2/Smagorinsky2.C b/src/turbulenceModels/incompressible/LES/Smagorinsky2/Smagorinsky2.C
index 95574f28ef7ba39f3ed64d40ae03a4c536557eb5..cfd9e5505825e1b2ad56930d81af62136236bb30 100644
--- a/src/turbulenceModels/incompressible/LES/Smagorinsky2/Smagorinsky2.C
+++ b/src/turbulenceModels/incompressible/LES/Smagorinsky2/Smagorinsky2.C
@@ -75,7 +75,7 @@ Smagorinsky2::Smagorinsky2
tmp<volSymmTensorField> Smagorinsky2::B() const
{
- volSymmTensorField D = dev(symm(fvc::grad(U())));
+ volSymmTensorField D(dev(symm(fvc::grad(U()))));
return (((2.0/3.0)*I)*k() - 2.0*nuSgs_*D - (2.0*cD2_)*delta()*(D&D));
}
@@ -86,7 +86,7 @@ tmp<fvVectorMatrix> Smagorinsky2::divDevBeff
volVectorField& U
) const
{
- volTensorField gradU = fvc::grad(U);
+ volTensorField gradU(fvc::grad(U));
volSymmTensorField aniNuEff
(