<html>
<head>
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<div class="moz-cite-prefix">Hi,<br>
I guess the problem is the line<br>
eg = gllel(e)<br>
<br>
You should use <br>
eg = lglel(e)<br>
It would translate your local element number to the global one.<br>
Regards<br>
Adam<br>
<br>
On 09/02/17 08:34, <a class="moz-txt-link-abbreviated" href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a> wrote:<br>
</div>
<blockquote
cite="mid:mailman.230.1486644024.2967.nek5000-users@lists.mcs.anl.gov"
type="cite">
<div dir="ltr">
<div>Hi, neks.<br>
</div>
I am trying to get the 1d profile. vx(y) in turbChannel, for or
example. Here the code below:<br>
<br>
real myvx(ny1,nely)<br>
<br>
do e=1,nelv<br>
<br>
eg = gllel(e)<br>
call get_exyz(ex,ey,ez,eg,nelx,nely,nelz)<br>
<br>
if (ex.eq.1 .and. ez.eq.1) then<br>
do i=1,ny1<br>
myvx(i,ey) = vx(1,i,1,e)<br>
enddo<br>
endif<br>
<br>
enddo<br>
<div><br>
</div>
<div>But after this, myvx contains incorrect value, that not
looks like really vx(y).<br>
</div>
<div>What's wrong?<br>
</div>
<div><br clear="all">
<div>
<div>
<div>
<div class="gmail_signature">
<div dir="ltr">
<div>
<div dir="ltr">Vatslav</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
Nek5000-users mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a>
<a class="moz-txt-link-freetext" href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a>
</pre>
</blockquote>
<p><br>
</p>
</body>
</html>