; $Id: pspec_comp.pro,v 1.2 2004/04/11 16:11:47 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=14,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
;  comparing decay laws for different runs
;  mv idl.ps ~/tex/hydro/hyper/fig/pspec_comp.ps
;
power,v1='u',spec1=spec,k=k,i=i,/noplot,tt=tt
ns=(size(spec))[2]
del_itt=10
;
!x.title='!8k!6 (!7m!6!d3!n!8t!6)!u1/6!n'
!y.title='!8E!6(!8k,t!6) !8k!6!u5/3!n!7e!6!u-2/3!n (!8t!6/!7s!6)!u5/4!n'
!p.charsize=1.8
;
;  the value of eps512 must be consistent with that in hydro.pro
;  in tex/hydro/hyper/idl
;
eps512=0.000271407/1.7
nu=5e-13
fac=k^1.667/eps512^.667
i1=40 & i2=75 & i3=5
dt=tt-tt(0)
dt_exponent=1.25
tau=24. ;(the same as in urms^2/urms0^2 = [(t-t0)/tau]^1.25
;
plot_oo,k*(nu*dt(i1))^0.1667,fac*spec(*,i1)*(dt(i1)/tau)^dt_exponent,xr=[.018,3.4],yr=[.1,5]
for i=i1,i2,i3 do begin
  print,i,tt(i)-tt(0)
  oplot,k*(nu*dt(i))^0.1667,fac*spec(*,i)*(dt(i)/tau)^dt_exponent
  wait,.04
end
;
siz=2.
xx=[.06,.4] & oplot,xx,2.0/xx^.03,li=2
xyouts,.34,1.4,'!9A!8k!6!u-0.03',siz=siz
;
END