###########################################################################
#
#	GEMPAK FOUR-PANEL SCRIPT FOR EXAMINING CSI 
#          IN CROSS SECTIONS AND HORIZONTAL MAPS
#
#                      WRITTEN BY:
#           PHIL SCHUMACHER AND DAVID SCHULTZ
#                     23 APRIL 1998
#
#           CONTACTS:   Philip.Schumacher@noaa.gov
#			schultz@nssl.noaa.gov
#
#            AVAILABLE FOR FREE DISTRIBUTION
#
#      NOTE:  This script will only work as is with GEMPAK version 5.4.
#      For earlier versions of GEMPAK, pvor(geo) and pvor(obs) should
#      replace the present versions of pvor.  There will be no simple
#      way to include the effects of theta-e* (THES) at this time.
#      Also, version 5.2 pvor does not multiply by g (gravity) to get
#      PVU units.  Modify the SCALE parameter accordingly (reduce it by
#      one).  Also thts replaces thes.   8/19/98  DS
#
###########################################################################
#
#  Select the following parameters:

#  Plot type:  
        set plot = "fourpanel1"

#  Set GEMPAK file name:
        set grid = /nssl/nsslsun/schultz/csi/gempak-phil/98010215_ruc.gem

#  Set Model name:
        set model = "RUC"

#  Set GEMPAK date of event:
        set gdat1 = "980102/1500F009"
        set date = "0000 UTC 3 JANUARY 1998"

#  Set Cross-section and horizontal-map areas:
        set area = "FAR>52;-98"
        set area = "ABR>51;-98"
	set maparea = "pir-"
	set maparea = "43;-105;53;-90"

#  Set Lower-tropospheric levels for assessment of stability:
	set lowlev = "600:650"
	set lowlev1 = 600
	set lowlev2 = 650

#  Set Lower-tropospheric level for assessment of frontogenesis:
	set frntlev = "800"

#  Set Upper-tropospheric levels for assessment of forcing:
	set uplev = "300:400"
	set uplev1 = 300
	set uplev2 = 400

#  Thermodynamic variable:
#     Use theta = thte (theta-e) to assess PSI
#     Use theta = thes (theta-e*) to assess CSI
#	set theta = thes
	set theta = thte

#  Plotting controls:
        set lift = -2

###########################################################################
#
#		No changes necessary below this point.
#
###########################################################################
if ($theta == "thes") then
	set thetaname = "THETA-E*"
	set mpvname = "MPVg*"
	set instab = "CONDITIONAL INST."
endif
if ($theta == "thte") then
	set thetaname = "THETA-E"
	set mpvname = "MPVg"
	set instab = "POTENTIAL INST."
endif
#--------------------------------------------------------------------------
if ($plot == "fourpanel1") then
gdcross << end1
CXSTNS  = $area
GDATTIM = $gdat1
GVCORD  = pres 
GFUNC   = pvor($theta;geo) 
GVECT   =
GDFILE  = $grid
WIND    = bk07
REFVEC  =
PTYPE   = LOG
YAXIS   = 1000/200   
CINT    = -1;-0.5;-0.25;0
SCALE   = 6 ! 0 
LINE    = 32/1 
BORDER  = 1  
TITLE   = 1/$lift/(a) $mpvname (shaded), $thetaname 
CLEAR   = y  
DEVICE  = ps
TEXT    = 1//a
PANEL   = ul/1//2
CLRBAR  = 1/h/UL/.64;.99/.2;.03
CLRBAR  = 1/V/UL/.94;.95/.005;.05
CLRBAR  = 1/V/UL/.927;.805/.2;.015
CONTUR  = 3
FINT    = -1;0
FLINE   = 16;18;0
CTYPE   = f
r

CLEAR   = n  
GFUNC   = $theta
CINT    = 4
SCALE   = 0
LINE    = 1/1/1/3
TITLE   = 1/2/$model : $gdat1 : $date 
TEXT    = 1
CTYPE   = c
r

GFUNC   = relh
WIND    = bk07/0.7
LINE    = 1/1
TITLE   = 1/$lift/(b) RH (shaded), FRONTOGENESIS
PANEL   = ur/1//2
FINT    = 80;90 
FLINE   = 0;18;16
CTYPE   = f       
r

GFUNC   = $theta
CINT    = 4
SCALE   = 0
LINE    = 1/1/1/3
TITLE   = 0
TITLE   = 1/2/$model : $gdat1 : $date 
TEXT    = 1
CTYPE   = c


GFUNC   = frnt(thta;obs)
WIND    = bk07 
CINT    = 2;5;10;15;25 
SCALE   = 1 
LINE    = 1/-12/4/1 
r

exit
end1

gdplot << end2
GLEVEL  = $lowlev ! $lowlev
PANEL   = ll/1//2 
SCALE   = 6 ! 0 
GFUNC   = pvor($theta;geo) ! ldf($theta)
CTYPE   = f ! c  
CONTUR  = 3 
CINT    = -1;-0.5;-0.25;0;0.25;0.5 ! 1//0
LINE    = 32/1 ! 1//3 
FINT    = -1;0
FLINE   = 16;18;0
HILO    = 
HLSYM   = 
GVECT   = wnd 
WIND    = bk07 
REFVEC  = 
TITLE   = 1/$lift/(c) $lowlev1 - $lowlev2 mb: $mpvname (shaded), $instab 
GAREA   = $maparea
PROJ    = lcc 
MAP     = 10 
LATLON  = 0  
r

GLEVEL  = $uplev ! $uplev ! $frntlev 
PANEL   = lr/1//2 
SKIP    = 0/1
SCALE   = 6 ! 11! 1
GFUNC   = pvor(thta;wnd) !adv(pvor(thta;wnd);vlav(wnd))!frnt(thta;wnd)
CTYPE   = f ! c 
CONTUR  = 3  
CINT    = 0.5;1;1.5;2;3;4;5;6 ! 3;5;7;9;11 ! 5;10;15;25;50 
LINE    = 32/1 ! 1//3 ! 1/10/2 
FINT    = 0.5;2  
FLINE   = 0;18;16  
WIND    = bk07/0.7 
TITLE   = 1/$lift/(d) $uplev1 - $uplev2 mb: PV, PV-ADV; $frntlev mb: FRONT.
r

exit
end2
gpend
endif