Compute the Greenwich Mean Sidereal Time (GMST) for a UT1 instant.
- Parameters
-
| jd_ut1 | The julian day number, on the UT1 scale — NOT the JDE (TT) used by most other algorithms in this library. Feeding a TT value here silently shifts the result by ΔT ≈ 69 s of time (≈ 0.29° of rotation, at epoch ~2026). |
- Returns
- The GMST (Greenwich hour angle of the mean vernal equinox), normalized to [0°, 360°).
- Note
- Meeus (12.3) gives GMST at 0h UT as a polynomial in julian centuries of universal time; its daily advance is only the ~0.9856° 0h-to-0h drift — the 360°/day rotation drops out mod 360° only on the 0h UT grid (Meeus's Example 12.a is such a 0h case). Formula (12.4), used here, folds the full sidereal rate back in and stays valid at any UT1 instant: θ₀ = 280.46061837° + 360.98564736629°·d + 0.000387933°·T² − T³/38710000 (mod 360°), where the constant picks up 180° because J2000.0 falls at noon UT. On the 0h grid this coincides with (12.3) exactly; julian_day::jde_to_jc is not reused since its contract is TT-based.
- See also
- Jean Meeus, "Astronomical Algorithms", Second Edition, Chapter 12, Formulas (12.2)-(12.4).