applications/solvers/electromagnetics/mhdFoam/createPhiB.H

   1     IOobject phiBHeader
   2     (
   3         "phiB",
   4         runTime.timeName(),
   5         mesh,
   6         IOobject::NO_READ
   7     );
   8
   9
  10     surfaceScalarField* phiBPtr;
  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;