diff --git a/src/finiteVolume/finiteVolume/ddtSchemes/backwardDdtScheme/backwardDdtScheme.H b/src/finiteVolume/finiteVolume/ddtSchemes/backwardDdtScheme/backwardDdtScheme.H index 24812aee91..0d7e9174d9 100644 --- a/src/finiteVolume/finiteVolume/ddtSchemes/backwardDdtScheme/backwardDdtScheme.H +++ b/src/finiteVolume/finiteVolume/ddtSchemes/backwardDdtScheme/backwardDdtScheme.H @@ -101,7 +101,7 @@ public: { if (is.good() && !is.eof()) { - this->ddtPhiCoeff_ = readScalar(is); + is >> this->ddtPhiCoeff_; } // Ensure the old-old-time cell volumes are available diff --git a/src/finiteVolume/finiteVolume/ddtSchemes/ddtScheme/ddtScheme.C b/src/finiteVolume/finiteVolume/ddtSchemes/ddtScheme/ddtScheme.C index a1ef450c2d..f5e11840b9 100644 --- a/src/finiteVolume/finiteVolume/ddtSchemes/ddtScheme/ddtScheme.C +++ b/src/finiteVolume/finiteVolume/ddtSchemes/ddtScheme/ddtScheme.C @@ -171,10 +171,11 @@ tmp ddtScheme::fvcDdtPhiCoeff ) ); + surfaceScalarField& ddtCouplingCoeff = tddtCouplingCoeff.ref(); if (ddtPhiCoeff_ < 0) { - tddtCouplingCoeff.ref() =- min + ddtCouplingCoeff -= min ( mag(phiCorr) /(mag(phi) + dimensionedScalar("small", phi.dimensions(), SMALL)), @@ -183,14 +184,11 @@ tmp ddtScheme::fvcDdtPhiCoeff } else { - tddtCouplingCoeff.ref() = + ddtCouplingCoeff = dimensionedScalar("ddtPhiCoeff", dimless, ddtPhiCoeff_); } - surfaceScalarField& ddtCouplingCoeff = tddtCouplingCoeff.ref(); - - surfaceScalarField::Boundary& ccbf = - ddtCouplingCoeff.boundaryFieldRef(); + surfaceScalarField::Boundary& ccbf = ddtCouplingCoeff.boundaryFieldRef(); forAll(U.boundaryField(), patchi) {