#!/bin/csh

#########################################################################
#									#
#  Paul Marinsky 			                                #
#  METR2021								#
#  LABORATORY EXERCISE #3 (problem 4)					#
# 									#
#  Date: 22 April 2009							#
#									#
#  INSTRUCTOR:  Kevin Goebbert						#
#									#
#  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: Modified Script Lab3_3.g to loop RUC236 times		#
#									#
#                							#
#########################################################################


source /usr/local/nawips/Gemenviron

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


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

########################################################################
# THE LOOP BELOW WILL RUN FOR THE SELECTED HOURS SPECIFIED IN THE TIME
# VARIABLE.  THAT MODEL RUN TIME WILL THEN BE OUTPUT TO A GRAPHIC AND
# THE NEXT TIME WILL BE RUN THROUGH THE LOOP.
########################################################################

foreach TIME(00 03 06 09 12)

set GDFILE  = "/data/gempak/model/ruc/2009033021_ruc236.gem"
set GDATTIM = "f${TIME}"
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    = "5000-6000"
set LINE    = "32"
set FINT    = "-50;-45;-40;-35;-30;-25;-20;-15;-10;-5;00;05;10;15"
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/H5 Temperatures and Heights ${TIME}UTC 090330  Paul Marinsky 2009"
set TEXT    = "0.75"
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|/home/marinsky/http/lab3_4_${TIME}.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}

run

exit

EOF

gpend # Closes Graphics Window and puts image on website

end

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