AverageSides_Scalar1D_t Subroutine

  subroutine AverageSides_Scalar1D_t(this)
    implicit none
    class(Scalar1D_t),intent(inout) :: this
    ! Local
    integer :: iel
    integer :: ivar

    do concurrent(iel=1:this%nElem,ivar=1:this%nVar)
      ! Left side - we account for the -\hat{x} normal
      this%avgboundary(1,iel,ivar) = 0.5_prec*( &
                                     this%boundary(1,iel,ivar)+ &
                                     this%extBoundary(1,iel,ivar))

      ! Right side - we account for the +\hat{x} normal
      this%avgboundary(2,iel,ivar) = 0.5_prec*( &
                                     this%boundary(2,iel,ivar)+ &
                                     this%extBoundary(2,iel,ivar))
    enddo

  endsubroutine AverageSides_Scalar1D_t