applications/solvers/compressible/rhopSonicFoam/resetPhiPatches.H

   1 surfaceScalarField::GeometricBoundaryField& phiPatches = phi.boundaryField();
   2 const volVectorField::GeometricBoundaryField& rhoUpatches = rhoU.boundaryField();
   3 const surfaceVectorField::GeometricBoundaryField& SfPatches = mesh.Sf().boundaryField();
   4
   5 forAll(phiPatches, patchI)
   6 {
   7     if (phi.boundaryField().types()[patchI] == "calculated")
   8     {
   9         calculatedFvsPatchScalarField& phiPatch =
  10             refCast<calculatedFvsPatchScalarField>(phiPatches[patchI]);
  11
  12         phiPatch == (rhoUpatches[patchI] & SfPatches[patchI]);
  13     }
  14 }