In the grand scheme of boosting performance by making the whole RHS of the Ribasim ODE problem C1 smooth, doing this for TabulatedRatingCurve is one of the steps. Currently the $Q(h)$ relation is given by piecewise linear interpolation of the input data. To make this more smooth, we can use a cubic polynomial that rounds out the corners. As proposed by @Huite, this is easiest to understand by hydrologist by taking some fixed small interval centered around each data point where the linear interpolation is replaced by a polynomial (some care must be taken half this interval is bigger than the distance between points in the $x$ direction).

For the nicest integration this could be incorporated into DataInterpolations.jl. For some inspiration see here, which uses an interval around the data points whose size is relative to the distance between points instead of absolute as proposed here.