<div dir="ltr">SL,<div><br></div><div>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:</div><div><br></div><div><a href="https://github.com/psakievich/Nek5000">https://github.com/psakievich/Nek5000</a></div><div><br></div><div>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 "<a href="http://ps_hpts.in">ps_hpts.in</a>". Then in nek call ps_hpts("MYF") where "MYF" is the three character prefix that will be added to the fld files. </div><div><br></div><div>The format of "<a href="http://ps_hpts.in">ps_hpts.in</a>" is also a little different. The first line should be:</div><div>NPoints, NElm, NPX, NPY, NPZ</div><div><br></div><div>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 <a href="http://hpts.in">hpts.in</a> files. </div><div><br></div><div>A couple of additional points:</div><div><br></div><div>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) </div><div><br></div><div>2) You need to ensure enough memory is allocated for the number of hpts you are using through the SIZE file like normal.</div><div><br></div><div>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...</div><div><br></div><div>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.</div><div><br></div><div>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.</div><div><br></div><div>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.</div><div><br></div><div>Cheers,</div><div><br></div><div>Phil<br><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Mar 3, 2016 at 11:00 AM, <span dir="ltr"><<a href="mailto:nek5000-users-request@lists.mcs.anl.gov" target="_blank">nek5000-users-request@lists.mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Send Nek5000-users mailing list submissions to<br>
<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:nek5000-users-request@lists.mcs.anl.gov">nek5000-users-request@lists.mcs.anl.gov</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:nek5000-users-owner@lists.mcs.anl.gov">nek5000-users-owner@lists.mcs.anl.gov</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of Nek5000-users digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. extracting data in a smaller computational domain and in a<br>
plane (<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>)<br>
2. Re: extracting data in a smaller computational domain and in<br>
a plane (<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>)<br>
3. Re: extracting data in a smaller computational domain and in<br>
a plane (<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Thu, 03 Mar 2016 04:29:05 +0100<br>
From: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a><br>
To: Nek5000 Users <<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>><br>
Subject: [Nek5000-users] extracting data in a smaller computational<br>
domain and in a plane<br>
Message-ID:<br>
<<a href="mailto:mailman.16887.1456975753.24727.nek5000-users@lists.mcs.anl.gov">mailman.16887.1456975753.24727.nek5000-users@lists.mcs.anl.gov</a>><br>
Content-Type: text/plain; charset=UTF-8; format=flowed<br>
<br>
Hi Neks,<br>
<br>
I am trying to do some postprocessing of the Nek5000 results.<br>
<br>
I have a 3D domain, however I will be interested in taking a smaller<br>
piece of this domain. The idea is to generate a data base of snapshots<br>
containing this smaller part of the domain. For example,<br>
<br>
if my computational domain is x in (-15, 20), y in (-20, 20), z in(0,2)<br>
<br>
I just want to take a data set of snapshots containing the velocity<br>
fields in the locations:<br>
x in (-2, 5), y in (-5, 5), z in(0,2)<br>
<br>
Is there any way to extract this data?<br>
<br>
In addition, is there any way to extract the same data but in a 2d<br>
plane? I mean:<br>
x in (-2, 5), y in (-5, 5), z=1<br>
<br>
Thank you in advance.<br>
Cheers<br>
SL<br>
<br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Thu, 3 Mar 2016 03:58:48 +0000<br>
From: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a><br>
To: "<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>"<br>
<<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>><br>
Subject: Re: [Nek5000-users] extracting data in a smaller<br>
computational domain and in a plane<br>
Message-ID:<br>
<<a href="mailto:mailman.16888.1456977536.24727.nek5000-users@lists.mcs.anl.gov">mailman.16888.1456977536.24727.nek5000-users@lists.mcs.anl.gov</a>><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
Hi SL,<br>
<br>
The easiest way for a profile is to do the Lineout in VisIt and choose 'curve' format in 'Set Save Options'.<br>
<br>
More accurate way is to setup history/coordinate points in <a href="http://hpts.in" rel="noreferrer" target="_blank">hpts.in</a> file and call hpts routine from userchk which will also work for 2D plane points.<br>
<br>
Aleks<br>
________________________________________<br>
From: <a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov">nek5000-users-bounces@lists.mcs.anl.gov</a> [<a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov">nek5000-users-bounces@lists.mcs.anl.gov</a>] on behalf of <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a> [<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>]<br>
Sent: Wednesday, March 02, 2016 9:29 PM<br>
To: Nek5000 Users<br>
Subject: [Nek5000-users] extracting data in a smaller computational domain and in a plane<br>
<br>
Hi Neks,<br>
<br>
I am trying to do some postprocessing of the Nek5000 results.<br>
<br>
I have a 3D domain, however I will be interested in taking a smaller<br>
piece of this domain. The idea is to generate a data base of snapshots<br>
containing this smaller part of the domain. For example,<br>
<br>
if my computational domain is x in (-15, 20), y in (-20, 20), z in(0,2)<br>
<br>
I just want to take a data set of snapshots containing the velocity<br>
fields in the locations:<br>
x in (-2, 5), y in (-5, 5), z in(0,2)<br>
<br>
Is there any way to extract this data?<br>
<br>
In addition, is there any way to extract the same data but in a 2d<br>
plane? I mean:<br>
x in (-2, 5), y in (-5, 5), z=1<br>
<br>
Thank you in advance.<br>
Cheers<br>
SL<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Thu, 03 Mar 2016 05:26:16 +0100<br>
From: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a><br>
To: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a><br>
Subject: Re: [Nek5000-users] extracting data in a smaller<br>
computational domain and in a plane<br>
Message-ID:<br>
<<a href="mailto:mailman.16889.1456979182.24727.nek5000-users@lists.mcs.anl.gov">mailman.16889.1456979182.24727.nek5000-users@lists.mcs.anl.gov</a>><br>
Content-Type: text/plain; charset=UTF-8; format=flowed<br>
<br>
Thank you Aleks,<br>
<br>
that is also very useful. However, it would be nice if I could extract a<br>
smaller 3d domain from my big computational domain. The idea is to work<br>
with the entire 3d small computational domain...<br>
<br>
I need to do the same for several time steps (velocity snapshots,<br>
changing in time), so if it would be any way to do it directly from<br>
Nek5000 it would make the things much easier...<br>
<br>
Thanks again.<br>
SL<br>
<br>
El 03-03-2016 04:58, <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a> escribi?:<br>
> Hi SL,<br>
><br>
> The easiest way for a profile is to do the Lineout in VisIt and choose<br>
> 'curve' format in 'Set Save Options'.<br>
><br>
> More accurate way is to setup history/coordinate points in <a href="http://hpts.in" rel="noreferrer" target="_blank">hpts.in</a><br>
> file and call hpts routine from userchk which will also work for 2D<br>
> plane points.<br>
><br>
> Aleks<br>
> ________________________________________<br>
> From: <a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov">nek5000-users-bounces@lists.mcs.anl.gov</a><br>
> [<a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov">nek5000-users-bounces@lists.mcs.anl.gov</a>] on behalf of<br>
> <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a> [<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>]<br>
> Sent: Wednesday, March 02, 2016 9:29 PM<br>
> To: Nek5000 Users<br>
> Subject: [Nek5000-users] extracting data in a smaller computational<br>
> domain and in a plane<br>
><br>
> Hi Neks,<br>
><br>
> I am trying to do some postprocessing of the Nek5000 results.<br>
><br>
> I have a 3D domain, however I will be interested in taking a smaller<br>
> piece of this domain. The idea is to generate a data base of snapshots<br>
> containing this smaller part of the domain. For example,<br>
><br>
> if my computational domain is x in (-15, 20), y in (-20, 20), z in(0,2)<br>
><br>
> I just want to take a data set of snapshots containing the velocity<br>
> fields in the locations:<br>
> x in (-2, 5), y in (-5, 5), z in(0,2)<br>
><br>
> Is there any way to extract this data?<br>
><br>
> In addition, is there any way to extract the same data but in a 2d<br>
> plane? I mean:<br>
> x in (-2, 5), y in (-5, 5), z=1<br>
><br>
> Thank you in advance.<br>
> Cheers<br>
> SL<br>
> _______________________________________________<br>
> Nek5000-users mailing list<br>
> <a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
> <a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
> _______________________________________________<br>
> Nek5000-users mailing list<br>
> <a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
> <a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br>
<br>
<br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br>
<br>
End of Nek5000-users Digest, Vol 85, Issue 5<br>
********************************************<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Phil Sakievich</div><div><br></div><div><div>PhD Candidate - Mechanical Engineering</div><div>Arizona State University - Ira A. Fulton School for Engineering of Matter Transport and Energy<br>Tempe, Arizona</div></div></div></div></div></div>
</div></div></div>