;$Id: pkt2.pro,v 1.6 2011/12/31 08:54: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/pkt_u256d2_noforce2.eps
;  mv idl.ps ~/tex/tina/Decay/pkt_128c4_pm10_noforce.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=[2e-10,5e-6]
;
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
for it=8,n-2,120 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-7*xx^2
;xx=[1,4] & oplot,xx,1.3e-9*xx^4
xx=[1,4] & oplot,xx,1.8e-9*xx^3
;
k1=k
k1(0)=1.
k1=1./k1
stop
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,sqrt(D/E)
plot_oo,t-t(0),(II/E)
END