#!/bin/sh

set -x
##################################################
#  Fetch the Global Boundaries for the SREF
#  Programmer: Jun Du
#  Date:       04/25/2001
#  change log:
#   10/16/2003 Jun Du, $TOTALPAIR was replaced by 
#                $GLBPAIR to diff from regional ensemble
#                pair numbers $REGPAIR 
#   06/14/2005 Jun Du, 3hr-old global ensemble forecasts are
#                used as LBCs instead of previously 9hr-old ones
#                since the global ensemble is now available 4 cycles 
#                per day 
#   03/22/2006 Jun Du, add 3z and 15z two off-time cycles
#   04/04/2006 Jun Du, modified to match the new glb ens data output format
#   06/25/2006 HUI-YA CHUANG, modified to run 32 km ETA MOS
#
##################################################

HR=$1
INDEX=`expr $HR / $INCHOUR`

JOBDATA=$DATA
export DATA=$JOBDATA/global.$INDEX
mkdir $DATA
cd $DATA

$utilscript/setup.sh

#--------------------------------------------------------------
# global gfs and ensemble forecasts as LBCs
#--------------------------------------------------------------
############################################
# ctl (gfs with current resolution):
############################################
HR=$1
INDEX=`expr $HR / $INCHOUR`

PAIR=ctl

# ETA MOS IS RUNNING ON TOP OF NAM AND BEFORE GFS SO NEED TO USE GFS FROM PREVIOUS CYCLE
TSHIFT=`expr $HR + 6`
if [ $TSHIFT -lt 10 ]
then
   TSHIFT=0$TSHIFT
fi

if [ $cyc -eq 00 ]
then
   cp $COMGFS/gfs.$PDYm1/gfs.$CYC4.sf${TSHIFT} $JOBDATA/bc_gfs.$INDEX.old
fi
if [ $cyc -eq 06 ]
then
   cp $COMGFS/gfs.$PDY/gfs.$CYC1.sf$TSHIFT $JOBDATA/bc_gfs.$INDEX.old
fi
if [ $cyc -eq 12 ]
then
   cp $COMGFS/gfs.$PDY/gfs.$CYC3.sf${TSHIFT} $JOBDATA/bc_gfs.$INDEX.old
fi
if [ $cyc -eq 18 ]
then
   cp $COMGFS/gfs.$PDY/gfs.$CYC2.sf$TSHIFT $JOBDATA/bc_gfs.$INDEX.old
fi

refTime=`/nwprod/util/exec/ndate +$HR $PDY$cyc`
$chtime $refTime $JOBDATA/bc_gfs.$INDEX.old $JOBDATA/bc_gfs.$INDEX
rm -f $JOBDATA/bc_gfs.$INDEX.old

# truncating current global gfs resolution to T126L28
export SIGINP=$JOBDATA/bc_gfs.$INDEX
export SFCINP=NULL
export SIGOUT=$JOBDATA/bc_$cycle.$PAIR.$INDEX
export SFCOUT=sfcout
$chres

cp $JOBDATA/bc_$cycle.$PAIR.$INDEX $GESDIR/.
rm -f $JOBDATA/bc_gfs.$INDEX


cd $JOBDATA
rm -rf $DATA

exit