Hi Paul,<br>
<br>
first of all thanks for the quick reply!<br>
<br>
<br>
<blockquote cite="mid:Pine.LNX.4.64.1203102028490.28269@v8.mcs.anl.gov" type="cite">I guess my first question is the following.  When you say you
  <br>
want to
  <br>

  <br>
"add or subtract this profile which is given on the global nodes
  <br>
to each point of the field"
  <br>

  <br>
it seems we need to be talking in the context of some software
  <br>
that is going to read the field in and then process the data.
  <br>

  <br>
At present, I know of only 3 codes to do that: nek, VisIt, and 
postnek (though there have been custom, one-off, codes written
  <br>
in the past).
  <br>
</blockquote>
I want to do the whole post-processing in nek. With VisIt I only 
visualize my data generated with nek (the post-processed data as well as
 the raw simulation data). With postnek I haven't worked until now and I
 think postnek is not what I'm looking for.<br>
<br>
<br>
<blockquote cite="mid:Pine.LNX.4.64.1203102028490.28269@v8.mcs.anl.gov" type="cite">Second question is --- What do you want to do with this field
  <br>
after you have subtracted the mean profile?   Do you want to
  <br>
write it to disk, visualize it, use it in a calculation, 
analyze it with more statistics, or....?
  <br>
</blockquote>
I want to analyze the statistics, calculate higher order moments of the 
fluctuation field (to compare it with the channel data of Robert D. 
Moser, John Kim and Nagi N. Mansour,<br>
Physics of Fluids, vol 11(4), for example the skewness or the flatness 
of the field) or the gradients of the fluctuation field du'_i/dx_j to 
get (fluctuation based) \omega'_k or the (fluctuation based) strain rate
 tensor s'_ij = 1/2*(du'_i/dx_j + du'_j/dx_i).<br>
These are only few examples on what I would like to do...<br>
<br>
But for all these things I need the fluctuation field and that means u' = u - <u>.<br>
In this context <u> would be (for example if I don't have the 
temporal avg_all() means) the spatial (y-) mean over the full channel 
domain.<br>
<br>
When I do my production run with the code I can obtain the temporal 
averages by calling avg_all() every timestep and then write it to disk, thats clear to me.<br>
<br>
But my question aims on how to get the fluctuation field without having 
the temporal averages. Based on only a single velocity filed (as 
described above).<br>
<br>
<br>
<blockquote cite="mid:Pine.LNX.4.64.1203102028490.28269@v8.mcs.anl.gov" type="cite">In any case, it seems likely that processing in nek is the
  <br>
way to go because it has the ability to read the data, write
  <br>
the data, and readily compute statistics.   In this mode, you
  <br>
simply are using nek as a post-processor and not for timestepping.
  <br>

  <br>
Please advise if that's the path you wish to take, or if you
  <br>
had something else in mind.
  <br>
</blockquote>
Yes, thats the way I want  to do my post-processing (and already do it right now...).<br>
<br>
Regards<br>
Jan F.