#!/bin/ksh -x
#
###################################################################
# This script gets 3-hourly Stage II precipitation analysis 
# i.e  00/03/06/09/12/15/18/21Z for a day (21z-21z)
# 
# Change Log:
# 10/01/10: Jun Du, initial script for precip bias correction job
#
###################################################################
set -x
export XLFRTEOPTS="unit_vars=yes"

vday=`echo $1 | cut -c 1-8`
vdaym1=`/nwprod/util/exec/ndate -24 ${vday}12 | cut -c 1-8`

YYYY=`echo $vday | cut -c 1-4`
MM=`echo $vday   | cut -c 5-6`

export CPGB=${CPGB:-/nwprod/util/exec/copygb}
export EXECnam=${EXECnam:-/nwprod/exec}

export ftpdir=/com/hourly/prod/nam_pcpn_anal.
# get Stage2 hourly accumulated precip data from /com directory 
# covering 21Z $vdaym1 - 21Z $vday
cp ${ftpdir}${vdaym1}/ST2ml${vdaym1}22.Grb.Z .
cp ${ftpdir}${vdaym1}/ST2ml${vdaym1}23.Grb.Z .
cp ${ftpdir}${vday}/ST2ml${vday}0*.Grb.Z .
cp ${ftpdir}${vday}/ST2ml${vday}1*.Grb.Z .
cp ${ftpdir}${vday}/ST2ml${vday}20.Grb.Z .
cp ${ftpdir}${vday}/ST2ml${vday}21.Grb.Z .

uncompress ST2ml*.Grb.Z

# process hourly data into 3hrly accumulated precip
hr=0
accdate=${vday}00
while [ $accdate -le ${vday}21 ]; do
  if [ $hr -lt 10 ];then hr=0$hr;fi
  time1=`/nwprod/util/exec/ndate -2 $accdate`  
  time2=`/nwprod/util/exec/ndate -1 $accdate`  
  time3=$accdate
  cat > input_acc3h << EOF
obs
ST2.
ST2ml$time1.Grb
ST2ml$time2.Grb
ST2ml$time3.Grb
EOF

$EXECnam/nam_stage4_acc < input_acc3h
export err=$?; err_chk

# convert to grid 212
#  for grid in 212 218 
 for grid in 212 
 do
   $CPGB -g${grid} -i3 -x ST2.$accdate.03h ppt$hr
 done

 accdate=`/nwprod/util/exec/ndate +3 $accdate`  
 hr=`expr $hr + 3`
done  # finished accumulating 3-hourly Stage II data.
 
ls -l ppt*
rm -f ST* input*

exit