Lines Matching refs:vorticity
97 real*8, allocatable, dimension(:,:) :: ybar, strain, vorticity local
177 allocate(vorticity(nshg,5))
531 & call computeVort( vorticity, GradV,strain)
556 call write_field(myrank,'a','vorticity',9,vorticity,
632 & vorticity,yphbar,rerr,irank2ybar,irank2yphbar)
757 deallocate(strain,vorticity)
1726 subroutine computeVort( vorticity, GradV,strain) argument
1729 real*8 gradV(nshg,nsdsq), strain(nshg,6), vorticity(nshg,5) local
1732 vorticity(:,1) = GradV(:,8)-GradV(:,6) !omega_x
1733 vorticity(:,2) = GradV(:,3)-GradV(:,7) !omega_y
1734 vorticity(:,3) = GradV(:,4)-GradV(:,2) !omega_z
1735 vorticity(:,4) = sqrt( vorticity(:,1)*vorticity(:,1)
1736 & + vorticity(:,2)*vorticity(:,2)
1737 & + vorticity(:,3)*vorticity(:,3) )
1746 vorticity(:,5) = 0.25*( vorticity(:,4)*vorticity(:,4) !Q
1842 & vorticity,yphbar,rerr,irank2ybar,irank2yphbar) argument
1844 real*8 ybar(nshg,irank2yphbar),yold(nshg,ndof),vorticity(nshg,5) local
1911 ybar(:,14) = tfact*vorticity(:,1) +
1913 ybar(:,15) = tfact*vorticity(:,2) +
1915 ybar(:,16) = tfact*vorticity(:,3) +
1917 ybar(:,17) = tfact*vorticity(:,4) +
2000 & tfactphase*vorticity(:,1)
2003 & tfactphase*vorticity(:,2)
2006 & tfactphase*vorticity(:,3)
2009 & tfactphase*vorticity(:,4)