[Nek5000-users] Problems with post-processing multiple-boxes in MATLAB and Python

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Thu Sep 14 10:33:07 CDT 2017


Hi All,


Genbox generates a tensor-product array of elements with layout:


   (ix , iy, iz)    (ix fastest moving, iy 2nd fastest, iz 3rd fastest)


If you want to have multiple boxes that continue in this ordering, you could

do the following, say:


box (1)

-5 -5 -5

-1 1 1.0

-1 1 1.0

0 2   1.0


box(2)

-5 -5 -10

-1 1 1

-1 1 1

1.0 9.0


etc


This would stack a 5x5x10 box on top of a 5x5x5 box


and topologically is the same as a 5x5 x (15) box


... i.e., things get concatenated in the z direction.



Also, note that you can use custom distributions in each direction, e.g.,


box

3  3  -5

0 0.1 0.3  10

0  .9   22 33

-5  9  1.5


which would give spacing [ 0 .1 .3 10] in the x direction


[0 .9 22 33] in y


and geometric partition of 5 elements in z on [-5:9]

with size ratio 1.5 between each adjacent z level.


hth


Paul


________________________________
From: Nek5000-users <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: Thursday, September 14, 2017 10:20:27 AM
To: shivang agarwal
Subject: Re: [Nek5000-users] Problems with post-processing multiple-boxes in MATLAB and Python

Hi Daniel,

I've had a similar problem before. Can you send your .box file so we can have a look?

Best regards,

Goran

On Sep 14, 2017, at 11:13, "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>> wrote:


Hi All,


I am using Nek5000 to simulate a buoyant plume rising from the centre of the floor of a square 2D domain.  I am trying to use MATLAB and/or Python to post-process my data.  Everything works fine when I post-process a simulation that was initialised with only one box in my .box file.  When, however, I run a simulation with multiple boxes (i.e. to cluster nodes), the local to global mapping of elements appears to go wrong and puts data into the wrong places.  When plotting a contour plot of coordinates for example, where I would expect a constantly changing colour gradient, I instead get back colour gradients in four quarters.


I haven't used this list before so if I need to provide examples of images or code I can do.  If anybody else has had this problem and has found a way to overcome it by finding the correct mapping some help would be very much appreciated.


Best regards,


Daniel Ward

_______________________________________________
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/20170914/4ab9769c/attachment.html>


More information about the Nek5000-users mailing list