code: def fix_repeat_dt(dt_list, offset_s=0.001):
      
      idx = (np.diff(dt_list) == timedelta(0))
      while np.any(idx):
          dt_list[idx.nonzero()[0] + 1] += timedelta(seconds=offset_s)
          idx = (np.diff(dt_list) == timedelta(0))
      return dt_list
docs: Add some small offset to remove duplicate times
      Needed for xarray interp, which expects monotonically increasing times