This repository has been archived on 2023-11-20. You can view files and clone it, but cannot push or open issues or pull requests.
foam-extend4.1-coherent-io/bin/foamGraphResKE
2011-01-13 17:21:49 +00:00

82 lines
2.3 KiB
Tcsh
Executable file

#!/bin/csh -f
#------------------------------------------------------------------------------
# ========= |
# \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
# \\ / O peration |
# \\ / A nd | Copyright held by original author
# \\/ M anipulation |
#------------------------------------------------------------------------------
# License
# This file is part of OpenFOAM.
#
# OpenFOAM is free software; you can redistribute it and/or modify it
# under the terms of the GNU General Public License as published by the
# Free Software Foundation; either version 2 of the License, or (at your
# option) any later version.
#
# OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
# FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
# for more details.
#
# You should have received a copy of the GNU General Public License
# along with OpenFOAM; if not, write to the Free Software Foundation,
# Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
#
# Script
# foamGraphResKE
#
# Description
#
#------------------------------------------------------------------------------
set tmpK = newK$$.dat
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
set tmpEpsilon = newEpsilon$$.dat
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
set linesK = `cat $tmpK | wc -l`
set linesEpsilon = `cat $tmpEpsilon | wc -l`
set ngraphs = 0
if ($linesK > 1) then
@ ngraphs = $ngraphs + 1
endif
if ($linesEpsilon > 1) then
@ ngraphs = $ngraphs + 1
endif
cat << EOF > residualKE.dat
Solver Performance
number of iterations
residual
0 0
0 0
$ngraphs
EOF
if ($linesK > 1) then
cat << EOF >> residualKE.dat
k
0
$linesK
EOF
cat $tmpK >> residualKE.dat
endif
if ($linesEpsilon > 1) then
cat << EOF >> residualKE.dat
epsilon
0
$linesEpsilon
EOF
cat $tmpEpsilon >> residualKE.dat
endif
rm $tmpK
rm $tmpEpsilon
#------------------------------------------------------------------------------