[Nek5000-users] extracting data in a smaller computational domain and in a plane
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Fri Mar 4 16:28:09 CST 2016
SL,
Yes that will work, but you will need to modify "makefile.template" in the
[..]/trunk/nek. There is a copy of the one I use in my github repo as
well. You can use it as a starting point, or just copy it (provided it
references all the files in the version of Nek you are using). My reference
template file is based on having a directory [..]/trunk/nek/ps/ where
PhilFunc.f is located (at least that's the way I do it). If you update
your makefile.template correctly then it will be included as a file that is
compiled and copied to your obj directory for each run.
Hope that helps,
Phil
On Fri, Mar 4, 2016 at 1:52 AM, <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. Re: Find out which cell a meshpoint belongs to
> (nek5000-users at lists.mcs.anl.gov)
> 2. Re: extracting data in a smaller computational domain and in
> a plane (nek5000-users at lists.mcs.anl.gov)
> 3. Re: extracting data in a smaller computational domain and in
> a plane (nek5000-users at lists.mcs.anl.gov)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Thu, 3 Mar 2016 19:52:18 +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] Find out which cell a meshpoint belongs
> to
> Message-ID:
> <mailman.16965.1457034739.24727.nek5000-users at lists.mcs.anl.gov>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Hi Lee,
>
> Take a look at nr, ns, nt, nel return arguments of findpts_setup() in
> nek/jl/findpts.c
>
> An example of usage of this routine is in intpts_setup() and intpts() of
> nek/postpro.f
>
> 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: Tuesday, February 02, 2016 4:31 AM
> To: nek5000-users at lists.mcs.anl.gov
> Subject: [Nek5000-users] Find out which cell a meshpoint belongs to
>
> Hi neks,
>
> I'm just wondering if there's currently a function which acts on a
> meshpoint and returns the cell it is closest to.
>
> I have also tried doing this manually but I need to know the corners of
> each cell, if that is stored somewhere?
>
> Otherwise, is there a way of generating a list of every gridpoint so I can
> calculate the corners using brute force.
>
> Thanks,
>
> Lee
>
>
>
>
> On Tue, Feb 2, 2016 at 12:10 AM -0800, "nek5000-users at lists.mcs.anl.gov" <
> nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov>>
> wrote:
>
> Hi
>
> The adjoint implementation in nek assumes change of the direction of
> time evolution. The goal is to use the standard nek solver with minimal
> modifications (only advection terms have to be changed). You get these
> equations from you standard one by changing the time and pressure sign
> (t -> -t; p-> -p). This of course changes your time integration limits,
> but it shouldn't be a problem. So when you have your linear direct
> simulation, to turn it into dual one it is enough to include in your
> userchk include files ADJOINT and TSTEP, and for ISTEP.eq.0 set
> IFADJOINT to true. Something like:
> subroutine userchk
>
> include 'SIZE'
> include 'TSTEP' ! ISTEP
> include 'ADJOINT' ! IFADJ
>
> if(ISTEP.eq.0) ADJOINT=.TRUE.
>
> There is no need to modify ###.rea or ###.map files. The only problem
> are boundary conditions for open flows, as they would require
> non-homogeneous bc. In this case to avoid problems we usually use
> sufficiently large domain together with zero Dirichlet bc. The last
> issue is sufficient resolution for both direct and adjoin runs. I hope
> this gives you some idea about nek implementation.
> Regards
> Adam
>
> On 2016-02-01 22:17, nek5000-users at lists.mcs.anl.gov wrote:
> > Hi Neks,
> >
> > To my knowledge, Nek5000 is able to solve adjoint equations of the form:
> > $\partial u_a/\partial t + V.\nabla u_a - u_a.(\nabla V)^T + \nabla
> > \p_a + 1/Re \nabla^2 u_a=0$ and $\nabla.u_a=0$
> > where u_a and p_a are adjoint velocity and pressure. (Slightly
> > different formulation may be seen in literature).
> >
> > Let's assume the forward problem is already solved so that $V$, i.e.
> > velocity, (and possibly $p$ as the pressure) is (are) already know.
> > Thus, the initial and boundary conditions for adjoint NS are also
> > known. How can we now solve the dual/adjoint NS problem with Nek5000?
> > Is there any example on how to modify .rea files? Specailly, the
> > "convection term of $V.\nabla u_a - u_a.(\nabla V)^T$ needs to be
> > modified compared to forward NS but am not sure how.
> > My search of mailing list entails some modifications in "perturb.f"
> > subroutine, but a little bit more clarification on that would be
> > really appreciated.
> >
> > Yours,
> > Saleh
> > _______________________________________________
> > 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/20160303/67eefdce/attachment-0001.html
> >
>
> ------------------------------
>
> Message: 2
> Date: Fri, 04 Mar 2016 19:44:17 +1100
> From: nek5000-users at lists.mcs.anl.gov
> To: nek5000-users at lists.mcs.anl.gov
> Subject: Re: [Nek5000-users] extracting data in a smaller
> computational domain and in a plane
> Message-ID:
> <mailman.16998.1457081065.24727.nek5000-users at lists.mcs.anl.gov>
> Content-Type: text/plain; charset=UTF-8; format=flowed
>
> Hi Phil,
>
> thank you very much. I sounds great!
> I suppose that I should include PhilFunc.f in trunk/tools and then
> compile it?
>
> Cheers
> SL
>
> El 04-03-2016 05:39, nek5000-users at lists.mcs.anl.gov escribi?:
> > SL,
> >
> > I have developed a set of routines for this express purpose which I am
> > happy to share. They are modified versions of the hpts routines that
> > dump out *fld files of subdomains at arbitrary resolutions.? The
> > routines (and a few others) can be found at the following link in the
> > file PhilFunc.f:
> >
> > https://github.com/psakievich/Nek5000 [3]
> >
> > To use you must define your own elements that follow the same format
> > as nek (npnts_in_x=npnts_in_y=npnts_in_z) and place them in an ascii
> > file named "ps_hpts.in [4]".? Then in nek call ps_hpts("MYF") where
> > "MYF" is the three character prefix that will be added to the fld
> > files. ??
> >
> > The format of "ps_hpts.in [4]" is also a little different.? The first
> > line should be:
> > NPoints, NElm, NPX, NPY, NPZ
> >
> > where NElm is the number of elements you are creating, and NPX is the
> > points in the x direction for the given element (akin to lx1), etc.
> > After that just specify the x,y,z coordinates like the normal hpts.in
> > [2] files.?
> >
> > A couple of additional points:
> >
> > 1) You can crank up the resolution, and use uniform grid spacing to
> > get much better iso surfaces inside the subvolume you define. (Very
> > nice way to get high quality visualizations that harness the power of
> > SEM without interpolating the whole gird onto a finer grid and
> > creating huge data sets)?
> >
> > 2) You need to ensure enough memory is allocated for the number of
> > hpts you are using through the SIZE file like normal.
> >
> > 3) You can use a resolution much higher than the lx1 specified in
> > nek.? However, if you go too high then it may create a seg fault.? I
> > haven't been able to track down the exact buffer sizes, but they are
> > dependent on LELT, and LX1 in your size file.? For example, I was to
> > dump out elements of with 32 points in each direction when my lx1 was
> > set to 10 in the SIZE file (LELT around 1000 I think?). However, when
> > I went higher I received a seg fault.? I hope to track this down at
> > some point...
> >
> > 4) ?If you want the elements to meet each other in VisIT you must
> > specify redundant boundaries between them (just like in nek) because I
> > am using the Nek5000 native output routines.? If the elements don't
> > have redundant boundaries then there will be gaps between them. I
> > would love to get around this but I haven't had time to write a
> > routine for a new output format.
> >
> > 5) I have only tested these routines in 3D, and I make no promises
> > that they are bug free. Theoretically, they should work for 2D as
> > well, but I haven't tried that.
> >
> > 6) DO NOT use the normal hpts/intpnts routines while using the ps_hpts
> > routines.? They use the same variables and it will create undefined
> > behavior.
> >
> > Cheers,
> >
> > Phil
> >
> > On Thu, Mar 3, 2016 at 11:00 AM,
> > <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 [1]
> >> 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. extracting data in a smaller computational domain and in a
> >> ? ? ? plane (nek5000-users at lists.mcs.anl.gov)
> >> ? ?2. Re: extracting data in a smaller computational domain? ?
> >> ?and in
> >> ? ? ? a plane (nek5000-users at lists.mcs.anl.gov)
> >> ? ?3. Re: extracting data in a smaller computational domain? ?
> >> ?and in
> >> ? ? ? a plane (nek5000-users at lists.mcs.anl.gov)
> >>
> >>
> > ----------------------------------------------------------------------
> >>
> >> Message: 1
> >> Date: Thu, 03 Mar 2016 04:29:05 +0100
> >> From: nek5000-users at lists.mcs.anl.gov
> >> To: Nek5000 Users <nek5000-users at lists.mcs.anl.gov>
> >> Subject: [Nek5000-users] extracting data in a smaller computational
> >> ? ? ? ? domain? and in a plane
> >> Message-ID:
> >> ? ? ? ?
> >> <mailman.16887.1456975753.24727.nek5000-users at lists.mcs.anl.gov>
> >> Content-Type: text/plain; charset=UTF-8; format=flowed
> >>
> >> Hi Neks,
> >>
> >> I am trying to do some postprocessing of the Nek5000 results.
> >>
> >> I have a 3D domain, however I will be interested in taking a
> >> smaller
> >> piece of this domain. The idea is to generate a data base of
> >> snapshots
> >> containing this smaller part of the domain. For example,
> >>
> >> if my computational domain is x in (-15, 20), y in (-20, 20), z
> >> in(0,2)
> >>
> >> I just want to take a data set of snapshots containing the velocity
> >> fields in the locations:
> >> ? ?x in (-2, 5), y in (-5, 5), z in(0,2)
> >>
> >> Is there any way to extract this data?
> >>
> >> In addition, is there any way to extract the same data but in a 2d
> >> plane? I mean:
> >> ? ?x in (-2, 5), y in (-5, 5), z=1
> >>
> >> Thank you in advance.
> >> Cheers
> >> SL
> >>
> >> ------------------------------
> >>
> >> Message: 2
> >> Date: Thu, 3 Mar 2016 03:58:48 +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] extracting data in a smaller
> >> ? ? ? ? computational domain? ? and in a plane
> >> Message-ID:
> >> ? ? ? ?
> >> <mailman.16888.1456977536.24727.nek5000-users at lists.mcs.anl.gov>
> >> Content-Type: text/plain; charset="us-ascii"
> >>
> >> Hi SL,
> >>
> >> The easiest way for a profile is to do the Lineout in VisIt and
> >> choose 'curve' format in 'Set Save Options'.
> >>
> >> More accurate way is to setup history/coordinate points in hpts.in
> >> [2] file and call hpts routine from userchk which will also work for
> >> 2D plane points.
> >>
> >> 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: Wednesday, March 02, 2016 9:29 PM
> >> To: Nek5000 Users
> >> Subject: [Nek5000-users] extracting data in a smaller computational
> >> domain? ? ? and in a plane
> >>
> >> Hi Neks,
> >>
> >> I am trying to do some postprocessing of the Nek5000 results.
> >>
> >> I have a 3D domain, however I will be interested in taking a
> >> smaller
> >> piece of this domain. The idea is to generate a data base of
> >> snapshots
> >> containing this smaller part of the domain. For example,
> >>
> >> if my computational domain is x in (-15, 20), y in (-20, 20), z
> >> in(0,2)
> >>
> >> I just want to take a data set of snapshots containing the velocity
> >> fields in the locations:
> >> ? ?x in (-2, 5), y in (-5, 5), z in(0,2)
> >>
> >> Is there any way to extract this data?
> >>
> >> In addition, is there any way to extract the same data but in a 2d
> >> plane? I mean:
> >> ? ?x in (-2, 5), y in (-5, 5), z=1
> >>
> >> Thank you in advance.
> >> Cheers
> >> SL
> >> _______________________________________________
> >> Nek5000-users mailing list
> >> Nek5000-users at lists.mcs.anl.gov
> >> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users [1]
> >>
> >> ------------------------------
> >>
> >> Message: 3
> >> Date: Thu, 03 Mar 2016 05:26:16 +0100
> >> From: nek5000-users at lists.mcs.anl.gov
> >> To: nek5000-users at lists.mcs.anl.gov
> >> Subject: Re: [Nek5000-users] extracting data in a smaller
> >> ? ? ? ? computational domain? ? and in a plane
> >> Message-ID:
> >> ? ? ? ?
> >> <mailman.16889.1456979182.24727.nek5000-users at lists.mcs.anl.gov>
> >> Content-Type: text/plain; charset=UTF-8; format=flowed
> >>
> >> Thank you Aleks,
> >>
> >> that is also very useful. However, it would be nice if I could
> >> extract a
> >> smaller 3d domain from my big computational domain. The idea is to
> >> work
> >> with the entire 3d small computational domain...
> >>
> >> I need to do the same for several time steps (velocity snapshots,
> >> changing in time), so if it would be any way to do it directly from
> >> Nek5000 it would make the things much easier...
> >>
> >> Thanks again.
> >> SL
> >>
> >> El 03-03-2016 04:58, nek5000-users at lists.mcs.anl.gov escribi?:
> >>> Hi SL,
> >>>
> >>> The easiest way for a profile is to do the Lineout in VisIt and
> >> choose
> >>> 'curve' format in 'Set Save Options'.
> >>>
> >>> More accurate way is to setup history/coordinate points in
> >> hpts.in [2]
> >>> file and call hpts routine from userchk which will also work for
> >> 2D
> >>> plane points.
> >>>
> >>> 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: Wednesday, March 02, 2016 9:29 PM
> >>> To: Nek5000 Users
> >>> Subject: [Nek5000-users] extracting data in a smaller
> >> computational
> >>> domain? ? ? and in a plane
> >>>
> >>> Hi Neks,
> >>>
> >>> I am trying to do some postprocessing of the Nek5000 results.
> >>>
> >>> I have a 3D domain, however I will be interested in taking a
> >> smaller
> >>> piece of this domain. The idea is to generate a data base of
> >> snapshots
> >>> containing this smaller part of the domain. For example,
> >>>
> >>> if my computational domain is x in (-15, 20), y in (-20, 20), z
> >> in(0,2)
> >>>
> >>> I just want to take a data set of snapshots containing the
> >> velocity
> >>> fields in the locations:
> >>> ? ? x in (-2, 5), y in (-5, 5), z in(0,2)
> >>>
> >>> Is there any way to extract this data?
> >>>
> >>> In addition, is there any way to extract the same data but in a
> >> 2d
> >>> plane? I mean:
> >>> ? ? x in (-2, 5), y in (-5, 5), z=1
> >>>
> >>> Thank you in advance.
> >>> Cheers
> >>> SL
> >>> _______________________________________________
> >>> Nek5000-users mailing list
> >>> Nek5000-users at lists.mcs.anl.gov
> >>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users [1]
> >>> _______________________________________________
> >>> Nek5000-users mailing list
> >>> Nek5000-users at lists.mcs.anl.gov
> >>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users [1]
> >>
> >> ------------------------------
> >>
> >> _______________________________________________
> >> Nek5000-users mailing list
> >> Nek5000-users at lists.mcs.anl.gov
> >> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users [1]
> >>
> >> End of Nek5000-users Digest, Vol 85, Issue 5
> >> ********************************************
> >
> > --
> >
> > Phil Sakievich
> >
> > PhD Candidate - Mechanical Engineering
> > Arizona State University - Ira A. Fulton School for Engineering of
> > Matter Transport and Energy
> > Tempe, Arizona
> >
> > Links:
> > ------
> > [1] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> > [2] http://hpts.in
> > [3] https://github.com/psakievich/Nek5000
> > [4] http://ps_hpts.in
> >
> > _______________________________________________
> > Nek5000-users mailing list
> > Nek5000-users at lists.mcs.anl.gov
> > https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
>
>
>
> ------------------------------
>
> Message: 3
> Date: Fri, 04 Mar 2016 09:52:36 +0100
> From: nek5000-users at lists.mcs.anl.gov
> To: <nek5000-users at lists.mcs.anl.gov>
> Subject: Re: [Nek5000-users] extracting data in a smaller
> computational domain and in a plane
> Message-ID:
> <mailman.16999.1457081578.24727.nek5000-users at lists.mcs.anl.gov>
> Content-Type: text/plain; charset="ISO-8859-1"
>
> The build-in function g2gi() does something similar and works with big
> data sets.
>
> Cheers,
> Stefan
>
>
> Am 3/4/16, 9:44 AM schrieb "nek5000-users-bounces at lists.mcs.anl.gov on
> behalf of nek5000-users at lists.mcs.anl.gov" unter
> <nek5000-users-bounces at lists.mcs.anl.gov on behalf of
> nek5000-users at lists.mcs.anl.gov>:
>
> >Hi Phil,
> >
> >thank you very much. I sounds great!
> >I suppose that I should include PhilFunc.f in trunk/tools and then
> >compile it?
> >
> >Cheers
> >SL
> >
> >El 04-03-2016 05:39, nek5000-users at lists.mcs.anl.gov escribi?:
> >> SL,
> >>
> >> I have developed a set of routines for this express purpose which I am
> >> happy to share. They are modified versions of the hpts routines that
> >> dump out *fld files of subdomains at arbitrary resolutions. The
> >> routines (and a few others) can be found at the following link in the
> >> file PhilFunc.f:
> >>
> >> https://github.com/psakievich/Nek5000 [3]
> >>
> >> To use you must define your own elements that follow the same format
> >> as nek (npnts_in_x=npnts_in_y=npnts_in_z) and place them in an ascii
> >> file named "ps_hpts.in [4]". Then in nek call ps_hpts("MYF") where
> >> "MYF" is the three character prefix that will be added to the fld
> >> files.
> >>
> >> The format of "ps_hpts.in [4]" is also a little different. The first
> >> line should be:
> >> NPoints, NElm, NPX, NPY, NPZ
> >>
> >> where NElm is the number of elements you are creating, and NPX is the
> >> points in the x direction for the given element (akin to lx1), etc.
> >> After that just specify the x,y,z coordinates like the normal hpts.in
> >> [2] files.
> >>
> >> A couple of additional points:
> >>
> >> 1) You can crank up the resolution, and use uniform grid spacing to
> >> get much better iso surfaces inside the subvolume you define. (Very
> >> nice way to get high quality visualizations that harness the power of
> >> SEM without interpolating the whole gird onto a finer grid and
> >> creating huge data sets)
> >>
> >> 2) You need to ensure enough memory is allocated for the number of
> >> hpts you are using through the SIZE file like normal.
> >>
> >> 3) You can use a resolution much higher than the lx1 specified in
> >> nek. However, if you go too high then it may create a seg fault. I
> >> haven't been able to track down the exact buffer sizes, but they are
> >> dependent on LELT, and LX1 in your size file. For example, I was to
> >> dump out elements of with 32 points in each direction when my lx1 was
> >> set to 10 in the SIZE file (LELT around 1000 I think?). However, when
> >> I went higher I received a seg fault. I hope to track this down at
> >> some point...
> >>
> >> 4) If you want the elements to meet each other in VisIT you must
> >> specify redundant boundaries between them (just like in nek) because I
> >> am using the Nek5000 native output routines. If the elements don't
> >> have redundant boundaries then there will be gaps between them. I
> >> would love to get around this but I haven't had time to write a
> >> routine for a new output format.
> >>
> >> 5) I have only tested these routines in 3D, and I make no promises
> >> that they are bug free. Theoretically, they should work for 2D as
> >> well, but I haven't tried that.
> >>
> >> 6) DO NOT use the normal hpts/intpnts routines while using the ps_hpts
> >> routines. They use the same variables and it will create undefined
> >> behavior.
> >>
> >> Cheers,
> >>
> >> Phil
> >>
> >> On Thu, Mar 3, 2016 at 11:00 AM,
> >> <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 [1]
> >>> 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. extracting data in a smaller computational domain and in a
> >>> plane (nek5000-users at lists.mcs.anl.gov)
> >>> 2. Re: extracting data in a smaller computational domain
> >>> and in
> >>> a plane (nek5000-users at lists.mcs.anl.gov)
> >>> 3. Re: extracting data in a smaller computational domain
> >>> and in
> >>> a plane (nek5000-users at lists.mcs.anl.gov)
> >>>
> >>>
> >> ----------------------------------------------------------------------
> >>>
> >>> Message: 1
> >>> Date: Thu, 03 Mar 2016 04:29:05 +0100
> >>> From: nek5000-users at lists.mcs.anl.gov
> >>> To: Nek5000 Users <nek5000-users at lists.mcs.anl.gov>
> >>> Subject: [Nek5000-users] extracting data in a smaller computational
> >>> domain and in a plane
> >>> Message-ID:
> >>>
> >>> <mailman.16887.1456975753.24727.nek5000-users at lists.mcs.anl.gov>
> >>> Content-Type: text/plain; charset=UTF-8; format=flowed
> >>>
> >>> Hi Neks,
> >>>
> >>> I am trying to do some postprocessing of the Nek5000 results.
> >>>
> >>> I have a 3D domain, however I will be interested in taking a
> >>> smaller
> >>> piece of this domain. The idea is to generate a data base of
> >>> snapshots
> >>> containing this smaller part of the domain. For example,
> >>>
> >>> if my computational domain is x in (-15, 20), y in (-20, 20), z
> >>> in(0,2)
> >>>
> >>> I just want to take a data set of snapshots containing the velocity
> >>> fields in the locations:
> >>> x in (-2, 5), y in (-5, 5), z in(0,2)
> >>>
> >>> Is there any way to extract this data?
> >>>
> >>> In addition, is there any way to extract the same data but in a 2d
> >>> plane? I mean:
> >>> x in (-2, 5), y in (-5, 5), z=1
> >>>
> >>> Thank you in advance.
> >>> Cheers
> >>> SL
> >>>
> >>> ------------------------------
> >>>
> >>> Message: 2
> >>> Date: Thu, 3 Mar 2016 03:58:48 +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] extracting data in a smaller
> >>> computational domain and in a plane
> >>> Message-ID:
> >>>
> >>> <mailman.16888.1456977536.24727.nek5000-users at lists.mcs.anl.gov>
> >>> Content-Type: text/plain; charset="us-ascii"
> >>>
> >>> Hi SL,
> >>>
> >>> The easiest way for a profile is to do the Lineout in VisIt and
> >>> choose 'curve' format in 'Set Save Options'.
> >>>
> >>> More accurate way is to setup history/coordinate points in hpts.in
> >>> [2] file and call hpts routine from userchk which will also work for
> >>> 2D plane points.
> >>>
> >>> 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: Wednesday, March 02, 2016 9:29 PM
> >>> To: Nek5000 Users
> >>> Subject: [Nek5000-users] extracting data in a smaller computational
> >>> domain and in a plane
> >>>
> >>> Hi Neks,
> >>>
> >>> I am trying to do some postprocessing of the Nek5000 results.
> >>>
> >>> I have a 3D domain, however I will be interested in taking a
> >>> smaller
> >>> piece of this domain. The idea is to generate a data base of
> >>> snapshots
> >>> containing this smaller part of the domain. For example,
> >>>
> >>> if my computational domain is x in (-15, 20), y in (-20, 20), z
> >>> in(0,2)
> >>>
> >>> I just want to take a data set of snapshots containing the velocity
> >>> fields in the locations:
> >>> x in (-2, 5), y in (-5, 5), z in(0,2)
> >>>
> >>> Is there any way to extract this data?
> >>>
> >>> In addition, is there any way to extract the same data but in a 2d
> >>> plane? I mean:
> >>> x in (-2, 5), y in (-5, 5), z=1
> >>>
> >>> Thank you in advance.
> >>> Cheers
> >>> SL
> >>> _______________________________________________
> >>> Nek5000-users mailing list
> >>> Nek5000-users at lists.mcs.anl.gov
> >>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users [1]
> >>>
> >>> ------------------------------
> >>>
> >>> Message: 3
> >>> Date: Thu, 03 Mar 2016 05:26:16 +0100
> >>> From: nek5000-users at lists.mcs.anl.gov
> >>> To: nek5000-users at lists.mcs.anl.gov
> >>> Subject: Re: [Nek5000-users] extracting data in a smaller
> >>> computational domain and in a plane
> >>> Message-ID:
> >>>
> >>> <mailman.16889.1456979182.24727.nek5000-users at lists.mcs.anl.gov>
> >>> Content-Type: text/plain; charset=UTF-8; format=flowed
> >>>
> >>> Thank you Aleks,
> >>>
> >>> that is also very useful. However, it would be nice if I could
> >>> extract a
> >>> smaller 3d domain from my big computational domain. The idea is to
> >>> work
> >>> with the entire 3d small computational domain...
> >>>
> >>> I need to do the same for several time steps (velocity snapshots,
> >>> changing in time), so if it would be any way to do it directly from
> >>> Nek5000 it would make the things much easier...
> >>>
> >>> Thanks again.
> >>> SL
> >>>
> >>> El 03-03-2016 04:58, nek5000-users at lists.mcs.anl.gov escribi?:
> >>>> Hi SL,
> >>>>
> >>>> The easiest way for a profile is to do the Lineout in VisIt and
> >>> choose
> >>>> 'curve' format in 'Set Save Options'.
> >>>>
> >>>> More accurate way is to setup history/coordinate points in
> >>> hpts.in [2]
> >>>> file and call hpts routine from userchk which will also work for
> >>> 2D
> >>>> plane points.
> >>>>
> >>>> 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: Wednesday, March 02, 2016 9:29 PM
> >>>> To: Nek5000 Users
> >>>> Subject: [Nek5000-users] extracting data in a smaller
> >>> computational
> >>>> domain and in a plane
> >>>>
> >>>> Hi Neks,
> >>>>
> >>>> I am trying to do some postprocessing of the Nek5000 results.
> >>>>
> >>>> I have a 3D domain, however I will be interested in taking a
> >>> smaller
> >>>> piece of this domain. The idea is to generate a data base of
> >>> snapshots
> >>>> containing this smaller part of the domain. For example,
> >>>>
> >>>> if my computational domain is x in (-15, 20), y in (-20, 20), z
> >>> in(0,2)
> >>>>
> >>>> I just want to take a data set of snapshots containing the
> >>> velocity
> >>>> fields in the locations:
> >>>> x in (-2, 5), y in (-5, 5), z in(0,2)
> >>>>
> >>>> Is there any way to extract this data?
> >>>>
> >>>> In addition, is there any way to extract the same data but in a
> >>> 2d
> >>>> plane? I mean:
> >>>> x in (-2, 5), y in (-5, 5), z=1
> >>>>
> >>>> Thank you in advance.
> >>>> Cheers
> >>>> SL
> >>>> _______________________________________________
> >>>> Nek5000-users mailing list
> >>>> Nek5000-users at lists.mcs.anl.gov
> >>>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users [1]
> >>>> _______________________________________________
> >>>> Nek5000-users mailing list
> >>>> Nek5000-users at lists.mcs.anl.gov
> >>>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users [1]
> >>>
> >>> ------------------------------
> >>>
> >>> _______________________________________________
> >>> Nek5000-users mailing list
> >>> Nek5000-users at lists.mcs.anl.gov
> >>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users [1]
> >>>
> >>> End of Nek5000-users Digest, Vol 85, Issue 5
> >>> ********************************************
> >>
> >> --
> >>
> >> Phil Sakievich
> >>
> >> PhD Candidate - Mechanical Engineering
> >> Arizona State University - Ira A. Fulton School for Engineering of
> >> Matter Transport and Energy
> >> Tempe, Arizona
> >>
> >> Links:
> >> ------
> >> [1] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> >> [2] http://hpts.in
> >> [3] https://github.com/psakievich/Nek5000
> >> [4] http://ps_hpts.in
> >>
> >> _______________________________________________
> >> 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
>
>
>
>
> ------------------------------
>
> _______________________________________________
> 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 85, Issue 7
> ********************************************
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20160304/7b3baaf8/attachment-0001.html>
More information about the Nek5000-users
mailing list