MoonPhases.py

#Kuun vaiheet
# Juhani Kaukoranta 27.4.2021
from skyfield import api
from skyfield import almanac
ts = api.load.timescale()
eph = api.load('de421.bsp') # JPL efemerididataa
t0 = ts.utc(2021, 4, 20) # mistä alkaen
t1 = ts.utc(2021, 6, 30) # mihin asti
t, y = almanac.find_discrete(t0, t1, almanac.moon_phases(eph))
# t-vektori sisältää päiväykset
# y-vektori  sisältää Kuun vaiheiden nimet

for i in range(0,len(y)):
    print(almanac.MOON_PHASES[y[i]] , t.utc_iso()[i])