applications/solvers/compressible/sonicLiquidFoam/readThermodynamicProperties.H

   1     Info<< "Reading thermodynamicProperties\n" << endl;
   2
   3     IOdictionary thermodynamicProperties
   4     (
   5         IOobject
   6         (
   7             "thermodynamicProperties",
   8             runTime.constant(),
   9             mesh,
  10             IOobject::MUST_READ,
  11             IOobject::NO_WRITE
  12         )
  13     );
  14
  15     dimensionedScalar rho0
  16     (
  17         thermodynamicProperties.lookup("rho0")
  18     );
  19
  20     dimensionedScalar p0
  21     (
  22         thermodynamicProperties.lookup("p0")
  23     );
  24
  25     dimensionedScalar psi
  26     (
  27         thermodynamicProperties.lookup("psi")
  28     );
  29
  30     // Density offset, i.e. the constant part of the density
  31     dimensionedScalar rhoO("rhoO", rho0 - psi*p0);