8 scalar phaseVolume = 0;
12 label celli = cellLabels[cli];
13 phaseVolume += alpha1[celli]*V[celli];
16 reduce(phaseVolume, sumOp<scalar>());
18 Info<<
"Phase volume in zone " <<
mesh.cellZones()[czi].name()
19 <<
" = " << phaseVolume*1e6 <<
" ml " <<
endl;