diff --git a/bin/foamProbe b/bin/foamProbe new file mode 100755 index 000000000..3f155579b --- /dev/null +++ b/bin/foamProbe @@ -0,0 +1,135 @@ +#!/bin/bash +#------------------------------------------------------------------------------ +# ========= | +# \\ / 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 +# foamProbe +# +# Description +# Proccess a scalar, vetorial or tensorial field of OpenFOAM +# probe file (for monitoring points). The original file will not be changed +# and the fields of each monitoring point will be stored in a new +# directory named probe (where is the field name) in the +# directory where are the original probes files. +# The fourth optional parameter can be used to create xmgrace sequentially +# graphs (case is 0) +# for monitoring points or automatically create eps figures (case is 1). + +# +# Author: +# Jovani L. Favero, J. F. Mitre (2009) +# +#------------------------------------------------------------------------------ + +if [ $# -lt 3 ]; then + echo "Usage: foamProbe [ ]" + exit 1 +fi + +# File name, NAME. +NAME=$(basename "$1") +DNAME="$1" +shift + +# Number of points, P. +P="$1" +shift + +NAME=$(basename "$1") +CNAME="$1" +shift +xmgrace=0 + +if [ $# -eq 1 ]; then +# Choose to save eps. + save="$1" + xmgrace=1 +fi + +if [ ! -d "$DNAME" ]; then + echo "Directory $DNAME does not exist." + exit 1 +fi + +if [ $P -lt 1 ]; then + echo "Can not be less than 1 single monitoring point." + exit 1 +fi + +cd $DNAME + +if [ ! -f "$CNAME" ]; then + echo "Field $CNAME does not exist." + exit 1 +fi + +# Screen information +echo "" +echo "Wait: Processing $P point(s) in $CNAME file ...." + +# Base directory to place processed files +DIR="./probe$NAME/" +if [ ! -d "$DIR" ]; then + mkdir "$DIR" +fi +sed -e "s/ *(/\t/g" -e 's/)//g' "$CNAME" >"$DIR/$NAME" +cd $DIR +column=1 +while [ "$column" -le "$P" ]; do + point=$column + column=`expr $column + 1` + cut -f 1,$column "$NAME" |sed -e "s/\t/ /g" -e '/^#/d' >$NAME.base + echo -e "# Time \t Point_$point" >$NAME\_$point + cat $NAME.base >>$NAME\_$point + rm -f $NAME.base +done +cd - &>/dev/null + +if [ "$xmgrace" = 1 ]; then + cd probe$CNAME/ + field=1 + underscore=_ + while [ "$field" -le "$P" ]; do + point=$field + field=`expr $field + 1` + echo " Opening file $CNAME$underscore$point" + if [ "$save" = 0 ]; then + echo " Opened $CNAME$underscore$point" + xmgrace -nxy $CNAME\_$point -noask + fi + + if [ "$save" -ne 0 ]; then + xmgrace -nxy $CNAME\_$point -hardcopy -printfile $CNAME\_$point.eps -hdevice EPS + echo " Saving file $CNAME$underscore$point.eps" + fi + echo " Closing file $CNAME$underscore$point" + echo "" + done +fi + +echo "" +echo "Done!!!" +echo "" + +exit 0