common cdat,x,y,z,nx,ny,nz,nw,ntmax,date0,time0
;
restore,'data/VAR3_16'
restore,'data/eVAR3_16'
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)
oi=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))
  oi(*,*,*,*,i)=fft_3d_filt_rangev(oo,ki(i),ki(i+1))
  print,'i=',i
end
;
SM_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
  SM_kpq(k,p,q)=-mean(dot(ui(*,*,*,*,p),cross(ui(*,*,*,*,k),oi(*,*,*,*,q))))
  print,k,p,q
endfor
endfor
endfor
;
save,file='SM_kpq.sav',SM_kpq
END