applications/solvers/incompressible/icoFoam/createFields.H

   1     Info<< "Reading transportProperties\n" << endl;
   2
   3     IOdictionary transportProperties
   4     (
   5         IOobject
   6         (
   7             "transportProperties",
   8             runTime.constant(),
   9             mesh,
  10             IOobject::MUST_READ_IF_MODIFIED,
  11             IOobject::NO_WRITE
  12         )
  13     );
  14
  15     dimensionedScalar nu
  16     (
  17         transportProperties.lookup("nu")
  18     );
  19
  20     Info<< "Reading field p\n" << endl;
  21     volScalarField p
  22     (
  23         IOobject
  24         (
  25             "p",
  26             runTime.timeName(),
  27             mesh,
  28             IOobject::MUST_READ,
  29             IOobject::AUTO_WRITE
  30         ),
  31         mesh
  32     );
  33
  34
  35     Info<< "Reading field U\n" << endl;
  36     volVectorField U
  37     (
  38         IOobject
  39         (
  40             "U",
  41             runTime.timeName(),
  42             mesh,
  43             IOobject::MUST_READ,
  44             IOobject::AUTO_WRITE
  45         ),
  46         mesh
  47     );
  48
  49
  50 #   include "createPhi.H"
  51
  52
  53     label pRefCell = 0;
  54     scalar pRefValue = 0.0;
  55     setRefCell(p, mesh.solutionDict().subDict("PISO"), pRefCell, pRefValue);