getMIIseacurrent.sh

Go to the documentation of this file.

#!/bin/zsh
script=${0##*/}

# Script to request and transform data from the MII instrumented line.
# 1 - a request is formed to the DB storing the data. The order of the fields is
#     fixed and must stay the same for the next steps :
#     time, current speed, error on current speed, current direction, error on current direction
#     A given period can be specified by providing two strings of the form 2019-05-01T12:00:00Z
#     The data are written into a ASC file (csv text file).
# 2 - The awk program splits the file in two separated text files with the header dropped for the
#     moment. The angles are converted from degrees with 0 at North and counting clockwise into
#     radians with 0 at East and counting anticlockwise. They are not further reduded to [-pi,+pi] 
#     to keep a continuity of values (the current is often directed toward West)

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  DIR                 $JPP_DIR/examples/JAcoustics/
set_variable  WORKDIR             ./

set_variable  MII_ASC        $WORKDIR/MII.asc
set_variable  CurSpeed_TXT   $WORKDIR/MII_CurrentSpeed.txt
set_variable  CurSpeed_ROOT  $WORKDIR/MII_CurrentSpeed.root
set_variable  CurSpeed_GIF   $WORKDIR/MII_CurrentSpeed.gif
set_variable  CurDirec_TXT   $WORKDIR/MII_CurrentDirection.txt
set_variable  CurDirec_ROOT  $WORKDIR/MII_CurrentDirection.root
set_variable  CurDirec_GIF   $WORKDIR/MII_CurrentDirection.gif

if ( do_usage $* ); then
    usage "$script <tmin> <tmax> or"\
          "\n$script clean" \
          "\nwhere <tmin> and <tmax> are strings of the form 2019-05-01T12:00:00Z" \
          "\nand clean will remove the output file" \
          "\n" \
          "\nScript to request sea current data from the erddap.osupytheas.fr database." \
          "\nthe validity of the datetimes is not checked by the script."
fi

case $# in
    1)  set_variable OPTION       $argv[1];&
    2)  set_variable TMAX         $argv[2];
        set_variable TMIN         $argv[1];;
    *)  fatal "Wrong number of arguments.";;
esac

function clean()
{
    rm -f $MII_ASC
    rm -f $CurSpeed_TXT
    rm -f $CurDirec_TXT
    rm -f $CurSpeed_ROOT
    rm -f $CurDirec_ROOT
    rm -f $CurSpeed_GIF
    rm -f $CurDirec_GIF
    exit
}
case "$OPTION" in
    clean) clean;;
    "")    ;;
    *)     fatal "Wrong option $OPTION.";;
esac


if [[ -n "$TMIN" && -n "$TMAX" ]]; then
  uTMIN=`date -d $TMIN +%s --utc`
  uTMAX=`date -d $TMAX +%s --utc`
  echo "Obtaining sea current speed and direction between" $TMIN "("$uTMIN") and" $TMAX "("$uTMAX")"
  wget -O $MII_ASC  "https://erddap.osupytheas.fr/erddap/tabledap/Emso_Aquadopp_64c7_8059_4af9.asc?time,Speed_mean,Speed_std,Direction_mean,Direction_std&time>=$uTMIN&time<=$uTMAX"
fi
# The ASC use the comma as separator. We split the file in two simple "space separated text files. 
# One for the strength of the sea current. One for the direction of the sea current.
# The direction value are converted from deg to rad and inverted from clockwise to anticlockwise
# The header is dropped for the moment since it has to be correctly handled using awk command. 
# It has currently 11 lines. A smarter way of handling the header and the fields must be found. 

if [[ -f $MII_ASC ]]; then
  awk -F',' '{if ( NR>11 ) {print $1 " " $2 " " $3;};}' $MII_ASC > $CurSpeed_TXT
  awk -F',' '{if ( NR>11 ) {print $1 " " (90.0 - $4)/57.2958 " " $5/57.2958 ;}}' $MII_ASC > $CurDirec_TXT
fi
if [[ -f $CurSpeed_TXT ]]; then
    JGraph \
  -f $CurSpeed_TXT    \
  -o $CurSpeed_ROOT

    JPlot1D \
  -f $CurSpeed_ROOT:\.\*               \
  -w 2000x500                          \
  -t "$TIMESTAMP"                      \
  -\> "Date"                           \
  -\^ "Current speed [m.s-1]"          \
  -T "Sea current speed at ORCA site " \
  -o $CurSpeed_GIF
fi
if [[ -f $CurDirec_TXT ]]; then
    JGraph \
  -f $CurDirec_TXT    \
  -o $CurDirec_ROOT

    JPlot1D \
  -f $CurDirec_ROOT:\.\*               \
  -w 2000x500                          \
  -t "$TIMESTAMP"                      \
  -\> "Date"                           \
  -\^ "Current direction [rad]"        \
  -T "Sea current direction at ORCA site " \
  -o $CurDirec_GIF
fi