kuunvaiheet.py

#Kuun vaiheet aikana t0 -> t1
# tulostaa UT-ajassa eli Greenwichin ajassa
#Kuun vaiheet aikana t0 -> t1
# tulostaa UT-ajassa eli Greenwichin ajassa
# Juhani Kaukoranta 5.1.2023 (korjattu 6.1.2023)
# latausivusto JPL:n paketille de421.bsp:
# paketti ladataan vain kerran, se jää hakemistoon
print(" sivusto, josta voi ladata datapaketin de421.bsp")
print("https://github.com/…/python…/blob/master/ci/de421.bsp")
from skyfield import api
from skyfield import almanac
ts = api.load.timescale()
eph = api.load('de421.bsp') # on ladattuna NASA:lta
print("Mistä alkaen, vuosi kk päivä, esim 2023 1 1")
line = [int(x) for x in input().split() ]
year1,month1,day1 = line[0],line[1],line[2]
print("Mihin asti, vuosikk päivä, esim 2023 6 1")
line = [int(x) for x in input().split()]
year2,month2,day2 = line[0],line[1],line[2]
#t = ts.utc(2023, 2, 28) # malli
t0 = ts.utc(year1,month1,day1)  #  alkupäiväys
t1 = ts.utc(year2,month2,day2) # loppupäiväys
t, y = almanac.find_discrete(t0, t1, almanac.moon_phases(eph))
print("kellonajat UTC (Greenwich) ajassa")
for i in range(0,len(y)): print(almanac.MOON_PHASES[y[i]] , t.utc_iso()[i]," UT")