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