applications/solvers/compressible/sonicLiquidFoam/compressibleContinuityErrs.H

   1 {
   2 #   include "rhoEqn.H"
   3 }
   4 {
   5     scalar sumLocalContErr =
   6         (sum(mag(rho - rho0 - psi*(p - p0)))/sum(rho)).value();
   7
   8     scalar globalContErr = (sum(rho - rho0 - psi*(p - p0))/sum(rho)).value();
   9
  10     cumulativeContErr += globalContErr;
  11
  12     Info<< "time step continuity errors : sum local = " << sumLocalContErr
  13          << ", global = " << globalContErr
  14          << ", cumulative = " << cumulativeContErr << endl;
  15 }