1 const volVectorField::GeometricBoundaryField& Ubf = U.boundaryField();
2 wordList rhoUboundaryTypes = Ubf.types();
4 forAll(rhoUboundaryTypes, patchi)
6 if (Ubf[patchi].fixesValue())
8 rhoUboundaryTypes[patchi] = fixedRhoUFvPatchVectorField::typeName;