[Nek5000-users] 1d profile

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

Hi, neks.
I am trying to get the 1d profile. vx(y) in turbChannel, for or example.
Here the code below:

real myvx(ny1,nely)

do e=1,nelv

  eg = gllel(e)
  call get_exyz(ex,ey,ez,eg,nelx,nely,nelz)

  if (ex.eq.1 .and. ez.eq.1) then
    do i=1,ny1
      myvx(i,ey) = vx(1,i,1,e)
    enddo
  endif

enddo

But after this, myvx contains incorrect value, that not looks like really
vx(y).
What's wrong?

Vatslav
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170209/4b3d29f8/attachment.html>