;$Id: pprossby.pro,v 1.2 2006/12/21 06:49:06 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=12,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
;  plot rms velocities for x,y,z components
;  using data from save file
;
;  mv idl.ps ~/tex/hydro/shearlayer/fig/pprossby_256_512a_Om02.ps
;  mv idl.ps ~/tex/procs/prague06b/fig/pprossby_256_512a_Om02.ps
;
Omega=.2
gray=210 & thick=20
!p.charsize=1.8
!x.title='!8z!6/!8d!6'
file='var.dat.uxyzrms'
restore,file
;
file='../256_512a_Om02_theta60/TTmz.sav'
restore,file
csm=sqrt(2./3.*TTm)
grav=1.
gam=5./3.
Hp=csm^2/(gam*grav)
Omega=.2
;
!y.range=[.3,30.]
plot_io,zzz,2.*Omega*Hp/sqrt(uxrms^2-uxm^2),li=1
;
;  plot range of unstable region
;
loadct,0
z1=min(zzz)
z2=max(zzz)
polyfill,[z1,0.,0.,z1],[[1,1]*!y.range(0),[1,1]*!y.range(1)],col=gray
polyfill,[z2,1.,1.,z2],[[1,1]*!y.range(0),[1,1]*!y.range(1)],col=gray
;
;  plot extra dash-dotted line and plot thick line at the bottom
;
oplot,[z1,0.,0.,z1],[[1,1]*!y.range(0),[1,1]*!y.range(1)],li=3
oplot,[z2,1.,1.,z2],[[1,1]*!y.range(0),[1,1]*!y.range(1)],li=3
oplot,[z1,0.],[1,1]*!y.range(0),thick=thick
oplot,[1.,z2],[1,1]*!y.range(0),thick=thick
;
plot_io,zzz,2.*Omega*Hp/sqrt(uxrms^2-uxm^2)/csm,li=1,/noerase
oplot,zzz,2.*Omega*Hp/sqrt(uyrms^2-uym^2)/csm
oplot,zzz,2.*Omega*Hp/sqrt(uzrms^2-uzm^2)/csm,li=2
;
END