Forward compatibility: flex
[foam-extend-3.2.git] / applications / utilities / postProcessing / scalarField / pPrime2 / pPrime2.C
bloba11c46c2832921752dcbe24676d0d5c68b71ac0d
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.
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.
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.
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/>.
24 Application
25     pPrime2
27 Description
28     Calculates and writes the scalar field of pPrime2 (sqr(p - pMean)) at
29     each time
31 \*---------------------------------------------------------------------------*/
33 #include "fvCFD.H"
36 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
38 int main(int argc, char *argv[])
40     timeSelector::addOptions();
42 #   include "setRootCase.H"
43 #   include "createTime.H"
45     instantList timeDirs = timeSelector::select0(runTime, args);
47 #   include "createMesh.H"
49     runTime.setTime(timeDirs[timeDirs.size()-1], timeDirs.size()-1);
51     volScalarField pMean
52     (
53         IOobject
54         (
55             "pMean",
56             runTime.timeName(),
57             mesh,
58             IOobject::MUST_READ
59         ),
60         mesh
61     );
63     forAll(timeDirs, timeI)
64     {
65         runTime.setTime(timeDirs[timeI], timeI);
67         Info<< "Time = " << runTime.timeName() << endl;
69         IOobject pheader
70         (
71             "p",
72             runTime.timeName(),
73             mesh,
74             IOobject::MUST_READ
75         );
77         // Check p exists
78         if (pheader.headerOk())
79         {
80             mesh.readUpdate();
82             Info<< "    Reading p" << endl;
83             volScalarField p(pheader, mesh);
85             Info<< "    Calculating pPrime2" << endl;
86             volScalarField pPrime2
87             (
88                 IOobject
89                 (
90                     "pPrime2",
91                     runTime.timeName(),
92                     mesh,
93                     IOobject::NO_READ
94                 ),
95                 sqr(p - pMean)
96             );
97             pPrime2.write();
98         }
99         else
100         {
101             Info<< "    No p" << endl;
102         }
104         Info<< endl;
105     }
107     return 0;
111 // ************************************************************************* //