732642856
42 lines
1.7 KiB
Python
42 lines
1.7 KiB
Python
"""Swiss Ephemeris sunrise/sunset evidence for Saham day/night formula selection."""
|
|
|
|
from __future__ import annotations
|
|
|
|
from datetime import datetime
|
|
from typing import Any
|
|
|
|
import swisseph as swe
|
|
|
|
|
|
class SahamDayNightError(ValueError):
|
|
pass
|
|
|
|
|
|
def determine_daytime(moment: datetime, *, lat: float, lon: float, tz: float) -> dict[str, Any]:
|
|
if not -90 <= float(lat) <= 90 or not -180 <= float(lon) <= 180:
|
|
raise SahamDayNightError("invalid WGS84 latitude/longitude")
|
|
local = moment.replace(tzinfo=None)
|
|
jd = swe.julday(local.year, local.month, local.day, local.hour + local.minute / 60 + local.second / 3600 - float(tz))
|
|
geopos = (float(lon), float(lat), 0.0)
|
|
rise_status, rise = swe.rise_trans(jd - 1.0, swe.SUN, swe.CALC_RISE, geopos)
|
|
set_status, sunset = swe.rise_trans(jd - 1.0, swe.SUN, swe.CALC_SET, geopos)
|
|
if rise_status != 0 or set_status != 0:
|
|
raise SahamDayNightError("sunrise_or_sunset_unavailable_for_location_date")
|
|
sunrise_jd, sunset_jd = rise[0], sunset[0]
|
|
# Normalize the next daily events around the queried instant.
|
|
while sunrise_jd > jd:
|
|
sunrise_jd -= 1.0
|
|
while sunset_jd > jd:
|
|
sunset_jd -= 1.0
|
|
is_day = sunrise_jd <= jd < sunset_jd if sunrise_jd < sunset_jd else not (sunset_jd <= jd < sunrise_jd)
|
|
return {
|
|
"scope": "saham_daynight_swiss",
|
|
"status": "computed",
|
|
"is_daytime": is_day,
|
|
"julian_day_ut": jd,
|
|
"sunrise_jd_ut": sunrise_jd,
|
|
"sunset_jd_ut": sunset_jd,
|
|
"method": "swisseph.rise_trans",
|
|
"boundary": "Formula-specific +30 degree exceptions must be applied by the Saham rule layer, not inferred from house placement.",
|
|
}
|