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

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