#!/bin/sh

PATH=/usr/local/bin:/bin:/usr/bin:/opt/grads-2.0.a9/bin/:/home/toyoda/bin:/sbin:/usr/sbin:/usr/local/sbin

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

rm -f A_*.grib

areas="
 -90.0_0.0_150.0_240.0 -90.0_0.0_240.0_330.0
 -90.0_0.0_330.0_60.0  -90.0_0.0_60.0_150.0
 0.0_90.0_150.0_240.0 0.0_90.0_240.0_330.0
 0.0_90.0_330.0_60.0  0.0_90.0_60.0_150.0
"
levs="925hPa 70hPa"

for lev in $levs
do
for area in $areas
do
  : area="${area}"
  base="http://www.wis-jma.go.jp/d/o/RJTD/GRIB/Global_Spectral_Model/Latitude_Longitude/1.25_1.25/${area}/${lev}/${yy}${mm}${dd}/${hh}0000/"
  wget -q -O"zlist.html" $base
  ruby -ne 'puts "'"${base}"'/#{$1}" if /a href="(A_H[HRTUV].[ABCE][^"]+grib)"/' \
    zlist.html > zurl.txt
  cat zurl.txt 
  wget -i"zurl.txt"
  rm -f zurl.txt zlist.html
done
done
rm -f zurl.txt zlist.html
tar cvf grib.tar A*grib
rm -f A*grib