We should add a domain to the Spline protocol in `spline.py`. After having add a domain for splines, we have to adjust the upper and lower bound in `add_cubic.py` (see TODOs in corresponding module). [Reference to comment in PR](https://github.com/climate-resource/gradient-aware-harmonisation/pull/20#discussion_r2046516406)