OpenFOAM-4.x/applications/solvers/electromagnetics/mhdFoam/createPhiB.H
2016-05-09 14:08:15 +01:00

49 lines
797 B
C

IOobject phiBHeader
(
"phiB",
runTime.timeName(),
mesh,
IOobject::NO_READ
);
surfaceScalarField* phiBPtr = NULL;
if (phiBHeader.headerOk())
{
Info<< "Reading face flux ";
phiBPtr = new surfaceScalarField
(
IOobject
(
"phiB",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);
}
else
{
Info<< "Calculating face flux ";
phiBPtr = new surfaceScalarField
(
IOobject
(
"phiB",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
fvc::flux(B)
);
}
surfaceScalarField& phiB = *phiBPtr;
Info<< phiB.name() << nl << endl;