AverageSides_Vector2D_t Subroutine

  subroutine AverageSides_Vector2D_t(this)
    implicit none
    class(Vector2D_t),intent(inout) :: this
    ! Local
    integer :: iel
    integer :: iside
    integer :: ivar
    integer :: i
    integer :: idir

    do concurrent(i=1:this%interp%N+1,iside=1:4,iel=1:this%nElem, &
                  ivar=1:this%nVar,idir=1:2)
      this%avgboundary(i,iside,iel,ivar,idir) = 0.5_prec*( &
                                                this%boundary(i,iside,iel,ivar,idir)+ &
                                                this%extBoundary(i,iside,iel,ivar,idir))
    enddo

  endsubroutine AverageSides_Vector2D_t