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 database 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 ASCII 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:  FORMAT         GRAPHICS_FORMAT        gif
set_variable+  BATCH          GRAPHICS_BATCH         -B

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.$FORMAT
set_variable   CurDirec_TXT   $WORKDIR/MII_CurrentDirection.txt   
set_variable   CurDirec_ROOT  $WORKDIR/MII_CurrentDirection.root  
set_variable   CurDirec_GIF   $WORKDIR/MII_CurrentDirection.$FORMAT

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 ASCII 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                          $BATCH
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                          $BATCH
fi