<div dir="ltr">Hi Kamal,<div><br></div><div>This doesn't work as you expect because it is very unlikely that you have node points that fall exactly at x = 0 and y = 0.  </div><div><br></div><div>The best way to go about this is to use the hpts() routine.  Instructions can be found at</div>
<div><a href="http://nek5000.mcs.anl.gov/index.php/Data_processing_example">http://nek5000.mcs.anl.gov/index.php/Data_processing_example</a><br></div><div><br></div><div>This will interpolate your solution to the points given in your "<a href="http://hpts.in">hpts.in</a>" file.</div>
<div><br></div><div>Hope this helps!</div><div><br></div><div>Josh</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Apr 2, 2014 at 4:46 AM,  <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello Neks<br>
<br>
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.<br>
<br>
<br>
I tried to loop on the elements to get the value<br>
------------------------------<u></u>------------------------------<u></u>------------------------------<u></u>----<br>
n =  lx1*ly1*lz1*lelt<br>
do i=1,n<br>
<br>
        x = xm1(i,1,1,1)<br>
        y = ym1(i,1,1,1)<br>
        z = zm1(i,1,1,1)<br>
<br>
<br>
<br>
           if (x.eq.0 .AND. y.eq.0 .AND. z.gt.0 .AND. z.lt.25) then<br>
           l  = z<br>
           up = vz(i,1,1,1)<br>
           write(13,*) z,up<br>
           endif<br>
<br>
 enddo<br>
<br>
------------------------------<u></u>------------------------------<u></u>------------------------------<u></u>-----<br>
<br>
<br>
but it doesn't seem to work. Could some tell me what is the correct variable to loop on to get the values.<br>
<br>
Thank you,<br>
<br>
Kamal.<br>
______________________________<u></u>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.<u></u>gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/<u></u>mailman/listinfo/nek5000-users</a><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>Josh Camp<br><br>"All that is necessary for the triumph of evil is that good men do nothing" -- Edmund Burke 
</div>