;$Id: pvort.pro,v 1.1 2012/04/30 19:58:55 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=18,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
;  mv idl.ps ~/tex/tina/Inflation/fig/pvort512d.ps
;
!p.charsize=1.7
!x.margin=[8.8,.5]
!y.margin=[3.2,.8]
!p.multi=[0,1,2]
pc_read_kf,kf
;pc_read_param,/param2,o=param
;kf=2./param.width_ff
@parameters
cwd,run
;
!p.multi=[0,1,2]
;
default,it1,400
urms=mean(ts.urms(it1:*))
;
!x.title='!8t!6 !8c!6!ds!n!8k!6!d1!n'
!y.title='!7x!6/!8uk!6!df!n  and  !13E!6!dM!n/!13E!6!dK!n'
plot_oo,ts.t,ts.orms/(urms*kf),xr=[.1,1000],yr=[.03,2.],li=0
oplot,ts.t,1*(ts.brms/urms)^2,li=2
;
!x.title='!13E!6!dM!n/!13E!6!dK!n'
!y.title='!7x!6/!8uk!6!df!n'
plot_oo,(ts.brms/urms)^2,(ts.orms/(urms*kf)),ps=1,xr=[.01,1.],yr=[.01,2.]
;
B0=0.08
xx=grange(.03,.31,30,/log) & oplot,xx,.1*(xx/B0)^1.2,col=122
xx=grange(.15,.5,30,/log) & oplot,xx,.04*(xx/B0)^2.2,col=122
END