;$Id: pkt2304_helt2_both.pro,v 1.3 2014/02/02 16:48:17 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=7,yoffset=3
  !p.charthick=2.5 & !p.thick=2.5 & !x.thick=2.5 & !y.thick=2.5
end
;
;  mv idl.ps ~/tex/procs/bete12/fig/pkt2304_helt2_short2304pm1_kf60b_noforce_both.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.5
!p.multi=[0,2,1]
;
k0=60.
sexp=1.
!x.title='!8k!6/!8k!6!dM!n(!60!6)'
ytit='!8k!6!dM!n(!8t!6) !8E!6!dM!n!6(!8k,t!6)'
ytit0='!8k!6!dM!n(!60) !8E!6!dM!n!6(!8k,t!6)'
!x.margin=[7.1,.8]
!y.margin=[3.1,.8]
!p.charsize=1.3
tilde='!9!s!aA!n!r!6'
@parameters
syr=[5e-10,3e-6]
t0=t(0)
plot_oo,k,spec2(*,0),xr=[.01,20],yr=syr,ytit=ytit0,/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:*))
  print,i,it,t(it),dt,kM
  oplot,k/60.,spec2(*,it)/k0,li=i,thick=.5
  wait,.05
endfor
;default,sampl4,.28
;default,sampl53,.008
;xx=[.12,.4] & oplot,xx,sampl4*xx^4
;xx=[2.3,25.] & oplot,xx,sampl53/xx^2
;
;-----------------------------------------------------------------------------
!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)'
tilde='!9!s!aA!n!r!6'
;iit=indgen(n_elements(t))
default,syr,[1e-6,2e-2]
t0=t(0)
plot_oo,k,spec2(*,0),xr=[.1,50],yr=syr,/nodata,ytit=ytit
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:*))
  print,i,it,dt,kM
  oplot,k/kM,spec2(*,it)/kM^sexp,li=i,thick=.5
  wait,.05
endfor
sampl4=8e-6
sampl53=2e-6
xx=[.12,.4] & oplot,xx,sampl4*xx^4,col=122
xx=[2.3,25.] & oplot,xx,sampl53/xx^2,col=122
xyouts,.3,3e-8,siz=siz,col=122,'!9A!8k!6!u4!n'
xyouts,5.,1e-7,siz=siz,col=122,'!9A!8k!6!u-2!n'
;
END