applications/solvers/multiphase/settlingFoam/compressibleContinuityErrs.H

   1     scalar sumLocalContErr =
   2         runTime.deltaT().value()*
   3         mag
   4         (
   5             fvc::ddt(rho)
   6           + fvc::div(phi)
   7         )().weightedAverage(rho*mesh.V()).value();
   8
   9     scalar globalContErr =
  10           runTime.deltaT().value()*
  11           (
  12               fvc::ddt(rho)
  13             + fvc::div(phi)
  14           )().weightedAverage(rho*mesh.V()).value();
  15
  16     cumulativeContErr += globalContErr;
  17
  18     Info<< "time step continuity errors : sum local = " << sumLocalContErr
  19         << ", global = " << globalContErr
  20         << ", cumulative = " << cumulativeContErr
  21         << endl;