#!/bin/sh
export member=airplt_redbook.sh
#
#################################################################
#####
#
# NAME     : airplt_redbook.sh
# Purpose  : This script creates 6 air pollution redbook        
#          : graphics using gempak tools
# History  : 07-03-22 Larry Sager  new script 
#
#################################################################
#####
#
cd $DATA
set -x

echo ADDSTN  YES  > gemglb.nts

# 
#   Make the Redbook graphic for awips
#

gdcntr <<EOF                    
 CINT     = .25
 DEVICE   = rbk | i0p || C                               
 GAREA    = awna
 GDATTIM  = ${PDY}/00${cyc}
 GDFILE   = table54a.grd
 GFUNC    = PWAT
 CLEAR    = YES
 CLRBAR   =
 CONTUR   = 0
 CTYPE    = c
 GLEVEL   = 0
 GVCORD   = none
 HILO     =
 HLSYM    =
 IJSKIP   =
 IMCBAR   =
 LATLON   = 1
 LINE     = 1
 LUTFIL   =
 MAP      = 1
 MSCALE   = 0
 PANEL    = 0
 PROJ     =
 RADFIL   =
 SATFIL   =
 SCALE    = 0
 SKIP     = 0
 STNPLT   =
 TEXT     = 1.25/21//hw
 TITLE    =

 run

EOF

echo 'PRECIPITABLE WATER ANALYSIS '   > outtext.txt
echo '                              ' >> outtext.txt

gptext << EOF3
PANEL  = 0
COLORS = 3
TEXT   = 1.25/21//hw
CLEAR  = no
DEVICE = rbk | i0p || C
TXTFIL = outtext.txt
TXTLOC = #23.0;-113.
COLUMN = 1
r
EOF3


echo $TITLEIT   > outtext.txt
echo '         ' >> outtext.txt

gptext << EOF4
PANEL  = 0
COLORS = 3
TEXT   = 1.00/21//hw
CLEAR  = no
DEVICE = rbk | i0p || C
TXTFIL = outtext.txt
TXTLOC = #70.0;-180.
COLUMN = 1
r
EOF4
$GEMEXE/gpend

#
#     Freezing Level 
#
echo ADDSTN  YES  > gemglb.nts

gdcntr <<EOFX
 CLRBAR   =
 CONTUR   = 0
 CLEAR    = YES
 CTYPE    = c
 DEVICE   = rbk| l0f  ||  C                                 
 GAREA    = awna
 GDATTIM  = ${PDY}/00${cyc}
 GDFILE   = table54b.grd
 GFUNC    = flv1
 GLEVEL   = 0
 GVCORD   = none
 HILO     =
 HLSYM    =
 IJSKIP   =
 IMCBAR   =
 LATLON   = 1
 LINE     = 3 
 LUTFIL   =
 MAP      = 1
 MSCALE   = 0
 SCALE    = -2
 CINT     = 40/00/240
 PANEL    = 0
 RADFIL   =
 SATFIL   =
 SKIP     = 0
 STNPLT   =
 TEXT     = 1.25/21//hw
 TITLE    =                           

 run

EOFX

echo '   FREEZING LEVEL ANALYSIS   '   > outtext.txt
echo '      (HUNDREDS OF FEET)     ' >> outtext.txt

gptext << EOF3
PANEL  = 0
COLORS = 3
TEXT   = 1.25/21//hw
CLEAR  = no
DEVICE = rbk | l0f || C
TXTFIL = outtext.txt
TXTLOC = #23.0;-113.
COLUMN = 1
r
EOF3


echo $TITLEIT   > outtext.txt
echo '         ' >> outtext.txt

gptext << EOF4
PANEL  = 0
COLORS = 3
TEXT   = 1.00/21//hw
CLEAR  = no
DEVICE = rbk | l0f || C
TXTFIL = outtext.txt
TXTLOC = #70.0;-180.
COLUMN = 1
r
EOF4

$GEMEXE/gpend
#
#  Plot the precipitable water graphic
#

$GEMEXE/sfmap << ENDFILE
AREA   = dset
GAREA  = awna 
MAP    = 1
DATTIM = all
TEXT   = 1.25/21//hw
COLORS = 1
FILTER = 0
LATLON = 1
TITLE  =                               
SFFILE = table52.sfc
SFPARM  = ;;;;;;;;;stid;pcpw;
DEVICE  = rbk | p03 || C
CLEAR   = YES

run
ENDFILE

$GEMEXE/sfmap << ENDRILE

AREA    = dset
GAREA   = awna
MAP     = 1
DATTIM  = all
COLORS  = 1
TEXT    = 1.25/21//hw
FILTER  = 0
LATLON  = 1
TITLE   =
 SFFILE = table54.sfc
 SFPARM  = stid;
 DEVICE  = rbk | p03  ||  C
 CLEAR   = NO

 run

ENDRILE

echo '   PRECIPITABLE WATER PLOT    '   > outtext.txt
echo '      PERCENT OF NORMAL       ' >> outtext.txt

gptext << EOF3
PANEL  = 0
COLORS = 3
TEXT   = 1.25/21//hw
CLEAR  = no
DEVICE = rbk | p03 || C
TXTFIL = outtext.txt
TXTLOC = #23.0;-113.
COLUMN = 1
r
EOF3


echo $TITLEIT   > outtext.txt
echo '         ' >> outtext.txt

gptext << EOF4
PANEL  = 0
COLORS = 3
TEXT   = 1.00/21//hw
CLEAR  = no
DEVICE = rbk | p03 || C
TXTFIL = outtext.txt
TXTLOC = #70.0;-180.
COLUMN = 1
r
EOF4

$GEMEXE/gpend
#
#  Plot the Redbook freezing level graphic

$GEMEXE/sfmap << ENDMILE

AREA    = dset
GAREA   = awna
MAP     = 1
DATTIM  = all
TEXT    = 1.25/21//hw
COLORS  = 1              
FILTER  = 0
LATLON  = 1
TITLE   =                          
SFFILE = table53.sfc
SFPARM  = ;;;;;;;;;stid;
DEVICE  = rbk | l0a || C                                 
CLEAR   = YES

run

ENDMILE

$GEMEXE/sfmap << ENDRILE

AREA    = dset
GAREA   = awna
MAP     = 1
DATTIM  = all
COLORS  = 1
TEXT    = 1.25/21//hw
FILTER  = 0
LATLON  = 1
TITLE   =
 SFFILE = table54.sfc
 SFPARM  = stid;
 DEVICE  = rbk | l0a  ||  C
 CLEAR   = NO

 run

ENDRILE
echo '   FREEZING LEVEL PLOT   '   > outtext.txt
echo '    (HUNDREDS OF FEET)   ' >> outtext.txt

gptext << EOF3
PANEL  = 0
COLORS = 3
TEXT   = 1.25/21//hw
CLEAR  = no
DEVICE = rbk | l0a || C
TXTFIL = outtext.txt
TXTLOC = #23.0;-113.
COLUMN = 1
r
EOF3


echo $TITLEIT   > outtext.txt
echo '         ' >> outtext.txt

gptext << EOF4
PANEL  = 0
COLORS = 3
TEXT   = 1.00/21//hw
CLEAR  = no
DEVICE = rbk | l0a || C
TXTFIL = outtext.txt
TXTLOC = #70.0;-180.
COLUMN = 1
r
EOF4
$GEMEXE/gpend

#
#  Plot the Lifted/K Index graphic
#

$GEMEXE/sfmap << ENDFILE
AREA    = dset
GAREA   = awna
MAP     = 1
DATTIM  = all
TEXT    = 1.25/21//hw
COLORS  = 1               
FILTER  = 0
LATLON  = 1
TITLE   =                          
 SFFILE = table56.sfc
 SFPARM  = ;;;;;;;;;lftd;stid
 DEVICE  = rbk | p01  || C                       
 CLEAR   = YES

 run

ENDFILE

$GEMEXE/sfmap << ENDFILE

AREA    = dset
GAREA   = awna
MAP     = 1
DATTIM  = all
COLORS  = 1
TEXT    = 1.25/21//hw
FILTER  = 0
LATLON  = 1
TITLE   =
 SFFILE = table54.sfc
 SFPARM  = stid;
 DEVICE  = rbk | p01  ||  C
 CLEAR   = NO

 run

ENDFILE

echo '   LIFTED INDEX PLOT    '   > outtext.txt
echo '       K INDEX PLOT        ' >> outtext.txt
echo '         ' >> outtext.txt

gptext << EOF3
PANEL  = 0
COLORS = 3
TEXT   = 1.25/21//hw
CLEAR  = no
DEVICE = rbk | p01 || C
TXTFIL = outtext.txt
TXTLOC = #23.0;-113.
COLUMN = 1
r
EOF3


echo $TITLEIT   > outtext.txt
echo '         ' >> outtext.txt

gptext << EOF4
PANEL  = 0
COLORS = 3
TEXT   = 1.00/21//hw
CLEAR  = no
DEVICE = rbk | p01 || C
TXTFIL = outtext.txt
TXTLOC = #70.0;-180.
COLUMN = 1
r
EOF4
$GEMEXE/gpend

#  Plot the Relative Humidity graphic
#

$GEMEXE/sfmap << ENDFILE

AREA    = dset
GAREA   = awna
MAP     = 1
DATTIM  = all
TEXT    = 1.25/21//hw
COLORS  = 1              
FILTER  = 0
TITLE   =                 
 SFFILE = table55.sfc    
 SFPARM  = ;;;;;;;;;stid;
 DEVICE  = rbk | p0d  ||  C                         
 CLEAR   = YES

 run
 
ENDFILE

$GEMEXE/sfmap << ENDFILE

AREA    = dset
GAREA   = awna
MAP     = 1
DATTIM  = all
COLORS  = 1              
TEXT    = 1.25/21//hw
FILTER  = 0
LATLON  = 1 
TITLE   =                            
 SFFILE = table54.sfc    
 SFPARM  = stid;
 DEVICE  = rbk | p0d  ||  C                         
 CLEAR   = NO  

 run
 
ENDFILE


echo '  RELATIVE HUMIDITY PLOT    '   > outtext.txt
echo 'SURFACE TO 500MB  (PERCENT) ' >> outtext.txt
echo '                            ' >> outtext.txt
echo '         ' >> outtext.txt

gptext << EOF3
PANEL  = 0
COLORS = 3
TEXT   = 1.25/21//hw 
CLEAR  = no
DEVICE = rbk | p0d || C                                 
TXTFIL = outtext.txt
TXTLOC = #23.0;-113.
COLUMN = 1
r
EOF3


echo $TITLEIT   > outtext.txt
echo '         ' >> outtext.txt

gptext << EOF4
PANEL  = 0
COLORS = 3
TEXT   = 1.00/21//hw
CLEAR  = no
DEVICE = rbk | p0d || C
TXTFIL = outtext.txt
TXTLOC = #70.0;-180.
COLUMN = 1
r
EOF4
$GEMEXE/gpend

#  Add the TOC header and send 

  export PDYHH=${PDY}${cyc}
  if test "$SENDCOM" = 'YES'
  then

    INPATH=$DATA/NMCGPHI0P
    OUTPATH=$pcom/i0p.redbook.pewa34.${job}
    WMO=PEWA34
    $USHutil/make_ntc_bull.pl redb NONE KWBC NONE $INPATH $OUTPATH
    
    INPATH=$DATA/NMCGPHL0F
    OUTPATH=$pcom/l0f.redbook.ptwa93.${job}
    WMO=PTWA93
    $USHutil/make_ntc_bull.pl redb NONE KWNO NONE $INPATH $OUTPATH

    INPATH=$DATA/NMCGPHP01
    OUTPATH=$pcom/p01.redbook.pxwa50.${job}
    WMO=PXWA50
    $USHutil/make_ntc_bull.pl redb NONE KWBC NONE $INPATH $OUTPATH

    INPATH=$DATA/NMCGPHP03
    OUTPATH=$pcom/p03.redbook.pewa33.${job}
    WMO=PEWA33
    $USHutil/make_ntc_bull.pl redb NONE KWBC NONE $INPATH $OUTPATH

    INPATH=$DATA/NMCGPHP0D
    OUTPATH=$pcom/p0d.redbook.prwa51.${job}
    WMO=PRWA51
    $USHutil/make_ntc_bull.pl redb NONE KWBC NONE $INPATH $OUTPATH

    INPATH=$DATA/NMCGPHL0A
    OUTPATH=$pcom/l0a.redbook.ptwa94.${job}
    WMO=PTWA94
    $USHutil/make_ntc_bull.pl redb NONE KWBC NONE  $INPATH $OUTPATH

  fi

exit
