FreeFOAM The Cross-Platform CFD Toolkit
MapConsistentVolFields.H
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 #ifndef MapConsistentVolFields_H
27 #define MapConsistentVolFields_H
28 
30 #include <sampling/meshToMesh.H>
31 #include <OpenFOAM/IOobjectList.H>
32 
33 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
34 
35 namespace Foam
36 {
37 
38 template<class Type>
40 (
41  const IOobjectList& objects,
42  const meshToMesh& meshToMeshInterp
43 )
44 {
45  const fvMesh& meshSource = meshToMeshInterp.fromMesh();
46  const fvMesh& meshTarget = meshToMeshInterp.toMesh();
47 
48  word fieldClassName
49  (
51  );
52 
53  IOobjectList fields = objects.lookupClass(fieldClassName);
54 
55  for
56  (
57  IOobjectList::iterator fieldIter = fields.begin();
58  fieldIter != fields.end();
59  ++fieldIter
60  )
61  {
62  Info<< " interpolating " << fieldIter()->name()
63  << endl;
64 
65  // Read field
67  (
68  *fieldIter(),
69  meshSource
70  );
71 
72  IOobject fieldTargetIOobject
73  (
74  fieldIter()->name(),
75  meshTarget.time().timeName(),
76  meshTarget,
79  );
80 
81  if (fieldTargetIOobject.headerOk())
82  {
83  // Read fieldTarget
85  (
86  fieldTargetIOobject,
87  meshTarget
88  );
89 
90  // Interpolate field
91  meshToMeshInterp.interpolate
92  (
93  fieldTarget,
94  fieldSource,
96  );
97 
98  // Write field
99  fieldTarget.write();
100  }
101  else
102  {
103  fieldTargetIOobject.readOpt() = IOobject::NO_READ;
104 
105  // Interpolate field
107  (
108  fieldTargetIOobject,
109  meshToMeshInterp.interpolate
110  (
111  fieldSource,
113  )
114  );
115 
116  // Write field
117  fieldTarget.write();
118  }
119  }
120 }
121 
122 
123 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
124 
125 } // End namespace Foam
126 
127 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
128 
129 #endif
130 
131 // ************************ vim: set sw=4 sts=4 et: ************************ //