docs: def at(self, t):
  		"""
  		Return the modelled tidal height at given times.
  		Arguments:
  		t -- array of times at which to evaluate the tidal height
  		"""
  		t0 = t[0]
  		hours = self._hours(t0, t)
  		partition = 240.0
  		t = self._partition(hours, partition)
  		times = self._times(t0, [(i + 0.5)*partition for i in range(len(t))])
  		speed, u, f, V0 = self.prepare(t0, times, radians = True)
  		H = self.model['amplitude'][:, np.newaxis]
  		p = d2r*self.model['phase'][:, np.newaxis]
  
  		return np.concatenate([
  			Tide._tidal_series(t_i, H, p, speed, u_i, f_i, V0)
  			for t_i, u_i, f_i in izip(t, u, f)
  		])