[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