diff --git a/src/postProcessing/functionObjects/field/fieldAverage/fieldAverage/fieldAverageTemplates.C b/src/postProcessing/functionObjects/field/fieldAverage/fieldAverage/fieldAverageTemplates.C index a69dce41b7..4e59e8a30b 100644 --- a/src/postProcessing/functionObjects/field/fieldAverage/fieldAverage/fieldAverageTemplates.C +++ b/src/postProcessing/functionObjects/field/fieldAverage/fieldAverage/fieldAverageTemplates.C @@ -63,21 +63,22 @@ void Foam::fieldAverage::addMeanField const fieldType& baseField = obr_.lookupObject(fieldName); - fieldType* fPtr = new fieldType - ( - IOobject - ( - meanFieldName, - obr_.time().timeName(), - obr_, - IOobject::READ_IF_PRESENT, - IOobject::NO_WRITE - ), - baseField - ); - // Store on registry - fPtr->store(); + obr_.store + ( + new fieldType + ( + IOobject + ( + meanFieldName, + obr_.time().timeName(), + obr_, + IOobject::READ_IF_PRESENT, + IOobject::NO_WRITE + ), + 1*baseField + ) + ); meanFieldList[fieldI] = meanFieldName; } @@ -121,22 +122,22 @@ void Foam::fieldAverage::addPrime2MeanField const fieldType1& meanField = obr_.lookupObject(meanFieldList[fieldI]); - fieldType2* fPtr = new fieldType2 + obr_.store ( - IOobject + new fieldType2 ( - meanFieldName, - obr_.time().timeName(), - obr_, - IOobject::READ_IF_PRESENT, - IOobject::NO_WRITE - ), - sqr(baseField) - sqr(meanField) + IOobject + ( + meanFieldName, + obr_.time().timeName(), + obr_, + IOobject::READ_IF_PRESENT, + IOobject::NO_WRITE + ), + sqr(baseField) - sqr(meanField) + ) ); - // Store on registry - fPtr->store(); - prime2MeanFieldList[fieldI] = meanFieldName; } } diff --git a/src/thermophysicalModels/basic/derivedFvPatchFields/fixedEnthalpy/fixedEnthalpyFvPatchScalarField.C b/src/thermophysicalModels/basic/derivedFvPatchFields/fixedEnthalpy/fixedEnthalpyFvPatchScalarField.C index e409d243c8..19e345abe0 100644 --- a/src/thermophysicalModels/basic/derivedFvPatchFields/fixedEnthalpy/fixedEnthalpyFvPatchScalarField.C +++ b/src/thermophysicalModels/basic/derivedFvPatchFields/fixedEnthalpy/fixedEnthalpyFvPatchScalarField.C @@ -102,12 +102,13 @@ void fixedEnthalpyFvPatchScalarField::updateCoeffs() ( "thermophysicalProperties" ); - + const label patchi = patch().index(); - fvPatchScalarField& Tw = + fvPatchScalarField& Tw = const_cast(thermo.T().boundaryField()[patchi]); Tw.evaluate(); + operator==(thermo.h(Tw, patchi)); fixedValueFvPatchScalarField::updateCoeffs();