add limit parameter to interpolation functions

The functions interpolate and linear currently don't expose the parameter limit which makes the unsuitable to do interpolations in a classical sense, i.e. fill all missing values. This MR introduces the parameter and sets them to None.

After we decided on !645 (closed), I would also like to expose extrapolate.