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

   1 label offset = Num_variables - nSecondMeshVariables - nFaMeshVariables;
   2
   3 if (nVar < offset || nVar >= (Num_variables-nFaMeshVariables))
   4 {
   5     return Z_UNDEF;
   6 }
   7
   8 IOobject fieldObjectPtr
   9 (
  10     secondMeshFieldNames[secondMeshVar2field[nVar-offset]],
  11     runTime.timeName(),
  12     *secondMeshPtr,
  13     IOobject::NO_READ
  14 );
  15
  16 if (!fieldObjectPtr.headerOk())
  17 {
  18     return Z_UNDEF;
  19 }
  20
  21 IOobject fieldObject
  22 (
  23     secondMeshFieldNames[secondMeshVar2field[nVar-offset]],
  24     runTime.timeName(),
  25     *secondMeshPtr,
  26     IOobject::MUST_READ,
  27     IOobject::NO_WRITE
  28 );
  29
  30 volScalarField scalarField
  31 (
  32     fieldObject,
  33     *secondMeshPtr
  34 );
  35
  36 const cellShapeList& cellShapes = secondMeshPtr->cellShapes();
  37
  38 // hexa's
  39 if (which_type == Z_HEX08)
  40 {
  41     const cellModel& hex = *(cellModeller::lookup("hex"));
  42     //const cellModel& wedge = *(cellModeller::lookup("wedge"));
  43
  44     label counter = 1;
  45     for (label celli=0; celli<nCells; celli++)
  46     {
  47         const cellShape& cellShape = cellShapes[celli];
  48         const cellModel& cellModel = cellShape.model();
  49
  50         if (cellModel == hex) // || (cellModel == wedge))
  51         {
  52             var_array[counter++] = scalarField[celli];
  53         }
  54     }
  55 }
  56
  57 // penta's
  58 if (which_type == Z_PEN06)
  59 {
  60     const cellModel& prism = *(cellModeller::lookup("prism"));
  61
  62     label counter = 1;
  63     for (label n=0; n<nCells; n++)
  64     {
  65         const cellShape& cellShape = cellShapes[n];
  66         const cellModel& cellModel = cellShape.model();
  67
  68         if (cellModel == prism)
  69         {
  70             var_array[counter++] = scalarField[n];
  71         }
  72     }
  73 }
  74
  75 // pyramids's
  76 if (which_type == Z_PYR05)
  77 {
  78     const cellModel& pyr = *(cellModeller::lookup("pyr"));
  79
  80     label counter = 1;
  81     for (label n=0; n<nCells; n++)
  82     {
  83         const cellShape& cellShape = cellShapes[n];
  84         const cellModel& cellModel = cellShape.model();
  85
  86         if (cellModel == pyr)
  87         {
  88             var_array[counter++] = scalarField[n];
  89         }
  90     }
  91 }
  92
  93 // tet's
  94 if (which_type == Z_TET04)
  95 {
  96     const cellModel& tet = *(cellModeller::lookup("tet"));
  97
  98     label counter = 1;
  99     for (label n=0; n<nCells; n++)
 100     {
 101         const cellShape& cellShape = cellShapes[n];
 102         const cellModel& cellModel = cellShape.model();
 103
 104         if (cellModel == tet)
 105         {
 106             var_array[counter++] = scalarField[n];
 107         }
 108     }
 109 }
 110
 111 if (which_type == Z_NFACED)
 112 {
 113     const cellList& cells = secondMeshPtr->cells();
 114
 115     label counter = 1;
 116     for (label n=0; n<nCells; n++)
 117     {
 118         const labelList& points = cellShapes[n];
 119         label nFacesInCell = cells[n].size();
 120
 121         if ((nFacesInCell == 6) && (points.size() == 8))
 122         {}
 123         else if ((nFacesInCell == 4) && (points.size() == 4))
 124         {}
 125         else if (nFacesInCell == 5)
 126         {
 127             if (points.size() == 6)
 128             {}
 129             else if (points.size() == 5)
 130             {}
 131             else
 132             {
 133                 var_array[counter++] = scalarField[n];
 134             }
 135         }
 136         else
 137         {
 138             var_array[counter++] = scalarField[n];
 139         }
 140     }
 141
 142 }