#!/bin/csh

####################################################################
#Paul Marinsky                                                     #
#METR 2021-Introduction to Meteorology II Laboratory               #
#                                                                  #
#Instructor: Kevin Goebbert                                        #
#                                                                  #
#Date: 03/25/2009                                                  #
#                                                                  #
#DESCRIPTION AND MODIFICATIONS                                     #
#  This script is designed to create contour plots of both         #
#  temperature and wind.                                           #
####################################################################

#First task is to copy the appropriate color table to our directory#
#To do this we must input a command line as follows                #

cp /home/metr2021/colortables/coltbl.xwp_white coltbl.xwp

#We need a FOREACH loop to process the data through some time series#
#Foreach loop for time indicated in UTC                             #

foreach UTC (00 12)
set LEVEL = "850"

####################################################################
#VARIABLES TO BE USED IN PLOTTING THE MAPS                         #
####################################################################
#Variables used in both SNMAP and GDPLOT2                          #
####################################################################

set DATE = "20060313"
set AREA = "15;-126;55;-40;"
set GAREA = "18;-123;50;-50;"
set PROJ = "lcc/55;-100;20/4;4;4;4"
set DEVICE = "gif | /home/marinsky/http/Lab7_850.gif | 1200;800"
set TITLE = "31/0/850mb Plot  Temp and Wind     Paul Marinsky"
set FILTER = "0.50"
set TEXT = "0.75"
set PANEL = "ALL"
set MAP = "31"
set MSCALE = "0"
set LATLON = "0"

##################################################################
#VARIABLES USED ONLY IN SNMAP                                    #
##################################################################
set SNFILE = "/archive/class/metr2021/${DATE}/upperair/${DATE}_upa.gem"
set DATTIM = "${UTC}"
set SNPARM = "brbk:0.75:0.75:112;tmpc;;stdz;;dpdc"
set COLORS = "31;31;31;31;31;31"
set VCOORD = "PRES"

###################################################################
#VARIABLES USED IN ONLY GDPLOT2                                   #
###################################################################
set GDFILE = "/archive/class/metr2021/${DATE}/models/${DATE}${UTC}_ruc211.gem"
set SKIP = "0"
set SCALE = "999"
set CONTUR = "4"
set GDPFUN1 = "TMPC"
set TYPE1 = "C/F"
set CINT1 = "2/-50/50"
set LINE1 = "2/3/2/1/"
set FINT1 = "${CINT1}"
set FLINE1 = "0"
set GVCORD1 = "PRES"
set GLEVEL1 = "${LEVEL}"
set GDPFUN2 = "HGHT"
set TYPE2 = "C"
set CINT2 = "30/0000/6000"
set LINE2 = "31/1/2/1/"
set GVCORD2 = "PRES"
set GLEVEL2 = "${LEVEL}"
set HILO = ""
set HLSYM = ""
set CLRBAR = "31/V/LL/0.005;0.05/0.5;0.01/1/|0.5/2"

##############################################################################
#GDPLOT2 must be run before SNMAP.  To run the program we have the following:#
##############################################################################

gdplot2 << EOF

GDFILE=${GDFILE}
GDATTIM=F00
GLEVEL=${GLEVEL1} !${GLEVEL2}
GVCORD=${GVCORD1} !${GVCORD2}
PANEL=${PANEL}
SKIP=${SKIP}
SCALE={SCALE}
GDPFUN=${GDPFUN1} !${GDPFUN2}
TYPE=${TYPE1} !${TYPE2}
CONTUR=${CONTUR}
CINT=${CINT1} !${CINT2}
LINE=${LINE1} !${LINE2}
FINT=${FINT1}
FLINE=${FLINE1}
HILO=${HILO}
HLSYM=${HLSYM}
CLRBAR=${CLRBAR}
WIND=
REFVEC=
TITLE=${TITLE}
TEXT=${TEXT}
CLEAR=no
GAREA=${GAREA}
IJSKIP=
PROJ=${PROJ}
MAP=${MAP}
MSCALE=${MSCALE}
LATLON=${LATLON}
DEVICE=${DEVICE}
STNPLOT=
SATFIL=
RADFIL=
IMCBAR=
LUTFIL=
STREAM=
POSN=
COLORS=
MARKER=
GRDLBL=
FILTER=${FILTER}


run


exit

EOF

###################################################
#After running GDPLOT2 we want to run SNMAP       #
###################################################

snmap << EOF

AREA=${AREA}
GAREA=${GAREA}
SATFIL=
RADFIL=
IMCBAR=
SNPARM=${SNPARM}
DATTIM=${DATTIM}
LEVELS=${LEVEL}
VCOORD=${VCOORD}
SNFILE=${SNFILE}
COLORS=${COLORS}
MAP=${MAP}
MSCALE=${MSCALE}
LATLON={LATLON}
TITLE=
CLEAR=NO
PANEL=${PANEL}
DEVICE=${DEVICE}
PROJ=${PROJ}
FILTER=${FILTER}
TEXT=${TEXT}
LUTFIL=
STNPLT=
CLRBAR=
LSTPRM=


run


exit

EOF


gpend

end

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