#!/bin/ksh
################################################################################
####  UNIX Script Documentation Block
#                      .                                             .
# Script name:         exdgex_real.sh.sms
# Script description:  Runs REAL to get 78-h NAM forecast in WRF-NMM input file
                       format to start WRF-NNM DGEX
#
# Author:        Eric Rogers       Org: NP22         Date: 2006-04-21
#

set -x

mod=ETA

export MP_BINDPROC=NO

# NOTE : for prod DOMNAM, CYC, DATA will be set in J-job

msg="JOB $job FOR WRF-NMM HAS BEGUN"
postmsg "$jlogfile" "$msg"

#
# Get needed variables from exdgex_prelim.sh.sms
#
. $COMOUT/${RUN}_${REGION}.t${cyc}z.envir.sh

cd $DATA

cp /com/date/t${cyc}z .

CYCer=`cat t${cyc}z | cut -c 7-16`
CYCstart=`/nwprod/util/exec/ndate -0 $CYCer`
CYCdgex=`/nwprod/util/exec/ndate 78 $CYCer`

echo DATEXX${CYCstart} > nmcdate
echo DATEXX${CYCdgex} > nmcdate.dgex

PDY=`cat nmcdate | cut -c7-14`
yy=`cat nmcdate | cut -c7-10`
mm=`cat nmcdate | cut -c11-12`
dd=`cat nmcdate | cut -c13-14`
cyc=`cat nmcdate | cut -c15-16`

PDYnam=`cat nmcdate | cut -c7-14`
yynam=`cat nmcdate | cut -c7-10`
mmnam=`cat nmcdate | cut -c11-12`
ddnam=`cat nmcdate | cut -c13-14`
cycnam=`cat nmcdate | cut -c15-16`

PDYdgex=`cat nmcdate.dgex | cut -c7-14`
yydgex=`cat nmcdate.dgex | cut -c7-10`
mmdgex=`cat nmcdate.dgex | cut -c11-12`
dddgex=`cat nmcdate.dgex | cut -c13-14`
cycdgex=`cat nmcdate.dgex | cut -c15-16`

### modify namelist file
ystart=`cat nmcdate.dgex | cut -c7-10`
mstart=`cat nmcdate.dgex | cut -c11-12`
dstart=`cat nmcdate.dgex | cut -c13-14`
hstart=`cat nmcdate.dgex | cut -c15-16`
cycstart=`cat nmcdate.dgex | cut -c7-16`
PDYdgex=`cat nmcdate.dgex | cut -c7-14`

start=$ystart$mstart$dstart

end=`/nwprod/util/exec/ndate 3 $cycstart`

yend=`echo $end | cut -c1-4`
mend=`echo $end | cut -c5-6`
dend=`echo $end | cut -c7-8`
hend=`echo $end | cut -c9-10`

cp $PARMdgex/dgex_${REGION}_wrfnmm_namelist.input.in namelist.input_in

cat namelist.input_in | sed s:YSTART:$ystart: | sed s:MSTART:$mstart: \
 | sed s:DSTART:$dstart: | sed s:HSTART:$hstart: | sed s:YEND:$yend: \
 | sed s:MEND:$mend:     | sed s:DEND:$dend: | sed s:HEND:$hend: > namelist.input

cp $DATA/nmcdate.dgex t${cyc}z 

if [ $? -ne 0 ] 
then
exit
fi

rm fort.*

cp $FIXdgex/dgex_LANDUSE.TBL LANDUSE.TBL
cp $PARMdgex/dgex_ETAMPNEW_DATA ETAMPNEW_DATA
cp $FIXdgex/dgex_VEGPARM.TBL VEGPARM.TBL
cp $FIXdgex/dgex_SOILPARM.TBL SOILPARM.TBL
cp $FIXdgex/dgex_GENPARM.TBL GENPARM.TBL
cp $FIXdgex/dgex_co2.60_hybrid_30mb co2.60_hybrid_30mbtop
cp $FIXdgex/dgex_co2.60_hybrid_30mb co2_trans

cp $FIXdgex/hiresw_RRTM_DATA RRTM_DATA
cp $FIXdgex/hiresw_tr49t67 tr49t67
cp $FIXdgex/hiresw_tr49t85 tr49t85
cp $FIXdgex/hiresw_tr67t85 tr67t85
cp $PARMdgex/nam_micro_lookup.dat eta_micro_lookup.dat

ln -sf co2.60_hybrid_30mbtop fort.14

cp $PARMdgex/nam_fcst_gribmap.txt gribmap.txt

export pgm=dgex_nmm_real
. prep_step

export XLFUNIT_14=co2.60_hybrid_30mbtop

startmsg
$EXECdgex/dgex_nmm_real >> $pgmout 2>errfile
export err=$?;err_chk

cp wrfinput_d01 $COMOUT/dgex_${REGION}.t${cyc}z.wrfinput_d01
cp wrfbdy_d01 $COMOUT/dgex_${REGION}.t${cycnam}z.wrfbdy_d01_fromreal

# write out 78-h NAM boundary made by real into a binary file that can be read by COMBC 

mv wrfbdy_d01 wrfbdy_d01_fromreal
echo "nothing" > gribmap.txt

export pgm=dgex_extract_wrfbdy_${REGION}
. prep_step

startmsg
$EXECdgex/dgex_extract_wrfbdy_${REGION} >> $pgmout 2> errfile
export err=$?;err_chk

mv wrfbdy_extract.bin $COMOUT/dgex_${REGION}.t${cycnam}z.wrfbdy_nam78h.bin

exit