common cdat,x,y,z,nx,ny,nz,nw,ntmax,date0,time0
;
;restore,'data/VAR3_16'
;restore,'data/VAR1_2'
restore,'data/VAR5_8'
s=size(bb) & nx=s[1] & ny=s[2] & nz=s[3]
x=2.*!pi*findgen(nx)/nx
y=2.*!pi*findgen(ny)/ny
z=2.*!pi*findgen(nz)/nz
;
kmax=nx/2
m=24
ki=0.5+findgen(m+1)
print,ki
;
ui=fltarr(nx,ny,nz,3,m)
bi=fltarr(nx,ny,nz,3,m)
ji=fltarr(nx,ny,nz,3,m)
;
for i=0,m-1 do begin
  ui(*,*,*,*,i)=fft_3d_filt_rangev(uu,ki(i),ki(i+1))
  bi(*,*,*,*,i)=fft_3d_filt_rangev(bb,ki(i),ki(i+1))
  ji(*,*,*,*,i)=fft_3d_filt_rangev(jj,ki(i),ki(i+1))
  print,'i=',i
end
;
TM_kpq=fltarr(m,m,m)
for q=0,m-1 do begin
for p=0,m-1 do begin
for k=0,m-1 do begin
  TM_kpq(k,p,q)=-mean(dot(ui(*,*,*,*,p),cross(Ji(*,*,*,*,k),Bi(*,*,*,*,q))))
  print,k,p,q
endfor
endfor
endfor
;
save,file='TM_kpq.sav',TM_kpq
END