[Nek5000-users] avg_all and visnek
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Mon Jul 25 00:38:18 CDT 2016
Hi Marco,
It seems like ifxyo overwrites itself to .false. after istep .eq. 0
conditional statement. If you have the flexibility of running the code
again (if it runs quickly), you can set ifxyo = .true. just before the
avg_all() call. Otherwise, you can run, nek5000 is post-processing mode as
follows.
subroutine userchk
include 'SIZE'
include 'TOTAL'
include 'ZPER' ! Needed to communicate nelx,nely,nelz to my_stats
include 'RESTART'
character*80 filename
filename = 'avg.....fld..' <original filename of average file>
call load_fld(filename)
ifxyo = .true.
call outpost(vx,vy,vz,pr,t,'avx')
call exitt
return
end
On Fri, Jul 22, 2016 at 12:48 PM, <nek5000-users at lists.mcs.anl.gov> wrote:
> 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 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 [
> 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:* Friday, July 22, 2016 1:52 PM
> *To:* 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 wrote:
>
> firsttimestep
>
>
> _______________________________________________
> 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
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20160724/637d861f/attachment.html>
More information about the Nek5000-users
mailing list