applications/solvers/surfaceTracking/include/freeSurfaceContinuityErrs.H

   1 // if(Pstream::master())
   2 // {
   3     Info<< "Free surface flux: sum local = "
   4         << gSum(mag(phi.boundaryField()[interface.aPatchID()]))
   5         << ", global = "
   6         << gSum(phi.boundaryField()[interface.aPatchID()]) << endl;
   7
   8     scalarField netPhiA =
   9         phi.boundaryField()[interface.aPatchID()]
  10       - fvc::meshPhi(rho,U)().boundaryField()[interface.aPatchID()];
  11
  12     Info<< "Free surface continuity error : sum local = "
  13         << gSum(mag(netPhiA)) << ", global = " << gSum(netPhiA)
  14         << endl << endl;
  15 // }