[Nek5000-users] Average solution along lines

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

Thank you Paul.

The routine should work for non-uniform cartesian mesh also right, if the
dx2 is calculated. This is done in y_average function but commented out for
dy2 ???

Regards
praveen

On Wed, Nov 13, 2013 at 5:10 PM, <nek5000-users at lists.mcs.anl.gov> wrote:

> Hi Praveen,
>
> Yes, you are correct in assuming that y_average computes
>
> ua(x) = 1/L int_0^L u(x,y) dy
>
> under the assumptions of: 2D, uniform element length in y,
> and tensor-product array of elements (e.g., as generated
> by genbox).
>
> I wrote an x_average() routine a few weeks back and have
> tested it.  I've just now added this to navier5.f --- hopefully it
> will work for you.
>
> Be certain you set lelx, lely, and lelz in SIZE, and likewise
> nelx,nely,nelz in (say) usrdat(), where you must also add:
>
>      include 'ZPER'
>
> Of course, nelz=lelz=1 in your case.
>
> I
>
> Paul
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20131113/e7a7d7c5/attachment.html>