src/postProcessing/foamCalcFunctions/field/componentsTurbo/componentsTurbo.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 Author
  25     Martin Beaudoin, Hydro-Quebec, 2012.  All rights reserved
  26
  27 \*---------------------------------------------------------------------------*/
  28
  29 #include "componentsTurbo.H"
  30 #include "addToRunTimeSelectionTable.H"
  31
  32 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  33
  34 namespace Foam
  35 {
  36     namespace calcTypes
  37     {
  38         defineTypeNameAndDebug(componentsTurbo, 0);
  39         addToRunTimeSelectionTable(calcType, componentsTurbo, dictionary);
  40     }
  41 }
  42
  43
  44 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  45
  46 Foam::calcTypes::componentsTurbo::componentsTurbo()
  47 :
  48     calcType()
  49 {}
  50
  51
  52 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
  53
  54 Foam::calcTypes::componentsTurbo::~componentsTurbo()
  55 {}
  56
  57
  58 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  59
  60 void Foam::calcTypes::componentsTurbo::init()
  61 {
  62     argList::validArgs.append("componentsTurbo");
  63     argList::validArgs.append("fieldName");
  64     //argList::validOptions.insert("rotationAxis", "rotation axis");
  65     //argList::validOptions.insert("rotationAxisOrigin", "rotation axis origin");
  66 }
  67
  68
  69 void Foam::calcTypes::componentsTurbo::preCalc
  70 (
  71     const argList& args,
  72     const Time& runTime,
  73     const fvMesh& mesh
  74 )
  75 {}
  76
  77
  78 void Foam::calcTypes::componentsTurbo::calc
  79 (
  80     const argList& args,
  81     const Time& runTime,
  82     const fvMesh& mesh
  83 )
  84 {
  85     const word& fieldName = args.additionalArgs()[1];
  86
  87     IOobject fieldHeader
  88     (
  89         fieldName,
  90         runTime.timeName(),
  91         mesh,
  92         IOobject::MUST_READ
  93     );
  94
  95     // Check field exists
  96     if (fieldHeader.headerOk())
  97     {
  98         bool processed = false;
  99
 100         writeComponentFields<vector>(fieldHeader, mesh, processed);
 101         //writeComponentFields<sphericalTensor>(fieldHeader, mesh, processed);
 102         //writeComponentFields<symmTensor>(fieldHeader, mesh, processed);
 103         //writeComponentFields<tensor>(fieldHeader, mesh, processed);
 104
 105         if (!processed)
 106         {
 107             FatalError
 108                 << "Unable to process " << fieldName << nl
 109                 << "No call to componentsTurbo for fields of type "
 110                 << fieldHeader.headerClassName() << nl << nl
 111                 << exit(FatalError);
 112         }
 113     }
 114     else
 115     {
 116         Info<< "    No " << fieldName << endl;
 117     }
 118 }
 119
 120
 121 // ************************************************************************* //