[Nek5000-users] Nodal points using face and element number

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

Hi Sijo,


You can use the function "facind" to obtain a set of index ranges on one face, denoted by the input argument "iface". Then, to obtain the x-coordinates you just access entries in array xm1.



      call facind (kx1,kx2,ky1,ky2,kz1,kz2,nx,ny,nz,iface)
      do 100 iz=kz1,kz2
      do 100 iy=ky1,ky2
      do 100 ix=kx1,kx2
         write(6,*) xm1(ix,iy,iz,ie)
  100 continue
      

Hope this helps,
Li



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: Wednesday, March 14, 2018 10:24

To: nek5000-users

Subject: [Nek5000-users] Nodal points using face and element number











Hi Neks,







Could you please tell me, if i give the face number and element number, can i get the nodal co ordinates? My problem is I have a body in fluid domain and I have the element number and face number of elements which is in contact with the surface of the body. I need to get the x cordinates of the Nodes which are on the body.











Regards







Sijo