{
  volScalarField contErr(fvc::div(phi));

  scalar sumLocalContErr = runTime.deltaTValue()*
  mag(contErr)().weightedAverage(mesh.V()).value();

  scalar globalContErr = runTime.deltaTValue()*
  contErr.weightedAverage(mesh.V()).value();
  cumulativeContErr[i] += globalContErr;

  Info<< "time step continuity errors : sum local = " << sumLocalContErr
  << ", global = " << globalContErr
  << ", cumulative = " << cumulativeContErr[i]
  << endl;
}
