2defba00a9
old "positions" file form
The change to barycentric-based tracking changed the contents of the
cloud "positions" file to a new format comprising the barycentric
co-ordinates and other cell position-based info. This broke
backwards compatibility, providing no option to restart old cases
(v1706 and earlier), and caused difficulties for dependent code, e.g.
for post-processing utilities that could only infer the contents only
after reading.
The barycentric position info is now written to a file called
"coordinates" with provision to restart old cases for which only the
"positions" file is available. Related utilities, e.g. for parallel
running and data conversion have been updated to be able to support both
file types.
To write the "positions" file by default, use set the following option
in the InfoSwitches section of the controlDict:
writeLagrangianPositions 1;
96 lines
3.0 KiB
C
96 lines
3.0 KiB
C
/*---------------------------------------------------------------------------*\
|
|
========= |
|
|
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
|
\\ / O peration |
|
|
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
|
|
\\/ M anipulation | Copyright (C) 2016 OpenCFD Ltd.
|
|
-------------------------------------------------------------------------------
|
|
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 <http://www.gnu.org/licenses/>.
|
|
|
|
\*---------------------------------------------------------------------------*/
|
|
|
|
#include "ensightOutputSerialCloud.H"
|
|
#include "ensightPTraits.H"
|
|
|
|
#include "passiveParticle.H"
|
|
#include "IOField.H"
|
|
#include "volFields.H"
|
|
#include "surfaceFields.H"
|
|
|
|
// * * * * * * * * * * * * * * * Global Functions * * * * * * * * * * * * * //
|
|
|
|
void Foam::ensightSerialCloud::writePositions
|
|
(
|
|
const polyMesh& mesh,
|
|
const word& cloudName,
|
|
autoPtr<ensightFile> output
|
|
)
|
|
{
|
|
label nTotParcels = 0;
|
|
autoPtr<Cloud<passiveParticle>> cloudPtr;
|
|
|
|
cloudPtr.reset(new Cloud<passiveParticle>(mesh, cloudName, false));
|
|
nTotParcels = cloudPtr().size();
|
|
|
|
Cloud<passiveParticle> parcels(mesh, cloudName, false);
|
|
|
|
if (Pstream::master())
|
|
{
|
|
ensightFile& os = output();
|
|
os.beginParticleCoordinates(nTotParcels);
|
|
|
|
// binary write is Ensight6 - first ids, then positions
|
|
if (os.format() == IOstream::BINARY)
|
|
{
|
|
// 1-index
|
|
for (label parcelId = 0; parcelId < nTotParcels; ++parcelId)
|
|
{
|
|
os.write(parcelId+1);
|
|
}
|
|
|
|
|
|
forAllConstIter(Cloud<passiveParticle>, cloudPtr(), elmnt)
|
|
{
|
|
const vector p(elmnt().position());
|
|
|
|
os.write(p.x());
|
|
os.write(p.y());
|
|
os.write(p.z());
|
|
}
|
|
}
|
|
else
|
|
{
|
|
// ASCII id + position together
|
|
|
|
label parcelId = 0;
|
|
forAllConstIter(Cloud<passiveParticle>, cloudPtr(), elmnt)
|
|
{
|
|
const vector p(elmnt().position());
|
|
|
|
os.write(++parcelId, 8); // unusual width
|
|
os.write(p.x());
|
|
os.write(p.y());
|
|
os.write(p.z());
|
|
os.newline();
|
|
}
|
|
}
|
|
}
|
|
}
|
|
|
|
|
|
// ************************************************************************* //
|