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

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