sound-velocity.sh

Go to the documentation of this file.

#!/usr/bin/env 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   WORKDIR     ${TMPDIR:-/tmp}/
set_variable   DETECTOR    49
set_variable:  FORMAT      GRAPHICS_FORMAT        gif
set_variable+  BATCH       GRAPHICS_BATCH         -B
 
if do_usage $*; then
    usage "$script <input file>"\
          "\nExample input file https://git.km3net.de/calibration/input_tables/-/blob/master/data/sound/SV_Profile_ORCA_250120.txt"
fi
 
if (( $# != 1 )); then
    fatal "Wrong number of arguments."
fi
 
DIR=`pwd`
 
cd $WORKDIR
 
rm -f $WORKDIR/sound_velocity.txt
 
JAcoustics.sh $DETECTOR
 
tail -1 sound_velocity.txt | read A B DEPTH
 
cd $DIR
 
echo $A $B $DEPTH
 
JF1 \
    -F "$A - 1.0 * $B * (x + $DEPTH)"  \
    -x "2000 2600"                     \
    -o $WORKDIR/f1.root
 
JGraph -f $1 -o $WORKDIR/graph.root
 
JPlot1D \
    -f $WORKDIR/graph.root:\.\*        \
    -f $WORKDIR/f1.root:\.\*           \
    -> "D [m]"                        \
    -\^ "V [m/s]"                      \
    -T ""                              \
    -o sound-velocity.$FORMAT          $BATCH