#!/bin/sh

PATH=/bin:/usr/bin:/usr/local/bin:/opt/grads-2.0.a9/bin
cd `dirname $0`

if [[ ! -t 2 ]]; then
  exec > zsurf.log 2>&1
fi

: ${ofs:=4000}
: ${verbose:=}

export TZ=UTC+0

bt=`ruby -e "puts Time.at((Time.now.utc - $ofs).to_i / 21600 * 21600
  ).strftime('%Y %m %d %H %a')"`
set -- $bt
bthh=$4
: ${hh:=${bthh}}
if [[ $bthh < $hh ]]; then
  bt=`ruby -e "puts (Time.gm(*ARGV) - 86400).strftime('%Y %m %d %H %a')" $bt`
  set -- $bt
fi
yy=$1
mm=$2
dd=$3
wd=$5

zip=/var/www/html/wisdata/arch/bfsf/${yy}${mm}/bfsf${yy}${mm}${dd}.zip
if [ ! -f $zip ]; then
  echo $zip not found
  exit 1
fi
export BUFR_TABLES=/usr/share/libecbufr

dhm="${dd}${hh}00"

unzip -p $zip \
  "A_IS?[ADNX]??E[BGHI]*${dhm}*_C_*" \
  "A_IS?[ADNX]??LF??*${dhm}*_C_*" \
  > zsynop3.bufr
bufr_decoder -dump -inbufr zsynop3.bufr | ruby bufr2synop.rb > zsynop3.txt
if [ -d ../fakesynop ]; then
  cp zsynop3.txt ../fakesynop/synop${wd}${hh}.txt
fi

ruby gradswx.rb -m:40..70,-30..20 \
  -y${dhm} -ym${yy}-${mm} -f"synop-sa.tsv" zsynop3.txt
stnmap -i zgrads.ctl > /dev/null
grads -b -l -c 'run gradswx-sa.gs'
convert zgradssa.png -crop 1915x1400+170+230 zplotsa${hh}.png
cp zplotsa${hh}.png zgrads.tmp
mv -f zgrads.tmp zgradssa.png
convert zgradssa.png zgradssa.jpg