[Nek5000-users] userbc to set BC from velocity field computed in userchk

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

Thank you for your quick answer!

Best wishes,

Arnold

Le 06/10/2016 18:43, nek5000-users at lists.mcs.anl.gov 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: nek5000-users-bounces at lists.mcs.anl.gov [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: Thursday, October 06, 2016 10:37 AM
> To: nek5000-users at lists.mcs.anl.gov
> 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
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users