#!/bin/csh
#
# This script creates a 4-panel graphic with important QG parameters.
#
# Author: Michael Morris
# Date: November 1, 2005
#
cd /webhost/mmorris/http
rm -f qg_4panel.gif
cd /home/mmorris/gempak
#
setenv DISPLAY :1.0
set FRONTS = `ls -t /data/gempak/nwx/hpc/fronts/*.front | head -1`
source /usr/local/nawips/Gemenviron
#
#
#
/usr/local/nawips/bin/linux/gdplot <<END1

pan=1
proj = lcc
garea= uslcc
satfil = 
gdattim = last
gvect=
gvc=pres
gle=500
gfun=hght!tmpc!adv(vor(wnd),wnd))
gdf=/data/gempak/hrs/cur_upa
dev=gif|/webhost/mmorris/http/qg_4panel.gif|1500;1200
cint=60!5!1/2
lin=6/1/2!2/2/2!5/1/2
map=31/9/1
\$MAPFIL=default
tit=31/-3/~ 500 mb Analysis!31/-2/HGHT (blue) TMPC (red) REL VORT ADV (yel)
clrb=
cle=y
text=0.8
scale = 999
latlon=0
hilo = 
hls=0
cont=2
skip=0
filter=y
fint=
fline=
ctype=
l
r

exit
END1

/usr/local/nawips/bin/linux/gdcntr << END1
r last
pan=1
cle=no
tit=31/-2/HGHT (blue) TMPC (red) REL VORT ADV (yel)
l
r

exit
END1
#
# Now make the second panel (Quadrant I) with a 700 mb analysis.
#
/usr/local/nawips/bin/linux/gdplot <<END6
r last
pan=2
ctype=c
gvc=pres
gfu=hght!tmpc!adv(tmpc, wnd)
lin=6/1/2!2/2/2!17/1/2
fli=
fint=
cint=30!3!2/10
gle=700
tit=31/-2/700 MB HGHT (blu) TEMP (red) WAA (orange)
l
r

exit
END6

#
# Now plot 850 on the 3rd panel.
#
/usr/local/nawips/bin/linux/gdplot << END10
r last
cle=no
pan=3
gle=850
gvc=pres
ctype=c
gfu=hght!adv(tmpc, wnd)!tmpc
tit=31/-1/ 850 mb HGHT (blue) WAA (orange) Temp (red)
cint=30!3/12!3
line=6/1/2!18/1/2!2/2/2
l
r

exit
END10

/usr/local/nawips/bin/linux/gdplot <<END7
r last
gdattim=last
gdfi=/data/gempak/hrs/cur_upa
pan=4
gvc=pres
gle=700
gfu=div(qvec(thta, wnd))
map=31/9/1
mapfil=default
ctype=c
tit=31/-1/SFC Frontogenesis (red)  700 mb Q-Vector Convergence (yel)
cint=2//0
lin=5/1/2
l
r

exit
END7

/usr/local/nawips/bin/linux/gdcntr <<END7

r last
gdattim=last
gdfi=/data/gempak/hrs/cur_sfc
pan=4
gvc=none
gle=0
gfu=frnt(stha, wnd))
map=31/9/1
mapfil=default
ctype=c
cint=3/3
lin=2/1/2
l
r

exit
END7

/usr/local/nawips/bin/linux/gpend

exit