diff --git a/src/functionObjects/field/pressure/pressure.C b/src/functionObjects/field/pressure/pressure.C index e923c83975fb184bcd4a3cde388edce5baba258d..df050772ae727de451f98f9c0e1152131e0bff40 100644 --- a/src/functionObjects/field/pressure/pressure.C +++ b/src/functionObjects/field/pressure/pressure.C @@ -157,10 +157,8 @@ Foam::tmp<Foam::volScalarField> Foam::functionObjects::pressure::calcPressure } case ISENTROPIC: { -DebugVar("getting thermo"); const basicThermo* thermoPtr = p.mesh().lookupObjectPtr<basicThermo>(basicThermo::dictName); -DebugVar("got thermo"); if (!thermoPtr) { @@ -169,16 +167,13 @@ DebugVar("got thermo"); << "thermodynamics package" << exit(FatalError); } -DebugVar("calc gamma"); const volScalarField gamma(thermoPtr->gamma()); -DebugVar(gamma); const volScalarField Mb ( mag(lookupObject<volVectorField>(UName_)) /sqrt(gamma*tp.ref()/thermoPtr->rho()) ); -DebugVar(Mb); return tp()*(pow(1 + (gamma - 1)/2*sqr(Mb), gamma/(gamma - 1))); }