Hi Aleks,<br>
<br>
thank you for posting your routines. They are certainly very useful!<br>
But my point is to find an easy and quick way to check a single velocity
 field. That can usually be done by taking the means over the full box 
and subtracting them (as described in the first mail) from the 
instantaneous field.<br>
I don't wanna do my full postprocessing that way. It's just to have a quick look on the data (in particular the fluctuations)...<br>
Do you think that is possible?<br>
<br>
Thanks!<br>
Jan F.<br><br><div class="gmail_quote">Am 9. März 2012 15:34 schrieb  <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Jan,<br>
<br>
Here are routines that I sometimes use that call avg_all for the list of a dump binary files given in a.list<br>
<br>
<a href="http://www.mcs.anl.gov/%7Eobabko/avg_fld.usr" target="_blank">http://www.mcs.anl.gov/~obabko/avg_fld.usr</a><br>
<br>
Use<br>
<br>
     call avg_fld_all<br>
<br>
in post-processing mode of Nek5000 (nsteps=0 in .rea) with a created a.list file of your binary .f files (for .fld files uncomment line with p66 & p67 = 4)<br>
<br>
Best.<br>
Aleks<br>
<div><div class="h5"><br>
<br>
<br>
----- Original Message -----<br>
From: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a><br>
To: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a><br>
Sent: Wednesday, March 7, 2012 4:10:44 AM<br>
Subject: [Nek5000-users] Add/subtract 1d-mean profile<br>
<br>
<br>
Hi NEKs,<br>
<br>
I'm new to NEK and as a first exercise I'm calculating DNS of turbulent channel flow. So far everything is fine and works as is should.<br>
<br>
Now i wanna do some postprocessing with the (binary) files I dumped out. Let's say I wanna look at a single file:<br>
<br>
Is it possible to calc the 1d mean of u over the whole channel domain (for example with "call planar_average_s") then take this 1d mean profile and subtract it from u at every point to get u' (=> u'(i,j,k,e) = u(i,j,k,e) - <u>(j))?<br>

<br>
As I didn't dump the averages with "avg_all()" the hint given in <a href="https://lists.mcs.anl.gov/mailman/htdig/nek5000-users/2010-May/000451.html" target="_blank">https://lists.mcs.anl.gov/mailman/htdig/nek5000-users/2010-May/000451.html</a> in this case is not applicable.<br>

<br>
Thanks!<br>
Jan Frielinghausen<br>
<br>
</div></div>_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">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">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>
</blockquote></div><br>