#!/bin/sh
set -x

. $HOME/hl_home/OPR/progressPP.log
prevdtg=`$HOME/bin/mandtg $DTGPP + -6`
todir=$HOME/public_html/$prevdtg

date=`$HOME/bin/mandtg -date $prevdtg`
hour=`$HOME/bin/mandtg -hour $prevdtg`

[ -d $todir ] || mkdir $todir               # Create directory if not there yet.
[ -r $todir/use_span$hour.png ] && exit     # Already processed.

wd=/data/hirlam/WD$$
mkdir $wd || exit
cd $wd || exit

# Get the obsa file with Surface Parameter Analysis information in readable form

gunzip -c $ARCHDIR/obsa${date}_$hour.gz > span_list.dat

# Surface Parameter analysis data; uses land/sea mask for grid definition from first boundary

$HOME/PlotsOfTheDay/use-span `ls -1t /data/hirlam/hl_home/OPR/mf2*f | head -1` span_list.dat
mv SPAN_COUNT.txt $todir

. $HOME/PlotsOfTheDay/LONLAT.def

ln -s $HOME/PlotsOfTheDay/coastline.txt .

for h in 00 -3 -2 -1 +0 +1 +2
do

# Process the resulting text files and move result to $HOME/public_html/$todir

   if [ $h == "00" ]
   then
      dtgoff=${date}00
      houroff=00
   else
      dtgoff=`$HOME/bin/mandtg $prevdtg + $h`
      houroff=`$HOME/bin/mandtg -hour $dtgoff`
   fi
   plotcom="splot "

   if [ -s SPAN$houroff.txt -o -s PSEUDO$houroff.txt ]
   then

      [ -r PSEUDO$houroff.txt ]  && plotcom="$plotcom 'PSEUDO$houroff.txt'  title 'pseudo ${h}H'  with points palette pointtype 7 pointsize 0.05,"
      [ -r SPAN$houroff.txt ]    && plotcom="$plotcom 'SPAN$houroff.txt' title 'span ${h}H' with points palette pointtype 1,"
                                    plotcom="$plotcom 'coastline.txt' notitle with points pointtype 7 pointsize 0.05 linecolor 'medium-blue'"

# Construct the plot in PNG format

      gnuplot << EOF
         set terminal png truecolor size $XMAP,$YMAP
         set output 'use_span$houroff.png'
         set view map
         $XRANGE
         $YRANGE
         set pointsize 0.5
         set palette defined (0.0 'dark-gray', 0.5 'dark-gray', 0.5 'green', 1.5 'green', 1.5 'yellow', 2.5 'yellow', 2.5 'orange', 3.5 'orange', 3.5 'red', 4.0 'red')
         unset colorbox
         set zrange [0:4]
         set cbrange [0:4]
         set title 'Gebruik Oppervlaktewaarnemingen rond $dtgoff UTC (groen: gebruikt; rood: verworpen).'
         $plotcom
EOF

      mv use_span$houroff.png $todir

# Construct the plot in Postscript format

      gnuplot << EOF
         set terminal postscript landscape color
         set output 'use_span$houroff.ps'
         set view map
         $XRANGE
         $YRANGE
         set pointsize 0.25
         set palette defined (0.0 'dark-gray', 0.5 'dark-gray', 0.5 'green', 1.5 'green', 1.5 'yellow', 2.5 'yellow', 2.5 'orange', 3.5 'orange', 3.5 'red', 4.0 'red')
         unset colorbox
         set zrange [0:4]
         set cbrange [0:4]
         set title 'Gebruik Oppervlaktewaarnemingen rond $dtgoff UTC (groen: gebruikt; rood: verworpen).'
         $plotcom
EOF

      mv use_span$houroff.ps $todir

   fi

# Remove the intermediate text files

   rm SPAN$houroff.txt PSEUDO$houroff.txt

done

cd

rm -rf $wd

exit 0