47 rhoUf().boundaryFieldRef()[patchI] =
rhoUfint.boundaryField()[patchI];
static tmp< GeometricField< Type, faePatchField, edgeMesh > > interpolate(const GeometricField< Type, faPatchField, areaMesh > &tvf, const edgeScalarField &faceFlux, Istream &schemeData)
Interpolate field onto faces using scheme given by Istream.
GeometricField< vector, fvsPatchField, surfaceMesh > surfaceVectorField
forAll(rhoUf().boundaryField(), patchI)
surfaceVectorField rhoUfint(fvc::interpolate(rho *U))
const surfaceScalarField faceMaskOld(localMin< scalar >(mesh).interpolate(cellMask.oldTime()))
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField