40 using reportFields = foamToVtkReportFields;
42 autoPtr<faMesh> faMeshPtr;
44 const label nAreaFields =
45 objects.count(stringListOps::foundOp<word>(fieldTypes::area));
47 if (nAreaFields || withMeshIds)
49 faMeshPtr = faMesh::TryNew(meshProxy.baseMesh());
52 if (faMeshPtr && (nAreaFields || withMeshIds))
54 const faMesh& areaMesh = faMeshPtr();
56 reportFields::area(Info, objects);
58 const auto& pp = faMeshPtr->patch();
60 vtkWriterType_areaMesh
writer
66 /
"finiteArea" + timeDesc
73 writer.beginFile(areaMesh.name());
75 writer.writeTimeValue(timeValue);
81 (withMeshIds ? 2 + (
writer.parallel() ? 1 : 0) : 0)
87 const globalIndex
procAddr(areaMesh.nFaces());
94 writer.writeCellData(
"faceLabels", areaMesh.faceLabels());
111 if (Pstream::master())
115 fileName seriesName(vtk::seriesWriter::base(
outputName));
117 vtk::seriesWriter& series = vtkSeries(seriesName);
124 series.load(seriesName,
true, timeValue);
128 series.write(seriesName);
vtk::internalMeshWriter writer(topoMesh, topoCells, vtk::formatType::INLINE_ASCII, runTime.path()/"blockTopology")
globalIndex procAddr(aMesh.nFaces())
fileName relativePath(const fileName &input, const bool caseTag=false) const
word outputName("finiteArea-edges.obj")
Foam::argList args(argc, argv)