def local_hour_angle (self, h0, lat, delta) :
"""Local hour angle corresponding to the time of rise or set of a
celestial body.
"""
### J. Meeus, eq. (15.1), p. 102
cos_H0 = (h0.sin - lat.sin * delta.sin) / (lat.cos * delta.cos)
if abs (cos_H0) <= 1 :
return Angle_D.acos (cos_H0)
def local_hour_angle(self, h0, lat, delta):
"""Local hour angle corresponding to the time of rise or set of a
celestial body.
"""
### J. Meeus, eq. (15.1), p. 102
cos_H0 = (h0.sin - lat.sin * delta.sin) / (lat.cos * delta.cos)
if abs(cos_H0) <= 1:
return Angle_D.acos(cos_H0)
