Jpp 19.3.0-rc.2
the software that should make you happy
Loading...
Searching...
No Matches
JGeometry.sh
Go to the documentation of this file.
1#!/usr/bin/env zsh
2script=${0##*/}
3
4if [ -z $JPP_DIR ]; then
5 echo "Variable JPP_DIR undefined."
6 exit
7fi
8
9source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
10
11set_variable DEBUG 2
12set_variable WORKDIR ${TMPDIR:-/tmp}
13set_variable DIR $JPP_DIR/examples/JAcoustics/
14set_array GRID 18 60 570
15set_variable: FORMAT GRAPHICS_FORMAT gif
16set_variable+ BATCH GRAPHICS_BATCH -B
17
18if do_usage $*; then
19 usage "$script <a> <b>"
20fi
21
22if (( $# == 2 )); then
23 set_variable A $argv[1]
24 set_variable B $argv[2]
25else
26 fatal "Wrong number of arguments."
27fi
28
29
30for TX in 0.01 0.05 0.10; do
31
32 $DIR/JGeometry \
33 -M "$A $B" \
34 -S "$TX 0.00" \
35 -G "$GRID[*]" \
36 -o $WORKDIR/geometry\[${TX}\].txt
37
38 JGraph \
39 -f $WORKDIR/geometry\[${TX}\].txt \
40 -o $WORKDIR/geometry\[${TX}\].root >& /dev/null
41
42 tail -1 $WORKDIR/geometry\[${TX}\].txt | read XMAX ZMAX
43done
44
45set_variable GRAPHS `ls $WORKDIR/geometry\[*\].root | sed 's/.root/.root:\.\*/g' `
46
47let "XMAX = $XMAX + 2.0"
48let "YMAX = $GRID[3] + 10.0"
49
50JPlot1D \
51 -f "$GRAPHS" \
52 -x "0 $XMAX" \
53 -y "0 $YMAX" \
54 -> "x [m]" \
55 -\^ "z [m]" \
56 -G Y \
57 -L BR \
58 -T "" \
59 -d $DEBUG \
60 -o geometry_${A}_${B}.$FORMAT $BATCH