FreeFOAM The Cross-Platform CFD Toolkit
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 
47  surfaceScalarField& phiB = *phiBPtr;
48 
49 // ************************ vim: set sw=4 sts=4 et: ************************ //