Skip to content
Snippets Groups Projects
Commit 9fae7f30 authored by Andrew Heather's avatar Andrew Heather
Browse files

BUG: Removed Debug statements

parent fec74af0
No related merge requests found
...@@ -157,10 +157,8 @@ Foam::tmp<Foam::volScalarField> Foam::functionObjects::pressure::calcPressure ...@@ -157,10 +157,8 @@ Foam::tmp<Foam::volScalarField> Foam::functionObjects::pressure::calcPressure
} }
case ISENTROPIC: case ISENTROPIC:
{ {
DebugVar("getting thermo");
const basicThermo* thermoPtr = const basicThermo* thermoPtr =
p.mesh().lookupObjectPtr<basicThermo>(basicThermo::dictName); p.mesh().lookupObjectPtr<basicThermo>(basicThermo::dictName);
DebugVar("got thermo");
if (!thermoPtr) if (!thermoPtr)
{ {
...@@ -169,16 +167,13 @@ DebugVar("got thermo"); ...@@ -169,16 +167,13 @@ DebugVar("got thermo");
<< "thermodynamics package" << "thermodynamics package"
<< exit(FatalError); << exit(FatalError);
} }
DebugVar("calc gamma");
const volScalarField gamma(thermoPtr->gamma()); const volScalarField gamma(thermoPtr->gamma());
DebugVar(gamma);
const volScalarField Mb const volScalarField Mb
( (
mag(lookupObject<volVectorField>(UName_)) mag(lookupObject<volVectorField>(UName_))
/sqrt(gamma*tp.ref()/thermoPtr->rho()) /sqrt(gamma*tp.ref()/thermoPtr->rho())
); );
DebugVar(Mb);
return tp()*(pow(1 + (gamma - 1)/2*sqr(Mb), gamma/(gamma - 1))); return tp()*(pow(1 + (gamma - 1)/2*sqr(Mb), gamma/(gamma - 1)));
} }
......
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