gfs_bufr_netcdf.sh

#! /usr/bin/env bash

#
#  UTILITY SCRIPT NAME :  gfsbufr.sh
#               AUTHOR :  Hua-Lu Pan
#         DATE WRITTEN :  02/03/97
#
#  Abstract:  This utility script produces BUFR file of
#             station forecasts from the GFS suite.
#
#     Input:  none
# Script History Log:
# 2016-10-30  H Chuang: Tranistion to read nems output.
#             Change to read flux file fields in gfs_bufr
#             so remove excution of gfs_flux
# 2018-03-22 Guang Ping Lou: Making it works for either 1 hourly or 3 hourly output
# 2018-05-22 Guang Ping Lou: Making it work for both GFS and FV3GFS 
# 2018-05-30  Guang Ping Lou: Make sure all files are available.
# 2019-10-10  Guang Ping Lou: Read in NetCDF files
# echo "History: February 2003 - First implementation of this utility script"
#
source "${USHgfs}/preamble.sh"

if test "$F00FLAG" = "YES"
then
   f00flag=".true."
else
   f00flag=".false."
fi

hh=$FSTART
while  test $hh -le $FEND
do  
   hh=$( expr $hh + $FINT )
   if test $hh -lt 10
   then
      hh=0$hh
   fi
done

export pgm="gfs_bufr.x"
#. prep_step

if test "$MAKEBUFR" = "YES"
then
   bufrflag=".true."
else
   bufrflag=".false."
fi

SFCF="sfc"
CLASS="class1fv3"
cat << EOF > gfsparm
 &NAMMET
  levs=$LEVS,makebufr=$bufrflag,
  dird="$COMOUT/bufr.${cycle}/bufr",
  nstart=$FSTART,nend=$FEND,nint=$FINT,
  nend1=$NEND1,nint1=$NINT1,nint3=$NINT3,
  nsfc=80,f00=$f00flag,
/
EOF

hh=$FSTART
   if test $hh -lt 100
   then
      hh1=$(echo "${hh#"${hh%??}"}")
      hh=$hh1
   fi

sleep_interval=10
max_tries=360   
while  test $hh -le $FEND
do  
   if test $hh -lt 100
   then
      hh2=0$hh
   else
      hh2=$hh
   fi

   filename="${COMIN}/${RUN}.${cycle}.logf${hh2}.txt"
   if ! wait_for_file "${filename}" "${sleep_interval}" "${max_tries}" ; then
     err_exit "FATAL ERROR COULD NOT LOCATE logf${hh2} file"
   fi

#------------------------------------------------------------------
   ${NLN} $COMIN/${RUN}.${cycle}.atmf${hh2}.nc sigf${hh}
   ${NLN} $COMIN/${RUN}.${cycle}.${SFCF}f${hh2}.nc flxf${hh}

   hh=$( expr $hh + $FINT )
   if test $hh -lt 10
   then
      hh=0$hh
   fi
done  

#  define input BUFR table file.
${NLN} ${PARMgfs}/product/bufr_gfs_${CLASS}.tbl fort.1
${NLN} ${STNLIST:-${PARMgfs}/product/bufr_stalist.meteo.gfs} fort.8
${NLN} ${PARMgfs}/product/bufr_ij13km.txt fort.7

${APRUN_POSTSND} "${EXECgfs}/${pgm}" < gfsparm > "out_gfs_bufr_${FEND}"
export err=$?

exit ${err}