R1i(dZddlmZmZddlmZmZddlmZm Z ddl m Z ddl m Z ddlmZdd lmZe d Zd Zgd Zgd ZgdZdZgdZdZdZddgZdZddgZdZddgZdZ dZ!dddd d!d"Z"d#Z# d&d$Z$y%)'zIRoutines to solve for circumstances like sunrise, sunset, and moon phase.)print_functiondivision)cos zeros_like)pitau)ecliptic_frame) find_discrete)iau2000b_radians)Anglecv|dj|j||}|j|dS)z~ .. deprecated:: 1.42 Use the :meth:`~skyfield.positionlib.ICRF.phase_angle()` position method instead. earthsun)atobserve phase_angle) ephemerisbodytps T/home/cursorai/projects/iching/venv/lib/python3.12/site-packages/skyfield/almanac.pyrrs< 'a ((49A ==5) **cRt|||j}ddt|zzS)z .. deprecated:: 1.42 Use the :meth:`~skyfield.positionlib.ICRF.fraction_illuminated()` position method instead. g?g?)rradiansr)rrras rfraction_illuminatedrs+ ItQ'//A #A, r)SpringSummerAutumnWinter)zVernal EquinoxzSummer SolsticezAutumnal EquinoxzWinter Solstice)z March Equinoxz June SolsticezSeptember EquinoxzDecember Solsticec8|d|dfd}d|_|S)aBuild a function of time that returns the quarter of the year. The function that this returns will expect a single argument that is a :class:`~skyfield.timelib.Time` and will return 0 through 3 for the seasons Spring, Summer, Autumn, and Winter. rrc t||_j|}|jj j t \}}}|jtdz zdzjtS)z8Return season 0 (Spring) through 3 (Winter) at time `t`. r _nutation_angles_radiansrrapparent frame_latlonr rr astypeint)re_slonrrs r season_atzseasons..season_atGsg%5a%8" HHQKYYs^,,.;;NK 4 q)A-55c::rgV@ step_days)rr.rrs @@rseasonsr1<s. g E E C;I r)zNew Moonz First Quarterz Full Moonz Last Quartercj|dj|}|d|d}}|j|jjt\}}}|j|jjt\}}}t |j |j z tzS)uhReturn the Moon phase 0°–360° at time ``t``, where 180° is Full Moon. More precisely: this returns an :class:`~skyfield.units.Angle` giving the difference between the geocentric apparent ecliptic longitudes of the Moon and Sun, constrained to the interval 0°–360° (0–𝜏 radians) where 0° is New Moon and 180° is Full Moon. rmoonr)r)rrr'r(r r rr )rrr+r3rr,mlonr-s r moon_phaser5Xs 'a A&!9U#3#D4))+88HJAtQ3((*77GJAtQ $,,5< ==rcF|d|d|dfd}d|_|S)aPBuild a function of time that returns the moon phase 0 through 3. The function that this returns will expect a single argument that is a :class:`~skyfield.timelib.Time` and will return the phase of the moon as an integer. See the accompanying array ``MOON_PHASES`` if you want to give string names to each phase. rr3rct||_j|}|jj j t \}}}|jj j t \}}}|j|jz tdz zdzjtS)5Return the phase of the moon 0 through 3 at time `t`.r$r%)rr+r,r4r-rr3rs r moon_phase_atz"moon_phases..moon_phase_atts%5a%8" HHQKYYt_--/<<^L 4YYs^,,.;;NK 4 ,#':Q>FFsKKrg@r/)rr9rr3rs @@@r moon_phasesr:gs9 g E V D E CL"M r descending ascendingc8|d|dfd}d|_|S)aBuild a function of time that identifies lunar nodes. This returns a function taking a :class:`~skyfield.timelib.Time` and returning ``True`` if the Moon is above the ecliptic else ``False``. See :ref:`lunar-nodes` for how to use this routine. rr3cj|}|jjjt\}}}|j dkDS)r8g)rrr'r(r r)rr+latr,rr3s r moon_node_atz moon_nodes..moon_node_atsF HHQKIIdO,,.;;NK Q{{S  rg(@r/)rr@rr3s @@r moon_nodesrAs. g E V D! "L r conjunction oppositionc|dj|dfd}jdk(r d|_|Sd|_|S)zBuild a function to find oppositions and conjunctions with the Sun. See :ref:`oppositions-conjunctions` for how to call this routine and interpret the results. rrcP|}|jjjt\}}}|jjjt\}}}|j|jz t z dzj dS)z5Return whether the target is east or west of the Sun.g@int8)rr'r(r rrr))rr+r,r-tlonearth_atrtargets rleading_or_trailingz5oppositions_conjunctions..leading_or_trailings QKYYs^,,.;;NK 4YYv&//1>>~N 4 ,2S8@@HHri-()rrIr0)rrIrJrHrs ` @@roppositions_conjunctionsrMsU!$$H E CI}}(*% )+% rzAntimeridian transitzMeridian transitcL|dzjfd}d|_|S)uBuild a function of time for finding when a body transits the meridian. The returned function accepts a :class:`~skyfield.timelib.Time` argument and returns ``True`` if the ``target`` body is west of the observer’s meridian at that time, and otherwise returns ``False.`` See :ref:`transits` for how to use this to search for a body’s meridian transits and antimeridian transits. rc.t||_j|jd\}}}|j j jd\}}}|j |j z tztkS)uAReturn `True` if the target is west of the observer’s meridian.date)epoch) r r&rradecrr'rr r)rra1r,ra2rItopostopos_ats rwest_of_meridian_atz.meridian_transits..west_of_meridian_ats%5a%8" HHQK%%F%3 QQK''/88:@@v@N Q ckk)S0255rg?rr0)rrIrUrWrVs `` @rmeridian_transitsrYs0'"U*..H 6%(! rcR|d|d|zjfd}d|_|S)ulBuild a function of time that returns whether the Sun is up. The function that is returned will expect a single argument that is a :class:`~skyfield.timelib.Time`, and will return ``True`` if the sun is up, else ``False``. Skyfield uses the same definition as the United States Naval Observatory: the Sun is up when its center is 0.8333 degrees below the horizon, which accounts for both its apparent radius of around 16 arcminutes and also for the 34 arcminutes by which atmospheric refraction on average lifts the image of the Sun. If you need to provide a custom value for refraction, adjust the estimate of the Sun’s radius, or account for a vantage point above the Earth’s surface, see :ref:`risings-and-settings` to learn about the more versatile :func:`~skyfield.almanac.risings_and_settings()` routine. rrct||_|jjj dj dk\S)zReturn `True` if the sun has risen by time `t`. The Sun has risen if its altitude above the horizon is greater than -0.8333 degrees. r7d꿩r r&rr'altazdegrees)rrrVs r is_sun_up_atz$sunrise_sunset..is_sun_up_atsI&6a%8"{""3'00288:1=EEPPr{Gz?rX)rrUr`rrVs @@rsunrise_sunsetrbs:( E C'"U*..HQ"L rNightzAstronomical twilightzNautical twilightzCivil twilightDay)rrr$cR|d|d|zjfd}d|_|S)u`Build a function of time returning whether it is dark, twilight, or day. The function that this returns will expect a single argument that is a :class:`~skyfield.timelib.Time` and will return: | 0 — Dark of night. | 1 — Astronomical twilight. | 2 — Nautical twilight. | 3 — Civil twilight. | 4 — Sun is up. rrct||_|jjj dj }t |t}d||dk\<d||dk\<d||dk\<d||d k\<|S) zAReturn whether the Sun is up, down, or whether there is twilight.rrg2reg(rfgr$r\)r r&rr'r^r_rr*)rr_rrrVs ris_it_dark_twilight_day_atz5dark_twilight_day..is_it_dark_twilight_day_at s%5a%8"1+%%c*335;;=a@HH w $'U 'U 'T/ !'W rrarX)rrUrjrrVs @@rdark_twilight_dayrks: E C'"U*..H ,0( %%rcV|d|zj||z fd}d|_|S)uBuild a function of time that returns whether a body is up. This returns a function taking a :class:`~skyfield.timelib.Time` argument returning ``True`` if the body’s altazimuth altitude angle plus ``radius_degrees`` is greater than ``horizon_degrees``, else ``False``. See :ref:`risings-and-settings` to learn about how to search for risings and settings, and to see more about using the parameters ``horizon_degrees`` and ``radius_degrees``. rct||_|jjj dj kDS)z2Return `True` if the target has risen by time `t`.rr])rhrIrVs r is_body_up_atz+risings_and_settings..is_body_up_at)sG%5a%8"{""6*335;;=a@HH1LLrg?rX)rrIrUhorizon_degreesradius_degreesrornrVs ` @@rrisings_and_settingsrrs:'"U*..H.(AM #M rN)g""""""r)%__doc__ __future__rrnumpyrr constantsrr framelibr searchlibr nutationlibr unitsr rrSEASONS SEASON_EVENTSSEASON_EVENTS_NEUTRALr1 MOON_PHASESr5r: MOON_NODESrA CONJUNCTIONSrMMERIDIAN_TRANSITSrYrb TWILIGHTSrkrrrrrsO/!$$)+   *  >2 * .,-?@4"J   &>EFr