;$Id: psij_proj.pro,v 1.6 2014/07/27 22:30:07 brandenb Exp $
;
;  prepare aa where the average is removed
;  data were read first with peig.pro
;
aa2=aa
x=findgen(144)/144.*2*!pi
;
vecs=reform(vecs,nx,ny,nz,3,4)
aa2=dot(aa,float(vecs(*,*,*,*,1)))^2
;
;  projections against sij
;
u=1./mean(dot2(uu))
o=1./mean(dot2(oo))
a=1./mean(dot2(aa))
a2=1./mean(aa2)
b=1./mean(dot2(bb))
j=1./mean(dot2(jj))
;
suu=fltarr(1L*nx*ny*nz)
soo=fltarr(1L*nx*ny*nz)
saa=fltarr(1L*nx*ny*nz)
;saa2=fltarr(1L*nx*ny*nz)
sbb=fltarr(1L*nx*ny*nz)
sjj=fltarr(1L*nx*ny*nz)
for i=0,2 do begin
for j=0,2 do begin
  suu=suu+u*sij(*,i,j)*uu(*,i)*uu(*,j)
  soo=soo+o*sij(*,i,j)*oo(*,i)*oo(*,j)
  saa=saa+a*sij(*,i,j)*aa(*,i)*aa(*,j)
  ;saa2=saa2+a*sij(*,i,j)*aa2(*,i)*aa2(*,j)
  sbb=sbb+b*sij(*,i,j)*bb(*,i)*bb(*,j)
  sjj=sjj+j*sij(*,i,j)*jj(*,i)*jj(*,j)
  print,i,j
end
end
saa2=a2*float(vecs(*,*,*,1,3))*aa2
;
suu_mean=mean(suu) & suu_rms=sqrt(mean(suu^2))
soo_mean=mean(soo) & soo_rms=sqrt(mean(soo^2))
saa_mean=mean(saa) & saa_rms=sqrt(mean(saa^2))
saa2_mean=mean(saa2) & saa2_rms=sqrt(mean(saa2^2))
sbb_mean=mean(sbb) & sbb_rms=sqrt(mean(sbb^2))
sjj_mean=mean(sjj) & sjj_rms=sqrt(mean(sjj^2))
;
save,file='psij_proj.sav',suu_mean,suu_rms,soo_mean,soo_rms,saa_mean,saa_rms,saa2_mean,saa2_rms,sbb_mean,sbb_rms,sjj_mean,sjj_rms
;
END