;$Id: peig.pro,v 1.6 2014/07/27 22:19:14 brandenb Exp $
;
;  Procedure to compute rate of strain tensor,
;  as well as its eigenvalues and eigenvectors.
;
file='VAR3_16'
restore,'data/'+file
restore,'data/e'+file
;
;  dimensions
;
s=size(uij)
nx=s[1]
ny=s[2]
nz=s[3]
;
bb=reform(bb,1L*nx*ny*nz,3)
oo=reform(oo,1L*nx*ny*nz,3)
uij=reform(uij,1L*nx*ny*nz,3,3)
sij=uij
;
;  compute sij
;
for i=0,2 do begin
for j=0,2 do begin
  sij(*,i,j)=.5*(uij(*,i,j)+uij(*,j,i))
end
end
;
;  the last row, vecs(*,*,3), contains the eigenvalues
;
;stop
o=1./sqrt(dot2(oo))
b=1./sqrt(dot2(bb))
vecs=eigvec3_arr(sij)
coso1=o*dot(oo,vecs(*,*,0)) & pdf,fmin=0,fmax=1,n=500,abs(coso1),xo1,yo1
coso2=o*dot(oo,vecs(*,*,1)) & pdf,fmin=0,fmax=1,n=500,abs(coso2),xo2,yo2
coso3=o*dot(oo,vecs(*,*,2)) & pdf,fmin=0,fmax=1,n=500,abs(coso3),xo3,yo3
cosb1=b*dot(bb,vecs(*,*,0)) & pdf,fmin=0,fmax=1,n=500,abs(cosb1),xb1,yb1
cosb2=b*dot(bb,vecs(*,*,1)) & pdf,fmin=0,fmax=1,n=500,abs(cosb2),xb2,yb2
cosb3=b*dot(bb,vecs(*,*,2)) & pdf,fmin=0,fmax=1,n=500,abs(cosb3),xb3,yb3
save,file='align_'+file+'.sav',xo1,yo1,xo2,yo2,xo3,yo3,xb1,yb1,xb2,yb2,xb3,yb3
;
;  pdfs of eigenvalues
;
pdf,n=500,reform(vecs(*,0,3)),xe1,ye1
pdf,n=500,reform(vecs(*,1,3)),xe2,ye2
pdf,n=500,reform(vecs(*,2,3)),xe3,ye3
save,file='eigenvals_'+file+'.sav',xe1,ye1,xe2,ye2,xe3,ye3
;
END