applications/utilities/mesh/manipulation/mergeMeshes/mergeMeshes.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 Description
  25     Merge two meshes.
  26
  27 \*---------------------------------------------------------------------------*/
  28
  29 #include "argList.H"
  30 #include "Time.H"
  31 #include "mergePolyMesh.H"
  32
  33 using namespace Foam;
  34
  35 void getRootCase(fileName& casePath)
  36 {
  37     casePath.clean();
  38
  39     if (casePath.empty() || casePath == ".")
  40     {
  41         // handle degenerate form and '.'
  42         casePath = cwd();
  43     }
  44     else if (casePath[0] != '/' && casePath.name() == "..")
  45     {
  46         // avoid relative cases ending in '..' - makes for very ugly names
  47         casePath = cwd()/casePath;
  48         casePath.clean();
  49     }
  50 }
  51
  52
  53 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  54 // Main program:
  55
  56 int main(int argc, char *argv[])
  57 {
  58     argList::addNote
  59     (
  60         "merge two meshes"
  61     );
  62
  63     argList::noParallel();
  64     argList::validArgs.append("masterCase");
  65     argList::addOption
  66     (
  67         "masterRegion",
  68         "name",
  69         "specify alternative mesh region for the master mesh"
  70     );
  71
  72     argList::validArgs.append("addCase");
  73     argList::addOption
  74     (
  75         "addRegion",
  76         "name",
  77         "specify alternative mesh region for the additional mesh"
  78     );
  79
  80     argList args(argc, argv);
  81     if (!args.check())
  82     {
  83          FatalError.exit();
  84     }
  85
  86     fileName masterCase = args[1];
  87     word masterRegion = polyMesh::defaultRegion;
  88     args.optionReadIfPresent("masterRegion", masterRegion);
  89
  90     fileName addCase = args[2];
  91     word addRegion = polyMesh::defaultRegion;
  92     args.optionReadIfPresent("addRegion", addRegion);
  93
  94     getRootCase(masterCase);
  95     getRootCase(addCase);
  96
  97     Info<< "Master:      " << masterCase << "  region " << masterRegion << nl
  98         << "mesh to add: " << addCase    << "  region " << addRegion << endl;
  99
 100     #include "createTimes.H"
 101
 102     Info<< "Reading master mesh for time = " << runTimeMaster.timeName() << nl;
 103
 104     Info<< "Create mesh\n" << endl;
 105     mergePolyMesh masterMesh
 106     (
 107         IOobject
 108         (
 109             masterRegion,
 110             runTimeMaster.timeName(),
 111             runTimeMaster
 112         )
 113     );
 114
 115
 116     Info<< "Reading mesh to add for time = " << runTimeToAdd.timeName() << nl;
 117
 118     Info<< "Create mesh\n" << endl;
 119     polyMesh meshToAdd
 120     (
 121         IOobject
 122         (
 123             addRegion,
 124             runTimeToAdd.timeName(),
 125             runTimeToAdd
 126         )
 127     );
 128
 129     runTimeMaster++;
 130
 131     Info<< "Writing combined mesh to " << runTimeMaster.timeName() << endl;
 132
 133     masterMesh.addMesh(meshToAdd);
 134     masterMesh.merge();
 135     masterMesh.polyMesh::write();
 136
 137     Info<< "\nEnd\n" << endl;
 138
 139     return 0;
 140 }
 141
 142
 143 // ************************************************************************* //