AverageSides_Scalar2D_t Subroutine

  subroutine AverageSides_Scalar2D_t(this)
    implicit none
    class(Scalar2D_t),intent(inout) :: this
    ! Local
    integer :: iel
    integer :: iside
    integer :: ivar
    integer :: i

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

  endsubroutine AverageSides_Scalar2D_t