;$Id: pkt1.pro,v 1.3 2010/04/13 21:20:02 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/pkt1.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.6
;
;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
yr=[1e-9,1e-1]
yr=[2e-8,2e-3]
iit=[1,2,5,10,20,50,100,200,500]
iit=[1,3,10,30,100,500]
nnt=n_elements(iit)
plot_oo,k,s*spec2(*,1),xr=[1,128],yr=yr
;oplot,k,k*0+1e-2
for jt=0,nnt-1 do begin
  it=iit(jt)
  oplot,k,s*spec2(*,it),li=0
  oplot,k,s*spec1(*,it),li=2
  print,t(it)
  wait,w
endfor
it=n-2 & oplot,k,s*spec2(*,it),li=2,thick=5
it=n-2 & oplot,k,s*spec1(*,it),li=1,thick=5
;
k1=k
k1(0)=1.
k1=1./k1
xx=[1,20] & oplot,xx,xx^4*1e-8,li=3
xx=[1,20] & oplot,xx,xx^2*1e-5,li=3
;
END