setboundarycondition_DGModel2D_t Subroutine

  subroutine setboundarycondition_DGModel2D_t(this)
    !! Boundary conditions for the solution are set to
    !! 0 for the external state to provide radiation type
    !! boundary conditions.
    implicit none
    class(DGModel2D_t),intent(inout) :: this
    ! local
    integer :: i,iEl,j,e2,bcid
    real(prec) :: nhat(1:2),x(1:2)

    do iEl = 1,this%solution%nElem ! Loop over all elements
      do j = 1,4 ! Loop over all sides

        bcid = this%mesh%sideInfo(5,j,iEl) ! Boundary Condition ID
        e2 = this%mesh%sideInfo(3,j,iEl) ! Neighboring Element ID

        if(e2 == 0) then
          if(bcid == SELF_BC_PRESCRIBED) then

            do i = 1,this%solution%interp%N+1 ! Loop over quadrature points
              x = this%geometry%x%boundary(i,j,iEl,1,1:2)

              this%solution%extBoundary(i,j,iEl,1:this%nvar) = &
                this%hbc2d_Prescribed(x,this%t)
            enddo

          elseif(bcid == SELF_BC_RADIATION) then

            do i = 1,this%solution%interp%N+1 ! Loop over quadrature points
              nhat = this%geometry%nhat%boundary(i,j,iEl,1,1:2)

              this%solution%extBoundary(i,j,iEl,1:this%nvar) = &
                this%hbc2d_Radiation(this%solution%boundary(i,j,iEl,1:this%nvar),nhat)
            enddo

          elseif(bcid == SELF_BC_NONORMALFLOW) then

            do i = 1,this%solution%interp%N+1 ! Loop over quadrature points
              nhat = this%geometry%nhat%boundary(i,j,iEl,1,1:2)

              this%solution%extBoundary(i,j,iEl,1:this%nvar) = &
                this%hbc2d_NoNormalFlow(this%solution%boundary(i,j,iEl,1:this%nvar),nhat)
            enddo

          endif
        endif

      enddo
    enddo

  endsubroutine setboundarycondition_DGModel2D_t