;$Id: pkt512.pro,v 1.2 2012/11/17 15:26:44 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/Decay/fig/pkt128_128c5_noforce.ps
;  mv idl.ps ~/tex/tina/Decay/fig/pkt128_128c5_noforce3.ps
;  mv idl.ps ~/tex/tina/Decay/fig/pkt128_128c5_noforce3_pm10.ps
;  mv idl.ps ~/tex/tina/Decay/fig/pkt128_128c5_noforce3_pm100.ps
;  mv idl.ps ~/tex/tina/Decay/fig/pkt256_256c5_pm100_noforce1.ps
;  mv idl.ps ~/tex/tina/Decay/fig/pkt256_256c5_pm100.ps
;  mv idl.ps ~/tex/tina/Decay/fig/pkt256_short256pm1000a_noforce.ps
;  mv idl.ps ~/tex/tina/Decay/fig/pkt256_short256pm1000b_noforce.ps
;  mv idl.ps ~/tex/tina/Decay/fig/short256pm1b_noforce.ps
;  mv idl.ps ~/tex/tina/Decay/fig/short256pm10b_noforce.ps
;  mv idl.ps ~/tex/tina/Decay/fig/short256pm100b_noforce.ps
;  mv idl.ps ~/tex/tina/Decay/fig/short256pm1000b_noforce.ps
;  mv idl.ps ~/tex/tina/Decay/fig/short256pm1a2_warp_noforce.ps
;  mv idl.ps ~/tex/tina/Decay/fig/short256pm1a4_noforce.ps
;  mv idl.ps ~/tex/tina/Decay/fig/short256pm10a_noforce.ps
;  mv idl.ps ~/tex/tina/Decay/fig/short256pm100a_noforce.ps
;
pc_read_kf,kf
default,iread,0
if iread eq 0 then begin
  power,'_kin','_mag',k=k,spec1=spec1,spec2=spec2,i=n,tt=t,/noplot & print,n
  iread=1
endif
;
siz=1.7
!p.multi=0
!x.title='!8k!6/!8k!6!d0!n'
!y.title='!8E!6!dM!n!6(!8k,!7g!6) and !8E!6!dK!n!6(!8k,!7g!6)'
!x.margin=[9.8,.8]
!y.margin=[3.3,.8]
!p.charsize=1.7
yr=[2e-5,2e-0]
tilde='!9!s!aA!n!r!6'
@parameters
;
;w=0.
default,w,.1
s=kf/total(spec2(*,0))
E=t
EK=t
D=t
II=t
IIK=t
plot_oo,k/kf,s*spec2(*,0),xr=[1,256]/kf,yr=yr,/nodata
default,it3,80
;oplot,k,k*0+1e-2
;for it=8,n-2,it3 do begin
nit=n_elements(iit)
for i=0,nit-1 do begin
  it=iit(i)
  oplot,k/kf,s*spec2(*,it),li=0,thick=.5
  oplot,k/kf,s*spec1(*,it),li=1,thick=.5
  print,it,t(it)-t(0)
  wait,w
endfor
;
;  overplot special time in thick
;
default,itspecial,n-2
it=itspecial
oplot,k/kf,s*spec2(*,it),li=0,thick=8
oplot,k/kf,s*spec1(*,it),li=2,thick=8
;
;  overplot slope
;
@postproc
default,fact4,4e+1
default,fact2,7e-2
xx=xx0/kf & oplot,xx,fact4*xx^4,col=122,li=3
xx=xx0/kf & oplot,xx,fact2*xx^2,col=122,li=3
;xx=[12.,34.] & oplot,xx,3e-2/xx^1.67,col=122
;xx=[1.7,6] & oplot,xx,2e-7/xx^.5,col=122
;xx=[1,6] & oplot,xx,2e-8/xx ;(for PrM=100)
;
END