JGeometry.sh

Go to the documentation of this file.

#!/bin/zsh
script=${0##*/}
 
if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi

source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null

set_variable   DEBUG       2                              
set_variable   WORKDIR     ${TMPDIR:-/tmp}                           
set_variable   DIR         $JPP_DIR/examples/JAcoustics/  
set_array      GRID        18 60 600
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B

if do_usage $*; then
    usage "$script <a> <b>"
fi

if (( $# == 2 )); then
    set_variable  A               $argv[1]
    set_variable  B               $argv[2]
else
    fatal "Wrong number of arguments."
fi


for TX in 0.01 0.05 0.10; do

    $DIR/JGeometry \
        -M "$A $B"                           \
        -S "$TX 0.00"                        \
        -G "$GRID[*]"                        \
        -o $WORKDIR/geometry\[${TX}\].txt    >& /dev/null

    JGraph \
        -f $WORKDIR/geometry\[${TX}\].txt    \
        -o $WORKDIR/geometry\[${TX}\].root   >& /dev/null

    tail -1 $WORKDIR/geometry\[${TX}\].txt   | read XMAX ZMAX
done

set_variable  GRAPHS  `ls $WORKDIR/geometry\[*\].root | sed 's/.root/.root:\.\*/g' `

let "XMAX = $XMAX    +  2.0"
let "YMAX = $GRID[3] + 10.0"

JPlot1D \
    -f "$GRAPHS"                     \
    -x "0 $XMAX"                     \
    -y "0 $YMAX"                     \
    -\> "x [m]"                      \
    -\^ "z [m]"                      \
    -L BR                            \
    -T ""                            \
    -d $DEBUG                        \
    -o geometry_${A}_${B}.$FORMAT    $BATCH