exgfs_atmos_nawips.sh

#! /usr/bin/env bash

###################################################################
# echo "exnawips - convert NCEP GRIB files into GEMPAK Grids"
###################################################################

#### If EMC GFS PARA runs hourly file are not available, The ILPOST
#### will set to 3 hour in EMC GFS PARA.
#### Note:  ILPOST default set to 1
export ILPOST=${ILPOST:-1}

cd "${DATA}" || exit 1
grid=$1
fhr3=$2
DBN_ALERT_TYPE=$3
destination=$4

DATA_RUN="${DATA}/${grid}"
mkdir -p "${DATA_RUN}"
cd "${DATA_RUN}" || exit 1

# "Import" functions used in this script
source "${USHglobal}/product_functions.sh"

NAGRIB="${GEMEXE}/nagrib2"

cpyfil=gds
garea=dset
gbtbls=
maxgrd=4999
kxky=
grdarea=
proj=
output=T
pdsext=no

sleep_interval=10
max_tries=360

mkdir -p "lock.${fhr3}"
cd "lock.${fhr3}" || exit 1

for table in g2varswmo2.tbl g2vcrdwmo2.tbl g2varsncep1.tbl g2vcrdncep1.tbl; do
    source_table="${HOMEglobal}/gempak/fix/${table}"
    if [[ ! -f "${source_table}" ]]; then
        err_exit "FATAL ERROR: ${table} is missing"
    fi
    cpreq "${source_table}" "${table}"
done

GEMGRD="${RUN}_${grid}_${PDY}${cyc}f${fhr3}"

# Set type of Interpolation for WGRIB2
export opt1=' -set_grib_type same -new_grid_winds earth '
export opt1uv=' -set_grib_type same -new_grid_winds grid '
export opt21=' -new_grid_interpolation bilinear -if '
export opt22=":(CSNOW|CRAIN|CFRZR|CICEP|ICSEV):"
export opt23=' -new_grid_interpolation neighbor -fi '
export opt24=' -set_bitmap 1 -set_grib_max_bits 16 -if '
export opt25=":(APCP|ACPCP|PRATE|CPRAT):"
export opt26=' -set_grib_max_bits 25 -fi -if '
export opt27=":(APCP|ACPCP|PRATE|CPRAT|DZDT):"
export opt28=' -new_grid_interpolation budget -fi '

case ${grid} in
    # TODO: Why aren't we interpolating from the 0p25 grids for 35-km and 40-km?
    '0p50' | '0p25') grid_in=${grid} ;;
    *) grid_in="1p00" ;;
esac

source_var="COMIN_ATMOS_GRIB_${grid_in}"
export GRIBIN="${!source_var}/${RUN}.${cycle}.pres_a.${grid_in}.f${fhr3}.grib2"
GRIBIN_chk="${!source_var}/${RUN}.${cycle}.pres_a.${grid_in}.f${fhr3}.grib2.idx"

if ! wait_for_file "${GRIBIN_chk}" "${sleep_interval}" "${max_tries}"; then
    export err=7
    err_exit "After 1 hour of waiting for ${GRIBIN_chk} file at F${fhr3} to end."
fi

case "${grid}" in
    35km_pac) grid_spec='latlon 130.0:416:0.312 75.125:186:-0.312' ;;
    35km_atl) grid_spec='latlon 230.0:480:0.312 75.125:242:-0.312' ;;
    40km) grid_spec='lambert:265.0:25.0:25.0 226.541:185:40635.0 12.19:129:40635.0' ;;
    *) grid_spec='' ;;
esac

if [[ "${grid_spec}" != "" ]]; then
    # shellcheck disable=SC2086,SC2248
    "${WGRIB2}" "${GRIBIN}" ${opt1uv} ${opt21} ${opt22} ${opt23} ${opt24} ${opt25} ${opt26} ${opt27} ${opt28} -new_grid ${grid_spec} "grib${fhr3}"
    trim_rh "grib${fhr3}"
else
    cpreq "${GRIBIN}" "grib${fhr3}"
fi

export pgm="nagrib2 F${fhr3}"
startmsg

${NAGRIB} << EOF
GBFILE   = grib${fhr3}
INDXFL   =
GDOUTF   = ${GEMGRD}
PROJ     = ${proj}
GRDAREA  = ${grdarea}
KXKY     = ${kxky}
MAXGRD   = ${maxgrd}
CPYFIL   = ${cpyfil}
GAREA    = ${garea}
OUTPUT   = ${output}
GBTBLS   = ${gbtbls}
GBDIAG   =
PDSEXT   = ${pdsext}
l
r
EOF

export err=$?
if [[ ${err} -ne 0 ]]; then
    err_exit
fi

"${GEMEXE}/gpend"
export err=$?
if [[ ${err} -ne 0 ]]; then
    err_exit "${GEMEXE}/gpend failed!"
fi

cpfs "${GEMGRD}" "${destination}/${GEMGRD}"

if [[ ${SENDDBN} == "YES" ]]; then
    "${DBNROOT}/bin/dbn_alert" MODEL "${DBN_ALERT_TYPE}" "${job}" \
        "${destination}/${GEMGRD}"
fi
cd "${DATA_RUN}" || exit 1

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