;$Id: pkt.pro,v 1.10 2012/01/28 21:40:34 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
;
;  mv idl.ps ~/tex/tina/Kazantsev/fig/pkt_256d_noforce2.eps
;  mv idl.ps ~/tex/tina/Decay/fig/128h2_noforce1_abisko.eps
;
power,'_kin','_mag',k=k,spec1=spec1,spec2=spec2,i=n,tt=t,/noplot & print,n
;
!p.multi=0
!x.title='!8k!6/!8k!6!d1!n'
!y.title='!8E!6!dM!n!6(!8k,t!6) and !8E!6!dK!n!6(!8k,t!6)'
!y.margin=[3.3,.3]
!p.charsize=1.7
yr=[1e-9,1e-1]
yr=[1e-7,1.6e-3]
@parameters
;
;w=0.
default,w,.1
s=1./k^2.
s=k^1.667
s=k^2.5
s=1.
E=t
EK=t
D=t
II=t
IIK=t
plot_oo,k,s*spec2(*,0),xr=[1,128],yr=yr
default,it3,80
;oplot,k,k*0+1e-2
;for it=8,n-2,it3 do begin
for it=8,it2,it3 do begin
  oplot,k,s*spec2(*,it),li=0,thick=.5
  oplot,k,s*spec1(*,it),li=1,thick=.5
  print,t(it),t(it)-t(0)
  wait,w
endfor
;--it=n-2 & oplot,k,s*spec2(*,it),li=0,thick=8
;--it=n-2 & oplot,k,s*spec1(*,it),li=2,thick=8
;
;  overplot slope
;
;--xx=[5,20] & oplot,xx,.01/xx^2.
;--xx=[4,14] & oplot,xx,.00005/xx^0.5
;xx=[1,4] & oplot,xx,3e-5*xx^2
;xx=[1,3] & oplot,xx,1e-5*xx^4
;--xx=[1,3] & oplot,xx,2e-5*xx^3
;
k1=k
k1(0)=1.
k1=1./k1
;
;wait,.5
;for it=0,n-2 do IIK(it)=total(k1*spec1(*,it))
;for it=0,n-2 do II(it)=total(k1*spec2(*,it))
;for it=0,n-2 do D(it)=total(k^2*spec2(*,it))
;for it=0,n-2 do E(it)=total(spec2(*,it))
;for it=0,n-2 do EK(it)=total(spec1(*,it))
;plot_oo,t-t(0),(II/E)
END