FreeFOAM The Cross-Platform CFD Toolkit
MapVolFields.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 MapVolFields_H
27 #define MapVolFields_H
28 
30 #include <sampling/meshToMesh.H>
31 #include <OpenFOAM/IOobjectList.H>
32 
33 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
34 
35 namespace Foam
36 {
37 
38 template<class Type>
39 void MapVolFields
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  IOobject fieldTargetIOobject
63  (
64  fieldIter()->name(),
65  meshTarget.time().timeName(),
66  meshTarget,
69  );
70 
71  if (fieldTargetIOobject.headerOk())
72  {
73  Info<< " interpolating " << fieldIter()->name()
74  << endl;
75 
76  // Read field fieldSource
78  (
79  *fieldIter(),
80  meshSource
81  );
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  }
102 }
103 
104 
105 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
106 
107 } // End namespace Foam
108 
109 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
110 
111 #endif
112 
113 // ************************ vim: set sw=4 sts=4 et: ************************ //