4 # Script to request and transform data from the MII instrumented line.
5 # 1 - a request is formed to the database storing the data. The order of the fields is
6 # fixed and must stay the same for the next steps :
7 # time, current speed, error on current speed, current direction, error on current direction
8 # A given period can be specified by providing two strings of the form 2019-05-01T12:00:00Z
9 # The data are written into a ASCII file (csv text file).
10 # 2 - The awk program splits the file in two separated text files with the header dropped for the
11 # moment. The angles are converted from degrees with 0 at North and counting clockwise into
12 # radians with 0 at East and counting anticlockwise. They are not further reduded to [-pi,+pi]
13 # to keep a continuity of values (the current is often directed toward West)
15 if [ -z $JPP_DIR ]; then
16 echo "Variable JPP_DIR undefined."
20 source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null
23 set_variable DIR $JPP_DIR/examples/JAcoustics/
24 set_variable WORKDIR ./
25 set_variable: FORMAT GRAPHICS_FORMAT gif
26 set_variable+ BATCH GRAPHICS_BATCH -B
28 set_variable MII_ASC $WORKDIR/MII.asc
29 set_variable CurSpeed_TXT $WORKDIR/MII_CurrentSpeed.txt
30 set_variable CurSpeed_ROOT $WORKDIR/MII_CurrentSpeed.root
31 set_variable CurSpeed_GIF $WORKDIR/MII_CurrentSpeed.$FORMAT
32 set_variable CurDirec_TXT $WORKDIR/MII_CurrentDirection.txt
33 set_variable CurDirec_ROOT $WORKDIR/MII_CurrentDirection.root
34 set_variable CurDirec_GIF $WORKDIR/MII_CurrentDirection.$FORMAT
37 usage "$script <tmin> <tmax> or"\
39 "\nwhere <tmin> and <tmax> are strings of the form 2019-05-01T12:00:00Z" \
40 "\nand clean will remove the output file" \
42 "\nScript to request sea current data from the erddap.osupytheas.fr database." \
43 "\nthe validity of the datetimes is not checked by the script."
47 1) set_variable OPTION $argv[1];&
48 2) set_variable TMAX $argv[2];
49 set_variable TMIN $argv[1];;
50 *) fatal "Wrong number of arguments.";;
68 *) fatal "Wrong option $OPTION.";;
72 if [[ -n "$TMIN" && -n "$TMAX" ]]; then
73 uTMIN=`date -d $TMIN +%s --utc`
74 uTMAX=`date -d $TMAX +%s --utc`
75 echo "Obtaining sea current speed and direction between" $TMIN "("$uTMIN") and" $TMAX "("$uTMAX")"
76 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"
79 # The ASCII use the comma as separator. We split the file in two simple space separated text files.
80 # One for the strength of the sea current. One for the direction of the sea current.
81 # The direction value are converted from deg to rad and inverted from clockwise to anticlockwise
82 # The header is dropped for the moment since it has to be correctly handled using awk command.
83 # It has currently 11 lines. A smarter way of handling the header and the fields must be found.
85 if [[ -f $MII_ASC ]]; then
86 awk -F',' '{if ( NR>11 ) {print $1 " " $2 " " $3;};}' $MII_ASC > $CurSpeed_TXT
87 awk -F',' '{if ( NR>11 ) {print $1 " " (90.0 - $4)/57.2958 " " $5/57.2958 ;}}' $MII_ASC > $CurDirec_TXT
90 if [[ -f $CurSpeed_TXT ]]; then
97 -f $CurSpeed_ROOT:\.\* \
101 -\^ "Current speed [m.s-1]" \
102 -T "Sea current speed at ORCA site " \
103 -o $CurSpeed_GIF $BATCH
106 if [[ -f $CurDirec_TXT ]]; then
113 -f $CurDirec_ROOT:\.\* \
117 -\^ "Current direction [rad]" \
118 -T "Sea current direction at ORCA site " \
119 -o $CurDirec_GIF $BATCH