BoundaryInterp_Tensor2D_t Subroutine

  subroutine BoundaryInterp_Tensor2D_t(this)
    implicit none
    class(Tensor2D_t),intent(inout) :: this
! Local
    integer :: i,ii,idir,jdir,iel,ivar
    real(prec) :: fbs,fbe,fbn,fbw

    do concurrent(i=1:this%N+1, &
                  iel=1:this%nelem,ivar=1:this%nvar, &
                  idir=1:2,jdir=1:2)

      fbs = 0.0_prec
      fbe = 0.0_prec
      fbn = 0.0_prec
      fbw = 0.0_prec
      do ii = 1,this%N+1
        fbs = fbs+this%interp%bMatrix(ii,1)*this%interior(i,ii,iel,ivar,idir,jdir) ! South
        fbe = fbe+this%interp%bMatrix(ii,2)*this%interior(ii,i,iel,ivar,idir,jdir) ! East
        fbn = fbn+this%interp%bMatrix(ii,2)*this%interior(i,ii,iel,ivar,idir,jdir) ! North
        fbw = fbw+this%interp%bMatrix(ii,1)*this%interior(ii,i,iel,ivar,idir,jdir) ! West
      enddo

      this%boundary(i,1,iel,ivar,idir,jdir) = fbs
      this%boundary(i,2,iel,ivar,idir,jdir) = fbe
      this%boundary(i,3,iel,ivar,idir,jdir) = fbn
      this%boundary(i,4,iel,ivar,idir,jdir) = fbw

    enddo

  endsubroutine BoundaryInterp_Tensor2D_t