[Nek5000-users] Velocity at the center line

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

Hello Neks

I would like to plot the value of axial-velocity at the center line of a 
pipe flow.  At present am using Visit to do that, But I would like to 
know If it is possible to get it directly from Nek.


I tried to loop on the elements to get the value
----------------------------------------------------------------------------------------------
n =  lx1*ly1*lz1*lelt
do i=1,n

         x = xm1(i,1,1,1)
         y = ym1(i,1,1,1)
         z = zm1(i,1,1,1)



            if (x.eq.0 .AND. y.eq.0 .AND. z.gt.0 .AND. z.lt.25) then
            l  = z
            up = vz(i,1,1,1)
            write(13,*) z,up
            endif

  enddo

-----------------------------------------------------------------------------------------------


but it doesn't seem to work. Could some tell me what is the correct 
variable to loop on to get the values.

Thank you,

Kamal.