applications/utilities/mesh/conversion/foamToStarMesh/foamToStarMesh.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
   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 Application
  25     foamToStarMesh
  26
  27 Description
  28     Reads an OpenFOAM mesh and writes a pro-STAR (v4) bnd/cel/vrt format.
  29
  30 Usage
  31     - foamToStarMesh [OPTION] \n
  32     Reads an OpenFOAM mesh and writes a pro-STAR (v4) bnd/cel/vrt format.
  33
  34     \param -noBnd \n
  35     Suppress writing the \c .bnd file
  36
  37     \param -scale \<factor\>\n
  38     Specify an alternative geometry scaling factor.
  39     The default is \b 1000 (scale \em [m] to \em [mm]).
  40
  41 Note
  42     The cellTable information available in the files
  43     \c constant/cellTable and \c constant/polyMesh/cellTableId
  44     will be used if available. Otherwise the cellZones are used when
  45     creating the cellTable information.
  46
  47 See Also
  48     Foam::cellTable, Foam::meshWriter and Foam::meshWriters::STARCD
  49
  50 \*---------------------------------------------------------------------------*/
  51
  52 #include "argList.H"
  53 #include "timeSelector.H"
  54 #include "Time.H"
  55 #include "polyMesh.H"
  56 #include "STARCDMeshWriter.H"
  57
  58 using namespace Foam;
  59
  60 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  61 // Main program:
  62
  63 int main(int argc, char *argv[])
  64 {
  65     argList::addNote
  66     (
  67         "read OpenFOAM mesh and write a pro-STAR (v4) bnd/cel/vrt format"
  68     );
  69     argList::noParallel();
  70     timeSelector::addOptions();
  71
  72     argList::addOption
  73     (
  74         "scale",
  75         "factor",
  76         "geometry scaling factor - default is 1000 ([m] to [mm])"
  77     );
  78     argList::addBoolOption
  79     (
  80         "noBnd",
  81         "suppress writing a boundary (.bnd) file"
  82     );
  83
  84 #   include "setRootCase.H"
  85 #   include "createTime.H"
  86
  87     instantList timeDirs = timeSelector::select0(runTime, args);
  88
  89     fileName exportName = meshWriter::defaultMeshName;
  90     if (args.optionFound("case"))
  91     {
  92         exportName += '-' + args.globalCaseName();
  93     }
  94
  95     // default: rescale from [m] to [mm]
  96     scalar scaleFactor = 1000;
  97     if (args.optionReadIfPresent("scale", scaleFactor))
  98     {
  99         if (scaleFactor <= 0)
 100         {
 101             scaleFactor = 1;
 102         }
 103     }
 104
 105 #   include "createPolyMesh.H"
 106
 107     forAll(timeDirs, timeI)
 108     {
 109         runTime.setTime(timeDirs[timeI], timeI);
 110
 111 #       include "getTimeIndex.H"
 112
 113         polyMesh::readUpdateState state = mesh.readUpdate();
 114
 115         if (!timeI || state != polyMesh::UNCHANGED)
 116         {
 117             meshWriters::STARCD writer(mesh, scaleFactor);
 118
 119             if (args.optionFound("noBnd"))
 120             {
 121                 writer.noBoundary();
 122             }
 123
 124             fileName meshName(exportName);
 125             if (state != polyMesh::UNCHANGED)
 126             {
 127                 meshName += '_' + runTime.timeName();
 128             }
 129
 130             writer.write(meshName);
 131         }
 132
 133         Info<< nl << endl;
 134     }
 135
 136     Info<< "End\n" << endl;
 137
 138     return 0;
 139 }
 140
 141 // ************************************************************************* //