applications/utilities/postProcessing/graphics/newEnsightFoamReader/getPatchFieldVector.H

   1 label patchi = which_part - 2;
   2
   3 if
   4 (
   5     nVar >= Num_variables - nSprayVariables
   6   - nSecondMeshVariables - nFaMeshVariables
   7 )
   8 {
   9     return Z_UNDEF;
  10 }
  11
  12 IOobject fieldObjectPtr
  13 (
  14     fieldNames[var2field[nVar]],
  15     runTime.timeName(),
  16     mesh,
  17     IOobject::NO_READ
  18 );
  19
  20 if (!fieldObjectPtr.headerOk())
  21 {
  22     return Z_UNDEF;
  23 }
  24
  25 IOobject fieldObject
  26 (
  27     fieldNames[var2field[nVar]],
  28     runTime.timeName(),
  29     mesh,
  30     IOobject::MUST_READ,
  31     IOobject::NO_WRITE
  32 );
  33
  34 volVectorField sf
  35 (
  36     fieldObject,
  37     mesh
  38 );
  39
  40 const vectorField& sfb = sf.boundaryField()[patchi];
  41 const polyBoundaryMesh& bMesh = meshPtr->boundaryMesh();
  42
  43 if (which_type == Z_TRI03)
  44 {
  45     label counter = 1;
  46     for (label facei=0; facei<sfb.size(); facei++)
  47     {
  48         label nPoints = bMesh[patchi][facei].size();
  49         if (nPoints == 3)
  50         {
  51             var_array[counter++] = sfb[facei][component];
  52         }
  53     }
  54 }
  55
  56 if (which_type == Z_QUA04)
  57 {
  58     label counter = 1;
  59     for (label facei=0; facei<sfb.size(); facei++)
  60     {
  61         label nPoints = bMesh[patchi][facei].size();
  62         if (nPoints == 4)
  63         {
  64             var_array[counter++] = sfb[facei][component];
  65         }
  66     }
  67 }
  68
  69 if (which_type == Z_NSIDED)
  70 {
  71     label counter = 1;
  72     for (label facei=0; facei<sfb.size(); facei++)
  73     {
  74         label nPoints = bMesh[patchi][facei].size();
  75         if ((nPoints != 3) && (nPoints != 4))
  76         {
  77             var_array[counter++] = sfb[facei][component];
  78         }
  79     }
  80 }