applications/utilities/postProcessing/sampling/probeLocations/probeLocations.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2010 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
  13     the Free Software Foundation, either version 3 of the License, or
  14     (at your 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, see <http://www.gnu.org/licenses/>.
  23
  24 Description
  25     Probe locations.
  26
  27 \*---------------------------------------------------------------------------*/
  28
  29 #include "argList.H"
  30 #include "timeSelector.H"
  31 #include "IOprobes.H"
  32
  33 using namespace Foam;
  34
  35 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  36 // Main program:
  37
  38 int main(int argc, char *argv[])
  39 {
  40     timeSelector::addOptions();
  41 #   include "addRegionOption.H"
  42
  43 #   include "setRootCase.H"
  44 #   include "createTime.H"
  45     instantList timeDirs = timeSelector::select0(runTime, args);
  46 #   include "createNamedMesh.H"
  47
  48     IOprobes sniff
  49     (
  50         probes::typeName,
  51         mesh,
  52         "probesDict",
  53         IOobject::MUST_READ,
  54         true
  55     );
  56
  57     forAll(timeDirs, timeI)
  58     {
  59         runTime.setTime(timeDirs[timeI], timeI);
  60         Info<< "Time = " << runTime.timeName() << endl;
  61
  62         // Handle geometry/topology changes
  63         polyMesh::readUpdateState state = mesh.readUpdate();
  64
  65         if
  66         (
  67             state == polyMesh::POINTS_MOVED
  68          || state == polyMesh::TOPO_CHANGE
  69         )
  70         {
  71             sniff.read();
  72         }
  73
  74         sniff.write();
  75
  76         Info<< endl;
  77     }
  78
  79     Info<< "End\n" << endl;
  80
  81     return 0;
  82 }
  83
  84
  85 // ************************************************************************* //