openfoam/bin/foamGraphResUVWP
2017-01-26 20:18:50 +00:00

65 lines
2.0 KiB
Bash
Executable File

#!/bin/sh
#------------------------------------------------------------------------------
# ========= |
# \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
# \\ / O peration |
# \\ / A nd | Copyright (C) 2011-2017 OpenFOAM Foundation
# \\/ 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 3 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, see <http://www.gnu.org/licenses/>.
#
# Script
# foamGraphResUVWP
#
# Description
# Deprecated script extract initial velocity and pressure residuals for each
# time-step from a log file for graphing.
#
# Superseded by the more general foamLog script.
#------------------------------------------------------------------------------
Script=${0##*/}
if [ ! -f "$1" ]; then
echo "$Script: $1 is not a valid filename"
exit 1
fi
scan()
{
tmpFile=new$1$$.dat
cat $2 | \
grep "Solving for $1" | \
grep -v "solution singularity" | \
sed s/,//g | \
awk 'BEGIN { NENTRIES = 0 }{NENTRIES++; printf("%d %e\n", NENTRIES, $8)}' - > $tmpFile
if [ "$(cat $tmpFile | wc -l)" -gt 1 ]; then
echo >> residualUVWP.dat
cat $tmpFile >> residualUVWP.dat
fi
rm $tmpFile
}
rm -f residualUVWP.dat
scan "Ux" $1
scan "Uy" $1
scan "Uz" $1
scan "p" $1
#------------------------------------------------------------------------------