[Nek5000-users] Nek5000-users Digest, Vol 105, Issue 32

Thu Nov 30 17:05:30 CST 2017

Dear Ketan,

Many thanks for your reply. Basically I want to get the averaged velocity,
rms velocities and shear stress where I got two questions.

1. Once the subroutines were writen in usr file, Is calling avg_all still
needed? In other words, the post-processing is done after longtime running,
and usr file is changed instead of calling avg_all, isn't it?

2. After trying the subroutines you provided, there isnt no files generated
in the directory. Can you give me more infomation about the averaging
process.

Best wishes,

ZJ

On 30 November 2017 at 18:00, <nek5000-users-request at lists.mcs.anl.gov>
wrote:

> Send Nek5000-users mailing list submissions to
>         nek5000-users at lists.mcs.anl.gov
>
> To subscribe or unsubscribe via the World Wide Web, visit
>         https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> or, via email, send a message with subject or body 'help' to
>         nek5000-users-request at lists.mcs.anl.gov
>
> You can reach the person managing the list at
>         nek5000-users-owner at lists.mcs.anl.gov
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of Nek5000-users digest..."
>
>
> Today's Topics:
>
>    1. tools in new version (nek5000-users at lists.mcs.anl.gov)
>    2. Re: avg_all averaging (nek5000-users at lists.mcs.anl.gov)
>    3. Re: tools in new version (nek5000-users at lists.mcs.anl.gov)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Thu, 30 Nov 2017 06:26:10 +0300
> From: nek5000-users at lists.mcs.anl.gov
> To: nek5000-users at lists.mcs.anl.gov
> Subject: [Nek5000-users] tools in new version
> Message-ID:
>         <mailman.3666.1512012374.19102.nek5000-users at lists.mcs.anl.gov>
> Content-Type: text/plain; charset="utf-8"
>
> Hi, Neks!
>
> I am trying to create a complex geometry and use nektools. I use nekmerge2
> script to merge two parts of my mesh, and then I want to use reatore2 and
> genmap. In previous version the same script worked find, but now
> (~/Nek5000/bin/) reatore2 gives a note:?
>
> Note: The following floating-point exceptions are signalling:
> IEEE_INVALID_FLAG IEEE_DENORMAL
>
> and after that genmap gives a error:
>
> ABORT: SELF-CHK 1 2 4 0
> Try to tighten the mesh tolerance!
>
> I have tested this script with old tools (on another computer) and it's
> still works fine. Do you have any ideas why does this problem occur?
>
> Best regards,
> Mark
>
>
>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/
> attachments/20171130/1450ceba/attachment-0001.html>
>
> ------------------------------
>
> Message: 2
> Date: Thu, 30 Nov 2017 05:01:23 +0000
> From: nek5000-users at lists.mcs.anl.gov
> To: "nek5000-users at lists.mcs.anl.gov"
>         <nek5000-users at lists.mcs.anl.gov>
> Subject: Re: [Nek5000-users] avg_all averaging
> Message-ID:
>         <mailman.3672.1512018086.19102.nek5000-users at lists.mcs.anl.gov>
> Content-Type: text/plain; charset="utf-8"
>
> Hi ZJ,
>
> A quick search of the archives with keyword ?avg_all? will give you answer
> to related questions:
> https://lists.mcs.anl.gov/mailman/mmsearch/nek5000-
> users?config=nek5000-users.htsearch&restrict=&exclude=&
> method=and&format=short&sort=score&words=avg_all
>
> avg_all is a collection of time-averaged information in the avg, rms and
> rm2 files. The avg files has time-averaged u,v and w; rms has time-averaged
> u^2, v^2 and w^2; and rm2 has uv, vw, and uw. The avg_all routines is
> located in navier5.f in case you want to get more information.
>
> Typically the first step in post processing is to use the routine
> ?auto_averager? to combine all the files output from avg_all. autoaverager
> takes as input the name of a file (example: fname1.list) which has a list
> of all the files that you want to combine (example: the contents of
> fname1.list will be: avgcasename0.f00001 avgcasename0.f00002 ? Note- there
> should be only 1 filename per line in the fname1.list file). You can use
> the auto_averager routine to combine all your avgcasename0.f00* files and
> output them to a single file or do more post-processing with it. The
> quantities available in avg, rms and rm2 files are sufficient to get the
> actual urms, vrms velocities etc? Additionally, once you time-average the
> data using auto_averager, you can use space averaging routines such as
> z_avg etc? to average quantities in homogeneous direction of your domain.
>
> It would also help if you specify what exactly are you trying to get from
> these avg_all files i.e. are you looking for rms velocities or something
> else?
>
> Here is a snippet of my code that I have used in the past to post-process
> these files:
>       call auto_averager(fname1)
>       call copy(uk,vx)
>       call copy(vk,vy)
>       call copy(wk,vz)
>       call col2(uk,uk,lt)
>       call col2(vk,vk,lt)
>       call col2(wk,wk,lt)
>       call auto_averager(fname2)
>       call sub2(vx,uk,lt)
>       call sub2(vy,vk,lt)
>       call sub2(vz,wk,lt)
>       call z_avg(vxa,vx,gs_avg_hndl,nelxy,ifld)
>       call z_avg(vya,vy,gs_avg_hndl,nelxy,ifld)
>       call z_avg(vza,vz,gs_avg_hndl,nelxy,ifld)
>       call outpost(vxa,vya,vza,vxa,vya,'   ')
>
> Here fname1 has the list of all avgcasename0.f00* files and fname2 has the
> list of all the rmscasename0.f00* files.
>
>
> Ketan
>
>
>
>
>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/
> attachments/20171130/b19f2925/attachment-0001.html>
>
> ------------------------------
>
> Message: 3
> Date: Thu, 30 Nov 2017 11:11:43 +0100
> From: nek5000-users at lists.mcs.anl.gov
> To: nek5000-users at lists.mcs.anl.gov     <nek5000-users at lists.mcs.anl.gov>
> Subject: Re: [Nek5000-users] tools in new version
> Message-ID:
>         <mailman.3678.1512036715.19102.nek5000-users at lists.mcs.anl.gov>
> Content-Type: text/plain; charset=utf-8
>
> Can you try again with the latest version on GitHub (just submitted a
> potential fix).
>
> -----Original message-----
> > From:nek5000-users at lists.mcs.anl.gov <nek5000-users at lists.mcs.anl.gov>
> > Sent: Thursday 30th November 2017 4:27
> > To: nek5000-users at lists.mcs.anl.gov
> > Subject: [Nek5000-users] tools in new version
> >
> > Hi, Neks!
> >
> > I am trying to create a complex geometry and use nektools. I use
> nekmerge2 script to merge two parts of my mesh, and then I want to use
> reatore2 and genmap. In previous version the same script worked find, but
> now (~/Nek5000/bin/) reatore2 gives a note:?
> >
> > Note: The following floating-point exceptions are signalling:
> IEEE_INVALID_FLAG IEEE_DENORMAL
> >
> > and after that genmap gives a error:
> >
> > ABORT: SELF-CHK 1 2 4 0
> > Try to tighten the mesh tolerance!
> >
> > I have tested this script with old tools (on another computer) and it's
> still works fine. Do you have any ideas why does this problem occur?
> >
> > Best regards,
> > Mark
> >
> >
> >
> > _______________________________________________
> > Nek5000-users mailing list
> > Nek5000-users at lists.mcs.anl.gov
> > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
>
> ------------------------------
>
> Subject: Digest Footer
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
>
> ------------------------------
>
> End of Nek5000-users Digest, Vol 105, Issue 32
> **********************************************
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20171130/dd7a68db/attachment.html>