docs: def get_unit(a):
      """Extract the time unit from array's dtype"""
      typestr = a.dtype.str
      i = typestr.find('[')
      if i == -1:
          raise TypeError("Expected a datetime64 array, not %s", a.dtype)
      return typestr[i + 1: -1]