docs: def sideral(date, longitude=0., model='mean', eop_correction=True, terms=106):  # pragma: no cover
      """Sideral time as a rotation matrix
      """
      theta = _sideral(date, longitude, model, eop_correction, terms)
      return rot3(np.deg2rad(-theta))