Hi Josh,<br><br>The routine that you want is attached to this email or you can find it here:<br><a href="http://www.mcs.anl.gov/~jwalker/interp_v.f">http://www.mcs.anl.gov/~jwalker/interp_v.f</a><br><br>You give it a packed list of the points that you want to interpolate  <br>
<br>x(1),y(1),z(1),....,x(n),y(n),z(n)<br><br>along with the number of points, n, and it gives you back the velocities at each point<br><br>u(1),v(1),w(1),...,u(n),v(n),w(n).<br><br>Make sure to change lpart in SIZE to the number of points you want.<br>
<br><br>Best,<br>Justin<br><br><br><div class="gmail_quote">On Thu, May 19, 2011 at 3:09 PM, Paul Fischer <span dir="ltr"><<a href="mailto:fischer@mcs.anl.gov">fischer@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
Justin,<br>
<br>
Can you post or send the latest interp_uvw() routine that<br>
you used in your hemi case?<br>
<br>
I guess we need to get the hemi example into the repo.<br>
<br>
Thanks,<br>
<br>
Paul<br>
<br>
<br>
---------- Forwarded message ----------<br>
Date: Thu, 19 May 2011 15:06:05 -0500 (CDT)<br>
From: <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a><br>
To: <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a><br>
Subject: Re: [Nek5000-users] Average in arbitrary direction<br>
<br>
Hi Josh,<br>
<br>
yes, there are 3 intpts* routines from postpro.f that you can use. I have an example but I need to update it for the current repo version of this routines -- I'll send it to you once I modify and test it.<br>
<br>
Best,<br>
Aleks<br>
<br>
<br>
----- Original Message -----<br>
From: <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a><br>
To: <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a><br>
Sent: Thursday, May 19, 2011 9:21:03 AM<br>
Subject: Re: [Nek5000-users] Average in arbitrary direction<br>
<br>
Hello Paul,<br>
<br>
This sound like exactly what I need.  I didn't find anything like<br>
interp() while using grep in the repo--are you referring to the<br>
routines found in postpro.f ??<br>
<br>
Thanks,<br>
<br>
Josh<br>
<br>
On Wed, May 18, 2011 at 3:59 PM,  <<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<br>
Hi Josh,<br>
<br>
You can use the interp() routines --- let me check if we have<br>
examples in the repo.<br>
<br>
These are somewhat expensive, so if you're also integrating in<br>
time (even over just 10 steps), it would pay to only interpolate<br>
the time-averaged data, rather than interpolating every timestep<br>
and then averaging the profiles.   However, I can understand that<br>
there are times you'd want the full time signal.<br>
<br>
The best integration formulae for periodic signals (e.g., in your theta<br>
direction and perhaps in z, if you have periodicity)<br>
is to simply sum the values, sampled uniformly in space (z or<br>
theta), and dived by n, the number of samples.   You do need<br>
to sample with sufficient density to capture all the wiggles ---<br>
another reason that time averaging has value.<br>
<br>
Paul<br>
<br>
<br>
<br>
On Wed, 18 May 2011, <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="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Hello Neks,<br>
<br>
I have looked at the turbChannel example for an example on planar<br>
averaging.<br>
<br>
However, I was wondering if there was any advice  on how to do this if<br>
the elements aren't exactly "nice"... (i.e., the elements may have<br>
arbitrary rotations)<br>
<br>
Specifically, I have a pipe flow, and I am wanting to do an average in<br>
the z (axial) direction.  My initial thoughts were to interpolate my<br>
values to set R, theta components, and then average this along z.  I'm<br>
stuck on the best way to do either one of these operations, however<br>
(especially because I would like to use an integration rule to compute<br>
the average in the z direction).<br>
<br>
If anyone could point me in the right direction, that would be a big<br>
help.  Thanks!<br>
<br>
--<br>
Josh Camp<br>
<br>
"All that is necessary for the triumph of evil is that good men do<br>
nothing" -- Edmund Burke<br>
_______________________________________________<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>
<br>
</blockquote>
_______________________________________________<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>
<br>
</blockquote>
<br>
<br>
<br>
-- <br>
Josh Camp<br>
<br>
"All that is necessary for the triumph of evil is that good men do<br>
nothing" -- Edmund Burke<br>
_______________________________________________<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>
_______________________________________________<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></blockquote></div><br>