applications/solvers/compressible/rhoSonicFoam/compressibleContinuityErrs.H

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