[Nek5000-users] How to get another element information that in the different processor?

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Mon Jul 24 21:06:50 CDT 2017


Hi all,
I want to get another element information when I handle one element, and the two elements maybe not in the same processor. How to get another element information that in the different processor?

Below is some of my code:
1)for element e, get global element number;
2)get the correspond element number at the inlet;
3)call the ym1 value at the inlet element
But I don't know how to get ym1(i,j,k,e_ele) in the e_nid processor.

 do e=1,nelv
	eg=lglel(e)
	call get_exyz(ex,ey,ez,eg,nelx,nely,nelz)
	e_in=(ez-1)*nelx*nely+(ey-1)*nelx+1
	e_nid=gllnid(e_in)
	e_ele=gllel(e_in)
c	if (nid.eq.e_nid) then	   
	   do k=1,nz1
	    do j=1,ny1
	     do i=1,nx1
		temp=-ym1(i,j,k,e)*u_t/vis/A_plus
		cdyn(i,j,k,e)=cs0*(1-exp(temp))
	     enddo
	    enddo
	   enddo
c	endif	
      enddo

Best wishes!
Hu


More information about the Nek5000-users mailing list