[Nek5000-users] auto_averager()

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

Hi Mohsen,

I would suggest checking the range of the temperature inputs as they are read
in auto_averager and then checking the range of the averages.

Something like:

      tmax = glmax(t,n)
      tmin  = glmin(t,n)
      if (nid.eq.0) write(6,*) k,tmin,tmax,' tmin'

would produce those values...  then you can do something
similar for the averaged result.

Perhaps the failure occurs during the read process ... or.. etc.

Paul
 
________________________________________
From: nek5000-users-bounces at lists.mcs.anl.gov [nek5000-users-bounces at lists.mcs.anl.gov] on behalf of nek5000-users at lists.mcs.anl.gov [nek5000-users at lists.mcs.anl.gov]
Sent: Thursday, December 04, 2014 8:03 PM
To: nek5000-users at lists.mcs.anl.gov
Subject: [Nek5000-users] auto_averager()

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
_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users