applications/utilities/preProcessing/mapFields/MapVolFields.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2004-2011 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
  29 #include "GeometricField.H"
  30 #include "meshToMesh.H"
  31 #include "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     const meshToMesh::order& mapOrder
  44 )
  45 {
  46     const fvMesh& meshSource = meshToMeshInterp.fromMesh();
  47     const fvMesh& meshTarget = meshToMeshInterp.toMesh();
  48
  49     word fieldClassName
  50     (
  51         GeometricField<Type, fvPatchField, volMesh>::typeName
  52     );
  53
  54     IOobjectList fields = objects.lookupClass(fieldClassName);
  55
  56     forAllIter(IOobjectList, fields, fieldIter)
  57     {
  58         IOobject fieldTargetIOobject
  59         (
  60             fieldIter()->name(),
  61             meshTarget.time().timeName(),
  62             meshTarget,
  63             IOobject::MUST_READ,
  64             IOobject::AUTO_WRITE
  65         );
  66
  67         if (fieldTargetIOobject.headerOk())
  68         {
  69             Info<< "    interpolating " << fieldIter()->name()
  70                 << endl;
  71
  72             // Read field fieldSource
  73             GeometricField<Type, fvPatchField, volMesh> fieldSource
  74             (
  75                 *fieldIter(),
  76                 meshSource
  77             );
  78
  79             // Read fieldTarget
  80             GeometricField<Type, fvPatchField, volMesh> fieldTarget
  81             (
  82                 fieldTargetIOobject,
  83                 meshTarget
  84             );
  85
  86             // Interpolate field
  87             meshToMeshInterp.interpolate(fieldTarget, fieldSource, mapOrder);
  88
  89             // Write field
  90             fieldTarget.write();
  91         }
  92     }
  93 }
  94
  95
  96 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  97
  98 } // End namespace Foam
  99
 100 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 101
 102 #endif
 103
 104 // ************************************************************************* //