restore,/ver,'vxt_96_2144-2689.save'
restore,/ver,'vyt_96_2144-2689.save'
restore,/ver,'vzt_96_2144-2689.save'
s=size(vx01t)
nx=s[1]
nz=s[2]
nt=s[3]
print,nx,nz,nt
pvx=fltarr((nx-1)/2,nt)
pvy=fltarr((nx-1)/2,nt)
pvz=fltarr((nx-1)/2,nt)
help,pvx,pvy,pvz

for t=0,nt-1 do begin $
   sizesp,reform(vx01t[*,*,t]),size=96,spectrum=tmp,wavenumbers=k & pvx[*,t]=tmp &$
   sizesp,reform(vy01t[*,*,t]),size=96,spectrum=tmp,wavenumbers=k & pvy[*,t]=tmp &$
   sizesp,reform(vz01t[*,*,t]),size=96,spectrum=tmp,wavenumbers=k & pvz[*,t]=tmp &$
   if mod(t,100) eq 0) then print,t &$
endfor

save,/ver,file='spectrumxyz.sav',k,pvx,pvy,pvz

pvxm=fltarr((nx-1)/2)
pvym=fltarr((nx-1)/2)
pvzm=fltarr((nx-1)/2)
for t=0,nt-1 do begin $
   pvxm=pvxm+pvx[*,t]/nt &$
   pvym=pvym+pvy[*,t]/nt &$
   pvzm=pvzm+pvz[*,t]/nt &$
endfor

save,/ver,file='spectrumxyz.sav',k,pvx,pvy,pvz,pvxm,pvym,pvzm