#!/bin/csh

########################################################################
#
#  Patrick Marsh
#  METR2021
#  Lab 3
# 
#  Date: 9 February 2009
#
#  Name:  500mb Temperature and Heights Map
#
#  Description:  This script creates a 500mb map with filled
#                temperatures and contoured heights.  The model used is
#                the RUC 236.
#
#
#
#  Modifications: Fixed the contour/fill switcheroo and added my name.
#		  Travis Darling, 18 February 2009
#
#                
########################################################################


source /usr/local/nawips/Gemenviron

# Copy special upperair colortable used on the hoot site
cp /home/metr2021/colortables/upperairtemp.xwp coltbl.xwp

########################################################################
# FOREACH LOOP
########################################################################

foreach hour (f00 f01 f02 f03 f06 f09 f12)

########################################################################
# SET VARIABLES FOR GEMPAK
########################################################################

set GDFILE  = "/data/gempak/model/ruc/2009021018_ruc236.gem"
set GDATTIM = "${hour}"
set GLEVEL  = "500"
set GVCORD  = "PRES"
set PANEL   = "ALL"
set SKIP    = "0"
set SCALE   = "999"
set GDPFUN2 = "HGHT"
set GDPFUN1 = "TMPC"
set TYPE1   = "F"
set TYPE2   = "C"
set CONTUR  = "5"
set CINT    = "60" 
set LINE    = "32"
set FINT    = "5"
set FLINE   = "03;04;05;06;07;08;09;10;11;12;13;14;15;16;17"
set HILO    = 
set HLSYM   = 
set CLRBAR  = "32/V/LL/0.005;0.05/0.5;0.01/1/|0.6/2"
set WIND    = 
set REFVEC  = 
set TITLE   = "32/0/Travis Darling's Modified Map at time ${hour}"
set TEXT    = "1"
set CLEAR   = "YES"
set GAREA   = "16.28;-126.14;55.48;-57.38"
set IJSKIP  = 
set PROJ    = "lcc/25.0;-95.0;25.0/4;4;4;4"
set MAP     = "32"
set MSCALE  = "0"
set LATLON  = 
set DEVICE  = "gif|/webhost/darlingt/http/2021L/4/500mb${hour}.gif|1200;800"
set STNPLT  = 
set SATFIL  = 
set RADFIL  = 
set IMCBAR  = 
set LUTFIL  = 
set STREAM  = 
set POSN    = "0"
set COLORS  = "1"
set MARKER  = "0"
set GRDLBL  = "0"
set FILTER  = "0"


########################################################################
# RUN GEMPAK
########################################################################

gdplot2 << EOF

    GDFILE  = ${GDFILE}
    GDATTIM = ${GDATTIM}
    GLEVEL  = ${GLEVEL}
    GVCORD  = ${GVCORD}
    PANEL   = ${PANEL}
    SKIP    = ${SKIP}
    SCALE   = ${SCALE}
    GDPFUN  = ${GDPFUN1} !${GDPFUN2}
    TYPE    = ${TYPE1} !${TYPE2}
    CONTUR  = ${CONTUR}
    CINT    = ${CINT}
    LINE    = ${LINE}
    FINT    = ${FINT}
    FLINE   = ${FLINE}
    HILO    = ${HILO}
    HLSYM   = ${HLSYM}
    CLRBAR  = ${CLRBAR}
    WIND    = ${WIND}
    REFVEC  = ${REFVEC}
    TITLE   = ${TITLE}
    TEXT    = ${TEXT}
    CLEAR   = ${CLEAR}
    GAREA   = ${GAREA}
    IJSKIP  = ${IJSKIP}
    PROJ    = ${PROJ}
    MAP     = ${MAP}
    MSCALE  = ${MSCALE}
    LATLON  = ${LATLON}
    DEVICE  = ${DEVICE}
    STNPLT  = ${STNPLT}
    SATFIL  = ${SATFIL}
    RADFIL  = ${RADFIL}
    IMCBAR  = ${IMCBAR}
    LUTFIL  = ${LUTFIL}
    STREAM  = ${STREAM}
    POSN    = ${POSN}
    COLORS  = ${COLORS}
    MARKER  = ${MARKER}
    GRDLBL  = ${GRDLBL}
    FILTER  = ${FILTER}

r

exit

EOF

# Closes Graphics Window and puts image on website
gpend

# Ends the foreach loop
end

rm -f *.nts
rm -f *.xwp