Chage default behavior of interpolate

`rdtools.interpolate()` default behavior is to use a `max_timedelta` based on the median timedelta of the input. This can lead to a high fraction of NaNs in the result for some irregular time series. Changing changing the default behavior to 2x the median timedelta would prevent this. It may also help to add a warning when a substantial fraction of returned points are NaN. 