;$Id: pkt2304_kalp.pro,v 1.4 2015/12/03 05:27:34 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/pkt2304_hel_short2304pm1_kf60b_kalp.ps
;  mv idl.ps ~/tex/tina/Classes/fig/pkt2304_hel_short2304pm1_kf60b_kalp.ps
;  convert idl.ps ~/tex/prop/NSF/2015/machine/collapse.png
;
default,iread,0
if iread eq 0 then begin
  ;power,'hel_mag','_mag',k=k,spec1=spec1,spec2=spec2,i=n,tt=t,/noplot & print,n
  ;power,'hel_kin','hel_mag',k=k,spec1=spec1h,spec2=spec2h,i=n,tt=t,/noplot & print,n
  power,'_kin','_mag',k=k,spec1=spec1,spec2=spec2,i=n,tt=t,/noplot & print,n
  iread=1
endif
;
siz=2.4
thick=.5
thick=3.
!p.multi=0
!x.title='!8k!6!7n!6(!8t!6)'
!y.title='!8E!6!dM!n!6(!8k,t!6) and !8E!6!dK!n!6(!8k,t!6)'
!x.margin=[7.9,.5]
!y.margin=[3.1,.5]
!p.charsize=2.4
tilde='!9!s!aA!n!r!6'
@parameters
pc_read_param,/param2,o=param
;
;w=0.
default,w,.1
s=1./k^2.
s=k^2.5
;s=k^1.667 & yr=[2e-5,2e-0]
vA0=.15 & kM0=60.
s=1./(vA0^2*kM0) & default,yr,[2e-9,2e-2]
k0=15.
k0=60.
E=t
EK=t
D=t
II=t
IIK=t
nu=param.nu
eta=param.eta
help,s
;
;  compute k1
;
k1=k
k1(0)=1.
k1=1./k1
;
;  normalize with respect to kM
;
plot_oo,k/k0,s*spec2(*,0),xr=xr2,yr=yr2,/nodata
default,it3,80
nit=n_elements(iit2)
for i=0,nit-1 do begin
  it=iit2(i)
  kM=total(spec2(*,it))/total(k1*spec2(*,it))
  epsK=total(2.*nu*k^2*spec1(*,it))
  epsM=total(2.*eta*k^2*spec2(*,it))
  urms=total(2.*spec1(*,it))^.5
  vA=total(2.*spec2(*,it))^.5
  s=1./kM
  print,t(it),epsM,epsK,vA
  oplot,k/kM,s*spec2(*,it),li=0,thick=thick
  oplot,k/kM,s*spec1(*,it),li=1,thick=thick
  print,it,t(it)-t(0)
  wait,w
endfor
help,s
;
;  overplot slope
;
if nopost eq 1 then stop
@postproc3
;
help,s
END