exgfs_wave_nawips.sh

#! /usr/bin/env bash

###################################################################
# echo "----------------------------------------------------"
# echo "exnawips - convert NCEP GRIB files into GEMPAK Grids"
# echo "----------------------------------------------------"
# echo "History: Mar 2000 - First implementation of this new script."
# echo "Sept 2011 - First implementation of this new script based on"
# echo "               /nwprod/scripts/exnawips.sh.sms"
# echo " March 2020- Modified for GEFSv12.0"
#  March-2020 Roberto.Padilla@noaa.gov                                   
#####################################################################

source "${USHgfs}/preamble.sh"
source "${USHgfs}/wave_domain_grid.sh"

#export grids=${GEMPAK_GRIDS:-'glo_30m at_10m ep_10m wc_10m ao_9km'} #Interpolated grids
export grids=${GEMPAK_GRIDS:-${waveinterpGRD:-'glo_30m'}}  #Native grids
export RUNwave=${RUNwave:-${RUN}wave}
export fstart=${fstart:-0}
export FHMAX_WAV=${FHMAX_WAV:-180}  #180 Total of hours to process
export FHMAX_HF_WAV=${FHMAX_HF_WAV:-72}
export FHOUT_WAV=${FHOUT_WAV:-6}
export FHOUT_HF_WAV=${FHOUT_HF_WAV:-3}
export maxtries=${maxtries:-720}
export cycle=${cycle:-t${cyc}z}
export GEMwave=${GEMwave:-${HOMEgfs}/gempak}
export DATA=${DATA:-${DATAROOT:?}/${jobid}}
if [ ! -d ${DATA} ];then
  mkdir -p ${DATA}
fi

cd ${DATA}
cp ${GEMwave}/fix/g2varswmo2.tbl .

cpyfil=gds
garea=dset
gbtbls=
maxgrd=4999
kxky=
grdarea=
proj=
output=T
pdsext=no
g2tbls=g2varswmo2.tbl
NAGRIB=nagrib2

sleep_interval=20
maxtries=15
fhcnt=${fstart}
while [ ${fhcnt} -le ${FHMAX_WAV} ]; do
  fhr=$(printf "%03d" ${fhcnt})
  for grid in ${grids};do
    case ${grid} in
      ao_9km)  grdIDin='arctic.9km'
               #grdIDout='gfswaveao9km' ;;
               grdIDout='gfswavearc' ;;
      at_10m)  grdIDin='atlocn.0p16'
               grdIDout='gfswaveat10m' ;;
      ep_10m)  grdIDin='epacif.0p16'
               grdIDout='gfswaveep10m' ;;
      wc_10m)  grdIDin='wcoast.0p16'
               grdIDout='gfswavewc10m' ;;
      glo_30m) grdIDin='global.0p25'
               grdIDout='gfswavegl30m' ;;
      glo_10m) grdIDin='global.0p16'   
               #grdIDout='gfswaveg16k' ;;
               grdIDout='gfswavenh' ;;
      gso_15m) grdIDin='gsouth.0p25' 
               #grdIDout='gfswaves25k' ;;
               grdIDout='gfswavesh' ;;
      glo_200) grdIDin='global.2p00'
               grdIDout='gfswaves200k' ;;
      *)       grdIDin= 
               grdIDout= ;;
    esac
    process_grdID "${grid}"
    com_varname="COMIN_WAVE_GRID_${GRDREGION}_${GRDRES}"
    com_dir=${!com_varname}
    GRIBIN="${com_dir}/${RUNwave}.${cycle}.${grdIDin}.f${fhr}.grib2"
    GRIBIN_chk=${GRIBIN}.idx
    if ! wait_for_file "${GRIBIN_chk}" "${sleep_interval}" "${maxtries}"; then
      echo "FATAL ERROR: ${GRIBIN_chk} not found after waiting $((sleep_interval * ( maxtries - 1))) secs"
      echo "${RUNwave} ${grdID} ${fhr} prdgen ${date} ${cycle} : GRIB file missing." >> "${wavelog}"
      err=1;export err;"${errchk}" || exit "${err}"
    fi

    #if [ "$grdIDin" = "global.0p25" && "$grid" = "glo_30m" ]; then
    if [ "${grdIDin}" = "global.0p25" ]; then
      ${WGRIB2} -lola 0:720:0.5 -90:361:0.5 gribfile.${grdIDout}.f${fhr}  grib \
                                          ${GRIBIN} 1> out 2>&1
      OK=$?
      if [ "${OK}" != '0' ]; then
        msg="ABNORMAL EXIT: ERROR IN interpolation the global grid"
        #set +x
        echo ' '
        echo '************************************************************* '
        echo '*** FATAL ERROR : ERROR IN making  gribfile.$grdID.f${fhr}*** '
        echo '************************************************************* '
        echo ' '
        echo ${msg}
        #set_trace
        echo "${RUNwave} ${grdID} prdgen ${date} ${cycle} : error in grbindex." >> ${wavelog}
        err=2;export err;err_chk
      else
        #cp $GRIBIN gribfile.$grdID.f${fhr}
        GRIBIN=gribfile.${grdIDout}.f${fhr}
      fi
    fi
    echo ${GRIBIN}

    GEMGRD=${grdIDout}_${PDY}${cyc}f${fhr}

    cp ${GRIBIN} grib_${grid}

    startmsg

	${NAGRIB} <<-EOF
	GBFILE   = grib_${grid}
	INDXFL   = 
	GDOUTF   = ${GEMGRD}
	PROJ     = ${proj}
	GRDAREA  = ${grdarea}
	KXKY     = ${kxky}
	MAXGRD   = ${maxgrd}
	CPYFIL   = ${cpyfil}
	GAREA    = ${garea}
	OUTPUT   = ${output}
	GBTBLS   = ${gbtbls}
	G2TBLS   = ${g2tbls}
	GBDIAG   = 
	PDSEXT   = ${pdsext}
	l
	r
	EOF
    export err=$?;pgm=${NAGRIB};err_chk
    #####################################################
    # GEMPAK DOES NOT ALWAYS HAVE A NON ZERO RETURN CODE
    # WHEN IT CAN NOT PRODUCE THE DESIRED GRID.  CHECK
    # FOR THIS CASE HERE.
    #####################################################
    ls -l ${GEMGRD}
    export err=$?;export pgm="GEMPAK CHECK FILE";err_chk

    if [ "${NAGRIB}" = "nagrib2" ] ; then
      gpend
    fi

    cpfs "${GEMGRD}" "${COMOUT_WAVE_GEMPAK}/${GEMGRD}"
    if [ ${SENDDBN} = "YES" ] ; then
        "${DBNROOT}/bin/dbn_alert" MODEL "${DBN_ALERT_TYPE}" "${job}" "${COMOUT_WAVE_GEMPAK}/${GEMGRD}"
    else
        echo "##### DBN_ALERT is: MODEL ${DBN_ALERT_TYPE} ${job} ${COMOUT_WAVE_GEMPAK}/${GEMGRD}#####"
    fi
    rm grib_${grid}
  done
  if [ ${fhcnt} -ge ${FHMAX_HF_WAV} ]; then
    inc=${FHOUT_WAV}
  else
    inc=${FHOUT_HF_WAV}
  fi
  let fhcnt=fhcnt+inc
done
#####################################################################


############################### END OF SCRIPT #######################