1 /*---------------------------------------------------------------------------*\
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 -------------------------------------------------------------------------------
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/>.
28 Calculates and writes the scalar field of pPrime2 (sqr(p - pMean)) at
31 \*---------------------------------------------------------------------------*/
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);
63 forAll(timeDirs, timeI)
65 runTime.setTime(timeDirs[timeI], timeI);
67 Info<< "Time = " << runTime.timeName() << endl;
78 if (pheader.headerOk())
82 Info<< " Reading p" << endl;
83 volScalarField p(pheader, mesh);
85 Info<< " Calculating pPrime2" << endl;
86 volScalarField pPrime2
101 Info<< " No p" << endl;
111 // ************************************************************************* //