9 IOobject::READ_IF_PRESENT,
IOdictionary porosityProperties(IOobject("porosityProperties", runTime.constant(), runTime, IOobject::READ_IF_PRESENT, IOobject::NO_WRITE))
const bool porosityEnabled(porosityProperties.getOrDefault< bool >("porosityEnabled", false))
tmp< volScalarField > tporosity