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) 1991-2008 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 the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     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, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 \*---------------------------------------------------------------------------*/
  26
  27 #ifndef MapVolFields_H
  28 #define MapVolFields_H
  29
  30 #include "GeometricField.H"
  31 #include "meshToMesh.H"
  32 #include "cuttingPlane.H"
  33 #include "IOobjectList.H"
  34
  35 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  36
  37 namespace Foam
  38 {
  39
  40 template<class Type>
  41 void MapVolFields
  42 (
  43     const IOobjectList& objects,
  44     const meshToMesh& meshToMeshInterp
  45 )
  46 {
  47     const fvMesh& meshSource = meshToMeshInterp.fromMesh();
  48     const fvMesh& meshTarget = meshToMeshInterp.toMesh();
  49
  50     word fieldClassName
  51     (
  52         GeometricField<Type, fvPatchField, volMesh>::typeName
  53     );
  54
  55     IOobjectList fields = objects.lookupClass(fieldClassName);
  56
  57     for
  58     (
  59         IOobjectList::iterator fieldIter = fields.begin();
  60         fieldIter != fields.end();
  61         ++fieldIter
  62     )
  63     {
  64         IOobject fieldTargetIOobject
  65         (
  66             fieldIter()->name(),
  67             meshTarget.time().timeName(),
  68             meshTarget,
  69             IOobject::MUST_READ,
  70             IOobject::AUTO_WRITE
  71         );
  72
  73         if (fieldTargetIOobject.headerOk())
  74         {
  75             Info<< "    interpolating " << fieldIter()->name()
  76                 << endl;
  77
  78             // Read field fieldSource
  79             GeometricField<Type, fvPatchField, volMesh> fieldSource
  80             (
  81                 *fieldIter(),
  82                 meshSource
  83             );
  84
  85             // Read fieldTarget
  86             GeometricField<Type, fvPatchField, volMesh> fieldTarget
  87             (
  88                 fieldTargetIOobject,
  89                 meshTarget
  90             );
  91
  92             // Interpolate field
  93             meshToMeshInterp.interpolate
  94             (
  95                 fieldTarget,
  96                 fieldSource,
  97                 meshToMesh::INTERPOLATE
  98             );
  99
 100             // Write field
 101             fieldTarget.write();
 102         }
 103     }
 104 }
 105
 106
 107 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 108
 109 } // End namespace Foam
 110
 111 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 112
 113 #endif
 114
 115 // ************************************************************************* //