#!/bin/sh
export member=airplt_gempak.sh
#
######################################################################     
#
# NAME     : airplt_gempak.sh
# Purpose  : This script creates 4 air pollution gifs using GEMPAK tools
#          : 
# History  : 02-15-2007   Larry Sager    
#            03-09-2009 C. Magee Change from ngm to hrly
#
######################################################################     
#
cd $DATA 

#
#  Convert air pollution tables into gempak formated files. 
#

set -x
cp $FIXgraph/hrlyairp_table5*.pack .
cp $FIXgraph/sffdsf* .
cp $FIXgempak/coltbl.xwp.wbg coltbl.xwp

cat <<\gemEOF >gem_namlst
 &INPARM npgiv=6/
gemEOF

export XLFUNIT_9="gem_namlst"

$EXECgraph/sffdsf << EOF > /dev/null     
 SFEFIL   = table54.dat
 SFOUTF   = table54.sfc  
 AREA     = -90;-180;90;180
 SFPRMF   = hrlyairp_table54.pack
 DATTIM   = ${PDY}/00${cyc} 
 TIMSTN   = 1/3999
 FORMAT   = mklein
r

EOF

cat <<\gemEOF >gem_namlst
 &INPARM npgiv=2/
gemEOF

$EXECgraph/sffdsf <<  EOF9                  
 SFEFIL   = table52.dat
 SFOUTF   = table52.sfc
 AREA     = -90;-180;90;180
 SFPRMF   = hrlyairp_table52.pack
 DATTIM   = ${PDY}/00${cyc} 
 TIMSTN   = 1/3999
 FORMAT   = mklein
r

EOF9

cat <<\gemEOF >gem_namlst
 &INPARM npgiv=3/
gemEOF

$EXECgraph/sffdsf <<  EOFILE                  
 SFEFIL   = table53.dat
 SFOUTF   = table53.sfc
 AREA     = -90;-180;90;180
 SFPRMF   = hrlyairp_table53.pack
 DATTIM   = ${PDY}/00${cyc}
 TIMSTN   = 1/3999
 FORMAT   = mklein
r

EOFILE


cat <<\gemEOF >gem_namlst
 &INPARM npgiv=2/
gemEOF

$EXECgraph/sffdsf <<  EOFILE                    
 SFEFIL   = table55.dat
 SFOUTF   = table55.sfc
 AREA     = -90;-180;90;180
 SFPRMF   = hrlyairp_table55.pack
 DATTIM   = ${PDY}/00${cyc}
 TIMSTN   = 1/3999
 FORMAT   = mklein
r

EOFILE

cat <<\gemEOF >gem_namlst
 &INPARM npgiv=3/
gemEOF

$EXECgraph/sffdsf <<  EOFILE                   
 SFEFIL   = table56.dat
 SFOUTF   = table56.sfc
 AREA     = -90;-180;90;180
 SFPRMF   = hrlyairp_table56.pack
 DATTIM   = ${PDY}/00${cyc}
 TIMSTN   = 1/3999
 FORMAT   = mklein
r

EOFILE

# Precipitable Water

echo ADDSTN  YES  > gemglb.nts

gdcfil <<EOF                 
 GDOUTF   = table54a.grd
 PROJ     = str/55;-100;0
 GRDAREA  = 105
 KXKY     = 61;41
 MAXGRD   = 100
 CPYFIL   =
 ANLYSS   =

 run

EOF

oabsfc <<EOFILE                   
 SFFILE   = table54.sfc
 GDFILE   = table54a.grd   
 SFPARM   = pwat  
 DATTIM   = ${PDY}/00${cyc}
 DTAAREA  = 105
 GUESS    =
 GAMMA    = 0.3
 SEARCH   = 50/ex
 NPASS    = 2
 QCNTL    =

 run

EOFILE

# Freezing Level

echo ADDSTN  YES  > gemglb.nts

gdcfil <<EOF                 
 GDOUTF   = table54b.grd
 PROJ     = str/55;-100;0
 GRDAREA  = 105
 KXKY     = 61;41
 MAXGRD   = 100
 CPYFIL   =
 ANLYSS   =

 run

EOF

oabsfc <<EOFILE                 
 SFFILE   = table54.sfc
 GDFILE   = table54b.grd   
 SFPARM   = flv1  
 DATTIM   = ${PDY}/00${cyc}
 DTAAREA  = 105
 GUESS    =
 GAMMA    = 0.3
 SEARCH   = 50/ex
 NPASS    = 2
 QCNTL    =

 run

EOFILE

# Lifted Index

echo ADDSTN  YES  > gemglb.nts

gdcfil <<EOF                   
 GDOUTF   = table56.grd
 PROJ     = str/55;-100;0
 GRDAREA  = 105
 KXKY     = 61;41
 MAXGRD   = 100
 CPYFIL   =
 ANLYSS   =

 run

EOF

oabsfc <<EOFILE                 
 SFFILE   = table56.sfc
 GDFILE   = table56.grd   
 SFPARM   = lftd  
 DATTIM   = ${PDY}/00${cyc}
 DTAAREA  = 105
 GUESS    =
 GAMMA    = 0.3
 SEARCH   = 50/ex
 NPASS    = 2
 QCNTL    =

 run

EOFILE

#
#  Relative Humidity
#

echo ADDSTN  YES  > gemglb.nts

gdcfil <<EOF                 
 GDOUTF   = table55.grd
 PROJ     = str/55;-100;0
 GRDAREA  = 105
 KXKY     = 61;41
 MAXGRD   = 100
 CPYFIL   =
 ANLYSS   =

 run

EOF

oabsfc <<EOFILE                
 SFFILE   = table55.sfc
 GDFILE   = table55.grd   
 SFPARM   = rhum  
 DATTIM   = ${PDY}/00${cyc}
 DTAAREA  = 105
 GUESS    =
 GAMMA    = 0.3
 SEARCH   = 50/ex
 NPASS    = 2
 QCNTL    =

 run

EOFILE

$GEMEXE/gpend

exit