[Nek5000-users] Data Structure

[nek5000-users at lists.mcs.anl.gov]

Hi Shriram,

{XM1,YM1,ZM1} stores the GLL mesh for the velocity and {XM2,YM2,ZM2} stores GL or GLL points for the pressure (defined in the GEOM CB). The mesh points are stored in a lexicographical order.

{XC,YC,ZC} stores the element vertices (defined in the INPUT CB). Note: we read the vertices into NEK using a preprocessor notation. However internally we use a different symmetric notation. For more details check the header of genxyzl().

In general you have access to all the arrays you like as long as you have a proper include statement for the CB in question.

Is there any specific reason why you want to have access to a particular point in the mesh?


Stefan

 
On Apr 19, 2010, at 6:24 AM, nek5000-users at lists.mcs.anl.gov wrote:

> Hello,
> 
> I am trying to understand the structure of the arrays that store velocities and co-ordinates in nek that is available under several subroutines like uskchk,usrdat2,usrdat etc. 
> Based on my understanding, I found that the velocities, say vx has the form vx(lx1,ly1,lz1,e) but I found that there are many arrays that address the co-ordinates namely xm0,xm1,xc etc though I am not sure what it denotes. Can all these be accessed from any subroutine in the usr file ? I have attached herewith a picture where I have mentioned how I thought the data structure might be. But, I would really like to know how one would access the co-ordinates of say ,a particular element and what xm0,xm1,xc denote for. Please correct me if I am wrong about the data structure. Thanks.
> 
> Regards
> Shriram Jagannathan
> <Array_Struct.pdf>_______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users