# Path: skyfield/framelib.py
# def build_matrix():
#     def rotation_at(t):
# def build_ecliptic_matrix(t):
#     def rotation_at(t):
#     def rotation_at(t):
#     def _dRdt_times_RT_at(t):
#     def rotation_at(t):
#     def _dRdt_times_RT_at(t):
#     def rotation_at(self, t):
#     def __init__(self, doc, matrix):
#     def rotation_at(self, t):
# class ICRS(object):
# class true_equator_and_equinox_of_date(object):
# class tirs(object):
# class itrs(object):
# class ecliptic_frame(object):
# class InertialFrame(object):
#         R = mxm(rot_z(-t.gast * tau / 24.0), t.M)
#             R = mxm(t.polar_motion_matrix(), R)
# 
# Path: skyfield/api.py
# N = E = +1.0
# S = W = -1.0
# 
# Path: skyfield/constants.py
# AU_M = 149597870700             # per IAU 2012 Resolution B2
# 
# ERAD = 6378136.6                # meters
# 
# Path: skyfield/positionlib.py
# class Geocentric(ICRF):
#     """An |xyz| position measured from the center of the Earth.
# 
#     A geocentric position is the difference between the position of the
#     Earth at a given instant and the position of a target body at the
#     same instant, without accounting for light-travel time or the effect
#     of relativity on the light itself.
# 
#     Its ``.position`` and ``.velocity`` vectors have |xyz| axes that
#     are those of the Geocentric Celestial Reference System (GCRS), an
#     inertial system that is an update to J2000 and that does not rotate
#     with the Earth itself.
# 
#     """
#     _default_center = 399
# 
#     def itrf_xyz(self):
#         """Deprecated; instead, call ``.frame_xyz(itrs)``. \
#         See `reference_frames`."""
#         return self.frame_xyz(framelib.itrs)
# 
#     def subpoint(self):
#         """Deprecated; instead, call either ``iers2010.subpoint(pos)`` or \
#         ``wgs84.subpoint(pos)``."""
#         from .toposlib import iers2010
#         return iers2010.subpoint(self)


from numpy import cos
from skyfield import framelib
from skyfield.api import Topos, load, wgs84
from skyfield.constants import AU_M, ERAD
from skyfield.positionlib import Geocentric

def test_radec_and_altaz_angles_and_rates():
    # HORIZONS test data in Skyfield repository: authorities/radec-altaz-rates
    ts = load.timescale()
    t = ts.utc(2021, 2, 3)

# Path: skyfield/framelib.py
# def build_matrix():
#     def rotation_at(t):
# def build_ecliptic_matrix(t):
#     def rotation_at(t):
#     def rotation_at(t):
#     def _dRdt_times_RT_at(t):
#     def rotation_at(t):
#     def _dRdt_times_RT_at(t):
#     def rotation_at(self, t):
#     def __init__(self, doc, matrix):
#     def rotation_at(self, t):
# class ICRS(object):
# class true_equator_and_equinox_of_date(object):
# class tirs(object):
# class itrs(object):
# class ecliptic_frame(object):
# class InertialFrame(object):
#         R = mxm(rot_z(-t.gast * tau / 24.0), t.M)
#             R = mxm(t.polar_motion_matrix(), R)
# 
# Path: skyfield/api.py
# N = E = +1.0
# S = W = -1.0
# 
# Path: skyfield/constants.py
# AU_M = 149597870700             # per IAU 2012 Resolution B2
# 
# ERAD = 6378136.6                # meters
# 
# Path: skyfield/positionlib.py
# class Geocentric(ICRF):
#     """An |xyz| position measured from the center of the Earth.
# 
#     A geocentric position is the difference between the position of the
#     Earth at a given instant and the position of a target body at the
#     same instant, without accounting for light-travel time or the effect
#     of relativity on the light itself.
# 
#     Its ``.position`` and ``.velocity`` vectors have |xyz| axes that
#     are those of the Geocentric Celestial Reference System (GCRS), an
#     inertial system that is an update to J2000 and that does not rotate
#     with the Earth itself.
# 
#     """
#     _default_center = 399
# 
#     def itrf_xyz(self):
#         """Deprecated; instead, call ``.frame_xyz(itrs)``. \
#         See `reference_frames`."""
#         return self.frame_xyz(framelib.itrs)
# 
#     def subpoint(self):
#         """Deprecated; instead, call either ``iers2010.subpoint(pos)`` or \
#         ``wgs84.subpoint(pos)``."""
#         from .toposlib import iers2010
#         return iers2010.subpoint(self)


# Path: skyfield/tests/test_frames.py
from numpy import cos
from skyfield import framelib
from skyfield.api import Topos, load, wgs84
from skyfield.constants import AU_M, ERAD
from skyfield.positionlib import Geocentric


def test_radec_and_altaz_angles_and_rates():
    # HORIZONS test data in Skyfield repository: authorities/radec-altaz-rates
    ts = load.timescale()
    t = ts.utc(2021, 2, 3)

    top = wgs84.latlon(35.1844866, 248.347300, elevation_m=2106.9128)