[Nek5000-users] Finding the maximum

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

Hi Hank and Paul ,

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 !

Regards
Shriram


On 29 October 2010 15:04, <nek5000-users at lists.mcs.anl.gov> wrote:

>
> Hi Shriram,
>
> You can do this in userchk as follows:
>
>
>      n = nx1*ny1*nz1*nelv
>      vmax = -9e20
>      do i=1,n
>         x=xm1(i,1,1,1)
>         y=ym1(i,1,1,1)
>         z=zm1(i,1,1,1)
>         if (x0.le.x.and.x.le.x1 .and.
>     $       y0.le.y.and.y.le.y1 .and.
>     $       z0.le.z.and.z.le.z1 )
>     $       vmax=max(vmax,eddyv(i,1,1,1))
>      enddo
>      vmax = glmax(vmax,1)
>
> where x0,x1,...,z1 are discrimators of your choice
> and eddyv, say, is the array in question.
>
> Paul
>
>
>
>
> On Fri, 29 Oct 2010, nek5000-users at lists.mcs.anl.gov wrote:
>
>  Hi,
>>
>> For a computation with eddy viscosity, I would like to find the maximum
>> eddy
>> viscosity within a particular region of my domain and also the location of
>> the max/min value.
>> Is there a routine in nek to give limits on x,y,z and then find max of a
>> particular variable ? glmax() finds the maximum in the entire domain . . .
>>
>> Regards
>> Shriram
>>
>>  _______________________________________________
> 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/20101029/4449a6d5/attachment.html>