Given a date, time, geographic location, and equatorial coordinates of a celestial body,
this function returns horizontal coordinates (azimuth and altitude) relative to the horizon.
The ra and dec must be equator-of-date coordinates. Equator-of-date coordinates can be
obtained by calling astro_equator() with equdate = "EQUATOR_OF_DATE" and
aberration = "ABERRATION".
Atmospheric refraction correction is recommended. Pass refraction = "REFRACTION_NORMAL"
to correct for optical lensing of the Earth's atmosphere that causes objects to appear
higher above the horizon than they actually are.