[Nek5000-users] auto_averager()

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

Hello all,
I am trying to take an average from previously averaged output 
"avgblah.f" files using the following calls in userchk in postprocessing 
mode:

...
call auto_averager(filename)
call outpost(vx,vy,vz,pr,t,'AVG')
call exitti('exiting AVG routine$',icount)
...

While I got the correct velocity field average, the average for the 
temperature field is not correct, as it is out of min/max range in the 
output files. I went through the "auto_averager" routine, but couldn't 
find anything wrong with that.

So I was wondering if anyone has any thought on this issue that can help me.

Thanks,
Mohsen