#
######################################################
echo "-----------------------------------------------"
echo "     exspcsref_xhail_gempak.sh.sms             "
echo "   SPC SREF XHAIL KF eta model and ct12        "
echo "-----------------------------------------------"
######################################################

#-----------------------------------------------------
#     This script is for KF eta model and ctl member,
#     do the eta experimental post processing 
#     including the hail model output
#-----------------------------------------------------

set -x
cd $DATA
  
cp $FIXspc/spc_fields_namX fields.etaX
cp $FIXspc/spc_nlist.txt nlist.txt
cp $FIXspc/spc_hail_model_flag_gem hail_model_flag_gem.out
chmod 775 *
#
# Check the next hour, current hour, and prev are available...otherwise wait for it.
#
export hour=$fcst_hr
export next_hour=$fcsthrp3
export prev_hour=$prev_hr

#
if [ ! -f ${COMIN}/${PDY}_sref_em.${cycle}.pgrb212.ctl.f${next_hour}.gem ]
then
    err=1;export err; err_chk
else
    cp ${COMIN}/${PDY}_sref_em.${cycle}.pgrb212.ctl.f${next_hour}.gem .
fi

echo "Next hour avbl...move on"
echo "Next hour avbl...move on" >> ./status.report
#

if [ ! -f ${COMIN}/${PDY}_sref_em.${cycle}.pgrb212.ctl.f${prev_hour}.gem ]
then
    err=1;export err; err_chk
else
    cp ${COMIN}/${PDY}_sref_em.${cycle}.pgrb212.ctl.f${prev_hour}.gem .
fi

echo "Prev hour avbl...move on"
echo "Prev hour avbl...move on" >> ./status.report
#

if [ ! -f ${COMIN}/${PDY}_sref_em.${cycle}.pgrb212.ctl.f${hour}.gem ]
then
    err=1;export err; err_chk
else
    cp ${COMIN}/${PDY}_sref_em.${cycle}.pgrb212.ctl.f${hour}.gem .
    echo "${PDY}_sref_em.${cycle}.pgrb212.ctl.f${hour}.gem" > cape_wx_model.input
fi

echo "Cur hour avbl...move on"
echo "Cur hour avbl...move on" >> ./status.report
#

# interpolate zagl levels now...
$GEMEXE/gdvint << EOF

 GDFILE   = ./${PDY}_sref_em.${cycle}.pgrb212.ctl.f${hour}.gem 
 GDOUTF   = ./${PDY}_sref_em.${cycle}.pgrb212.ctl.f${hour}.gem 
 GDATTIM  = LAST
 GVCORD   = pres/zagl
 GLEVEL   = 500-8000-500
 MAXGRD   = 5000
 GAREA    = grid
 VCOORD   = mslv;esfc

r

EOF

#$GEMEXE/gpend

${EXECspc}/spc_get_etaX_data_scape_sars_mcs full

if [ "$SENDCOM" = 'YES' ]
then
   cp ${PDY}_sref_em.${cycle}.pgrb212.ctl.f${hour}.gem ${COMOUT}/namkf40_${PDY}${cyc}f0${hour}
fi

if [ "$SENDDBN" = 'YES' ]
then
         $DBNROOT/bin/dbn_alert MODEL SREF_GEMPAK SPC ${COMOUT}/namkf40_${PDY}${cyc}f0${hour}
fi

msg="$job completed normally"
postmsg "$jlogfile" "$msg"

# End of SPC SREF XHAIL GEMPAK script -------------------------------------- #