src/postProcessing/foamCalcFunctions/basic/addSubtract/writeAddSubtractField.C

   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 \*---------------------------------------------------------------------------*/
  25
  26 template<class Type>
  27 void Foam::calcTypes::addSubtract::writeAddSubtractField
  28 (
  29     const IOobject& baseHeader,
  30     const IOobject& addHeader,
  31     const fvMesh& mesh,
  32     bool& processed
  33 )
  34 {
  35     typedef GeometricField<Type, fvPatchField, volMesh> fieldType;
  36
  37     if
  38     (
  39         baseHeader.headerClassName() == fieldType::typeName
  40      && baseHeader.headerClassName() == addHeader.headerClassName()
  41     )
  42     {
  43         if (resultName_ == "")
  44         {
  45             if (calcMode_ == ADD)
  46             {
  47                 resultName_ = baseHeader.name() + "_add_" + addHeader.name();
  48             }
  49             else
  50             {
  51                 resultName_ = baseHeader.name() + "_subtract_"
  52                     + addHeader.name();
  53             }
  54         }
  55
  56         Info<< "    Reading " << baseHeader.name() << endl;
  57         fieldType baseField(baseHeader, mesh);
  58
  59         Info<< "    Reading " << addHeader.name() << endl;
  60         fieldType addField(addHeader, mesh);
  61
  62         if (baseField.dimensions() == addField.dimensions())
  63         {
  64             Info<< "    Calculating " << resultName_ << endl;
  65
  66             fieldType newField
  67             (
  68                 IOobject
  69                 (
  70                     resultName_,
  71                     mesh.time().timeName(),
  72                     mesh,
  73                     IOobject::NO_READ
  74                 ),
  75                 calcMode_ == ADD ? baseField + addField : baseField - addField
  76             );
  77             newField.write();
  78         }
  79         else
  80         {
  81             Info<< "    Cannot calculate " << resultName_ << nl
  82                 << "    - inconsistent dimensions: "
  83                 << baseField.dimensions() << " - " << addField.dimensions()
  84                 << endl;
  85         }
  86
  87         processed = true;
  88     }
  89 }
  90
  91
  92 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //