applications/solvers/compressible/rhoSimpleFoam/rhoPorousMRFSimpleFoam/createZones.H

   1     MRFZones mrfZones(mesh);
   2     mrfZones.correctBoundaryVelocity(U);
   3
   4     thermalPorousZones pZones(mesh);
   5     Switch pressureImplicitPorosity(false);
   6
   7     // nUCorrectors used for pressureImplicitPorosity
   8     int nUCorr = 0;
   9     if (pZones.size())
  10     {
  11         // nUCorrectors for pressureImplicitPorosity
  12         if (mesh.solutionDict().subDict("SIMPLE").found("nUCorrectors"))
  13         {
  14             nUCorr = readInt
  15             (
  16                 mesh.solutionDict().subDict("SIMPLE").lookup("nUCorrectors")
  17             );
  18         }
  19
  20         if (nUCorr > 0)
  21         {
  22             pressureImplicitPorosity = true;
  23             Info<< "Using pressure implicit porosity" << endl;
  24         }
  25         else
  26         {
  27             Info<< "Using pressure explicit porosity" << endl;
  28         }
  29     }