docs: def _times(t0, hours):
  		"""
  		Return a (list of) datetime(s) given an initial time and an (list of) hourly offset(s).
  		Arguments:
  		t0 -- initial time
  		hours -- hourly offsets from t0
  		"""
  		if not isinstance(hours, Iterable):
  			return Tide._times(t0, [hours])[0]
  		elif not isinstance(hours[0], datetime):
  			return np.array([t0 + timedelta(hours=h) for h in hours])
  		else:
  			return np.array(hours)