38 int main(
int argc,
char *argv[])
40 argList::noParallel();
41 argList::validArgs.append(
"SMAP fileName");
56 fieldNameDict.add(fieldNames[i],
word(fieldNames[i]));
60 if (fieldNameDict.found(
"U")) nameMap.
add(
"SU",
word(
"U"));
61 if (fieldNameDict.found(
"p")) nameMap.
add(
"P",
word(
"p"));
62 if (fieldNameDict.found(
"T")) nameMap.
add(
"T",
word(
"T"));
63 if (fieldNameDict.found(
"rho")) nameMap.
add(
"DENS",
word(
"rho"));
64 if (fieldNameDict.found(
"k")) nameMap.
add(
"TE",
word(
"k"));
65 if (fieldNameDict.found(
"epsilon")) nameMap.
add(
"ED",
word(
"epsilon"));
66 if (fieldNameDict.found(
"nuEff")) nameMap.
add(
"VIS",
word(
"nuEff"));
75 <<
"Cannot open SMAP file " << smapFile.name()
79 while (!smapFile.eof())
83 token fieldName(smapFile);
92 fieldName.type() != token::WORD
93 && fieldName.wordToken() !=
"CELL"
97 <<
"Expected first CELL, found "
103 smapFile >> fieldName;
104 while (fieldName.type() == token::WORD)
106 starFieldNames[nCols++] = fieldName.wordToken();
107 smapFile >> fieldName;
113 reinterpret_cast<volScalarField*>(0)
119 reinterpret_cast<volVectorField*>(0)
125 if (nameMap.
found(starFieldNames[i]))
127 if (starFieldNames[i] ==
"SU")
134 nameMap.
lookup(starFieldNames[i]),
152 nameMap.
lookup(starFieldNames[i]),
185 smapFile >> (*sFields[i])[celli];
190 smapFile >> (*vFields[i])[celli].x();
191 smapFile >> (*vFields[i])[celli].
y();
192 smapFile >> (*vFields[i])[celli].z();
203 for (label i=0; i<nCols; i++)
207 sFields[i]->correctBoundaryConditions();
214 vFields[i]->correctBoundaryConditions();
227 <<
"Expected first SMAP dummy entry to be cell 0, found "
232 for (label i=0; i<nCols; i++)