initial commit for version 1.6.x patch release
[OpenFOAM-1.6.x.git] / bin / foamGraphResKE
blobae9d713ef2561c4183eb70521a4cb41355d72275
1 #!/bin/csh -f
2 #------------------------------------------------------------------------------
3 # ========= |
4 # \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
5 # \\ / O peration |
6 # \\ / A nd | Copyright (C) 1991-2009 OpenCFD Ltd.
7 # \\/ M anipulation |
8 #-------------------------------------------------------------------------------
9 # License
10 # This file is part of OpenFOAM.
12 # OpenFOAM is free software; you can redistribute it and/or modify it
13 # under the terms of the GNU General Public License as published by the
14 # Free Software Foundation; either version 2 of the License, or (at your
15 # option) any later version.
17 # OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
18 # ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
19 # FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
20 # for more details.
22 # You should have received a copy of the GNU General Public License
23 # along with OpenFOAM; if not, write to the Free Software Foundation,
24 # Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
26 # Script
27 # foamGraphResKE
29 # Description
31 #------------------------------------------------------------------------------
33 set tmpK = newK$$.dat
34 cat $1 | grep "Solving for k" | grep -v "solution singularity" | sed s/,//g | awk 'BEGIN { NENTRIES = 0 }{NENTRIES++; printf("%f %e\n", NENTRIES, $8)}' - > $tmpK
36 set tmpEpsilon = newEpsilon$$.dat
37 cat $1 | grep "Solving for epsilon" | grep -v "solution singularity" | sed s/,//g | awk 'BEGIN { NENTRIES = 0 }{NENTRIES++; printf("%f %e\n", NENTRIES, $8)}' - > $tmpEpsilon
39 set linesK = `cat $tmpK | wc -l`
40 set linesEpsilon = `cat $tmpEpsilon | wc -l`
42 set ngraphs = 0
44 if ($linesK > 1) then
45 @ ngraphs = $ngraphs + 1
46 endif
48 if ($linesEpsilon > 1) then
49 @ ngraphs = $ngraphs + 1
50 endif
52 cat << EOF > residualKE.dat
53 Solver Performance
54 number of iterations
55 residual
56 0 0
57 0 0
58 $ngraphs
59 EOF
61 if ($linesK > 1) then
62 cat << EOF >> residualKE.dat
65 $linesK
66 EOF
67 cat $tmpK >> residualKE.dat
68 endif
70 if ($linesEpsilon > 1) then
71 cat << EOF >> residualKE.dat
72 epsilon
74 $linesEpsilon
75 EOF
76 cat $tmpEpsilon >> residualKE.dat
77 endif
79 rm $tmpK
80 rm $tmpEpsilon
82 #------------------------------------------------------------------------------