1 Info<< "Reading field p\n" << endl;
15 p = dimensionedScalar("zero", p.dimensions(), 0.0);
18 Info<< "Reading field U\n" << endl;
32 U = dimensionedVector("0", U.dimensions(), vector::zero);
34 surfaceScalarField phi
44 fvc::interpolate(U) & mesh.Sf()
49 scalar pRefValue = 0.0;