fill_series.py

def fill_series(gridded, 
                max_contiguous_na=3, 
                periods_for_extrapolation=5,
                method='slinear'):
    # Applies the fill_values function over every grid box
    filled = np.apply_along_axis(func1d=lambda s: fill_values(s, 
                                                              max_contiguous_na, 
                                                              periods_for_extrapolation, 
                                                              method),
                                 arr=gridded.reshape(gridded.shape[0], -1).T,
                                 axis=1)
    filled = filled.T.reshape(gridded.shape)
    return filled