<div dir="ltr">Hi Hank and Paul ,<br><br>Both the method works. I am looking at time dependent max values of eddy viscosity and so have this routine by Paul implemented in nek. The one on VisIt helps me in getting the location of max values. Thanks !<br>

<br clear="all">Regards<br>Shriram<br>
<br><br><div class="gmail_quote">On 29 October 2010 15:04,  <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">

<br>
Hi Shriram,<br>
<br>
You can do this in userchk as follows:<br>
<br>
<br>
      n = nx1*ny1*nz1*nelv<br>
      vmax = -9e20<br>
      do i=1,n<br>
         x=xm1(i,1,1,1)<br>
         y=ym1(i,1,1,1)<br>
         z=zm1(i,1,1,1)<br>
         if (x0.le.x.and.x.le.x1 .and.<br>
     $       y0.le.y.and.y.le.y1 .and.<br>
     $       z0.le.z.and.z.le.z1 )<br>
     $       vmax=max(vmax,eddyv(i,1,1,1))<br>
      enddo<br>
      vmax = glmax(vmax,1)<br>
<br>
where x0,x1,...,z1 are discrimators of your choice<br>
and eddyv, say, is the array in question.<br>
<br>
Paul<div><div></div><div class="h5"><br>
<br>
<br>
<br>
On Fri, 29 Oct 2010, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a> wrote:<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi,<br>
<br>
For a computation with eddy viscosity, I would like to find the maximum eddy<br>
viscosity within a particular region of my domain and also the location of<br>
the max/min value.<br>
Is there a routine in nek to give limits on x,y,z and then find max of a<br>
particular variable ? glmax() finds the maximum in the entire domain . . .<br>
<br>
Regards<br>
Shriram<br>
<br>
</blockquote></div></div><div><div></div><div class="h5">
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</div></div></blockquote></div><br></div>