applications/utilities/mesh/conversion/mshToFoam/mshToFoam.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 Description
  25     Converts .msh file generated by the Adventure system.
  26
  27     Note: the .msh format does not contain any boundary information. It is
  28     purely a description of the internal mesh.
  29
  30     Can read both linear-tet format (i.e. 4 verts per tet) and linear-hex
  31     format (8 verts per hex) (if provided with the -hex option)
  32     (Note: will bomb out if not supplied with the correct option for the
  33      file format)
  34
  35     Not extensively tested.
  36
  37 \*---------------------------------------------------------------------------*/
  38
  39 #include "argList.H"
  40 #include "objectRegistry.H"
  41 #include "foamTime.H"
  42 #include "polyMesh.H"
  43 #include "IFstream.H"
  44 #include "polyPatch.H"
  45 #include "ListOps.H"
  46 #include "cellModeller.H"
  47
  48 #include <fstream>
  49
  50 using namespace Foam;
  51
  52 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  53
  54 // Main program:
  55
  56 int main(int argc, char *argv[])
  57 {
  58     argList::noParallel();
  59     argList::validArgs.append(".msh file");
  60     argList::validOptions.insert("hex", "");
  61
  62 #   include "setRootCase.H"
  63 #   include "createTime.H"
  64
  65     bool readHex = args.optionFound("hex");
  66
  67     fileName mshFile(args.additionalArgs()[0]);
  68     IFstream mshStream(mshFile);
  69
  70     label nCells;
  71     mshStream >> nCells;
  72
  73     if (readHex)
  74     {
  75         Info<< "Trying to read " << nCells << " hexes." << endl << endl;
  76     }
  77     else
  78     {
  79         Info<< "Trying to read " << nCells << " tets." << endl << endl;
  80     }
  81
  82     cellShapeList cells(nCells);
  83
  84     const cellModel& tet = *(cellModeller::lookup("tet"));
  85     const cellModel& hex = *(cellModeller::lookup("hex"));
  86
  87     labelList tetPoints(4);
  88     labelList hexPoints(8);
  89
  90     if (readHex)
  91     {
  92         for (label cellI = 0; cellI < nCells; cellI++)
  93         {
  94             for (label cp = 0; cp < 8; cp++)
  95             {
  96                 mshStream >> hexPoints[cp];
  97             }
  98             cells[cellI] = cellShape(hex, hexPoints);
  99         }
 100     }
 101     else
 102     {
 103         for (label cellI = 0; cellI < nCells; cellI++)
 104         {
 105             for (label cp = 0; cp < 4; cp++)
 106             {
 107                 mshStream >> tetPoints[cp];
 108             }
 109             cells[cellI] = cellShape(tet, tetPoints);
 110         }
 111     }
 112
 113
 114     label nPoints;
 115
 116     mshStream >> nPoints;
 117
 118     Info<< "Trying to read " << nPoints << " points." << endl << endl;
 119
 120     pointField points(nPoints);
 121
 122
 123     for (label pointI = 0; pointI < nPoints; pointI++)
 124     {
 125         scalar x, y, z;
 126
 127         mshStream >> x >> y >> z;
 128
 129         points[pointI] = point(x, y, z);
 130     }
 131
 132
 133     polyMesh mesh
 134     (
 135         IOobject
 136         (
 137             polyMesh::defaultRegion,
 138             runTime.constant(),
 139             runTime
 140         ),
 141         xferMove(points),
 142         cells,
 143         faceListList(0),
 144         wordList(0),
 145         wordList(0),
 146         "defaultFaces",
 147         polyPatch::typeName,
 148         wordList(0)
 149     );
 150
 151     Info<< "Writing mesh ..." << endl;
 152
 153     mesh.write();
 154
 155
 156     Info<< "End\n" << endl;
 157
 158     return 0;
 159 }
 160
 161
 162 // ************************************************************************* //