#!/bin/sh


#########################################################################
#
#   Script:  gempak_gfs_f00_gif.sh
#
#   This scripts creates GEMPAK .gif images of 00HR/Analysis fields from
#   GFS model output for archiving at NCDC.
#
#
#   History:   Ralph Jones     02/16/2005   JIF original version.
#   History:   Steve Lilly     04/30/2008   Change font size of the Titles
#                                           from .8 to a larger size (1 or 2)
#
#
#########################################################################




   msg=" Make GEMPAK GIFS utility"
   postmsg "$jlogfile" "$msg"


  set -x


  MAPAREA="normal"

  LATVAL="1/1/1/1/5;5"

  pixels="1728;1472"

  cp $FIXgempak/coltbl.spc coltbl.xwp


#################################################################
#                       ANALYSIS CHARTS                         # 
#################################################################


# Create time stamp (bottom) label 

  echo 0000${PDY}${cyc} > dates
  export XLFUNIT_55="title.output"
  $utilexec/webtitle < dates
  export TITLE=`cat fort.55`
  echo "\n\n TITLE = $TITLE \n"


# Define labels and file names for analysis charts

  hgttmp700lab="700MB ANALYSIS  HEIGHTS/TEMPERATURE"
  hgttmp700dev="gfs_700_hgt_tmp_nh_anl_${cyc}.gif"

  hgttmp500lab="500MB ANALYSIS  HEIGHTS/TEMPERATURE"
  hgttmp500dev="gfs_500_hgt_tmp_nh_anl_${cyc}.gif"

  hgtiso500lab="500MB ANALYSIS  HEIGHTS/ISOTACHS"
  hgtiso500dev="gfs_500_hgt_iso_nh_anl_${cyc}.gif"

  hgtiso300lab="300MB ANALYSIS  HEIGHTS/ISOTACHS"
  hgtiso300dev="gfs_300_hgt_iso_nh_anl_${cyc}.gif"

  hgtiso250lab="250MB ANALYSIS  HEIGHTS/ISOTACHS"
  hgtiso250dev="gfs_250_hgt_iso_nh_anl_${cyc}.gif"

  hgttmp250lab="250MB ANALYSIS  HEIGHTS/TEMPERATURE"
  hgttmp250dev="gfs_250_hgt_tmp_nh_anl_${cyc}.gif"

  hgtiso200lab="200MB ANALYSIS  HEIGHTS/ISOTACHS"
  hgtiso200dev="gfs_200_hgt_iso_nh_anl_${cyc}.gif"

  hgttmp200lab="200MB ANALYSIS  HEIGHTS/TEMPERATURE"
  hgttmp200dev="gfs_200_hgt_tmp_nh_anl_${cyc}.gif"

  hgtiso100lab="100MB ANALYSIS  HEIGHTS/ISOTACHS"
  hgtiso100dev="gfs_100_hgt_iso_nh_anl_${cyc}.gif"

  hgttmp100lab="100MB ANALYSIS  HEIGHTS/TEMPERATURE"
  hgttmp100dev="gfs_100_hgt_tmp_nh_anl_${cyc}.gif"

  hgtvor500lab="500MB ANALYSIS  HEIGHTS/VORTICITY"
  hgtvor500dev="gfs_500_hgt_vor_nh_anl_${cyc}.gif"

  hgtvor500usdev="gfs_500_hgt_vor_uscan_anl_${cyc}.gif"

  mslpthksfclab="ANALYSIS  MEAN SEA LEVEL PRESSURE/1000-500MB THICKNESS"
  mslpthksfcdev="gfs_sfc_mslp_thk_nh_anl_${cyc}.gif"
 
  mslpthksfcusdev="gfs_sfc_mslp_thk_uscan_anl_${cyc}.gif"

  rhvvel700lab="700MB ANALYSIS  RH/VERT VEL"
  rhvvel700dev="gfs_700_rh_vvel_nh_anl_${cyc}.gif"

  liftlab="ANALYSIS  LIFTED INDEX"
  liftdev="gfs_lift_nh_anl_${cyc}.gif"

  prswshtroplab="TROPOPAUSE PRESSURE/WIND SHEAR"
  prswshtropdev="gfs_trop_prs_wsh_nh_anl_${cyc}.gif"


# Set grid date and input file name

  gdattim=`echo ${PDY} | cut -c3-8`/${cyc}00F000
  gdfile=gem_grids${fhr}.gem



#  Execute the GEMPAK program

  $GEMEXE/gdplot2_gif << EOF


! 700MB HEIGHTS/TEMPERATURES

  restore $NTS/base_nh
  restore $NTS/700_hgt_tmp

  CLEAR   = yes
  GDFILE  = $gdfile
  GDATTIM = $gdattim
  MAP     = 1
  DEVICE  = gif | ${hgttmp700dev} | $pixels
  TITLE   =
  TEXT    = 1/2/2/c/sw
  LATLON  = $LATVAL
  l
  r

  CLEAR   = no
  GDPFUN  =
  TITLE   = 1/-4/$TITLE
  TEXT    = 2/2/2/c/sw 
  LATLON  = 0
  l
  r

  TITLE   = 1/3/${hgttmp700lab}
  l
  r


! 500MB HEIGHTS/TEMPERATURES

  restore $NTS/base_nh
  restore $NTS/500_hgt_tmp

  CLEAR   = yes
  GDFILE  = $gdfile
  GDATTIM = $gdattim
  MAP     = 1
  DEVICE  = gif | ${hgttmp500dev} | $pixels
  TITLE   =
  TEXT    = 1/2/2/c/sw
  LATLON  = $LATVAL
  l
  r

  CLEAR   = no
  GDPFUN  =
  TITLE   = 1/-4/$TITLE
  TEXT    = 2/2/2/c/sw 
  LATLON  = 0
  l
  r

  TITLE   = 1/3/${hgttmp500lab}
  l
  r


! 300MB HEIGHTS/ISOTACHS

  restore $NTS/base_nh
  restore $NTS/300_hgt_iso

  CLEAR   = yes
  GDFILE  = $gdfile
  GDATTIM = $gdattim
  MAP     = 1
  DEVICE  = gif | ${hgtiso300dev} | $pixels
  TITLE   =
  TEXT    = 1/2/2/c/sw 
  LATLON  = 1/1/1/1/5;5             !
  l
  r

  CLEAR   = no
  GDPFUN  =
  TITLE   = 1/-4/$TITLE
  TEXT    = 2/2/2/c/sw 
  LATLON  = 0
  l
  r

  TITLE   = 1/3/${hgtiso300lab}
  l
  r


! 250MB HEIGHTS/TEMPERATURES

  restore $NTS/base_nh
  restore $NTS/250_hgt_tmp

  CLEAR   = yes
  GDFILE  = $gdfile
  GDATTIM = $gdattim
  MAP     = 1
  DEVICE  = gif | ${hgttmp250dev} | $pixels
  TITLE   =
  TEXT    = 1/2/2/c/sw 
  LATLON  = $LATVAL
  l
  r

  CLEAR   = no
  GDPFUN  =
  TITLE   = 1/-4/$TITLE
  TEXT    = 2/2/2/c/sw 
  LATLON  = 0
  l
  r

  TITLE   = 1/3/${hgttmp250lab}
  l
  r


! 250MB ANALYSIS HEIGHTS/ISOTACHS

  restore $NTS/base_nh
  restore $NTS/250_hgt_iso

  CLEAR   = yes 
  GDFILE  = ${gdfile}
  GDATTIM = ${gdattim}
  MAP     = 1                       
  DEVICE  = gif | ${hgtiso250dev} | $pixels 
  TITLE   = 
  TEXT    = 1/2/2/c/sw 
  LATLON  = $LATVAL 
  l
  r

  CLEAR   = no
  MAP     = 0
  GDPFUN  =
  TITLE   = 1/-4/${TITLE}
  TEXT    = 2/2/2/c/sw 
  LATLON  = 0
  l
  r

  TITLE   = 1/3/${hgtiso250lab}
  l
  r


! 200MB HEIGHTS/ISOTACHS

  restore $NTS/base_nh
  restore $NTS/200_hgt_iso

  CLEAR   = yes
  GDFILE  = $gdfile
  GDATTIM = $gdattim
  MAP     = 1
  DEVICE  = gif | ${hgtiso200dev} | $pixels
  TITLE   =
  TEXT    = 1/2/2/c/sw 
  LATLON  = 1/1/1/1/5;5             !
  l
  r

  CLEAR   = no
  GDPFUN  =
  TITLE   = 1/-4/$TITLE
  TEXT    = 2/2/2/c/sw 
  LATLON  = 0
  l
  r

  TITLE   = 1/3/${hgtiso200lab}
  l
  r


! 100MB HEIGHTS/TEMPERATURES

  restore $NTS/base_nh
  restore $NTS/100_hgt_tmp

  CLEAR   = yes
  GDFILE  = $gdfile
  GDATTIM = $gdattim
  MAP     = 1
  DEVICE  = gif | ${hgttmp100dev} | $pixels
  TITLE   =
  TEXT    = 1/2/2/c/sw 
  LATLON  = $LATVAL
  l
  r

  CLEAR   = no
  GDPFUN  =
  TITLE   = 1/-4/$TITLE
  TEXT    = 2/2/2/c/sw 
  LATLON  = 0
  l
  r

  TITLE   = 1/3/${hgttmp100lab}
  l
  r


! 100MB HEIGHTS/ISOTACHS

  restore $NTS/base_nh
  restore $NTS/100_hgt_iso

  CLEAR   = yes
  GDFILE  = $gdfile
  GDATTIM = $gdattim
  MAP     = 1
  DEVICE  = gif | ${hgtiso100dev} | $pixels
  TITLE   =
  TEXT    = 1/2/2/c/sw 
  LATLON  = 1/1/1/1/5;5             !
  l
  r

  CLEAR   = no
  GDPFUN  =
  TITLE   = 1/-4/$TITLE
  TEXT    = 2/2/2/c/sw 
  LATLON  = 0
  l
  r

  TITLE   = 1/3/${hgtiso100lab}
  l
  r


! ANALYSIS MSLP/1000-500 THICKNESS

  restore $NTS/base_nh
  restore $NTS/sfc_mslp_thk

  CLEAR   = yes
  GDFILE  = ${gdfile}
  GDATTIM = ${gdattim}
  MAP     = 1
  DEVICE  = gif | ${mslpthksfcdev} | $pixels
  TITLE   =
  TEXT    = 1/2/2/c/sw 
  LATLON  = $LATVAL
  l
  r

  CLEAR   = no
  MAP     = 0
  GDPFUN  =
  TITLE   = 1/-4/${TITLE}
  TEXT    = 2/2/2/c/sw 
  LATLON  = 0
  l
  r

  TITLE   = 1/3/${mslpthksfclab}
  l
  r


! ANALYSIS MSLP/1000-500 THICKNESS (US/CANADA)

  restore $NTS/base_uscan
  restore $NTS/sfc_mslp_thk

  CLEAR   = yes
  GDFILE  = ${gdfile}
  GDATTIM = ${gdattim}
  MAP     = 1
  DEVICE  = gif | ${mslpthksfcusdev} | $pixels
  TITLE   =
  TEXT    = 1/2/2/c/sw 
  LATLON  = $LATVAL
  l
  r

  CLEAR   = no
  MAP     = 0
  GDPFUN  =
  TITLE   = 1/-4/${TITLE}
  TEXT    = 2/2/2/c/sw 
  LATLON  = 0
  l
  r

  TITLE   = 1/3/${mslpthksfclab}
  l
  r



! 500MB ANALYSIS  HEIGHTS/VORTICITY

  restore $NTS/base_nh
  restore $NTS/500_hgt_vor

  CLEAR   = yes
  GDFILE  = ${gdfile}
  GDATTIM = ${gdattim}
  MAP     = 1
  DEVICE  = gif | ${hgtvor500dev} | $pixels
  TITLE   =
  TEXT    = 1/2/2/c/sw 
  LATLON  = $LATVAL
  l
  r

  CLEAR   = no
  MAP     = 0
  GDPFUN  =
  TITLE   = 1/-4/${TITLE}
  TEXT    = 2/2/2/c/sw 
  LATLON  = 0
  l
  r

  TITLE   = 1/3/${hgtvor500lab}
  l
  r



! 500MB ANALYSIS  HEIGHTS/VORTICITY (US/CANADA)

  restore $NTS/base_uscan
  restore $NTS/500_hgt_vor

  CLEAR   = yes
  GDFILE  = ${gdfile}
  GDATTIM = ${gdattim}
  MAP     = 1
  DEVICE  = gif | ${hgtvor500usdev} | $pixels
  TITLE   =
  TEXT    = 1/2/2/c/sw 
  LATLON  = $LATVAL
  l
  r

  CLEAR   = no
  GDPFUN  =
  TITLE   = 1/-4/${TITLE}
  TEXT    = 2/2/2/c/sw 
  LATLON  = 0
  l

  TITLE   = 1/3/${hgtvor500lab}
  l
  r


! ANALYSIS  LIFTED INDEX

  restore $NTS/base_nh
  restore $NTS/100_lift

  CLEAR   = yes
  GDFILE  = ${gdfile}
  GDATTIM = ${gdattim}
  MAP     = 1
  DEVICE  = gif | ${liftdev} | $pixels
  TITLE   =
  TEXT    = 1/2/2/c/sw 
  LATLON  = $LATVAL
  l
  r

  CLEAR   = no
  GDPFUN  =
  TITLE   = 1/-4/${TITLE}
  TEXT    = 2/2/2/c/sw 
  LATLON  = 0
  l
  r

  TITLE   = 1/3/${liftlab}
  l
  r


! ANALYSIS  TROPOPAUSE PRESSURE/WIND SHEAR

  restore $NTS/base_nh
  restore $NTS/trop_pres_wshr

  CLEAR   = yes
  GDFILE  = ${gdfile}
  GDATTIM = ${gdattim}
  MAP     = 1
  DEVICE  = gif | ${prswshtropdev} | $pixels
  TITLE   =
  TEXT    = 1/2/2/c/sw 
  LATLON  = $LATVAL
  l
  r

  CLEAR   = no
  GDPFUN  =
  TITLE   = 1/-4/${TITLE}
  TEXT    = 2/2/2/c/sw 
  LATLON  = 0
  l
  r

  TITLE   = 1/3/${prswshtroplab}
  l
  r


! ANALYSIS 700MB RELATIVE HUMIDITY AND VERTICAL VELOCITY

  restore $NTS/base_nh
  restore $NTS/700_rel_vvel

  CLEAR   = yes
  GDFILE  = ${gdfile}
  GDATTIM = ${gdattim}
  MAP     = 1
  DEVICE  = gif | ${rhvvel700dev} | $pixels
  TITLE   =
  TEXT    = 1/2/2/c/sw 
  LATLON  = $LATVAL
  l
  r

  CLEAR   = no
  GDPFUN  =
  TITLE   = 1/-4/${TITLE}
  TEXT    = 2/2/2/c/sw 
  LATLON  = 0
  l
  r

  TITLE   = 1/3/${rhvvel700lab}
  l
  r

  exit
EOF


$GEMEXE/gpend


if [ $SENDCOM = YES ]; then

# Copy the GIF images into my area

  cp ${hgttmp700dev}    $COMNCDC
  cp ${hgttmp500dev}    $COMNCDC
  cp ${hgtiso300dev}    $COMNCDC
  cp ${hgtiso250dev}    $COMNCDC
  cp ${hgttmp250dev}    $COMNCDC
  cp ${hgtiso200dev}    $COMNCDC
  cp ${hgtiso100dev}    $COMNCDC
  cp ${hgttmp100dev}    $COMNCDC
  cp ${mslpthksfcdev}   $COMNCDC
  cp ${mslpthksfcusdev} $COMNCDC
  cp ${hgtvor500dev}    $COMNCDC
  cp ${hgtvor500usdev}  $COMNCDC
  cp ${liftdev}         $COMNCDC
  cp ${prswshtropdev}   $COMNCDC
  cp ${rhvvel700dev}    $COMNCDC



# Copy the GIF images onto the NCDC area on the public ftp server

 if [ $SENDDBN = YES ]; then


  $DBNROOT/bin/dbn_alert MODEL NCDCGIF ${job} $COMNCDC/${hgttmp700dev}
  $DBNROOT/bin/dbn_alert MODEL NCDCGIF ${job} $COMNCDC/${hgttmp500dev}
  $DBNROOT/bin/dbn_alert MODEL NCDCGIF ${job} $COMNCDC/${hgtiso300dev}
  $DBNROOT/bin/dbn_alert MODEL NCDCGIF ${job} $COMNCDC/${hgtiso250dev}
  $DBNROOT/bin/dbn_alert MODEL NCDCGIF ${job} $COMNCDC/${hgttmp250dev}
  $DBNROOT/bin/dbn_alert MODEL NCDCGIF ${job} $COMNCDC/${hgtiso200dev}
  $DBNROOT/bin/dbn_alert MODEL NCDCGIF ${job} $COMNCDC/${hgttmp200dev}
  $DBNROOT/bin/dbn_alert MODEL NCDCGIF ${job} $COMNCDC/${hgtiso100dev}
  $DBNROOT/bin/dbn_alert MODEL NCDCGIF ${job} $COMNCDC/${hgttmp100dev}
  $DBNROOT/bin/dbn_alert MODEL NCDCGIF ${job} $COMNCDC/${mslpthksfcdev}
  $DBNROOT/bin/dbn_alert MODEL NCDCGIF ${job} $COMNCDC/${mslpthksfcusdev}
  $DBNROOT/bin/dbn_alert MODEL NCDCGIF ${job} $COMNCDC/${hgtvor500dev}
  $DBNROOT/bin/dbn_alert MODEL NCDCGIF ${job} $COMNCDC/${hgtvor500usdev}
  $DBNROOT/bin/dbn_alert MODEL NCDCGIF ${job} $COMNCDC/${liftdev}
  $DBNROOT/bin/dbn_alert MODEL NCDCGIF ${job} $COMNCDC/${prswshtropdev}
  $DBNROOT/bin/dbn_alert MODEL NCDCGIF ${job} $COMNCDC/${rhvvel700dev}


# Convert the 500mb NH Hgts/Temps chart to tif, attach a heading and 
#   send to TOC via the NTC

  export input=$COMNCDC/${hgttmp500dev}
  export HEADER=YES
  export OUTPATH=$DATA/gfs_500_hgt_tmp_nh_anl_${cyc}.tif
  /nwprod/util/ush/make_tif.sh

 fi

fi 



   msg=" GEMPAK_GIF ${fhr} hour completed normally"
   postmsg "$jlogfile" "$msg"

   exit