#!/bin/ksh

#------------------------------------------------------------------
# smoke_grib.sh
# Conversion of HYSPLIT Smoke Forecast Output to GRIB and ASCII
#------------------------------------------------------------------
# Last Revised: 29 Mar 2007 (RRD) - parallel version
#------------------------------------------------------------------

set -xa

  RUNDATE=$1

echo "start grib conversion: ${DATA}"

# concsfc - output file for ground-level concentrations
# concpbl - output file for column integrated concentrations
# concgrb - concpbl outfile file converted to GRIB

  concsfc="bin_sfc.1hr"
  concpbl="bin_pbl.1hr"
  gribpbl="grib_pbl.1hr"
  gribsfc="grib_sfc.1hr"

#----------------------------------------------------------------------
  echo "Create GRIB and ASCII files"       

# ascii file transferred to web to be picked up by EPA BlueRains
  export pgm=hysplit_con2asc  

  ${EXECsmoke}/hysplit_con2asc -s -i${DATA}/${concpbl} -onoaasmk_${PDY} \ 
              >> $pgmout 2>errfile

#----------------------------------------------------------------------
# create grib1 file from hysplit binary
  export pgm=hysplit_con2grib 

  ${EXECsmoke}/hysplit_con2grib -s1 -g0 -z1 -r${RUNDATE}${cyc} \
     -i${DATA}/${concpbl} -o${DATA}/${gribpbl} -m4 -n25:72 -x1.0E+09 \
      >> $pgmout 2>errfile
#     >> ${DATA}/con2gribpbl.out

  ${EXECsmoke}/hysplit_con2grib -s1 -g0 -z1 -r${RUNDATE}${cyc} \
     -i${DATA}/${concsfc} -o${DATA}/${gribsfc} -m4 -n25:72 -x1.0E+09 \
      >> $pgmout 2>errfile
#     >> ${DATA}/con2gribsfc.out


#---------------------------------------------------------------------
# GRIB1 to GRIB2 conversion

  cd ${DATA}  

  ${EXECUTILsmoke}/smoke_copygb -g "227" -x ${gribpbl} grib_pbl.1hr_227
  ${EXECUTILsmoke}/smoke_copygb -g "227" -x ${gribsfc} grib_sfc.1hr_227

  ${EXECUTILsmoke}/cnvgrib -g12 -m -p32 grib_pbl.1hr_227 grib2_pbl.1hr
  ${EXECUTILsmoke}/cnvgrib -g12 -m -p32 grib_sfc.1hr_227 grib2_sfc.1hr

exit