#!/bin/sh

set -x
##################################################
#  Create Nam LBC and Analysis from Global Guess
#  Programmer: Jun Du
#  Date:       04/25/2001
#  change log:
#  10/17/2003: Jun Du,  (1) resolution change from 48kmL45
#                32kmL60; (2) included 4 Ferrier's new cloud 
#                fields
#  06/14/2005: Jun Du, extended 63hr LBC to 87hr
#  06/25/2006  HUI-YA CHUANG, modified to run 32 km ETA MOS
#
##################################################

PAIR=$1

JOBDATA=$DATA
export DATA=$JOBDATA/convglb2eta.$PAIR
mkdir $DATA
cd $DATA

$utilscript/setup.sh

#-----------------------------------------
# Transform global analyses to eta format
#-----------------------------------------
############################################
# PAIRS ARE:
# ctl      - control
# p(1-2)   - positively perturbed ones
# n(1-2)   - negatively perturbed ones
############################################

#------------------------------------------
# process lateral boundaries
#------------------------------------------
#----------------------------------------------------------
# READ THE FORECAST FIELDS FROM THE NCEP
# GLOBAL MODEL, TRANSFORM THEM ONTO AN INTERMEDIATE
# LATITUDE/LONGITUDE GRID, AND EXTRACT LATERAL BOUNDARY
# VALUES FOR THE MESOSCALE ETA MODEL FORECAST
#
# The program is independent to global model's horizontal
# resolution but dependent to vertical resolution
#----------------------------------------------------------
export pgm=sref_eta_mkbnd
. prep_step
export XLFUNIT_11=$JOBDATA/bc_$cycle.$PAIR.0
export XLFUNIT_12=$JOBDATA/bc_$cycle.$PAIR.1
export XLFUNIT_13=$JOBDATA/bc_$cycle.$PAIR.2
export XLFUNIT_14=$JOBDATA/bc_$cycle.$PAIR.3
export XLFUNIT_15=$JOBDATA/bc_$cycle.$PAIR.4
export XLFUNIT_16=$JOBDATA/bc_$cycle.$PAIR.5
export XLFUNIT_17=$JOBDATA/bc_$cycle.$PAIR.6
export XLFUNIT_18=$JOBDATA/bc_$cycle.$PAIR.7
export XLFUNIT_19=$JOBDATA/bc_$cycle.$PAIR.8
export XLFUNIT_20=$JOBDATA/bc_$cycle.$PAIR.9
export XLFUNIT_21=$JOBDATA/bc_$cycle.$PAIR.10
export XLFUNIT_22=$JOBDATA/bc_$cycle.$PAIR.11
export XLFUNIT_23=$JOBDATA/bc_$cycle.$PAIR.12
export XLFUNIT_24=$JOBDATA/bc_$cycle.$PAIR.13
export XLFUNIT_25=$JOBDATA/bc_$cycle.$PAIR.14
export XLFUNIT_26=$JOBDATA/bc_$cycle.$PAIR.15
export XLFUNIT_27=$JOBDATA/bc_$cycle.$PAIR.16
export XLFUNIT_28=$JOBDATA/bc_$cycle.$PAIR.17
export XLFUNIT_29=$JOBDATA/bc_$cycle.$PAIR.18
export XLFUNIT_30=$JOBDATA/bc_$cycle.$PAIR.19
export XLFUNIT_31=$JOBDATA/bc_$cycle.$PAIR.20
export XLFUNIT_32=$JOBDATA/bc_$cycle.$PAIR.21
export XLFUNIT_33=$JOBDATA/bc_$cycle.$PAIR.22
export XLFUNIT_34=$JOBDATA/bc_$cycle.$PAIR.23
export XLFUNIT_35=$JOBDATA/bc_$cycle.$PAIR.24
export XLFUNIT_36=$JOBDATA/bc_$cycle.$PAIR.25
export XLFUNIT_37=$JOBDATA/bc_$cycle.$PAIR.26
export XLFUNIT_38=$JOBDATA/bc_$cycle.$PAIR.27
export XLFUNIT_39=$JOBDATA/bc_$cycle.$PAIR.28
export XLFUNIT_40=$JOBDATA/bc_$cycle.$PAIR.29
export XLFUNIT_41=$PARMeta/sref_modtop.parm
export XLFUNIT_42=$FIXeta/nam_deta_ldt1.60.25mb

export XLFUNIT_51=rvetalbc
export XLFUNIT_52=etabcs00

startmsg
$EXECsref/sref_eta_mkbnd < $PARMsref/eta32k_mkbnd.parm_3h.f84 >> $pgmout 2> errfile
export err=$?;err_chk

#--------------------------------------------------------
# COMPUTE LATERAL BOUNDARY TENDENCIES FOR THE
# MESOSCALE ETA MODEL FORECAST
#--------------------------------------------------------
export pgm=sref_eta_bctend_3h
. prep_step
export XLFUNIT_17=etabcs00
#export XLFUNIT_18=$JOBDATA/nbc${RES}_$PAIR
export XLFUNIT_18=$JOBDATA/ndas.$cycle.nbcout

export XLFUNIT_51=$JOBDATA/nbcout.${cycle}_$PAIR

startmsg
$EXECsref/sref_eta_bctend_3h >> $pgmout 2> errfile
export err=$?; err_chk

cd $JOBDATA
#rm -rf $DATA

exit