FreeFOAM The Cross-Platform CFD Toolkit
reconstructLagrangianFields.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration |
5  \\ / A nd | Copyright (C) 1991-2010 OpenCFD Ltd.
6  \\/ M anipulation |
7 -------------------------------------------------------------------------------
8 License
9  This file is part of OpenFOAM.
10 
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
13  the Free Software Foundation, either version 3 of the License, or
14  (at your option) any later version.
15 
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
19  for more details.
20 
21  You should have received a copy of the GNU General Public License
22  along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
23 
24 \*---------------------------------------------------------------------------*/
25 
26 #include <OpenFOAM/IOField.H>
27 #include <OpenFOAM/Time.H>
28 
29 // * * * * * * * * * * * * * * * Global Functions * * * * * * * * * * * * * //
30 
31 template<class Type>
33 (
34  const word& cloudName,
35  const polyMesh& mesh,
36  const PtrList<fvMesh>& meshes,
37  const word& fieldName
38 )
39 {
40  // Construct empty field on mesh
41  tmp<IOField<Type> > tfield
42  (
43  new IOField<Type>
44  (
45  IOobject
46  (
47  fieldName,
48  mesh.time().timeName(),
49  cloud::prefix/cloudName,
50  mesh,
51  IOobject::NO_READ,
52  IOobject::NO_WRITE
53  ),
54  Field<Type>(0)
55  )
56  );
57  Field<Type>& field = tfield();
58 
59  forAll(meshes, i)
60  {
61  // Check object on local mesh
62  IOobject localIOobject
63  (
64  fieldName,
65  meshes[i].time().timeName(),
66  cloud::prefix/cloudName,
67  meshes[i],
68  IOobject::MUST_READ,
69  IOobject::NO_WRITE
70  );
71 
72  if (localIOobject.headerOk())
73  {
74  IOField<Type> fieldi(localIOobject);
75 
76  label offset = field.size();
77  field.setSize(offset + fieldi.size());
78 
79  forAll(fieldi, j)
80  {
81  field[offset + j] = fieldi[j];
82  }
83  }
84  }
85 
86  return tfield;
87 }
88 
89 
90 template<class Type>
92 (
93  const word& cloudName,
94  const polyMesh& mesh,
95  const PtrList<fvMesh>& meshes,
96  const IOobjectList& objects
97 )
98 {
99  word fieldClassName(IOField<Type>::typeName);
100 
101  IOobjectList fields = objects.lookupClass(fieldClassName);
102 
103  if (fields.size())
104  {
105  Info<< " Reconstructing lagrangian "
106  << fieldClassName << "s\n" << endl;
107 
108  forAllIter(IOobjectList, fields, fieldIter)
109  {
110  Info<< " " << fieldIter()->name() << endl;
111  reconstructLagrangianField<Type>
112  (
113  cloudName,
114  mesh,
115  meshes,
116  fieldIter()->name()
117  )().write();
118  }
119 
120  Info<< endl;
121  }
122 }
123 
124 
125 // ************************ vim: set sw=4 sts=4 et: ************************ //