Thank you for your quick answer! Best wishes, Arnold Le 06/10/2016 18:43, [email protected] a �crit :
Hi Arnold,
Correct syntax would be:
subroutine userbc(i,j,k,f,eg) include 'SIZE' include 'TOTAL' YOUR COMMON BLOCK HERE (to pass boundx...)
integer e,eg,f
e = gllel(eg) ! To get the local element number
ux = boundx(i,j,k,e) uy = boundy(i,j,k,e) uz = boundz(i,j,k,e)
return end
Paul
________________________________________ From: [email protected] [[email protected]] on behalf of [email protected] [[email protected]] Sent: Thursday, October 06, 2016 10:37 AM To: [email protected] Subject: [Nek5000-users] userbc to set BC from velocity field computed in userchk
Good afternoon Neks,
I am trying to set boundary conditions from a velocity field computed in the userchk subroutine. The nature of the velocity field is defined as follows :
real boundx(lx1,ly1,lz1,lelt) real boundy(lx1,ly1,lz1,lelt) real boundz(lx1,ly1,lz1,lelt)
I was wondering what was the proper syntax to set the user defined Dirichlet velocity; is the following syntax in the userbc correct ?
ux=boundx(ix,iy,iz,iside,ie) uy=boundy(ix,iy,iz,iside,ie) uz=boundz(ix,iy,iz,iside,ie)
Thank you for your time,
Best wishes,
Arnold _______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users _______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users