; $Id: pslice.pro,v 1.1 2016/05/15 02:27:06 brandenb Exp $
if !d.name eq 'PS' then begin
  device,filename='ko.ps',xsize=18,ysize=14,yoffset=3,/color
  ;device,xsize=18,ysize=14,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
; run cpar.pro after commenting the write command before using this
;
;  mv ko.ps ~/pencil-code/nishant/Helioseismology/Notes-1/figures/ko_e256a.ps
;  mv ko.ps ~/pencil-code/nishant/Helioseismology/Notes-1/figures/ko_e256b.ps
;  mv ko.ps ~/pencil-code/nishant/Helioseismology/Notes-1/figures/ko_e256e.ps
;  mv ko.ps ~/pencil-code/nishant/Helioseismology/Notes-1/figures/ko_GSMa1-1.ps
;  mv ko.ps ~/pencil-code/nishant/Helioseismology/Notes-1/figures/ko_GSMa1-2.ps
;  mv ko.ps ~/pencil-code/nishant/Helioseismology/Notes-1/figures/ko_GSMa2.ps
;  mv ko.ps ~/pencil-code/nishant/Helioseismology/Notes-1/figures/ko_GSMa2-vert.ps
;  mv ko.ps ~/pencil-code/nishant/Helioseismology/Notes-1/figures/ko_GSMa3.ps
;  mv ko.ps ~/pencil-code/nishant/Helioseismology/Notes-1/figures/ko_GSMa4.ps
;  mv ko.ps ~/pencil-code/nishant/Helioseismology/Notes-1/figures/ko_GSMa6.ps
;  mv ko.ps ~/pencil-code/nishant/Helioseismology/Notes-1/figures/ko_GSMb1oblB-1.ps
;  mv ko.ps ~/pencil-code/nishant/Helioseismology/Notes-1/figures/ko_GSMb1oblB-2.ps
;
!p.charsize=3.0
!x.margin=[5.5,1.5]
!y.margin=[3.2,.5]
@parameters
;
restore,'ko_diagram.sav'
print,minmax(fff)
lev=grange(-11,-5,14)
lev=grange(lev1,lev2,14) ;(7)
;
om_th=sqrt(g*kkk*del_rho)
; Eqns (33) and (34) of Miles and Roberts (1992)
;ck2=(rho_up/(rho_up+rho_down))*(va_down*va_down)
;ck2=((rho_up*va_max*va_max)+(rho_down*va_down*va_down))/(rho_up+rho_down)
;ck2=(rho_up/(rho_up+rho_down))*(va_max*va_max)
;ck2=2.*((rho_up/(rho_up+rho_down))*(va_up*va_up))
ck2=2.*((rho_up/(rho_up+rho_down))*(vax_up*vax_up))
om_th_mag=sqrt(ck2*kkk*kkk+g*kkk*del_rho)
!x.title='!8k!dx!n!6'
!y.title='!7x!6'
;loadct,5 ; our default color table
contour,/fil,lev=lev,fff,kkk,ooo,xr=[0,15],yr=[0,10]
oplot,kkk,kkk,li=5,col=255
oplot,kkk,csmax*kkk,li=1,col=255
;oplot,kkk,sqrt(g*kkk),li=1,col=255
oplot,kkk,om_th,li=4,col=255
oplot,kkk,om_th_mag,li=2,col=255
;oplot,kkk,va_down*kkk,li=3,col=255 
;oplot,kkk,va_max*kkk,li=4,col=255 
;loadct,2  ;loads the green color table
;oplot,kkk,om_th_mag,li=2,col=30  ; produces green for loadct=2
;loadct,5
;oplot,kkk,va_down*kkk,li=3,col=255 ; not present
;oplot,kkk,va_max*kkk,li=4,col=255 ; seems to be present
END