[Nek5000-users] 1d profile

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

Hi,
I guess the problem is the line
eg = gllel(e)

You should use
eg = lglel(e)
It would translate your local element number to the global one.
Regards
Adam

On 09/02/17 08:34, nek5000-users at lists.mcs.anl.gov wrote:
> 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
>
>
> _______________________________________________
> 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/20170209/23c6d786/attachment.html>