1 /*---------------------------------------------------------------------------*\
3 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
5 \\ / A nd | Copyright held by original author
7 -------------------------------------------------------------------------------
9 This file is part of OpenFOAM.
11 OpenFOAM is free software; you can redistribute it and/or modify it
12 under the terms of the GNU General Public License as published by the
13 Free Software Foundation; either version 2 of the License, or (at your
14 option) any later version.
16 OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
17 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
18 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
21 You should have received a copy of the GNU General Public License
22 along with OpenFOAM; if not, write to the Free Software Foundation,
23 Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
25 \*---------------------------------------------------------------------------*/
27 #include "lagrangianWriter.H"
28 #include "writeFuns.H"
30 #include "passiveParticle.H"
32 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
34 // Construct from components
35 Foam::lagrangianWriter::lagrangianWriter
39 const fileName& fName,
40 const word& cloudName,
47 cloudName_(cloudName),
50 const fvMesh& mesh = vMesh_.mesh();
53 writeFuns::writeHeader(os_, binary_, mesh.time().caseName());
54 os_ << "DATASET POLYDATA" << std::endl;
60 os_ << "POINTS " << nParcels_ << " float" << std::endl;
64 Cloud<passiveParticle> parcels(mesh, cloudName_, false);
66 nParcels_ = parcels.size();
68 os_ << "POINTS " << nParcels_ << " float" << std::endl;
70 DynamicList<floatScalar> partField(3*parcels.size());
72 forAllConstIter(Cloud<passiveParticle>, parcels, elmnt)
74 writeFuns::insert(elmnt().position(), partField);
76 writeFuns::write(os_, binary_, partField);
81 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
83 void Foam::lagrangianWriter::writeParcelHeader(const label nFields)
85 os_ << "POINT_DATA " << nParcels_ << std::endl
86 << "FIELD attributes " << nFields
91 // ************************************************************************* //