ó šßÈ[c@`sFdZddlmZmZmZmZddlmZyddlmZWn!e k roddl mZnXddl Z ddl mZddl mZdd lmZmZdd lmZmZdd l mZdd lmZd ddddddddddddddddddgZd Zd!Zed"ƒZed#ƒZd$Z ed%ƒZ!dZ"d&Z#d'Z$d(Z%d)Z&d*Z'd+Z(d,„Z)d-efd.„ƒYZ*d/ej+fd0„ƒYZ,e,ƒZ-d1e.fd2„ƒYZ/d3efd4„ƒYZ0d5e0fd6„ƒYZ1d7e0fd8„ƒYZ2d9e1fd:„ƒYZ3e2j4e0_4dS(;u¬ The astropy.utils.iers package provides access to the tables provided by the International Earth Rotation and Reference Systems Service, in particular allowing interpolation of published UT1-UTC values for given times. These are used in `astropy.time` to provide UT1 values. The polar motions are also used for determining earth orientation for celestial-to-terrestrial coordinate transformations (in `astropy.coordinates`). i(tabsolute_importtdivisiontprint_functiontunicode_literals(twarn(turlparseNi(tconfig(tunits(tTabletQTable(tget_pkg_data_filenametclear_download_cache(tutils(tAstropyWarninguConfuconfuIERSuIERS_BuIERS_Au IERS_Autou FROM_IERS_Bu FROM_IERS_AuFROM_IERS_A_PREDICTIONuTIME_BEFORE_IERS_RANGEuTIME_BEYOND_IERS_RANGEu IERS_A_FILEu IERS_A_URLu IERS_A_READMEu IERS_B_FILEu IERS_B_URLu IERS_B_READMEuIERSRangeErroruIERSStaleWarningufinals2000A.allu.http://maia.usno.navy.mil/ser7/finals2000A.alludata/ReadMe.finals2000Audata/eopc04_IAU2000.62-nowu<http://hpiers.obspm.fr/iers/eop/eopc04/eopc04_IAU2000.62-nowudata/ReadMe.eopc04_IAU2000iiiÿÿÿÿiþÿÿÿg@€OBAu¦interpolating from IERS_Auto using predictive values that are more than {0} days old. Normally you should not see this error because this class automatically downloads the latest IERS-A table. Perhaps you are offline? If you understand what you are doing then this error can be suppressed by setting the auto_max_age configuration variable to ``None``: from astropy.utils.iers import conf conf.auto_max_age = None cO`s9tjjjdtjƒtjj||ŽSWdQXdS(u Overload astropy.utils.data.download_file within iers module to use a custom (longer) wait time. This just passes through ``*args`` and ``**kwargs`` after temporarily setting the download_file remote timeout to the local ``iers.conf.remote_timeout`` value. uremote_timeoutN(R tdatatconftset_temptremote_timeoutt download_file(targstkwargs((s6lib/python2.7/site-packages/astropy/utils/iers/iers.pyRHstIERSStaleWarningcB`seZRS((t__name__t __module__(((s6lib/python2.7/site-packages/astropy/utils/iers/iers.pyRSstConfcB`sVeZdZejedƒZejddƒZejedƒZ ejddƒZ RS(u< Configuration parameters for `astropy.utils.iers`. u†Enable auto-downloading of the latest IERS data. If set to False then the local IERS-B file will be used by default. Default is True.g>@uMMaximum age (days) of predictive data before auto-downloading. Default is 30.u(URL for auto-downloading IERS file data.g$@u4Remote timeout downloading IERS file data (seconds).( RRt__doc__t_configt ConfigItemtTruet auto_downloadt auto_max_aget IERS_A_URLt iers_auto_urlR(((s6lib/python2.7/site-packages/astropy/utils/iers/iers.pyRWs   tIERSRangeErrorcB`seZdZRS(uJ Any error for when dates are outside of the valid range for IERS (RRR(((s6lib/python2.7/site-packages/astropy/utils/iers/iers.pyR!mstIERScB`s¡eZdZd Zed ed„ƒZed„ƒZdd„Z ded„Z ded„Z d„Z d d„Z d „Zd „Zd „Zed „ƒZRS(u¶Generic IERS table class, defining interpolation functions. Sub-classed from `astropy.table.QTable`. The table should hold columns 'MJD', 'UT1_UTC', and 'PM_x'/'PM_y'. cK`s„|dk s|jdkr}|dk rht|ƒjrU|jdt|d|ƒƒqh|jd|ƒn|j||_n|jS(u¶Open an IERS table, reading it from a file if not loaded before. Parameters ---------- file : str or None full local or network path to the ascii file holding IERS data, for passing on to the ``read`` class methods (further optional arguments that are available for some IERS subclasses can be added). If None, use the default location from the ``read`` class method. cache : bool Whether to use cache. Defaults to False, since IERS files are regularly updated. Returns ------- An IERS table class instance Notes ----- On the first call in a session, the table will be memoized (in the ``iers_table`` class attribute), and further calls to ``open`` will return this stored table if ``file=None`` (the default). If a table needs to be re-read from disk, pass on an explicit file location or use the (sub-class) close method and re-open. If the location is a network location it is first downloaded via download_file. For the IERS class itself, an IERS_B sub-class instance is opened. tfiletcacheN(tNonet iers_tableRtnetloctupdateRtread(tclsR#R$R((s6lib/python2.7/site-packages/astropy/utils/iers/iers.pytopen|s" cC`s d|_dS(uÃRemove the IERS table from the class. This allows the table to be re-read from disk during one's session (e.g., if one finds it is out of date and has updated the file). N(R%R&(R*((s6lib/python2.7/site-packages/astropy/utils/iers/iers.pytclose§sgcC`sdy|jj|jj}}Wntk r0nXtj|t|ƒ}|t||}||fS(uôTurn a time to MJD, returning integer and fractional parts. Parameters ---------- jd1 : float, array, or Time first part of two-part JD, or Time object jd2 : float or array, optional second part of two-part JD. Default is 0., ignored if jd1 is `~astropy.time.Time`. Returns ------- mjd : float or array integer part of MJD utc : float or array fractional part of MJD (tutctjd1tjd2t ExceptiontnptfloortMJD_ZERO(tselfR.R/tmjdR-((s6lib/python2.7/site-packages/astropy/utils/iers/iers.pytmjd_utc°s cC`s(|j||dg|r!|jndƒS(uÔInterpolate UT1-UTC corrections in IERS Table for given dates. Parameters ---------- jd1 : float, float array, or Time object first part of two-part JD, or Time object jd2 : float or float array, optional second part of two-part JD. Default is 0., ignored if jd1 is `~astropy.time.Time`. return_status : bool Whether to return status values. If False (default), raise ``IERSRangeError`` if any time is out of the range covered by the IERS table. Returns ------- ut1_utc : float or float array UT1-UTC, interpolated in IERS Table status : int or int array Status values (if ``return_status``=``True``):: ``iers.FROM_IERS_B`` ``iers.FROM_IERS_A`` ``iers.FROM_IERS_A_PREDICTION`` ``iers.TIME_BEFORE_IERS_RANGE`` ``iers.TIME_BEYOND_IERS_RANGE`` uUT1_UTCN(t _interpolatetut1_utc_sourceR%(R4R.R/t return_status((s6lib/python2.7/site-packages/astropy/utils/iers/iers.pytut1_utcÉscC`s+|j||ddg|r$|jndƒS(uKInterpolate polar motions from IERS Table for given dates. Parameters ---------- jd1 : float, float array, or Time object first part of two-part JD, or Time object jd2 : float or float array, optional second part of two-part JD. Default is 0., ignored if jd1 is `~astropy.time.Time`. return_status : bool Whether to return status values. If False (default), raise ``IERSRangeError`` if any time is out of the range covered by the IERS table. Returns ------- PM_x : Quantity with angle units x component of polar motion for the requested times PM_y : Quantity with angle units y component of polar motion for the requested times status : int or int array Status values (if ``return_status``=``True``):: ``iers.FROM_IERS_B`` ``iers.FROM_IERS_A`` ``iers.FROM_IERS_A_PREDICTION`` ``iers.TIME_BEFORE_IERS_RANGE`` ``iers.TIME_BEYOND_IERS_RANGE`` uPM_xuPM_yN(R7t pm_sourceR%(R4R.R/R9((s6lib/python2.7/site-packages/astropy/utils/iers/iers.pytpm_xyçscC`s(tj||kƒr$tdƒ‚ndS(uÛ Check that the indices from interpolation match those after clipping to the valid table range. This method gets overridden in the IERS_Auto class because it has different requirements. u8(some) times are outside of range covered by IERS table.N(R1tanyR!(R4t indices_origtindices_clippedt max_input_mjd((s6lib/python2.7/site-packages/astropy/utils/iers/iers.pyt_check_interpolate_indicesscC`sG|j||ƒ\}}t|dƒ p4|jdk}|rdtj|gƒ}tj|gƒ}n|j|ƒtj|dj|ddƒ}tj|dt |ƒdƒ} | d} |d| j|d| j} } g} xÈ|D]À}||| ||| }}||}|dkr5||j ƒ8}n||| || | |}||d||dk<||d||t |ƒk<|rœ|d}n| j |ƒqéW|r || ƒ}t ||dkeZdZdZedddd„ƒZd„Zd„ZRS(uIERS Table class targeted to IERS B, provided by IERS itself. These are final values; see http://www.iers.org/ Notes ----- If the package IERS B file (```iers.IERS_B_FILE``) is out of date, a new version can be downloaded from ``iers.IERS_B_URL``. icC`s[|dkrt}n|dkr*t}ntj|ddd|d|ƒ}||jƒƒS(u9Read IERS-B table from a eopc04_iau2000.* file provided by IERS. Parameters ---------- file : str full path to ascii file holding IERS-B data. Defaults to package version, ``iers.IERS_B_FILE``. readme : str full path to ascii file holding CDS-style readme. Defaults to package version, ``iers.IERS_B_README``. data_start : int starting row. Default is 14, appropriate for standard IERS files. Returns ------- ``IERS_B`` class instance RtucdsRut data_startN(R%t IERS_B_FILEt IERS_B_READMERR)Rx(R*R#RuR€tiers_b((s6lib/python2.7/site-packages/astropy/utils/iers/iers.pyR)ìs     cC`stj|ƒtS(u1Set UT1-UTC source flag for entries in IERS table(R1RyRz(R4RS((s6lib/python2.7/site-packages/astropy/utils/iers/iers.pyR8 scC`stj|ƒtS(u,Set PM source flag for entries in IERS table(R1RyRz(R4RS((s6lib/python2.7/site-packages/astropy/utils/iers/iers.pyR;sN( RRRR%R&ReR)R8R;(((s6lib/python2.7/site-packages/astropy/utils/iers/iers.pyRßs   t IERS_AutocB`sDeZdZdZed„ƒZd„Zd„Zed„ƒZ RS(up Provide most-recent IERS data and automatically handle downloading of updated values as necessary. cC`sòtjstjƒ|_|jS|jdk rV|jjjdƒtjkrV|jSnyt tjdt ƒ}WnKt k r¼}t t djtjt|ƒƒƒƒtjƒ|_|jSX|jd|ƒ|_ttjƒ|jjd<|jS(uÖIf the configuration setting ``astropy.utils.iers.conf.auto_download`` is set to True (default), then open a recent version of the IERS-A table with predictions for UT1-UTC and polar motion out to approximately one year from now. If the available version of this file is older than ``astropy.utils.iers.conf.auto_max_age`` days old (or non-existent) then it will be downloaded over the network and cached. If the configuration setting ``astropy.utils.iers.conf.auto_download`` is set to False then ``astropy.utils.iers.IERS()`` is returned. This is normally the IERS-B table that is supplied with astropy. On the first call in a session, the table will be memoized (in the ``iers_table`` class attribute), and further calls to ``open`` will return this stored table. Returns ------- `~astropy.table.QTable` instance with IERS (Earth rotation) data columns udata_urlR$u-failed to download {}, using local IERS-B: {}R#N(RRR"R+R&R%RotgetR RRR0RR RttstrR)(R*tfilenameterr((s6lib/python2.7/site-packages/astropy/utils/iers/iers.pyR+s    cC`sx|jd}tjdk r%tjntjtjƒj}||krt|jj ||krtt t j |ƒƒ‚ndS(uFCheck that the indices from interpolation match those after clipping to the valid table range. The IERS_Auto class is exempted as long as it has sufficiently recent available data so the clipped interpolation is always within the confidence bounds of current Earth rotation knowledge. upredictive_mjdN( RoRRR%R1tfinfotfloatRORdR5t ValueErrortINTERPOLATE_ERRORRt(R4R>R?R@tpredictive_mjdR((s6lib/python2.7/site-packages/astropy/utils/iers/iers.pyRALs   c C`s tj|ƒ}|jj}|jd}|jd}tjdk rMtjntjtj ƒj}|dkr}t dƒ‚n||kr|||kry&t tj ƒt tj dtƒ}Wn9tk rú}ttdjtj t|ƒƒƒƒdSX|jjd|ƒ} | d d |d d krætj| d j|d d ƒ} t|ƒ|} | | | | !||)| d | | |d d d tjkrªt dƒ‚nx#| | | D]} |j| ƒq¹W|jj| jƒqttdjtjƒƒƒndS(u—Potentially update the IERS table in place depending on the requested time values in ``mjd`` and the time span of the table. For IERS_Auto the behavior is that the table is refreshed from the IERS server if both the following apply: - Any of the requested IERS values are predictive. The IERS-A table contains predictive data out for a year after the available definitive values. - The first predictive values are at least ``conf.auto_max_age days`` old. In other words the IERS-A table was created by IERS long enough ago that it can be considered stale for predictions. upredictive_indexupredictive_mjdi uAIERS auto_max_age configuration value must be larger than 10 daysR$ufailed to download {}: {}. A coordinate or time-related calculation might be compromised or fail because the dates are not covered by the available IERS file. See the "IERS data access" section of the astropy documentation for additional information on working offline.NR#uMJDiÿÿÿÿRBurightgð?u0unexpected gap in MJD when refreshing IERS tableumIERS_Auto predictive values are older than {} days but downloading the latest table did not find newer values(R1RORdR5RoRRR%R‰RŠR‹R R RRR0RR RtR†t __class__R)RGRHRJtutdtadd_rowR(R( R4R5R@tnow_mjdtfpiRRR‡Rˆt new_tabletnew_fpit n_replacetrow((s6lib/python2.7/site-packages/astropy/utils/iers/iers.pyRF\s<       +cC`stjƒ}|d|dj}tj|dj|dddƒ}tj|dj|dddƒ}|||!}t|ƒ}|dkrÿtj|d| |djƒs½tdƒ‚n|d j|d|*|d j|d |*|d j|d |*n|S(uMSubstitute IERS B values with those from a real IERS B table. IERS-A has IERS-B values included, but for reasons unknown these do not match the latest IERS-B values (see comments in #4436). Here, we use the bundled astropy IERS-B table to overwrite the values in the downloaded IERS-A table. uMJDu UT1_UTC_BiRBuleftiÿÿÿÿurightuAunexpected mismatch when copying IERS-B values into IERS-A table.uUT1_UTCuPM_xuPM_X_BuPM_yuPM_Y_B( RR+RiR1RGRHRJtallcloseR‹(R*RqRƒtmjd_bRURTtn_iers_b((s6lib/python2.7/site-packages/astropy/utils/iers/iers.pyRj§s ##   !N( RRRR%R&ReR+RARFRj(((s6lib/python2.7/site-packages/astropy/utils/iers/iers.pyR„s 2  K(5Rt __future__RRRRtwarningsRRt ImportErrort urllib.parsetnumpyR1tRRRRRqRR t utils.dataR R R tutils.exceptionsR t__all__RvRRwRt IERS_B_URLR‚RzR{R|RMRNR3RŒRRtConfigNamespaceRRt IndexErrorR!R"RhRR„R)(((s6lib/python2.7/site-packages/astropy/utils/iers/iers.pyt sV"            ôx4³