1 if (nVar >= Num_variables - nSprayVariables)
7 IOobject fieldObjectPtr
9 fieldNames[var2field[nVar]],
15 if (!fieldObjectPtr.headerOk())
22 fieldNames[var2field[nVar]],
35 const cellShapeList& cellShapes = meshPtr->cellShapes();
38 if (which_type == Z_HEX08)
40 const cellModel& hex = *(cellModeller::lookup("hex"));
41 //const cellModel& wedge = *(cellModeller::lookup("wedge"));
44 for (label n=0; n<nCells; n++)
46 const cellShape& cellShape = cellShapes[n];
47 const cellModel& cellModel = cellShape.model();
49 if (cellModel == hex) // || (cellModel == wedge))
51 # include "tensorConversion.H"
57 if (which_type == Z_PEN06)
59 const cellModel& prism = *(cellModeller::lookup("prism"));
62 for (label n=0; n<nCells; n++)
64 const cellShape& cellShape = cellShapes[n];
65 const cellModel& cellModel = cellShape.model();
67 if (cellModel == prism)
69 # include "tensorConversion.H"
75 if (which_type == Z_PYR05)
77 const cellModel& pyr = *(cellModeller::lookup("pyr"));
80 for (label n=0; n<nCells; n++)
82 const cellShape& cellShape = cellShapes[n];
83 const cellModel& cellModel = cellShape.model();
87 # include "tensorConversion.H"
94 if (which_type == Z_TET04)
96 const cellModel& tet = *(cellModeller::lookup("tet"));
100 for (label n=0; n<nCells; n++)
102 const cellShape& cellShape = cellShapes[n];
103 const cellModel& cellModel = cellShape.model();
105 if (cellModel == tet)
107 # include "tensorConversion.H"
112 if (which_type == Z_NFACED)
114 const cellList& cells = meshPtr->cells();
117 for (label n=0; n<nCells; n++)
119 const labelList& points = cellShapes[n];
120 label nFacesInCell = cells[n].size();
122 if ((nFacesInCell == 6) && (points.size() == 8))
124 else if ((nFacesInCell == 4) && (points.size() == 4))
126 else if (nFacesInCell == 5)
128 if (points.size() == 6)
130 else if (points.size() == 5)
134 # include "tensorConversion.H"
139 # include "tensorConversion.H"