¿Cómo calculo el acimut del sol según el cenit, el ángulo horario, la declinación y la latitud?

Wikipedia cita esta fórmula para calcular el azimut solar:

Dónde:

• $\phi_s$es el ángulo azimutal solar
• $\theta_s$es el ángulo cenital solar
• $h$es el ángulo horario, en la hora solar local
• $\delta$es la declinación actual del sol
• $\Phi$es la latitud local

Estoy implementando esta ecuación en JavaScript, pero no parece correcta. ¿Puede alguien aquí ayudarme a confirmar si esta ecuación es correcta o no?

Aquí está mi implementación de JS. Lo único extraño está al final: limitar el dominio acos a [-1, 1] para evitar valores de NaN.

const getSolarAzimuth = function(
    zenithAngle, hourAngle, declination, latitude
) {
    const cos = Math.cos;
    const sin = Math.sin;

    const cos_phi = (
        sin(declination) * cos(latitude) - cos(hourAngle) * cos(declination) *
            sin(latitude)) /
          sin(zenithAngle);

    return Math.acos(Math.min(Math.max(cos_phi, -1), 1));
};