;$Id: pkt1_15t.pro,v 1.1 2012/04/30 19:58:08 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/Inflation/fig/pkt1_15t_u512_km1_a.eps
;  mv idl.ps ~/tex/tina/Inflation/fig/pkt1_15t_u512_km1_b.eps
;  mv idl.ps ~/tex/tina/Inflation/fig/pkt1_15t_u512_km1_c.eps
;  mv idl.ps ~/tex/tina/Inflation/fig/pkt1_15t_u512_km1_d.eps
;  mv idl.ps ~/tex/tina/Inflation/fig/pkt1_15t_u512_km1_e.eps
;
@parameters
power,'_kin','_mag',k=k,spec1=spec1,spec2=spec2,i=n,tt=t,/noplot & print,n
;
pc_read_kf,kf
!p.multi=0
!x.title='!8k!6/!8k!6!d0!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
ik=15
yr=[2e-7,5e-5]
;iit=[1,2,5,10,20,50,100,200,500]
;it=[1,3,10,30,100,200]
nnt=n_elements(iit)
plot_io,t,s*spec2(ik,*),yr=yr
oplot,t,s*spec1(ik,*),li=2
;
k1=k
k1(0)=1.
k1=1./k1
siz=1.6
;xx=[5,70] & oplot,xx/kf,xx^1*5e-9,li=3
;xyouts,1,6e-8,'!9A!8k!6',siz=siz
;xx=[1,3] & oplot,xx/kf,xx^2*2e-6,li=3
;xyouts,.082,6e-6,'!9A!8k!6!u2!n',siz=siz
;
END