def METHOD_NAME():
    da = make_array().rename_dims({'xx': 'time'})
    x = sc.arange(
        dim='time', start=0, stop=da.sizes['time'], step=1, unit='s', dtype='datetime64'
    )
    da.coords['time'] = x
    out = interp1d(da, 'time')(da.coords['time'], midpoints=True)
    int_x = (x - sc.epoch(unit=x.unit)).values
    midpoints = int_x[:-1] + 0.5 * (int_x[1:] - int_x[:-1])
    assert np.array_equal(
        out.values,
        theirs.interp1d(
            x=da.coords['time'].values.astype('int64'), y=da.values, axis=0
        )(midpoints),
    )