applications/utilities/postProcessing/miscellaneous/postChannel/readFields.H

   1     IOobject UMeanHeader
   2     (
   3         "UMean",
   4         runTime.timeName(),
   5         mesh,
   6         IOobject::MUST_READ
   7     );
   8
   9     if (!UMeanHeader.headerOk())
  10     {
  11         Info<< "    No UMean field" << endl;
  12         continue;
  13     }
  14
  15     volVectorField UMean
  16     (
  17         UMeanHeader,
  18         mesh
  19     );
  20
  21     volSymmTensorField UPrime2Mean
  22     (
  23         IOobject
  24         (
  25             "UPrime2Mean",
  26             runTime.timeName(),
  27             mesh,
  28             IOobject::MUST_READ
  29         ),
  30         mesh
  31     );
  32     volScalarField Rxx(UPrime2Mean.component(symmTensor::XX));
  33     volScalarField Ryy(UPrime2Mean.component(symmTensor::YY));
  34     volScalarField Rzz(UPrime2Mean.component(symmTensor::ZZ));
  35     volScalarField Rxy(UPrime2Mean.component(symmTensor::XY));
  36
  37     volScalarField pPrime2Mean
  38     (
  39         IOobject
  40         (
  41             "pPrime2Mean",
  42             runTime.timeName(),
  43             mesh,
  44             IOobject::MUST_READ
  45         ),
  46         mesh
  47     );
  48
  49     /*
  50     volScalarField epsilonMean
  51     (
  52         IOobject
  53         (
  54             "epsilonMean",
  55             runTime.timeName(),
  56             mesh,
  57             IOobject::MUST_READ
  58         ),
  59         mesh
  60     );
  61
  62     volScalarField nuMean
  63     (
  64         IOobject
  65         (
  66             "nuMean",
  67             runTime.timeName(),
  68             mesh,
  69             IOobject::MUST_READ
  70         ),
  71         mesh
  72     );
  73
  74     volScalarField gammaDotMean
  75     (
  76         IOobject
  77         (
  78             "gammaDotMean",
  79             runTime.timeName(),
  80             mesh,
  81             IOobject::MUST_READ
  82         ),
  83         mesh
  84     );
  85
  86     volScalarField nuPrime2
  87     (
  88         IOobject
  89         (
  90             "nuPrime",
  91             runTime.timeName(),
  92             mesh,
  93             IOobject::MUST_READ
  94         ),
  95         mesh
  96     );
  97     volScalarField nuPrime = sqrt(mag(nuPrime2 - sqr(nuMean)));
  98
  99
 100     volScalarField gammaDotPrime2
 101     (
 102         IOobject
 103         (
 104             "gammaDotPrime",
 105             runTime.timeName(),
 106             mesh,
 107             IOobject::MUST_READ
 108         ),
 109         mesh
 110     );
 111     volScalarField gammaDotPrime = sqrt(mag(gammaDotPrime2 -sqr(gammaDotMean)));
 112     */