Commit 1202a245 authored by henry's avatar henry
Browse files

Changed the constructor for pd to MUST_READ.

Calculate p from pd (note p is still read in by the termo package) and correct thermo.
parent cb08f660
......@@ -67,10 +67,14 @@
(
"pd",
runTime.timeName(),
mesh
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
p - rho*gh - pRef,
p.boundaryField().types()
mesh
);
p = pd + rho*gh + pRef;
thermo->correct();
dimensionedScalar initialMass = fvc::domainIntegrate(rho);
......@@ -64,12 +64,16 @@
(
"pd",
runTime.timeName(),
mesh
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
p - rho*gh - pRef,
p.boundaryField().types()
mesh
);
p = pd + rho*gh + pRef;
thermo->correct();
label pdRefCell = 0;
scalar pdRefValue = 0.0;
......
......@@ -65,12 +65,16 @@
(
"pd",
runTime.timeName(),
mesh
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
p - rho*gh - pRef,
p.boundaryField().types()
mesh
);
p = pd + rho*gh + pRef;
thermo->correct();
label pdRefCell = 0;
scalar pdRefValue = 0.0;
......
......@@ -173,7 +173,7 @@
Info<< " Updating p from pd\n" << endl;
thermof[i].p() == pdf[i] + rhof[i]*ghf[i] + pRef;
thermof[i].correct();
initialMassf[i] = fvc::domainIntegrate(rhof[i]).value();
}
......
......@@ -37,7 +37,7 @@
mesh
);
# include "compressibleCreatePhi.H"
#include "compressibleCreatePhi.H"
Info<< "Creating turbulence model\n" << endl;
......@@ -69,7 +69,7 @@
mesh
);
//pd = p - rho*gh - pRef;
p = pd + rho*gh + pRef;
thermo->correct();
dimensionedScalar initialMass = fvc::domainIntegrate(rho);
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment