def lunar_altitude(self, tee):
"""Return the geocentric altitude of moon at moment, tee,
at location, location, as a small positive/negative angle in degrees,
ignoring parallax and refraction. Adapted from 'Astronomical
Algorithms' by Jean Meeus, Willmann_Bell, Inc., 1998."""
lamb = Lunar.lunar_longitude(tee)
beta = Lunar.lunar_latitude(tee)
alpha = Astro.right_ascension(tee, beta, lamb)
delta = Astro.declination(tee, beta, lamb)
theta0 = Astro.sidereal_from_moment(tee)
cap_H = mod(theta0 + self.longitude - alpha, 360)
altitude = arcsin_degrees(
(sin_degrees(self.latitude) * sin_degrees(delta)) +
(cos_degrees(self.latitude) * cos_degrees(delta) * cos_degrees(cap_H)))
return mod(altitude + 180, 360) - 180
def sidereal_lunar_longitude(tee):
"""Return sidereal lunar longitude at moment, tee."""
return mod(Lunar.lunar_longitude(tee) - Astro.precession(tee) + SIDEREAL_START, 360)
