boundaryflux_advection_diffusion_1d Subroutine

  subroutine boundaryflux_advection_diffusion_1d(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_1d),intent(inout) :: this
    ! Local
    integer :: ndof

    ndof = this%solution%nelem*this%solution%nvar*2
    call boundaryflux_advection_diffusion_1d_gpu(this%solution%boundary_gpu, &
                                                 this%solution%extBoundary_gpu,this%solutionGradient%avgBoundary_gpu, &
                                                 this%flux%boundarynormal_gpu,this%u,this%nu,ndof)

  endsubroutine boundaryflux_advection_diffusion_1d