Sunrise Sunset

Rising and setting times of the Sun and the Moon and twilight times

Year, month, day, lambda (Geographic east longitude of the observer in [deg]), phi (Geographic latitude of the observer in [deg]), zone (local time - universal time in [hour]) and twilight (Indicates civil, nautical or astronomical twilight) are received, then rising and setting times of the Moon and the Sun and twilight times are computed.

In this version, the positions of the Sun and the Moon are computed using Jet Propulsion Laboratory Development Ephemeris (DE440). Moreover, Earth orientation parameters are not used. So, MJD_ET is utilized as an approximation of MJD_TT and MJD_TDB in SinAlt.m, and MJD_UTC is applied instead of MJD_UT1.

east longitude and latitude of the observer can be found at:

https://www.gaisma.com/en/dir/001-continent.html

https://eclipse.gsfc.nasa.gov/LEcat5/deltatpoly.html

References:

Montenbruck O., Pfleger T.; Astronomy on the Personal Computer; Springer Verlag, Heidelberg; 4th edition (2000).

https://ssd.jpl.nasa.gov/planets/eph_export.html