subroutine boundaryflux_advection_diffusion_3d(this)
! this method uses an linear upwind solver for the
! advective flux and the bassi-rebay method for the
! diffusive fluxes
implicit none
class(advection_diffusion_3d),intent(inout) :: this
call boundaryflux_advection_diffusion_3d_gpu(this%solution%boundary_gpu, &
this%solution%extBoundary_gpu,this%solutionGradient%avgBoundary_gpu, &
this%geometry%nhat%boundary_gpu,this%geometry%nscale%boundary_gpu, &
this%flux%boundarynormal_gpu,this%u,this%v,this%w, &
this%nu,this%solution%interp%N, &
this%solution%nelem,this%solution%nvar)
endsubroutine boundaryflux_advection_diffusion_3d