diff --git a/src/submodules/Line/src/LineInterpolationUtility@InterpolationMethods.F90 b/src/submodules/Line/src/LineInterpolationUtility@InterpolationMethods.F90
index 7e88f4c4..db12306a 100644
--- a/src/submodules/Line/src/LineInterpolationUtility@InterpolationMethods.F90
+++ b/src/submodules/Line/src/LineInterpolationUtility@InterpolationMethods.F90
@@ -391,7 +391,7 @@ SUBROUTINE handle_increasing
   !! layout INCREASING
   IF (layout(1:1) .EQ. "I") THEN
     DO ii = 1, nrow
-      CALL HeapSort(ans(ii, :))
+      CALL HeapSort(ans(ii, 1:ncol))
     END DO
   END IF
 END SUBROUTINE
diff --git a/src/submodules/Line/src/LineInterpolationUtility@QuadratureMethods.F90 b/src/submodules/Line/src/LineInterpolationUtility@QuadratureMethods.F90
index a2a42db6..21c8daaf 100644
--- a/src/submodules/Line/src/LineInterpolationUtility@QuadratureMethods.F90
+++ b/src/submodules/Line/src/LineInterpolationUtility@QuadratureMethods.F90
@@ -42,6 +42,11 @@
 CASE (qpopt%GaussLegendre, qpopt%GaussChebyshev, &
       qpopt%GaussJacobi, qpopt%GaussUltraspherical)
   ans = 1_I4B + INT(order / 2, kind=I4B)
+CASE (qpopt%GaussLegendreRadauRight, qpopt%GaussLegendreRadauLeft, &
+      qpopt%GaussChebyshevRadauLeft, qpopt%GaussChebyshevRadauRight, &
+      qpopt%GaussJacobiRadauLeft, qpopt%GaussJacobiRadauRight, &
+      qpopt%GaussUltraSphericalRadauLeft, qpopt%GaussUltraSphericalRadauRight)
+  ans = 2_I4B + INT((order - 1) / 2, kind=I4B)
 CASE DEFAULT
   ans = 2_I4B + INT(order / 2, kind=I4B)
 END SELECT