[Nek5000-users] avg_all and visnek

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

Here is what I have in my usrchk():

      subroutine userchk
      include 'SIZE'
      include 'TOTAL'
      include 'ZPER' ! Needed to communicate nelx,nely,nelz to my_stats

      common /my_recycle/ nslab,n_in_slab

      ifxyo = .true.
      if (istep.eq.0) call outpost(vx,vy,vz,pr,t,'   ')
      ifxyo = .false.

c      call set_inflow(nslab,n_in_slab)   ! set up inflow BCs

      nelx = n_in_slab
      nely = 1
      nelz = nelgv/(nelx*nely)

c     if (istep.gt.0) call my_stats

      call avg_all()

      return
      end

Marco
On Jul 22, 2016, at 3:34 PM, nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov> wrote:

Interesting... Check whether there is a statement

ifxyo = .false. in usrchk() which overrides .rea specification (typically after the first dump to save the disk space)

Aleks
________________________________
From: nek5000-users-bounces at lists.mcs.anl.gov<mailto:nek5000-users-bounces at lists.mcs.anl.gov> [nek5000-users-bounces at lists.mcs.anl.gov<mailto:nek5000-users-bounces at lists.mcs.anl.gov>] on behalf of nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov> [nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov>]
Sent: Friday, July 22, 2016 1:52 PM
To: nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov>
Subject: Re: [Nek5000-users] avg_all and visnek

Hi Aleks,

thanks for your help.

I did set CORDINATES to T in the block ‘SPECIFICATIONS FOLLOW’. It was originally set to F.

I copied my instantaneous dump file into avg*0.f00000 and then substituted firsttimestep:1 to firsttimestep:0. It did work.

This is the first time I have to do such a thing when using avg_all() even though I have used avg_all() in the past. Is there anything that has changed in the code?

Thanks,
Marco
On Jul 22, 2016, at 2:17 PM, nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov> wrote:

firsttimestep

_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov<mailto:Nek5000-users at lists.mcs.anl.gov>
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20160722/3cca4869/attachment-0001.html>