#!/bin/sh

. /etc/profile.local

#select LINE_ID, min(DATE_ACQ),max(DATE_ACQ),min(LON_ACQ),max(LAT_ACQ),min(TIME_ACQ),max(TIME_ACQ)  from NAV where CAM_ID='CA96' group by LINE_ID; 
# select LINE_ID, min(DATE_ACQ),max(DATE_ACQ),min(LON_ACQ),max(LON_ACQ),min(LAT_ACQ), max(LAT_ACQ),min(TIME_ACQ),max(TIME_ACQ)  from NAV where CAM_ID='VEMA00' group by LINE_ID; 

gmtset ANOT_FONT_SIZE 10 LABEL_FONT_SIZE 10

# GPL90 MISSING LAT-LON ONE LINE 
CRUISES="BOUVET94 BOUVET96"  # spostare NW-NE INSET
CRUISES="PRIMAR92 PRIMAR93 PRIMAR94 PRIMAR96"    # CAMBIARE -B troppo esteso lo span lon lat ..
CRUISES="VEMA00 VEMA92 VEMA98"  # spostare NW-NE INSET
CRUISES="GS90 GB90  CM91 PA91 PG91 PL91 GP92 ALB92 CM94 ADS94  AP96 CA96 BS98"
CRUISES="BAL84 AC85 ACC86  PS87 VP87 MATS87 PS88 INVO89"  # MATS87 spostare NW INSET
CRUISES="MESC2001 VOL2002 M2001 M2005  PANA07"
CRUISES="AZ99 M2000 V2002"  ######   NO DATA INSERIRE LE LINEE 
CRUISES="SCPM04 RS05"   ##########   LON LAT INVERTITE !!!
CRUISES="AB06"  #### MANCANO DATI ???
CRUISES="MRS06" ##### TOGLIERER LINEE A N !
CRUISES="CALAMARE08 ADR0208 MNG0109 MNG0209 M2009"  #########  INSERIRE
CRUISES="TY"
CRUISES="ETNA93"
CRUISES="RIFANO92"
CRUISES="CNR"

HOST="-h tornado.bo.ismar.cnr.it"
HOST=

for c in $CRUISES ; do
  mysql IGMNAV  $HOST  << EOF > IMG/NAV_DATA.${c}
    select LON_ACQ,LAT_ACQ,NFIX,LINE_ID,CAM_ID,DATE_ACQ,TIME_ACQ \
         from NAV where CAM_ID='${c}' order by CAM_ID,LINE_ID,NFIX;
EOF

  mysql IGMNAV  $HOST << EOF > IMG/LINES_DATA.${c}
    select *  from NAV_LINES where CAM_ID='${c}' order by DATE_S,TIME_S;
EOF

  mysql IGMNAV $HOST <<EOF > IMG/MINMAX.$c
    select min(LON_ACQ),max(LON_ACQ),min(LAT_ACQ),max(LAT_ACQ), min(DATE_ACQ), max(DATE_ACQ) \
      from NAV where CAM_ID='${c}';
EOF

  B=-Bg15ma30mWSen;  XOFF=0.5 ; YOFF=0.5; st_bound="lonw=10 lone=10 lats=10 latn=10"

  case $c in

  VEMA00) B=-Bg2a4WSen ; st_bound="lonw=10 lone=10 lats=5 latn=10";;
  VEMA92) YOFF=8.5 ;;
  VEMA98) B=-Bg30ma1f5mWSen ; st_bound="lonw=10 lone=10 lats=1.5 latn=10" ;;
  PRIMAR92) B=-Bg30ma1f5mWSen ; st_bound="lonw=10 lone=10 lats=2 latn=10" ;;
  PRIMAR93) B=-Bg1a4f10mWSen ; st_bound="lonw=10 lone=10 lats=5 latn=10" ;;
  PRIMAR94) B=-Bg2a4f10mWSen  ;;
  PRIMAR96) B=-Bg1a2f15mWSen ; st_bound="lonw=10 lone=10 lats=1 latn=10" ;;
  MATS87)  st_bound="lonw=10 lone=10 lats=2 latn=10" ;;
  BOUVET94) B=-Bg1a2f15mWSen ; st_bound="lonw=10 lone=10 lats=5 latn=10" ;;
  BOUVET96) B=-Bg1a2f15mWSen ; st_bound="lonw=10 lone=10 lats=5 latn=10" ;;
  M2001) B=-Bg10ma30mf1mWSen ; st_bound="lonw=10 lone=10 lats=1 latn=10" ;;
  VOL2002) B=-Bg1a5f10mWSen ; st_bound="lonw=10 lone=10 lats=1.5 latn=10" ;;
  PANA07)  st_bound="lonw=10 lone=10 lats=1.5 latn=10" ;;
  RIFANO92)  B=-Bg1a2f10mWSen ; st_bound="lonw=10 lone=10 lats=5 latn=10" ;;

  esac

  JJ=`perl -ane 'next if(/ACQ/);$slat=($F[2]+$F[3])/2;printf"M0/%d",$slat;' IMG/MINMAX.$c`

#  RR=( `perl -ane 'next if(/ACQ/);$slon=($F[1]-$F[0])/10;$slat=($F[3]-$F[2])/10; 
#    print $F[0]-$slon," ",$F[1]+$slon," ",$F[2]-$slat," ",$F[3]+$slat,"\n";' MINMAX.$c` )
  RR=( `awk 'NR>1{slon=($2-$1);slat=($4-$3);print $1-slon/lonw,$2+slon/lone,$3-slat/lats,$4+slat/latn}' \
    $st_bound IMG/MINMAX.$c` )

  echo $JJ ${RR[0]} ${RR[1]} ${RR[2]} ${RR[3]}


  PS=IMG/${c}
  DIM=18.5
  #echo $JJ |tr [A-Z] [a-z]`
  J="-J$JJ/$DIM" 
  R="-R${RR[0]}/${RR[1]}/${RR[2]}/${RR[3]}"
  echo "Done with $c .. $R .. $J .. $PS "
  pscoast $B $R $J -Df -W1 -X1.90 -Y2.5 -UISMAR-CNR -K -P > $PS
  psxy -H1 IMG/NAV_DATA.${c} $R $J -Sc0.015 -K -O >> $PS
  RR=`perl -ane 'next if /ACQ/;printf "%s/%s/%s/%s", $F[0]-5,$F[1]+5,$F[2]-3.5,$F[3]+3.5;' IMG/MINMAX.$c`
  R=-R$RR
  J=-J$JJ/3.5
  echo $R $J
  gmtset ANOT_FONT_SIZE 8 LABEL_FONT_SIZE 8
  pscoast -Dh $R $J -Bg1a5wsEN -X$XOFF -Y$YOFF -W1 -Gkhaki -Slightblue -K -O >> $PS
  RR=( `perl -ane 'next if(/ACQ/); print "$F[0] $F[1] $F[2] $F[3]\n";' IMG/MINMAX.$c` )
  psxy $R $J -m -K -O -W5,red <<EOF >> $PS
  ${RR[0]} ${RR[2]}
  ${RR[0]} ${RR[3]}
  ${RR[1]} ${RR[3]}
  ${RR[1]} ${RR[2]}
  ${RR[0]} ${RR[2]}
EOF
  R=-R0/$DIM/0/$DIM; J=-JX$DIM
  DATES=`awk 'NR>1{print "START: ",$5," END : ",$6}' IMG/MINMAX.$c`
  pstext $R $J -K -O <<EOF >> $PS
  4.5 0.25 14 0 1 1 $c $DATES
EOF

  echo "showpage" >> $PS
  ps2epsi $PS $PS.eps
  convert $PS.eps $PS.jpg
  rm $PS

done