#!/bin/ksh

#------------------------------------------------------------------
# dust_grib_conus.sh
# Conversion of HYSPLIT Conus Dust Forecast Output to GRIB and ASCII
#------------------------------------------------------------------
# Last Revised:  02 Apr 2010 (RRD,GSM) - initial version
#------------------------------------------------------------------

set -xa

echo "start dust grib conversion: ${DATA}"

# dustsfc - output file for ground-level concentrations
# dustpbl - output file for column integrated concentrations

  concsfc="dustsfc"
  concpbl="dustpbl"
  gribpbl="grib_pbl.t${cyc}z.1hr"
  gribsfc="grib_sfc.t${cyc}z.1hr"

#----------------------------------------------------------------------
# create GRIB1 file from hysplit binary

  export pgm=hysplit_con2grib

  ${EXECdust}/${pgm} -s1 -g0 -z1 -r${PDY}${cyc} -i${DATA}/${concpbl} \
                     -o${DATA}/${gribpbl} -m6 -x1.0E+06

  ${EXECdust}/${pgm} -s1 -g0 -z1 -r${PDY}${cyc} -i${DATA}/${concsfc} \
                     -o${DATA}/${gribsfc} -m6 -x1.0E+06

#---------------------------------------------------------------------
# make output on 227 grid 

  echo ' &NLCOPYGB IDS(164)=1, IBS(164)=-1 /' > dust_scale
  ${EXECUTILhysplt}/smoke_copygb -N dust_scale -g "227" -x ${gribpbl} grib_pbl.t${cyc}z.1hr_227
  ${EXECUTILhysplt}/smoke_copygb -N dust_scale -g "227" -x ${gribsfc} grib_sfc.t${cyc}z.1hr_227

# GRIB1 to GRIB2 conversion
  ${EXECUTILhysplt}/cnvgrib -g12 -m0 -p32 grib_pbl.t${cyc}z.1hr_227 grib2_pbl.t${cyc}z.1hr_227
  ${EXECUTILhysplt}/cnvgrib -g12 -m0 -p32 grib_sfc.t${cyc}z.1hr_227 grib2_sfc.t${cyc}z.1hr_227
  ${EXECUTILhysplt}/cnvgrib -g12 -m0 -p32 grib_pbl.t${cyc}z.1hr grib2_pbl.t${cyc}z.1hr
  ${EXECUTILhysplt}/cnvgrib -g12 -m0 -p32 grib_sfc.t${cyc}z.1hr grib2_sfc.t${cyc}z.1hr

  cp grib*pbl* ${COMOUT}/.
  cp grib*sfc* ${COMOUT}/.
exit