#!/bin/ksh

###############################################################
# Runs SPC gempak from SREF pgrb212 for requested forecast hour
###############################################################

set -xa
export PS4='+ $SECONDS + '
date

export job=${job:-${PBS_JOBNAME:?}}
export jobid=${jobid:-$job.${PBS_JOBID:?}}

export DATA=${DATA:-${DATAROOT:?}/${jobid:?}}
# Remove old directory if it exists
rm -rf $DATA
mkdir -p $DATA
cd $DATA

export pgmout="OUTPUT.$$"
export pgmerr=errfile
export cycle=t${cyc}z

###############################
# Specify NET and RUN name
##############################
export NET=spcsref
export RUN=spcsref
export fend=00
export finc=03
export fstart=00
export GRIB=
export EXT=""

####################################################
# SENDCOM  - Copy Files From TMPDIR to $COMOUTgempak
####################################################
export SENDCOM=${SENDCOM:-YES}

export HOMEspc=${HOMEspc:-${PACKAGEROOT}/spcsref.${spcsref_ver}}
export EXECspc=${EXECspc:-$HOMEspc/exec}
export FIXspc=${FIXspc:-$HOMEspc/fix}
export PARMspc=${PARMspc:-$HOMEspc/parm}
export USHspc=${USHspc:-$HOMEspc/ush}
export GEMPAKspc=${GEMPAKspc:-$HOMEspc/gempak}

# Run setpdy and initialize PDY variables
setpdy.sh
. ./PDY

#
# Copy special table wmogrib2.tbl, ncepgrib2.tbl
# and vcrdgrib1.tbl into $DATA
#
cp $GEMPAKspc/fix/sref_wmogrib2.tbl wmogrib2.tbl
cp $GEMPAKspc/fix/sref_ncepgrib2.tbl ncepgrib2.tbl
cp $GEMPAKspc/fix/sref_vcrdgrib1.tbl vcrdgrib1.tbl
cp $GEMPAKspc/fix/g2varsncep1.tbl g2varsncep1.tbl
cp $GEMPAKspc/fix/g2varswmo2.tbl g2varswmo2.tbl
cp $GEMPAKspc/fix/g2vcrdncep1.tbl g2vcrdncep1.tbl
cp $GEMPAKspc/fix/g2vcrdwmo2.tbl g2vcrdwmo2.tbl

############################################################
# Determine COMIN and COMOUT (this job gets input from SREF)
############################################################
export COMINsref=${COMINsref:-$(compath.py ${envir}/sref/${sref_ver})/sref.${PDY}/${cyc}/pgrb}
export COMOUTgempak=${COMOUTgempak:-$(compath.py -o ${NET}/${spcsref_ver}/${RUN}.${PDY}/gempak)}
export COMOUT=${COMOUT:-$(compath.py -o ${NET}/${spcsref_ver}/${RUN}.${PDY})}

if [ ! -f $COMOUTgempak ] ; then
  mkdir -p -m 775 $COMOUTgempak
fi
 
env

filecount=$(ls $COMINsref/sref_*.${cycle}.pgrb212.*.f${gempak_time}.grib2.idx |wc -l)
if [ $filecount -lt 1 ]; then
   msg="FATAL ERROR: Aborting with no sref members for f${gempak_time}"
   err_exit $msg
elif [ $filecount -ne 26 ]; then
   msg="WARNING: $filecount members available for processing. (26 expected)"
   postmsg "$msg"
fi

########################################################
# Execute the script for each available member.
for file in $(ls $COMINsref/sref_*.${cycle}.pgrb212.*.f${gempak_time}.grib2.idx)
do
    export model=$(basename $file | cut -d. -f1 | cut -c6- )
    export ens_mem=$(basename $file | cut -d. -f4)
    for member in $ens_mem
    do
      export member

#     produce the sref gempak for this member
      $HOMEspc/scripts/exspcsref_pgrb212_gempak.sh

    done
done
########################################################
cat $pgmout

msg="JOB COMPLETED NORMALLY"
postmsg "$msg"

#######################################################
# Remove the Working Directory if KEEPDATA is not "YES"
#######################################################
if [ ${KEEPDATA:-NO} != YES ]; then
  [[ -n "$DATAROOT" ]] && cd $DATAROOT
  rm -rf $DATA
fi

date