Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- # Lunar eclipses
- # Juhani Kaukoranta 18.1.2022
- from skyfield import api
- ts = api.load.timescale()
- eph = api.load('de406.bsp') # kattaa vuodet -3000 -> 3000
- from skyfield import eclipselib
- t0 = ts.utc(2022, 1, 1) # mistä alkaen lasketaan
- t1 = ts.utc(2030, 1, 1) # mihin asti lasketaan
- t, y, details = eclipselib.lunar_eclipses(t0, t1, eph)
- # tulostetaan kuunpimennykset
- # Total = täydellinen pimennys (Kuu kokonaan Maan varjossa)
- # Partial = osittainen pimennys (Osa Kuusta Maan varjossa)
- # Penumbral = puolivarjopimennys (Kuun Maan puolivarjossa)
- for ti, yi in zip(t, y):
- print(ti.utc_strftime('%Y-%m-%d %H:%M'),
- 'y={}'.format(yi),
- eclipselib.LUNAR_ECLIPSES[yi])
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement