#!/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
# 04/22/11: Jun Du, added an alternative option of using CCPA data
#
###################################################################
set -x
export XLFRTEOPTS="unit_vars=yes"

option=${option:-CCPA3}

vday=`echo $1 | cut -c 1-8`
vdaym1=`/nwprod/util/exec/ndate -24 ${vday}12 | cut -c 1-8`
vdayp1=`/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}

if [ $option = ST2 ]; then
# get Stage2 hourly data
 export OBSPPT_DIR=/com/hourly/prod/nam_pcpn_anal.
else
# get 3hrly CCPA data
 export OBSPPT_DIR=/com/gens/prod
fi

# Get Stage II data:
if [ $option = ST2 ]; then
# get Stage2 hourly accumulated precip data from /com directory 
# covering 21Z $vdaym1 - 21Z $vday
cp ${OBSPPT_DIR}${vdaym1}/ST2ml${vdaym1}22.Grb.Z .
cp ${OBSPPT_DIR}${vdaym1}/ST2ml${vdaym1}23.Grb.Z .
cp ${OBSPPT_DIR}${vday}/ST2ml${vday}0*.Grb.Z .
cp ${OBSPPT_DIR}${vday}/ST2ml${vday}1*.Grb.Z .
cp ${OBSPPT_DIR}${vday}/ST2ml${vday}20.Grb.Z .
cp ${OBSPPT_DIR}${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 
fi    # finished accumulating 3-hourly Stage II data.

# If use CCPA data:
grid=212
if [ $option = CCPA3 ]; then
 hr=0
 while [ $hr -le 21 ]; do
  if [ $hr -lt 10 ];then hr=0$hr;fi
# $CPGB -g${grid} -i3 -x ${OBSPPT_DIR}/gefs.$vday/$hr/ccpa/ccpa_conus_hrap_t${hr}z_03h ppt$hr

  if [ $hr -eq 00 ];then
   $CPGB -g${grid} -i3 -x ${OBSPPT_DIR}/gefs.$vday/00/ccpa/ccpa_conus_hrap_t${hr}z_03h ppt$hr
  fi
  if [ $hr -eq 03 -o $hr -eq 06 ];then
   $CPGB -g${grid} -i3 -x ${OBSPPT_DIR}/gefs.$vday/06/ccpa/ccpa_conus_hrap_t${hr}z_03h ppt$hr
  fi
  if [ $hr -eq 09 -o $hr -eq 12 ];then
   $CPGB -g${grid} -i3 -x ${OBSPPT_DIR}/gefs.$vday/12/ccpa/ccpa_conus_hrap_t${hr}z_03h ppt$hr
  fi
  if [ $hr -eq 15 -o $hr -eq 18 ];then
   $CPGB -g${grid} -i3 -x ${OBSPPT_DIR}/gefs.$vday/18/ccpa/ccpa_conus_hrap_t${hr}z_03h ppt$hr
  fi
  if [ $hr -eq 21 ];then
   $CPGB -g${grid} -i3 -x ${OBSPPT_DIR}/gefs.$vdayp1/00/ccpa/ccpa_conus_hrap_t${hr}z_03h ppt$hr
  fi

  hr=`expr $hr + 3`
 done
fi #finished getting CCPA data
 
ls -l ppt*
rm -f ST* input*

exit