Home
Downloads
Documentation
Installation
User Guide
man-pages
API Documentation
README
Release Notes
Changes
License
Support
SourceForge Project
Main Page
Related Pages
Namespaces
Classes
Files
File List
File Members
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Friends
Macros
Pages
applications
utilities
preProcessing
mapFields
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
29
#include <
OpenFOAM/GeometricField.H
>
30
#include <
sampling/meshToMesh.H
>
31
#include <
OpenFOAM/IOobjectList.H
>
32
33
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
34
35
namespace
Foam
36
{
37
38
template
<
class
Type>
39
void
MapConsistentVolFields
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
(
50
GeometricField<Type, fvPatchField, volMesh>::typeName
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
66
GeometricField<Type, fvPatchField, volMesh>
fieldSource
67
(
68
*fieldIter(),
69
meshSource
70
);
71
72
IOobject
fieldTargetIOobject
73
(
74
fieldIter()->
name
(),
75
meshTarget.
time
().
timeName
(),
76
meshTarget,
77
IOobject::MUST_READ
,
78
IOobject::AUTO_WRITE
79
);
80
81
if
(fieldTargetIOobject.headerOk())
82
{
83
// Read fieldTarget
84
GeometricField<Type, fvPatchField, volMesh>
fieldTarget
85
(
86
fieldTargetIOobject,
87
meshTarget
88
);
89
90
// Interpolate field
91
meshToMeshInterp.
interpolate
92
(
93
fieldTarget,
94
fieldSource,
95
meshToMesh::INTERPOLATE
96
);
97
98
// Write field
99
fieldTarget.
write
();
100
}
101
else
102
{
103
fieldTargetIOobject.readOpt() =
IOobject::NO_READ
;
104
105
// Interpolate field
106
GeometricField<Type, fvPatchField, volMesh>
fieldTarget
107
(
108
fieldTargetIOobject,
109
meshToMeshInterp.
interpolate
110
(
111
fieldSource,
112
meshToMesh::INTERPOLATE
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: ************************ //