applications/solvers/incompressible/nonNewtonianIcoFoam/createFields.H

   1     Info<< "Reading field p\n" << endl;
   2     volScalarField p
   3     (
   4         IOobject
   5         (
   6             "p",
   7             runTime.timeName(),
   8             mesh,
   9             IOobject::MUST_READ,
  10             IOobject::AUTO_WRITE
  11         ),
  12         mesh
  13     );
  14
  15
  16     Info<< "Reading field U\n" << endl;
  17     volVectorField U
  18     (
  19         IOobject
  20         (
  21             "U",
  22             runTime.timeName(),
  23             mesh,
  24             IOobject::MUST_READ,
  25             IOobject::AUTO_WRITE
  26         ),
  27         mesh
  28     );
  29
  30
  31 #   include "createPhi.H"
  32
  33
  34     singlePhaseTransportModel fluid(U, phi);
  35
  36
  37     label pRefCell = 0;
  38     scalar pRefValue = 0.0;
  39     setRefCell(p, mesh.solutionDict().subDict("PISO"), pRefCell, pRefValue);