applications/solvers/compressible/sonicDyMFoam/compressibleContinuityErrs.H

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