sound-velocity.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   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