;$Id: pkt512_helt2.pro,v 1.4 2012/03/01 22:32:46 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/short512pm1a3_noforce_comp.ps
;  mv idl.ps ~/tex/tina/ObsSigna/short512pm1a3_noforce_tdep.ps
;  mv idl.ps ~/tex/tina/ObsSigna/pkt512_helt2_short512pm1a3_noforce.ps
;  mv idl.ps ~/tex/tina/ObsSigna/pkt512_helt2_short512pm1b3_noforce.ps
;
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!dM!n(!8t!6)'
default,ytit,'!8k!6!dM!n(!8t!6) !8E!6!dM!n!6(!8k,t!6)'
!x.margin=[9.4,.8]
!y.margin=[3.3,.8]
!p.charsize=1.7
tilde='!9!s!aA!n!r!6'
@parameters
;iit=indgen(n_elements(t))
default,syr,[1e-6,2e-2]
;
t0=t(0)
plot_oo,k,spec2(*,0),xr=[.1,50],yr=syr,ytit=ytit,/nodata
nit=n_elements(iit)
for i=0,nit-1 do begin
  it=iit(i)
  dt=t(it)-t0
  kM=total(spec2(1:*,it))/total(spec2(1:*,it)/k(1:*))
  ;kM=38./dt^.6667
  print,i,it,dt,kM
  oplot,k/kM,spec2(*,it)/kM^sexp,li=i,thick=.5
  wait,.05
endfor
;
default,sampl4,.2
default,sampl53,.008
xx=[.12,.4] & oplot,xx,sampl4*xx^4
xx=[.8,3.] & oplot,xx,sampl53/xx^1.6667
;
END