#!/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:=3600}
: ${verbose:=}

export TZ=UTC+0

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

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

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

unzip -c $zip \
  "A_S?????RJTD*${dhm}*_C_*.txt" \
  "A_S?PA??KWBC*${dhm}*_C_*.txt" \
  "A_S?KO??RKSL*${dhm}*_C_*.txt" \
  "A_S?KR*${dhm}*_C_*.txt" \
  "A_S?PH*${dhm}*_C_*.txt" \
  "A_S?RA??"{RUHB,KWBC}"*${dhm}*_C_*.txt" \
  "A_S?MO*${dhm}*_C_*.txt" \
  "A_S?HK*${dhm}*_C_*.txt" \
  "A_S?CI*${dhm}*_C_*.txt" \
  "A_S?VB*${dhm}*_C_*.txt" \
  > zsynop.txt

unzip -c $zip \
  "A_S?RA??RUNW*${dhm}*_C_*.txt" \
  "A_S?NP*${dhm}*_C_*.txt" \
  "A_SMIN02DEMS*${dhm}*_C_*.txt" \
  "A_SIIN21DEMS*${dhm}*_C_*.txt" \
  "A_S?IN04KWBC*${dhm}*_C_*.txt" \
  "A_S?BW*${dhm}*_C_*.txt" \
  "A_S?BM*${dhm}*_C_*.txt" \
  "A_S?TH*${dhm}*_C_*.txt" \
  "A_S?KP*${dhm}*_C_*.txt" \
  "A_S?VS*${dhm}*_C_*.txt" \
  "A_S?LT*${dhm}*_C_*.txt" \
  "A_S?AK*${dhm}*_C_*.txt" \
  > zsynop2.txt

sh getgrib.sh $yy $mm $dd $hh

ruby gradswx.rb -f"synop-japan.tsv" -f"synop-fe.tsv" -m:0..55,85..180 \
  -y${dhm} -ym${yy}-${mm} zsynop.txt zsynop2.txt zbogus.txt
stnmap -i zgrads.ctl > /dev/null
grads -b -l -c 'run gradswx-fe.gs'
convert zgrads.png -crop 1915x1400+170+230 zplotfe${hh}.png
cp zplotfe${hh}.png zgrads.tmp
mv -f zgrads.tmp zgrads.png
convert zgrads.png zgrads.jpg


ruby ${verbose} jpreport.rb -f"synop-japan.tsv" -y${dhm} \
	zsynop.txt zbogus.txt > zreport.txt
ruby txt2htm.rb zreport.txt > zreport.xhtml