[Nek5000-users] auto_averager()

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

Thanks Paul for your suggestion.

The failure occurs due to common block /scrns/ with array "ta" in  
auto_averager. This block is also shared with the "restart" routine,  
in which new values are apparently reassigned to it after each call.  
It is solved by either defining a new block or moving "ta" to block / 
scruz/.

Mohsen

On 4-Dec-14, at 10:31 PM, <nek5000-users at lists.mcs.anl.gov> wrote:

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