1 surfaceScalarField::GeometricBoundaryField& phivPatches = phiv.boundaryField();
2 const volVectorField::GeometricBoundaryField& Upatches = U.boundaryField();
3 const surfaceVectorField::GeometricBoundaryField& SfPatches = mesh.Sf().boundaryField();
5 forAll(phivPatches, patchI)
7 if (phiv.boundaryField().types()[patchI] == "calculated")
9 calculatedFvsPatchScalarField& phivPatch =
10 refCast<calculatedFvsPatchScalarField>(phivPatches[patchI]);
12 phivPatch == (Upatches[patchI] & SfPatches[patchI]);