[Nek5000-users] Vorticity in entire domain

[nek5000-users at lists.mcs.anl.gov]

Hi Jurriaan,

Are you trying to compute a norm of the vorticity?

Otherwise, vorticity is a vector field and can be viewed
just like the velocity field, e.g.

       if (mod(istep,iostep).eq.0)
      $  call outpost(t(1,1,1,1,3),t(1,1,1,1,4),t(1,1,1,1,5),pr,t,'vrt')

would create a set of files with vorticity in place of velocity.

Paul

On Tue, 3 Dec 2013, nek5000-users at lists.mcs.anl.gov wrote:

> Dear Nek's
>
> I'm running the turbChannel case and trying to find the vorticity in the
> entire domain. With
>
> call comp_vort3(t(1,1,1,1,3),wo1,wo2,vx,vy,vz)
>
> I however get the 3 columns with length lx1*ly1*lz1*lelv. How do I make
> sure to gather the information from all processors?
> Any help would be greatly appreciated,
>
> Jurriaan
>