2023-12-02 17:01:06 +00:00
|
|
|
|
# sun_moon_test.py Test script for sun_moon.py
|
2023-12-01 16:41:40 +00:00
|
|
|
|
|
2023-12-02 17:01:06 +00:00
|
|
|
|
# Copyright (c) Peter Hinch 2023
|
|
|
|
|
|
# Released under the MIT license (see LICENSE)
|
2023-12-01 16:41:40 +00:00
|
|
|
|
|
2023-12-02 17:01:06 +00:00
|
|
|
|
# On mip-installed host:
|
|
|
|
|
|
# import sched.sun_moon_test
|
|
|
|
|
|
# On PC in astronomy directory:
|
|
|
|
|
|
# import sun_moon_test
|
|
|
|
|
|
|
|
|
|
|
|
try:
|
|
|
|
|
|
from .sun_moon import RiSet, abs_to_rel_days
|
|
|
|
|
|
except ImportError: # Running on PC in astronomy directory
|
|
|
|
|
|
from sun_moon import RiSet, abs_to_rel_days
|
|
|
|
|
|
|
|
|
|
|
|
nresults = [] # Times in seconds from local midnight
|
2023-12-01 16:41:40 +00:00
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
def show(rs):
|
|
|
|
|
|
print(f"Sun rise {rs.sunrise(3)} set {rs.sunset(3)}")
|
|
|
|
|
|
print(f"Moon rise {rs.moonrise(3)} set {rs.moonset(3)}")
|
2023-12-02 17:01:06 +00:00
|
|
|
|
nresults.extend([rs.sunrise(), rs.sunset(), rs.moonrise(), rs.moonset()])
|
2023-12-01 16:41:40 +00:00
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
print("4th Dec 2023: Seattle UTC-8")
|
|
|
|
|
|
rs = RiSet(lat=47.61, long=-122.35, lto=-8) # Seattle 47°36′35″N 122°19′59″W
|
|
|
|
|
|
rs.set_day(abs_to_rel_days(19695)) # 4th Dec 2023
|
|
|
|
|
|
show(rs)
|
|
|
|
|
|
print()
|
|
|
|
|
|
|
|
|
|
|
|
print("4th Dec 2023: Sydney UTC+11")
|
|
|
|
|
|
rs = RiSet(lat=-33.86, long=151.21, lto=11) # Sydney 33°52′04″S 151°12′36″E
|
|
|
|
|
|
rs.set_day(abs_to_rel_days(19695)) # 4th Dec 2023
|
|
|
|
|
|
show(rs)
|
|
|
|
|
|
print()
|
|
|
|
|
|
|
|
|
|
|
|
print("From 4th Dec 2023: UK, UTC")
|
|
|
|
|
|
rs = RiSet()
|
|
|
|
|
|
for day in range(7):
|
|
|
|
|
|
rs.set_day(abs_to_rel_days(19695 + day)) # Start 4th Dec 2023
|
|
|
|
|
|
print(f"Day: {day}")
|
|
|
|
|
|
show(rs)
|
|
|
|
|
|
|
2023-12-02 17:01:06 +00:00
|
|
|
|
# Expected results as computed on Unix build (64-bit FPU)
|
|
|
|
|
|
exp = [
|
|
|
|
|
|
27628,
|
|
|
|
|
|
58714,
|
|
|
|
|
|
85091,
|
|
|
|
|
|
46417,
|
|
|
|
|
|
20212,
|
|
|
|
|
|
71598,
|
|
|
|
|
|
2747,
|
|
|
|
|
|
41257,
|
|
|
|
|
|
29049,
|
|
|
|
|
|
57158,
|
|
|
|
|
|
82965,
|
|
|
|
|
|
46892,
|
|
|
|
|
|
29130,
|
|
|
|
|
|
57126,
|
|
|
|
|
|
None,
|
|
|
|
|
|
47460,
|
|
|
|
|
|
29209,
|
|
|
|
|
|
57097,
|
|
|
|
|
|
892,
|
|
|
|
|
|
47958,
|
|
|
|
|
|
29285,
|
|
|
|
|
|
57072,
|
|
|
|
|
|
5244,
|
|
|
|
|
|
48441,
|
|
|
|
|
|
29359,
|
|
|
|
|
|
57051,
|
|
|
|
|
|
9625,
|
|
|
|
|
|
48960,
|
|
|
|
|
|
29430,
|
|
|
|
|
|
57033,
|
|
|
|
|
|
14228,
|
|
|
|
|
|
49577,
|
|
|
|
|
|
29499,
|
|
|
|
|
|
57019,
|
|
|
|
|
|
19082,
|
|
|
|
|
|
50384,
|
|
|
|
|
|
]
|
|
|
|
|
|
print()
|
|
|
|
|
|
max_error = 0
|
|
|
|
|
|
for act, requirement in zip(nresults, exp):
|
|
|
|
|
|
if act is not None and requirement is not None:
|
|
|
|
|
|
err = abs(requirement - act)
|
|
|
|
|
|
max_error = max(max_error, err)
|
|
|
|
|
|
if err > 30:
|
|
|
|
|
|
print(f"Error {requirement - act}")
|
|
|
|
|
|
print(f"Maximum error {max_error}. Expect 0 on 64-bit platform, 30s on 32-bit")
|
|
|
|
|
|
|
2023-12-01 16:41:40 +00:00
|
|
|
|
# Times from timeanddate.com
|
|
|
|
|
|
# Seattle
|
|
|
|
|
|
# Sunrise 7:40 sunset 16:18 Moonrise 23:37 Moonset 12:53
|
|
|
|
|
|
# Sydney
|
|
|
|
|
|
# Sunrise 5:37 sunset 19:53 Moonrise 00:45 Moonset 11:28
|
|
|
|
|
|
# UK
|
|
|
|
|
|
# Sunrise 8:04 sunset 15:52 Moonrise 23:02 Moonset 13:01
|
