createPhiB.H
Go to the documentation of this file.
1  IOobject phiBHeader
2  (
3  "phiB",
4  runTime.timeName(),
5  mesh,
6  IOobject::NO_READ
7  );
8 
9 
11 
12  if (phiBHeader.headerOk())
13  {
14  Info<< "Reading field phiB\n" << endl;
15 
16  phiBPtr = new surfaceScalarField
17  (
18  IOobject
19  (
20  "phiB",
21  runTime.timeName(),
22  mesh,
23  IOobject::MUST_READ,
24  IOobject::AUTO_WRITE
25  ),
26  mesh
27  );
28  }
29  else
30  {
31  Info<< "Calculating face flux\n" << endl;
32 
33  phiBPtr = new surfaceScalarField
34  (
35  IOobject
36  (
37  "phiB",
38  runTime.timeName(),
39  mesh,
40  IOobject::NO_READ,
41  IOobject::AUTO_WRITE
42  ),
43  (fvc::interpolate(B) & mesh.Sf())
44  );
45  }
46 
messageStream Info
tmp< surfaceScalarField > interpolate(const RhoType &rho)
dynamicFvMesh & mesh
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
surfaceScalarField * phiBPtr
Definition: createPhiB.H:10
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
IOobject phiBHeader("phiB", runTime.timeName(), mesh, IOobject::NO_READ)