32 using reportFields = foamToVtkReportFields;
46 sScalars = readFields<surfaceScalarField>(meshProxy, objects);
48 reportFields::print(
" surfScalar :", Info,
sScalars);
58 sVectors = readFields<surfaceVectorField>(meshProxy, objects);
60 reportFields::print(
" surfVector :", Info,
sVectors);
74 surfaceVectorField unitNorm(
mesh.Sf()/
mesh.magSf());
81 surfaceVectorField* tsvfPtr = (ssf * unitNorm).ptr();
82 tsvfPtr->
rename(ssf.name());
90 vtk::surfaceFieldWriter
writer
96 /
"surface-fields"/
"surfaceFields" + timeDesc
105 writer.writeTimeValue(timeValue);
119 if (Pstream::master())
123 fileName seriesName(vtk::seriesWriter::base(
outputName));
125 vtk::seriesWriter& series = vtkSeries(seriesName);
132 series.load(seriesName,
true, timeValue);
136 series.write(seriesName);
144 if (!selectedFaceZones.empty() && !
mesh.faceZones().empty())
148 sScalars = readFields<surfaceScalarField>(meshProxy, objects);
151 reportFields::print(
" surfScalar :", Info,
sScalars);
160 sVectors = readFields<surfaceVectorField>(meshProxy, objects);
163 reportFields::print(
" surfVector :", Info,
sVectors);
170 for (
const faceZone& fz :
mesh.faceZones())
172 if (!selectedFaceZones.match(fz.name()))
178 const primitiveFacePatch& pp = fz();
180 vtkWriterType_faceZone
writer
186 / (meshProxy.useSubMesh() ? meshProxy.name() : fz.name())
192 Info<<
" FaceZone : "
196 writer.beginFile(fz.name());
197 writer.writeTimeValue(timeValue);
215 if (Pstream::master())
219 fileName seriesName(vtk::seriesWriter::base(
outputName));
221 vtk::seriesWriter& series = vtkSeries(seriesName);
228 series.load(seriesName,
true, timeValue);
232 series.write(seriesName);
vtk::internalMeshWriter writer(topoMesh, topoCells, vtk::formatType::INLINE_ASCII, runTime.path()/"blockTopology")
Info<< nl<< "Wrote faMesh in vtk format: "<< writer.output().name()<< nl;}{ vtk::lineWriter writer(aMesh.points(), aMesh.edges(), fileName(aMesh.mesh().time().globalPath()/"finiteArea-edges"));writer.writeGeometry();writer.beginCellData(4);writer.writeProcIDs();{ Field< scalar > fld(faMeshTools::flattenEdgeField(aMesh.magLe(), true))
fileName relativePath(const fileName &input, const bool caseTag=false) const
virtual void rename(const word &newName)
Rename.
label nSurfaceVectorField
PtrList< const surfaceVectorField > sVectors
label nSurfaceScalarField
PtrList< const surfaceScalarField > sScalars
word outputName("finiteArea-edges.obj")
Foam::argList args(argc, argv)