void lcsMassSourcePluginFunction::doEvaluation()
{
typedef DimensionedField<scalar,volMesh> dimScalarField;
autoPtr<dimScalarField> pSrho;
// pick up the first fitting class
castAndCall(pSrho,dimScalarField,basicReactingCloud,reactingCloud,Srho());
castAndCall(pSrho,dimScalarField,basicReactingMultiphaseCloud,reactingMultiphaseCloud,Srho());
noCloudFound(pSrho);
const dimScalarField &Srho=pSrho();
autoPtr<volScalarField> pSource(
new volScalarField(
IOobject(
cloudName()+"MassSource",
mesh().time().timeName(),
mesh(),
IOobject::NO_READ,
IOobject::NO_WRITE
),
mesh(),
Srho.dimensions(),
"zeroGradient"
)
);
pSource->internalField()=Srho.field();
result().setObjectResult(pSource);
}
void lcsSpeciesSourcePluginFunction::doEvaluation()
{
wordList speciesList;
label speciesIndex=getIndex(speciesList);
if(speciesList.size()==0) {
listAvailableClouds(Info);
FatalErrorIn("lcsSpeciesSourcePluginFunction::doEvaluation()")
<< "No species list found (probably wrong cloud type)"
<< endl
<< "Supported cloud types (typeid): " \
<< supportedTypes() << endl \
<< exit(FatalError);
}
if(speciesIndex<0) {
listAvailableClouds(Info);
FatalErrorIn("lcsSpeciesSourcePluginFunction::doEvaluation()")
<< "Species " << speciesName_ << " not found in gas composition "
<< speciesList
<< endl
<< "Supported cloud types (typeid): " \
<< supportedTypes() << endl \
<< exit(FatalError);
}
autoPtr<dimScalarField> pSrho=internalEvaluate(speciesIndex);;
noCloudFound(pSrho);
const dimScalarField &Srho=pSrho();
autoPtr<volScalarField> pSource(
new volScalarField(
IOobject(
cloudName()+"SpeciesSource",
mesh().time().timeName(),
mesh(),
IOobject::NO_READ,
IOobject::NO_WRITE
),
mesh(),
Srho.dimensions(),
"zeroGradient"
)
);
pSource->internalField()=Srho.field();
result().setObjectResult(pSource);
}