#!/bin/csh

###################################################################
#Paul Marinsky                                                    #
#METR 2021  Introduction to Meteorology II Laboratory             #
#                                                                 #
#Instructor: Kevin Goebbert                                       #
#                                                                 #
#Date: 03-25-2009                                                 #
#                                                                 #
#DESCRIPTION or MODIFICATIONS                                     #
#  This script is designed to run GDPLOT2 and SNMAP to create     #
#  an image of 850 temperature advection                          #
###################################################################

#############################################################
#use the -cp command line to copy the designated color table#
#############################################################

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

######################################################################
#we need a FOREACH loop to process the data over the appropriate time#
######################################################################

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

#############################
#Constant Variables         #
#############################

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_850advection.gif | 1200;800"
set TITLE = "31/0/850mb Plot  Temperature Advection     Paul Marinsky"
set FILTER = "0.50"
set TEXT = "1"
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 = "4 \!999"
set CONTUR = "5"
set GDPFUN1 = "ADV(TMPC,WND)"
set TYPE1 = "F/C"
set CINT1 = "2/-12/12"
set LINE1 = "31/1/1/1/"
set FINT1 = "2/-12/12/"
set FLINE1 = "24;25;26;27;6;32;32;11;12;13;2;8;32"
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=YES
GAREA=${GAREA}
IJSKIP=
PROJ=${PROJ}
MAP=${MAP}
MSCALE=${MSCALE}
LATLON=${LATLON}
DEVICE=${DEVICE}
STNPLT=
SATFIL=
RADFIL=
IMCBAR=
LUTFIL=
STREAM=
POSN=
COLORS=
MARKER=
GRDLBL=
FILTER=${FILTER}


run


exit

EOF

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=

exit

EOF

gpend

end

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

###########################################################################
#If this script has run correctly, two images should have plotted to your #
#web directory that shows 850mb temperature plotted with wind barbs.  If  #
#the script was not successful then there is obviously a problem within   #
#the code and you should take this time to debug.  Once you think you have#
#corrected any errors in the command code, try running the script again.  #
#If problems still exist, then you should contact your instructor.        #
###########################################################################