From d8d6030ab688b7d6abae158608ace0c35c3a5485 Mon Sep 17 00:00:00 2001 From: Andrew Heather Date: Fri, 7 Jul 2017 11:39:56 +0100 Subject: [PATCH] INT: Integration of Mattijs' collocated parallel IO additions Original commit message: ------------------------ Parallel IO: New collated file format When an OpenFOAM simulation runs in parallel, the data for decomposed fields and mesh(es) has historically been stored in multiple files within separate directories for each processor. Processor directories are named 'processorN', where N is the processor number. This commit introduces an alternative "collated" file format where the data for each decomposed field (and mesh) is collated into a single file, which is written and read on the master processor. The files are stored in a single directory named 'processors'. The new format produces significantly fewer files - one per field, instead of N per field. For large parallel cases, this avoids the restriction on the number of open files imposed by the operating system limits. The file writing can be threaded allowing the simulation to continue running while the data is being written to file. NFS (Network File System) is not needed when using the the collated format and additionally, there is an option to run without NFS with the original uncollated approach, known as "masterUncollated". The controls for the file handling are in the OptimisationSwitches of etc/controlDict: OptimisationSwitches { ... //- Parallel IO file handler // uncollated (default), collated or masterUncollated fileHandler uncollated; //- collated: thread buffer size for queued file writes. // If set to 0 or not sufficient for the file size threading is not used. // Default: 2e9 maxThreadFileBufferSize 2e9; //- masterUncollated: non-blocking buffer size. // If the file exceeds this buffer size scheduled transfer is used. // Default: 2e9 maxMasterFileBufferSize 2e9; } When using the collated file handling, memory is allocated for the data in the thread. maxThreadFileBufferSize sets the maximum size of memory in bytes that is allocated. If the data exceeds this size, the write does not use threading. When using the masterUncollated file handling, non-blocking MPI communication requires a sufficiently large memory buffer on the master node. maxMasterFileBufferSize sets the maximum size in bytes of the buffer. If the data exceeds this size, the system uses scheduled communication. The installation defaults for the fileHandler choice, maxThreadFileBufferSize and maxMasterFileBufferSize (set in etc/controlDict) can be over-ridden within the case controlDict file, like other parameters. Additionally the fileHandler can be set by: - the "-fileHandler" command line argument; - a FOAM_FILEHANDLER environment variable. A foamFormatConvert utility allows users to convert files between the collated and uncollated formats, e.g. mpirun -np 2 foamFormatConvert -parallel -fileHandler uncollated An example case demonstrating the file handling methods is provided in: $FOAM_TUTORIALS/IO/fileHandling The work was undertaken by Mattijs Janssens, in collaboration with Henry Weller. --- .../electromagnetics/mhdFoam/createPhiB.H | 2 +- .../heatTransfer/thermoFoam/setAlphaEff.H | 2 +- .../test/CompactIOList/Test-CompactIOList.C | 8 +- applications/test/IOField/Make/files | 3 + applications/test/IOField/Make/options | 2 + applications/test/IOField/Test-IOField.C | 189 ++ .../test/decomposedBlockData/Make/files | 3 + .../test/decomposedBlockData/Make/options | 0 .../Test-decomposedBlockData.C | 98 + .../Test-surfaceMeshConvert.C | 6 +- .../mesh/advanced/modifyMesh/modifyMesh.C | 3 +- .../mesh/advanced/selectCells/edgeStats.C | 15 +- .../foamToStarMesh/foamToStarMesh.C | 4 +- .../conversion/foamToSurface/foamToSurface.C | 2 +- .../DelaunayMesh/DelaunayMeshIO.C | 6 +- .../generation/snappyHexMesh/snappyHexMesh.C | 1 + .../mesh/manipulation/checkMesh/checkMesh.C | 3 +- .../mesh/manipulation/checkMesh/checkTools.C | 9 +- .../manipulation/createPatch/createPatch.C | 1 + .../mesh/manipulation/refineMesh/refineMesh.C | 3 +- .../mesh/manipulation/setSet/setSet.C | 8 +- .../mesh/manipulation/stitchMesh/stitchMesh.C | 5 +- .../mesh/manipulation/topoSet/topoSet.C | 9 +- .../foamFormatConvert/foamFormatConvert.C | 289 ++- .../foamFormatConvert/writeMeshObject.H | 14 +- .../decomposePar/decomposePar.C | 104 +- .../decomposePar/domainDecomposition.C | 3 - ...lagrangianFieldDecomposerDecomposeFields.C | 12 +- .../reconstructPar/reconstructPar.C | 117 +- .../reconstructParMesh/reconstructParMesh.C | 35 +- .../redistributePar/loadOrCreateMesh.C | 9 +- .../foamToEnsight/checkMeshMoving.H | 2 +- .../foamToEnsightParts/getTimeIndex.H | 2 +- .../utilities/postProcessing/noise/noise.C | 3 +- .../changeDictionary/changeDictionary.C | 455 +++-- .../mapFields/MapConsistentVolFields.H | 6 +- .../preProcessing/mapFields/MapVolFields.H | 15 +- .../surface/surfaceHookUp/surfaceHookUp.C | 3 +- .../surfaceMeshTriangulate.C | 1 + .../equilibriumCO/equilibriumCO.C | 1 + etc/controlDict | 14 + src/OSspecific/POSIX/POSIX.C | 300 ++- src/OSspecific/POSIX/regExp.H | 2 +- src/OpenFOAM/Make/files | 14 +- src/OpenFOAM/db/IOobject/IOobject.C | 187 +- src/OpenFOAM/db/IOobject/IOobject.H | 33 +- src/OpenFOAM/db/IOobject/IOobjectTemplates.C | 56 +- .../db/IOobject/IOobjectWriteHeader.C | 4 +- src/OpenFOAM/db/IOobjectList/IOobjectList.C | 27 +- .../IOobjects/CompactIOField/CompactIOField.C | 75 +- .../IOobjects/CompactIOField/CompactIOField.H | 10 +- .../IOobjects/CompactIOList/CompactIOList.C | 28 +- .../IOobjects/CompactIOList/CompactIOList.H | 5 +- .../IOobjects/GlobalIOField/GlobalIOField.H | 4 +- .../IOobjects/GlobalIOField/GlobalIOFields.C} | 41 +- .../db/IOobjects/GlobalIOList/GlobalIOList.H | 4 +- .../db/IOobjects/GlobalIOList/globalIOLists.C | 1 + .../db/IOobjects/GlobalIOList/globalIOLists.H | 9 +- src/OpenFOAM/db/IOobjects/IOField/IOField.C | 39 +- src/OpenFOAM/db/IOobjects/IOField/IOField.H | 5 +- src/OpenFOAM/db/IOobjects/IOMap/IOMap.C | 2 +- src/OpenFOAM/db/IOobjects/IOMap/IOMap.H | 4 +- src/OpenFOAM/db/IOobjects/IOMap/IOMapName.C | 10 +- .../db/IOobjects/IOPtrList/IOPtrList.C | 2 +- .../db/IOobjects/IOdictionary/IOdictionary.C | 2 +- .../db/IOobjects/IOdictionary/IOdictionary.H | 5 +- .../IOobjects/IOdictionary/baseIOdictionary.H | 8 +- .../IOdictionary/localIOdictionary.C | 17 +- .../IOdictionary/localIOdictionary.H | 10 +- .../IOdictionary/unwatchedIOdictionary.H | 4 +- .../decomposedBlockData/decomposedBlockData.C | 921 +++++++++ .../decomposedBlockData/decomposedBlockData.H | 197 ++ src/OpenFOAM/db/IOstreams/Fstreams/OFstream.C | 53 +- src/OpenFOAM/db/IOstreams/Fstreams/OFstream.H | 8 +- .../db/IOstreams/Fstreams/masterOFstream.C | 168 ++ .../db/IOstreams/Fstreams/masterOFstream.H | 100 + .../IOstreams/dummyISstream/dummyISstream.H | 160 ++ .../db/IOstreams/dummyIstream/dummyIstream.H | 161 ++ src/OpenFOAM/db/Time/Time.C | 133 +- src/OpenFOAM/db/Time/Time.H | 52 +- src/OpenFOAM/db/Time/TimeIO.C | 56 +- src/OpenFOAM/db/Time/findInstance.C | 149 +- src/OpenFOAM/db/Time/findTimes.C | 60 +- src/OpenFOAM/db/Time/instant/instant.C | 14 +- src/OpenFOAM/db/Time/instant/instant.H | 6 +- src/OpenFOAM/db/Time/timeSelector.C | 7 +- src/OpenFOAM/db/dictionary/dictionaryIO.C | 3 +- .../functionEntries/codeStream/codeStream.C | 2 +- .../functionEntries/codeStream/codeStream.H | 2 +- .../includeEntry/includeEntry.C | 15 +- .../includeEtcEntry/includeEtcEntry.C | 15 +- .../includeIfPresentEntry.C | 9 +- .../dlLibraryTable/dlLibraryTable.C | 8 +- .../functionObjectList/functionObjectList.C | 14 +- .../db/functionObjects/writeFile/writeFile.C | 5 +- .../db/objectRegistry/objectRegistry.C | 7 +- .../db/objectRegistry/objectRegistry.H | 303 +-- src/OpenFOAM/db/regIOobject/regIOobject.C | 82 +- src/OpenFOAM/db/regIOobject/regIOobject.H | 33 +- src/OpenFOAM/db/regIOobject/regIOobjectRead.C | 191 +- .../db/regIOobject/regIOobjectWrite.C | 67 +- .../GeometricField/GeometricField.C | 8 +- .../UniformDimensionedField.H | 4 +- ...eVaryingUniformFixedValuePointPatchField.C | 3 +- src/OpenFOAM/global/argList/argList.C | 62 +- .../collatedFileOperation/OFstreamCollator.C | 335 +++ .../collatedFileOperation/OFstreamCollator.H | 167 ++ .../collatedFileOperation.C | 415 ++++ .../collatedFileOperation.H | 138 ++ .../threadedCollatedOFstream.C | 65 + .../threadedCollatedOFstream.H | 92 + .../fileOperation/fileOperation.C | 666 ++++++ .../fileOperation/fileOperation.H | 431 ++++ .../masterUncollatedFileOperation.C | 1804 +++++++++++++++++ .../masterUncollatedFileOperation.H | 716 +++++++ .../masterUncollatedFileOperationTemplates.C | 148 ++ .../uncollatedFileOperation.C | 673 ++++++ .../uncollatedFileOperation.H | 283 +++ src/OpenFOAM/include/OSspecific.H | 31 +- .../interpolation2DTable.C | 2 +- .../interpolationTable/interpolationTable.C | 1 - .../tableReaders/csv/csvTableReader.C | 8 +- .../openFoam/openFoamTableReader.C | 8 +- .../uniformInterpolationTable.C | 3 +- .../polyBoundaryMesh/polyBoundaryMesh.C | 5 +- .../polyBoundaryMesh/polyBoundaryMesh.H | 5 +- src/OpenFOAM/meshes/polyMesh/polyMesh.C | 2 +- src/OpenFOAM/meshes/polyMesh/polyMeshIO.C | 2 +- .../preservePatchTypes/preservePatchTypes.C | 2 +- .../primitives/functions/Function1/CSV/CSV.C | 6 +- .../functions/Function1/TableFile/TableFile.C | 4 +- .../primitives/strings/fileName/fileName.C | 7 +- .../primitives/strings/fileName/fileName.H | 2 +- .../primitives/strings/string/string.H | 3 + .../primitives/strings/string/stringI.H | 6 + src/Pstream/mpi/UPstream.C | 19 +- src/conversion/common/reader/meshReader.C | 3 +- src/conversion/common/reader/meshReaderAux.C | 5 +- .../dynamicRefineFvMesh/dynamicRefineFvMesh.C | 7 +- .../dynamicRefineFvMesh/dynamicRefineFvMesh.H | 5 +- .../componentDisplacementMotionSolver.C | 2 +- .../points0/points0MotionSolver.C | 2 +- .../motionSolvers/motionSolver/motionSolver.C | 3 +- .../motionSolvers/motionSolver/motionSolver.H | 5 +- .../polyTopoChange/hexRef8/hexRef8.C | 10 +- .../polyTopoChange/hexRef8/hexRef8.H | 4 +- .../polyTopoChange/hexRef8/hexRef8Data.C | 2 +- src/engine/engineTime/engineTime.H | 4 +- .../cfdTools/general/MRF/IOMRFZoneList.C | 2 +- .../cfdTools/general/fvOptions/fvOptions.C | 2 +- src/finiteVolume/fvMesh/fvMesh.C | 17 +- src/finiteVolume/fvMesh/fvMesh.H | 7 +- .../fieldCoordinateSystemTransformTemplates.C | 8 +- .../field/readFields/readFieldsTemplates.C | 6 +- .../field/streamLine/streamLineParticle.C | 21 +- .../writeDictionary/writeDictionary.C | 5 +- .../displacementLaplacianFvMotionSolver.C | 2 +- .../Templates/DSMCParcel/DSMCParcelIO.C | 25 +- src/lagrangian/basic/Cloud/Cloud.H | 3 +- src/lagrangian/basic/Cloud/CloudIO.C | 41 +- src/lagrangian/basic/IOPosition/IOPosition.C | 26 +- src/lagrangian/basic/IOPosition/IOPosition.H | 6 +- .../basic/particle/particleTemplates.C | 46 +- .../Templates/KinematicCloud/KinematicCloud.C | 3 +- .../Templates/ReactingCloud/ReactingCloud.C | 2 +- .../ReactingMultiphaseCloud.C | 2 +- .../CollidingParcel/CollidingParcelIO.C | 63 +- .../KinematicParcel/KinematicParcelIO.C | 90 +- .../Templates/MPPICParcel/MPPICParcelIO.C | 15 +- .../ReactingMultiphaseParcelIO.C | 28 +- .../ReactingParcel/ReactingParcelIO.C | 26 +- .../Templates/ThermoParcel/ThermoParcelIO.C | 15 +- .../CloudFunctionObject/CloudFunctionObject.C | 1 + .../AveragingMethod/AveragingMethod.C | 10 +- .../AveragingMethod/AveragingMethod.H | 2 +- .../molecule/molecule/moleculeIO.C | 56 +- .../pairPotentialList/pairPotentialList.C | 28 +- .../solidParticle/solidParticleIO.C | 15 +- .../Templates/SprayParcel/SprayParcelIO.C | 102 +- .../extendedFeatureEdgeMesh.H | 3 +- .../featureEdgeMesh/featureEdgeMesh.H | 2 +- .../triSurfaceMesh/triSurfaceMesh.C | 14 +- .../triSurfaceMesh/triSurfaceMesh.H | 5 +- .../cellSources/fieldToCell/fieldToCell.C | 2 +- src/meshTools/sets/topoSets/cellZoneSet.C | 9 +- src/meshTools/sets/topoSets/cellZoneSet.H | 5 +- src/meshTools/sets/topoSets/faceZoneSet.C | 9 +- src/meshTools/sets/topoSets/faceZoneSet.H | 5 +- src/meshTools/sets/topoSets/pointZoneSet.C | 9 +- src/meshTools/sets/topoSets/pointZoneSet.H | 5 +- .../distributedTriSurfaceMesh.C | 5 +- .../distributedTriSurfaceMesh.H | 25 +- .../reconstruct/reconstructLagrangianFields.C | 2 +- .../regionModel/regionModel/regionModel.C | 3 +- .../rigidBodyMeshMotion/rigidBodyMeshMotion.C | 3 +- .../rigidBodyMeshMotion/rigidBodyMeshMotion.H | 3 +- .../rigidBodyMeshMotionSolver.C | 7 +- .../rigidBodyMeshMotionSolver.H | 5 +- src/sampling/probes/probes.C | 4 +- .../sampledSet/sampledSets/sampledSets.C | 6 +- .../isoSurface/sampledIsoSurface.C | 2 +- .../sampledSurfaces/sampledSurfaces.C | 6 +- .../sixDoFRigidBodyMotionSolver.C | 5 +- .../sixDoFRigidBodyMotionSolver.H | 5 +- .../radiationModel/radiationModelNew.C | 2 +- .../basicMultiComponentMixture.C | 63 +- .../pyrolysisChemistryModel.C | 2 +- tutorials/IO/fileHandler/0/U | 38 + tutorials/IO/fileHandler/Allrun | 28 + tutorials/IO/fileHandler/constant/g | 22 + .../constant/kinematicCloudPositions | 92 + .../constant/kinematicCloudProperties | 138 ++ .../fileHandler/constant/transportProperties | 24 + .../fileHandler/constant/turbulenceProperties | 21 + tutorials/IO/fileHandler/system/blockMeshDict | 72 + tutorials/IO/fileHandler/system/controlDict | 53 + .../IO/fileHandler/system/decomposeParDict | 45 + tutorials/IO/fileHandler/system/fvSchemes | 43 + tutorials/IO/fileHandler/system/fvSolution | 22 + .../BernardCells/system/blockMeshDict | 5 +- .../iglooWithFridges/Allrun | 16 + 221 files changed, 11118 insertions(+), 1946 deletions(-) create mode 100644 applications/test/IOField/Make/files create mode 100644 applications/test/IOField/Make/options create mode 100644 applications/test/IOField/Test-IOField.C create mode 100644 applications/test/decomposedBlockData/Make/files create mode 100644 applications/test/decomposedBlockData/Make/options create mode 100644 applications/test/decomposedBlockData/Test-decomposedBlockData.C rename src/OpenFOAM/{fields/Fields/globalFields/globalIOFields.C => db/IOobjects/GlobalIOField/GlobalIOFields.C} (64%) create mode 100644 src/OpenFOAM/db/IOobjects/decomposedBlockData/decomposedBlockData.C create mode 100644 src/OpenFOAM/db/IOobjects/decomposedBlockData/decomposedBlockData.H create mode 100644 src/OpenFOAM/db/IOstreams/Fstreams/masterOFstream.C create mode 100644 src/OpenFOAM/db/IOstreams/Fstreams/masterOFstream.H create mode 100644 src/OpenFOAM/db/IOstreams/dummyISstream/dummyISstream.H create mode 100644 src/OpenFOAM/db/IOstreams/dummyIstream/dummyIstream.H create mode 100644 src/OpenFOAM/global/fileOperations/collatedFileOperation/OFstreamCollator.C create mode 100644 src/OpenFOAM/global/fileOperations/collatedFileOperation/OFstreamCollator.H create mode 100644 src/OpenFOAM/global/fileOperations/collatedFileOperation/collatedFileOperation.C create mode 100644 src/OpenFOAM/global/fileOperations/collatedFileOperation/collatedFileOperation.H create mode 100644 src/OpenFOAM/global/fileOperations/collatedFileOperation/threadedCollatedOFstream.C create mode 100644 src/OpenFOAM/global/fileOperations/collatedFileOperation/threadedCollatedOFstream.H create mode 100644 src/OpenFOAM/global/fileOperations/fileOperation/fileOperation.C create mode 100644 src/OpenFOAM/global/fileOperations/fileOperation/fileOperation.H create mode 100644 src/OpenFOAM/global/fileOperations/masterUncollatedFileOperation/masterUncollatedFileOperation.C create mode 100644 src/OpenFOAM/global/fileOperations/masterUncollatedFileOperation/masterUncollatedFileOperation.H create mode 100644 src/OpenFOAM/global/fileOperations/masterUncollatedFileOperation/masterUncollatedFileOperationTemplates.C create mode 100644 src/OpenFOAM/global/fileOperations/uncollatedFileOperation/uncollatedFileOperation.C create mode 100644 src/OpenFOAM/global/fileOperations/uncollatedFileOperation/uncollatedFileOperation.H create mode 100644 tutorials/IO/fileHandler/0/U create mode 100755 tutorials/IO/fileHandler/Allrun create mode 100644 tutorials/IO/fileHandler/constant/g create mode 100644 tutorials/IO/fileHandler/constant/kinematicCloudPositions create mode 100644 tutorials/IO/fileHandler/constant/kinematicCloudProperties create mode 100644 tutorials/IO/fileHandler/constant/transportProperties create mode 100644 tutorials/IO/fileHandler/constant/turbulenceProperties create mode 100644 tutorials/IO/fileHandler/system/blockMeshDict create mode 100644 tutorials/IO/fileHandler/system/controlDict create mode 100644 tutorials/IO/fileHandler/system/decomposeParDict create mode 100644 tutorials/IO/fileHandler/system/fvSchemes create mode 100644 tutorials/IO/fileHandler/system/fvSolution diff --git a/applications/solvers/electromagnetics/mhdFoam/createPhiB.H b/applications/solvers/electromagnetics/mhdFoam/createPhiB.H index 21076a51a5..02b48f8d28 100644 --- a/applications/solvers/electromagnetics/mhdFoam/createPhiB.H +++ b/applications/solvers/electromagnetics/mhdFoam/createPhiB.H @@ -9,7 +9,7 @@ IOobject phiBHeader surfaceScalarField* phiBPtr = nullptr; - if (phiBHeader.typeHeaderOk(true)) +if (phiBHeader.typeHeaderOk(true)) { Info<< "Reading face flux "; diff --git a/applications/solvers/heatTransfer/thermoFoam/setAlphaEff.H b/applications/solvers/heatTransfer/thermoFoam/setAlphaEff.H index c760349a6e..e7da9c6af6 100644 --- a/applications/solvers/heatTransfer/thermoFoam/setAlphaEff.H +++ b/applications/solvers/heatTransfer/thermoFoam/setAlphaEff.H @@ -11,7 +11,7 @@ IOobject turbulencePropertiesHeader false ); -if (turbulencePropertiesHeader.typeHeaderOk(false)) +if (turbulencePropertiesHeader.typeHeaderOk(true)) { autoPtr turbulence ( diff --git a/applications/test/CompactIOList/Test-CompactIOList.C b/applications/test/CompactIOList/Test-CompactIOList.C index 12d64c6bde..2024e663e0 100644 --- a/applications/test/CompactIOList/Test-CompactIOList.C +++ b/applications/test/CompactIOList/Test-CompactIOList.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2011-2015 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -87,7 +87,8 @@ int main(int argc, char *argv[]) ( format, IOstream::currentVersion, - IOstream::UNCOMPRESSED + IOstream::UNCOMPRESSED, + true ); Info<< "Written old format faceList in = " @@ -149,7 +150,8 @@ int main(int argc, char *argv[]) ( format, IOstream::currentVersion, - IOstream::UNCOMPRESSED + IOstream::UNCOMPRESSED, + true ); Info<< "Written new format faceList in = " diff --git a/applications/test/IOField/Make/files b/applications/test/IOField/Make/files new file mode 100644 index 0000000000..74e614deae --- /dev/null +++ b/applications/test/IOField/Make/files @@ -0,0 +1,3 @@ +Test-IOField.C + +EXE = $(FOAM_USER_APPBIN)/Test-IOField diff --git a/applications/test/IOField/Make/options b/applications/test/IOField/Make/options new file mode 100644 index 0000000000..6a9e9810b3 --- /dev/null +++ b/applications/test/IOField/Make/options @@ -0,0 +1,2 @@ +/* EXE_INC = -I$(LIB_SRC)/cfdTools/include */ +/* EXE_LIBS = -lfiniteVolume */ diff --git a/applications/test/IOField/Test-IOField.C b/applications/test/IOField/Test-IOField.C new file mode 100644 index 0000000000..d49acaf8f5 --- /dev/null +++ b/applications/test/IOField/Test-IOField.C @@ -0,0 +1,189 @@ +/*---------------------------------------------------------------------------*\ + ========= | + \\ / F ield | OpenFOAM: The Open Source CFD Toolbox + \\ / O peration | + \\ / A nd | Copyright (C) 2017 OpenFOAM Foundation + \\/ M anipulation | +------------------------------------------------------------------------------- +License + This file is part of OpenFOAM. + + OpenFOAM is free software: you can redistribute it and/or modify it + under the terms of the GNU General Public License as published by + the Free Software Foundation, either version 3 of the License, or + (at your option) any later version. + + OpenFOAM is distributed in the hope that it will be useful, but WITHOUT + ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or + FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License + for more details. + + You should have received a copy of the GNU General Public License + along with OpenFOAM. If not, see . + +Application + Test-IOField + +Description + Test the processor-local reading of IOField (used in the lagrangian libs) + +\*---------------------------------------------------------------------------*/ + +#include "IOField.H" +#include "argList.H" +#include "polyMesh.H" +#include "Time.H" + +using namespace Foam; + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +void write(const IOobject& io, const label sz) +{ + IOField