common cdat,x,y,z,nx,ny,nz,nw,ntmax,date0,time0
;
restore,'data/VAR3_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=12
ki=0.5+findgen(m+1)
print,ki
;
ui=fltarr(nx,ny,nz,3,m)
bi=fltarr(nx,ny,nz,3,m)
;
j1=fft_3d_filt_rangev(jj,ki(0),ki(0+1))
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))
  print,'i=',i
end
;
END