applications/utilities/mesh/conversion/mshToFoam/mshToFoam.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2004-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     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 (at your 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 "Time.H"
  41 #include "polyMesh.H"
  42 #include "IFstream.H"
  43 #include "polyPatch.H"
  44 #include "ListOps.H"
  45 #include "cellModeller.H"
  46
  47 #include <fstream>
  48
  49 using namespace Foam;
  50
  51 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  52
  53 // Main program:
  54
  55 int main(int argc, char *argv[])
  56 {
  57     argList::noParallel();
  58     argList::validArgs.append(".msh file");
  59     argList::addBoolOption
  60     (
  61         "hex",
  62         "treat input as containing hex instead of tet cells"
  63     );
  64
  65 #   include "setRootCase.H"
  66 #   include "createTime.H"
  67
  68     const bool readHex = args.optionFound("hex");
  69     IFstream mshStream(args[1]);
  70
  71     label nCells;
  72     mshStream >> nCells;
  73
  74     if (readHex)
  75     {
  76         Info<< "Trying to read " << nCells << " hexes." << nl << endl;
  77     }
  78     else
  79     {
  80         Info<< "Trying to read " << nCells << " tets." << nl << endl;
  81     }
  82
  83     cellShapeList cells(nCells);
  84
  85     const cellModel& tet = *(cellModeller::lookup("tet"));
  86     const cellModel& hex = *(cellModeller::lookup("hex"));
  87
  88     labelList tetPoints(4);
  89     labelList hexPoints(8);
  90
  91     if (readHex)
  92     {
  93         for (label cellI = 0; cellI < nCells; cellI++)
  94         {
  95             for (label cp = 0; cp < 8; cp++)
  96             {
  97                 mshStream >> hexPoints[cp];
  98             }
  99             cells[cellI] = cellShape(hex, hexPoints);
 100         }
 101     }
 102     else
 103     {
 104         for (label cellI = 0; cellI < nCells; cellI++)
 105         {
 106             for (label cp = 0; cp < 4; cp++)
 107             {
 108                 mshStream >> tetPoints[cp];
 109             }
 110             cells[cellI] = cellShape(tet, tetPoints);
 111         }
 112     }
 113
 114
 115     label nPoints;
 116
 117     mshStream >> nPoints;
 118
 119     Info<< "Trying to read " << nPoints << " points." << endl << endl;
 120
 121     pointField points(nPoints);
 122
 123
 124     for (label pointI = 0; pointI < nPoints; pointI++)
 125     {
 126         scalar x, y, z;
 127
 128         mshStream >> x >> y >> z;
 129
 130         points[pointI] = point(x, y, z);
 131     }
 132
 133
 134     polyMesh mesh
 135     (
 136         IOobject
 137         (
 138             polyMesh::defaultRegion,
 139             runTime.constant(),
 140             runTime
 141         ),
 142         xferMove(points),
 143         cells,
 144         faceListList(0),
 145         wordList(0),
 146         wordList(0),
 147         "defaultFaces",
 148         polyPatch::typeName,
 149         wordList(0)
 150     );
 151
 152     Info<< "Writing mesh ..." << endl;
 153
 154     mesh.write();
 155
 156
 157     Info<< "End\n" << endl;
 158
 159     return 0;
 160 }
 161
 162
 163 // ************************************************************************* //