6 const labelList& cellLabels =
mesh.cellZones()[czi];
8 for (
const volScalarField&
alpha :
fluid.phases())
10 scalar phaseVolume = 0;
12 for (
const label celli : cellLabels)
14 phaseVolume +=
alpha[celli]*V[celli];
17 reduce(phaseVolume, sumOp<scalar>());
20 <<
" phase volume in zone " <<
mesh.cellZones()[czi].name()
21 <<
" = " << phaseVolume*1e6 <<
" ml " << endl;
reduce(hasMovingMesh, orOp< bool >())
Field< scalar > scalarField
Specialisation of Field<T> for scalar.
#define forAll(list, i)
Loop across all elements in list.