applications/utilities/preProcessing/mapFields/MapLagrangianFields.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | foam-extend: Open Source CFD
   4    \\    /   O peration     | Version:     3.2
   5     \\  /    A nd           | Web:         http://www.foam-extend.org
   6      \\/     M anipulation  | For copyright notice see file Copyright
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of foam-extend.
  10
  11     foam-extend 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 3 of the License, or (at your
  14     option) any later version.
  15
  16     foam-extend is distributed in the hope that it will be useful, but
  17     WITHOUT ANY WARRANTY; without even the implied warranty of
  18     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  19     General Public License for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 InNamespace
  25     Foam
  26
  27 Description
  28     Gets the indices of (source)particles that have been appended to the
  29     target cloud and maps the lagrangian fields accordingly.
  30
  31 \*---------------------------------------------------------------------------*/
  32
  33 #ifndef MapLagrangianFields_H
  34 #define MapLagrangianFields_H
  35
  36 #include "cloud.H"
  37 #include "GeometricField.H"
  38 #include "meshToMesh.H"
  39 #include "IOobjectList.H"
  40
  41 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  42
  43 namespace Foam
  44 {
  45
  46 //- Gets the indices of (source)particles that have been appended to the
  47 //  target cloud and maps the lagrangian fields accordingly.
  48 template<class Type>
  49 void MapLagrangianFields
  50 (
  51     const string& cloudName,
  52     const IOobjectList& objects,
  53     const meshToMesh& meshToMeshInterp,
  54     const labelList& addParticles
  55 )
  56 {
  57     const fvMesh& meshTarget = meshToMeshInterp.toMesh();
  58
  59     IOobjectList fields = objects.lookupClass(IOField<Type>::typeName);
  60
  61     forAllIter(IOobjectList, fields, fieldIter)
  62     {
  63         Info<< "    mapping lagrangian field "
  64             << fieldIter()->name() << endl;
  65
  66         // Read field (does not need mesh)
  67         IOField<Type> fieldSource(*fieldIter());
  68
  69         // Map
  70         IOField<Type> fieldTarget
  71         (
  72             IOobject
  73             (
  74                 fieldIter()->name(),
  75                 meshTarget.time().timeName(),
  76                 cloud::prefix/cloudName,
  77                 meshTarget,
  78                 IOobject::NO_READ,
  79                 IOobject::NO_WRITE,
  80                 false
  81             ),
  82             addParticles.size()
  83         );
  84         forAll(addParticles, i)
  85         {
  86             fieldTarget[i] = fieldSource[addParticles[i]];
  87         }
  88
  89         // Write field
  90         fieldTarget.write();
  91     }
  92 }
  93
  94
  95 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  96
  97 } // End namespace Foam
  98
  99 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 100
 101 #endif
 102
 103 // ************************************************************************* //