src/postProcessing/functionObjects/field/maxFieldCell/maxFieldCell.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     Hrvoje Jasak, Wikki Ltd.  All rights reserved
  26
  27 \*---------------------------------------------------------------------------*/
  28
  29 #include "maxFieldCell.H"
  30 #include "addToRunTimeSelectionTable.H"
  31 #include "volFields.H"
  32
  33 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  34
  35 namespace Foam
  36 {
  37     defineTypeNameAndDebug(maxFieldCell, 0);
  38
  39     addToRunTimeSelectionTable
  40     (
  41         functionObject,
  42         maxFieldCell,
  43         dictionary
  44     );
  45 }
  46
  47
  48 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  49
  50
  51 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  52
  53 Foam::maxFieldCell::maxFieldCell
  54 (
  55     const word& name,
  56     const Time& t,
  57     const dictionary& dict
  58 )
  59 :
  60     functionObject(name),
  61     time_(t),
  62     regionName_(polyMesh::defaultRegion),
  63     fieldName_(dict.lookup("name"))
  64 {
  65     if (dict.found("region"))
  66     {
  67         dict.lookup("region") >> regionName_;
  68     }
  69
  70     Info<< "Creating maxFieldCell for field "
  71         << fieldName_ << endl;
  72 }
  73
  74
  75 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  76
  77 bool Foam::maxFieldCell::start()
  78 {
  79     return true;
  80 }
  81
  82
  83 bool Foam::maxFieldCell::execute()
  84 {
  85     const fvMesh& mesh =
  86         time_.lookupObject<fvMesh>(regionName_);
  87
  88     if (mesh.foundObject<volScalarField>(fieldName_))
  89     {
  90         const volScalarField& f = mesh.lookupObject<volScalarField>
  91         (
  92             fieldName_
  93         );
  94         const scalarField& fIn = f.internalField();
  95
  96         if (!fIn.empty())
  97         {
  98             label maxCellID = 0;
  99             scalar maxValue = mag(fIn[maxCellID]);
 100
 101             forAll (fIn, cellI)
 102             {
 103                 if (mag(fIn[cellI]) > maxValue)
 104                 {
 105                     maxCellID = cellI;
 106                     maxValue = mag(fIn[maxCellID]);
 107                 }
 108             }
 109
 110             Info<< "Field " << fieldName_ << " max magnitude = " << maxValue
 111                 << " in cell " << maxCellID << endl;
 112
 113             return true;
 114         }
 115         else
 116         {
 117             Info<< "Field "  << fieldName_ << " empty.  Skipping." << endl;
 118
 119             return false;
 120         }
 121     }
 122     else if (mesh.foundObject<volVectorField>(fieldName_))
 123     {
 124         const volVectorField& f =
 125             mesh.lookupObject<volVectorField>(fieldName_);
 126
 127         const vectorField& fIn = f.internalField();
 128
 129         if (!fIn.empty())
 130         {
 131             label maxCellID = 0;
 132             scalar maxValue = mag(fIn[maxCellID]);
 133
 134             forAll (fIn, cellI)
 135             {
 136                 if (mag(fIn[cellI]) > maxValue)
 137                 {
 138                     maxCellID = cellI;
 139                     maxValue = mag(fIn[maxCellID]);
 140                 }
 141             }
 142
 143             Info<< "Field " << fieldName_ << " max magnitude = " << maxValue
 144                 << " in cell " << maxCellID << endl;
 145
 146             return true;
 147         }
 148         else
 149         {
 150             Info<< "Field "  << fieldName_ << " empty.  Skipping." << endl;
 151
 152             return false;
 153         }
 154     }
 155     else
 156     {
 157         Info<< "Field "  << fieldName_ << " not found.  Skipping." << endl;
 158
 159         return false;
 160     }
 161 }
 162
 163
 164 bool Foam::maxFieldCell::read(const dictionary& dict)
 165 {
 166     fieldName_ = word(dict.lookup("name"));
 167
 168     if (dict.found("region"))
 169     {
 170         dict.lookup("region") >> regionName_;
 171     }
 172
 173     return false;
 174 }
 175
 176 // ************************************************************************* //