#!/bin/ksh

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

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="binak_sfc.1hr"
  concpbl="binak_pbl.1hr"
  gribpbl="gribak_pbl.1hr"
  gribsfc="gribak_sfc.1hr"

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

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

  ${EXEChysplt}/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 \
     -i${DATA}/${concpbl} -o${DATA}/${gribpbl} -m4 -n25:72 -x1.0E+09 -r${RUNDATE}${cyc} \
      >> $pgmout 2>errfile
#     >> ${DATA}/con2gribpbl.out

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

exit