11 scalar phaseVolume = 0;
15 label celli = cellLabels[cli];
16 phaseVolume += alpha[celli]*V[celli];
19 reduce(phaseVolume, sumOp<scalar>());
22 <<
" phase volume in zone " <<
mesh.cellZones()[czi].name()
23 <<
" = " << phaseVolume*1e6 <<
" ml " <<
endl;