[Nek5000-users] lambda2 in turbChannel

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

You can run nek5000 in postprocessing mode (set NSTEPS = 0). Then in
userchk, you can do the following

     INCLUDE ’SIZE’
     INCLUDE ’TOTAL’
     INCLUDE ’RESTART’

     character*80 filename(9999)
     real l2(lx1,ly1,lz1,lelv)

     ntot   = nx1*ny1*nz1*nelv

        call load_ fld(filename)  ! select a filename fld file of your
choice.
        call lambdat2(l2)

        call outpost(.vx,vy,vz,pr,l2,'lam')
     enddo

If your file is say turbChannel.fld015, the postprocessing will generate,
lamturbChannel.fld01. You can type in visnek lamturbChannel in the terminal
and then open visit and visualize the temperature variable. That will your
l2.

Best Regards,
Tanmoy

On Thu, Jun 9, 2016 at 5:13 PM, <nek5000-users at lists.mcs.anl.gov> wrote:

> Hello,
>
>
> In the userchk subroutine for the turbChannel example, lambda2 was
> computed. How can I visualize the isosurfaces (streamwise structures) in
> visit?
>
>
>       if(mod(istep,iostep).eq.0 .and. istep.gt.0) then
>         ! compute lambda2 vortex
>         call lambda2(t(1,1,1,1,2))
>         ! compute vorticity --- 3 components returned in t(...)
>         if(ldimt.ge.5) call comp_vort3(t(1,1,1,1,3),wo1,wo2,vx,vy,vz)
>       endif
>
>
> Thank you,
>
> Virgilio A. Gomez
>
>
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20160613/06de4b04/attachment.html>