boundaryflux_advection_diffusion_2d Subroutine

  subroutine boundaryflux_advection_diffusion_2d(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_2d),intent(inout) :: this

    call boundaryflux_advection_diffusion_2d_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%nu,this%solution%interp%N, &
                                                 this%solution%nelem,this%solution%nvar)

  endsubroutine boundaryflux_advection_diffusion_2d