twoPhaseEulerFoam:frictionalStressModel/Schaeffer: Correct mut on processor boundaries
[OpenFOAM-1.7.x.git] / bin / foamGraphResKE
blob8cdd621481355ab05116307ba6b1696f707f9eae
1 #!/bin/csh -f
2 #---------------------------------*- sh -*-------------------------------------
3 # ========= |
4 # \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
5 # \\ / O peration |
6 # \\ / A nd | Copyright (C) 1991-2010 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
14 # the Free Software Foundation, either version 3 of the License, or
15 # (at your 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, see <http://www.gnu.org/licenses/>.
25 # Script
26 # foamGraphResKE
28 # Description
30 #------------------------------------------------------------------------------
32 set tmpK = newK$$.dat
33 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
35 set tmpEpsilon = newEpsilon$$.dat
36 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
38 set linesK = `cat $tmpK | wc -l`
39 set linesEpsilon = `cat $tmpEpsilon | wc -l`
41 set ngraphs = 0
43 if ($linesK > 1) then
44 @ ngraphs = $ngraphs + 1
45 endif
47 if ($linesEpsilon > 1) then
48 @ ngraphs = $ngraphs + 1
49 endif
51 cat << EOF > residualKE.dat
52 Solver Performance
53 number of iterations
54 residual
55 0 0
56 0 0
57 $ngraphs
58 EOF
60 if ($linesK > 1) then
61 cat << EOF >> residualKE.dat
64 $linesK
65 EOF
66 cat $tmpK >> residualKE.dat
67 endif
69 if ($linesEpsilon > 1) then
70 cat << EOF >> residualKE.dat
71 epsilon
73 $linesEpsilon
74 EOF
75 cat $tmpEpsilon >> residualKE.dat
76 endif
78 rm $tmpK
79 rm $tmpEpsilon
81 #------------------------------------------------------------------------------