A Community Devoted to the Preservation and Practice of Celestial Navigation and Other Methods of Traditional Wayfinding
From: Tony Oz
Date: 2021 Jul 20, 09:42 -0700
Usually I use the following formula for the azimuth:
cos(Az) = (sin(Dec) - sin(H)·sin(Lat)) / (cos(H)·cos(Lat))
It works good except on poles (Lat == +90° or -90°) where it fails with the division by zero.
Is there a formula that would work everywhere - poles included - and cleanly return the GHA of a body in question when the calculation is done for an observer on a pole?
Thank you in advance.