applications/utilities/mesh/manipulation/rotateMesh/rotateMesh.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 anispulation  |
   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 Application
  25     rotateMesh
  26
  27 Description
  28     Rotates the mesh and fields from the direcion n1 to the direction n2.
  29
  30 \*---------------------------------------------------------------------------*/
  31
  32 #include "argList.H"
  33 #include "timeSelector.H"
  34 #include "objectRegistry.H"
  35 #include "foamTime.H"
  36 #include "fvMesh.H"
  37 #include "volFields.H"
  38 #include "surfaceFields.H"
  39 #include "transformGeometricField.H"
  40 #include "IOobjectList.H"
  41
  42 using namespace Foam;
  43
  44 template<class GeometricField>
  45 void RotateFields
  46 (
  47     const fvMesh& mesh,
  48     const IOobjectList& objects,
  49     const tensor& T
  50 )
  51 {
  52     // Search list of objects for volScalarFields
  53     IOobjectList fields(objects.lookupClass(GeometricField::typeName));
  54
  55     forAllIter(IOobjectList, fields, fieldIter)
  56     {
  57         Info<< "    Rotating " << fieldIter()->name() << endl;
  58
  59         GeometricField theta(*fieldIter(), mesh);
  60         transform(theta, dimensionedTensor(T), theta);
  61         theta.write();
  62     }
  63 }
  64
  65
  66 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  67
  68 int main(int argc, char *argv[])
  69 {
  70     timeSelector::addOptions();
  71
  72     argList::validArgs.append("n1");
  73     argList::validArgs.append("n2");
  74
  75 #   include "setRootCase.H"
  76 #   include "createTime.H"
  77
  78     vector n1(IStringStream(args.additionalArgs()[0])());
  79     n1 /= mag(n1);
  80
  81     vector n2(IStringStream(args.additionalArgs()[1])());
  82     n2 /= mag(n2);
  83
  84     tensor T = rotationTensor(n1, n2);
  85
  86     {
  87         pointIOField points
  88         (
  89             IOobject
  90             (
  91                 "points",
  92                 runTime.findInstance(polyMesh::meshSubDir, "points"),
  93                 polyMesh::meshSubDir,
  94                 runTime,
  95                 IOobject::MUST_READ,
  96                 IOobject::NO_WRITE,
  97                 false
  98             )
  99         );
 100
 101         points = transform(T, points);
 102
 103         // Set the precision of the points data to 10
 104         IOstream::defaultPrecision(10);
 105
 106         Info << "Writing points into directory " << points.path() << nl << endl;
 107         points.write();
 108     }
 109
 110
 111     instantList timeDirs = timeSelector::select0(runTime, args);
 112
 113 #   include "createMesh.H"
 114
 115
 116     forAll(timeDirs, timeI)
 117     {
 118         runTime.setTime(timeDirs[timeI], timeI);
 119
 120         Info<< "Time = " << runTime.timeName() << endl;
 121
 122         // Search for list of objects for this time
 123         IOobjectList objects(mesh, runTime.timeName());
 124
 125         RotateFields<volVectorField>(mesh, objects, T);
 126         RotateFields<volSphericalTensorField>(mesh, objects, T);
 127         RotateFields<volSymmTensorField>(mesh, objects, T);
 128         RotateFields<volTensorField>(mesh, objects, T);
 129
 130         RotateFields<surfaceVectorField>(mesh, objects, T);
 131         RotateFields<surfaceSphericalTensorField>(mesh, objects, T);
 132         RotateFields<surfaceSymmTensorField>(mesh, objects, T);
 133         RotateFields<surfaceTensorField>(mesh, objects, T);
 134     }
 135
 136     Info<< "\nEnd.\n" << endl;
 137
 138     return 0;
 139 }
 140
 141
 142 // ************************************************************************* //