;$Id: pkft_new2.pro,v 1.3 2017/02/09 00:19:47 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=16,ysize=20,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
;  compute various length scales from the spectra and
;  calculate the instantaneous scaling exponents.
;
default,texdir3,'lengthscale'
default,ihydro,0
@parameters
default,iread,0
if iread eq 0 then begin
  if ihydro eq 1 then begin
    power,'_kin','hel_kin',k=k,spec1=spec1,spec2=spec2,i=n,tt=t,/noplot
    print,n
  endif else begin
    power,'_kin','_mag',k=k,spec1=spec1,spec2=spec2,i=n,tt=t,/noplot
    print,n
  endelse
  pc_read_param,o=param,/param2
  iread=1
endif
;
!p.charsize=1.8
!x.margin=[4.8,.5]
!y.margin=[3.2,.5]
!p.multi=[0,3,3]
fo2='(f5.2)'
cwd,run
si2=1.4
;
if ihydro eq 0 then begin
  eta=param.eta
  kIM=t
  kTM=t
  k4M=t
  kM=t
  EM=t
endif
nu=param.nu
kIK=t
kTK=t
k4K=t
kK=t
EK=t
nt=n_elements(t)
;
;  compute 1/k
;
k1=k & k1(0)=1. & k1=1./k1 & k1(0)=0.
;
;  loop over all times
;
E=t & for it=0,nt-1 do E(it)=total(spec1(*,it))
epsK=t & for it=0,nt-1 do epsK(it)=2.*nu*total(k^2*spec1(*,it))
for it=0,nt-1 do begin
  if ihydro eq 0 then begin
    kIM(it)=1./(total(spec2(*,it)*k1)/total(spec2(*,it)))
    kTM(it)=(total(spec2(*,it)*k^2)/total(spec2(*,it)))^.5
    k4M(it)=(total(spec2(*,it)*k^4)/total(spec2(*,it)))^.25
    kM(it)=(total(2*eta*spec2(*,it)*k^2)/eta^3)^.25
    EM(it)=total(spec2(*,it))
  endif
  kIK(it)=1./(total(spec1(*,it)*k1)/total(spec1(*,it)))
  kTK(it)=(total(spec1(*,it)*k^2)/total(spec1(*,it)))^.5
  k4K(it)=(total(spec1(*,it)*k^4)/total(spec1(*,it)))^.25
  kK(it)=(total(2*nu*spec1(*,it)*k^2)/nu^3)^.25
  EK(it)=total(spec1(*,it))
endfor
t0=t(0)
dt=t-t0+.4
cwd,run
;
!p.multi=[0,1,2]
!x.title='!6'
plot_oo,t,kTK,/nodata,yr=[1.,3e3]
oplot,t,kIM,col=122
oplot,t,kIK,col=55
;
oplot,t,kTM,col=122,li=2
oplot,t,kTK,col=55,li=2
;
oplot,t,k4M,col=122,li=1
oplot,t,k4K,col=55,li=1
;
oplot,t,kM,col=122,li=3
oplot,t,kK,col=55,li=3
;
!x.title='!8t!6/!8t!6!d0!n'
plot_oo,t,kTK,/nodata,yr=[1,400]
oplot,t,kIM/kIM,col=122
oplot,t,kIK/kIM,col=55
;
oplot,t,kTM/kIM,col=122,li=2
oplot,t,kTK/kIM,col=55,li=2
;
oplot,t,k4M/kIM,col=122,li=1
oplot,t,k4K/kIM,col=55,li=1
;
oplot,t,kM/kIM,col=122,li=3
oplot,t,kK/kIM,col=55,li=3
;
print,'mv idl.ps ~/tex/tina/'+texdir3+'/fig/pkft_new2_'+run+'.ps'
!p.multi=0
END