B îq\özã@s–dZddlmZyddlmZWn ek r@ddlmZYnXddlZddlm Z ddlm Z ddl mZmZddlmZmZdd lmZdd lmZd d d dddddddddddddddddgZdZd Zd!Zed"ƒZed#ƒZd$Zed%ƒZdZd&Zd'Z d(Z!d)Z"d*Z#d+Z$d,d-„Z%Gd.d„deƒZ&Gd/d „d e j'ƒZ(e(ƒZ)Gd0d„de*ƒZ+Gd1d „d eƒZ,Gd2d„de,ƒZ-Gd3d„de,ƒZ.Gd4d„de-ƒZ/e.j0e,_0dS)5a¬ 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`). é)Úwarn)ÚurlparseN)Úconfig)Úunits)ÚTableÚQTable)Úget_pkg_data_filenameÚclear_download_cache)Úutils)ÚAstropyWarningÚConfÚconfÚIERSÚIERS_BÚIERS_AÚ IERS_AutoÚ FROM_IERS_BÚ FROM_IERS_AÚFROM_IERS_A_PREDICTIONÚTIME_BEFORE_IERS_RANGEÚTIME_BEYOND_IERS_RANGEÚ IERS_A_FILEÚ IERS_A_URLÚIERS_A_URL_MIRRORÚ IERS_A_READMEÚ IERS_B_FILEÚ IERS_B_URLÚ IERS_B_READMEÚIERSRangeErrorÚIERSStaleWarningzfinals2000A.allz.http://maia.usno.navy.mil/ser7/finals2000A.allz/http://toshi.nofs.navy.mil/ser7/finals2000A.allzdata/ReadMe.finals2000Azdata/eopc04_IAU2000.62-nowz@zMMaximum age (days) of predictive data before auto-downloading. Default is 30.z(URL for auto-downloading IERS file data.z/Mirror URL for auto-downloading IERS file data.g$@z4Remote timeout downloading IERS file data (seconds).N) r+r,r-Ú__doc__Ú_configZ ConfigItemÚ auto_downloadÚ auto_max_agerÚ iers_auto_urlrÚiers_auto_url_mirrorr$r)r)r)r*r Vs c@seZdZdZdS)rzJ Any error for when dates are outside of the valid range for IERS N)r+r,r-r.r)r)r)r*rosc@s”eZdZdZdZeddd„ƒZedd„ƒZd d d „Zd!d d „Z d"d d„Z d#dd„Z dd„Z d$dd„Z dd„Zdd„Zdd„Zdd„Zedd„ƒZdS)%rzÍGeneric IERS table class, defining interpolation functions. Sub-classed from `astropy.table.QTable`. The table should hold columns 'MJD', 'UT1_UTC', 'dX_2000A'/'dY_2000A', and 'PM_x'/'PM_y'. NFcKsZ|dk s|jdkrT|dk rFt|ƒjr:|jt||ddn |j|d|jf|Ž|_|jS)a¶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. N)Úcache)Úfile)Ú iers_tablerZnetlocÚupdater&Úread)Úclsr5r4r(r)r)r*Úopen~s"  z IERS.opencCs d|_dS)zÃ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)r6)r9r)r)r*Úclose©sz IERS.closeçcCsVy|jj|jj}}Wntk r*YnXt |t|¡}|t||}||fS)aô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 )ÚutcÚjd1Újd2Ú ExceptionÚnpZfloorÚMJD_ZERO)Úselfr>r?Úmjdr=r)r)r*Úmjd_utc²sz IERS.mjd_utccCs| ||dg|r|jnd¡S)aÔ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`` ÚUT1_UTCN)Ú _interpolateÚut1_utc_source)rCr>r?Ú return_statusr)r)r*Úut1_utcÌs z IERS.ut1_utccCs| ||ddg|r|jnd¡S)a.Interpolate CIP 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 0., ignored if jd1 is 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 ------- D_x : Quantity with angle units x component of CIP correction for the requested times D_y : Quantity with angle units y component of CIP correction 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`` ÚdX_2000AÚdY_2000AN)rGÚ dcip_source)rCr>r?rIr)r)r*Údcip_xyêsz IERS.dcip_xycCs| ||ddg|r|jnd¡S)aKInterpolate 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`` ÚPM_xÚPM_yN)rGÚ pm_source)rCr>r?rIr)r)r*Úpm_xy sz IERS.pm_xycCst ||k¡rtdƒ‚dS)zÛ 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. z8(some) times are outside of range covered by IERS table.N)rAÚanyr)rCÚ indices_origÚindices_clippedÚ max_input_mjdr)r)r*Ú_check_interpolate_indices)szIERS._check_interpolate_indicescCs¾| ||¡\}}t|dƒ p$|jdk}|rBt |g¡}t |g¡}| |¡tj|dj|dd}t |dt |ƒd¡} | d} |d| j|d| j} } g} x¢|D]š}||| ||| }}||}|dkrâ||  ¡8}||| || | |}||d||dk<||d||t |ƒk<|r8|d}|   |¡q¨W|rŒ|| ƒ}t ||dk<t ||t |ƒk<|r~|d}|   |¡| S| | |t |¡¡t | ƒdkr¶| dS| SdS) NZ __array__rÚMJDÚright)Úsider rFr")rEÚhasattrÚndimrAZarrayÚ_refresh_table_as_neededÚ searchsortedÚvalueZclipÚlenÚroundÚappendrrrWÚmax)rCr>r?ÚcolumnsÚsourcerDr=Z is_scalarÚiÚi1Úi0Zmjd_0Zmjd_1ZresultsÚcolumnZval_0Zval_1Zd_valÚvalZstatusr)r)r*rG3s@       zIERS._interpolatecCsdS)zó Potentially update the IERS table in place depending on the requested time values in ``mdj`` and the time span of the table. The base behavior is not to update the table. ``IERS_Auto`` overrides this method. Nr))rCrDr)r)r*r]lszIERS._refresh_table_as_neededcCs t |¡S)z2Source for UT1-UTC. To be overridden by subclass.)rAÚ zeros_like)rCrfr)r)r*rHtszIERS.ut1_utc_sourcecCs t |¡S)z9Source for CIP correction. To be overridden by subclass.)rArk)rCrfr)r)r*rMxszIERS.dcip_sourcecCs t |¡S)z7Source for polar motion. To be overridden by subclass.)rArk)rCrfr)r)r*rQ|szIERS.pm_sourcecCs0ddlm}y|jStk r*| ¡SXdS)z“ Property to provide the current time, but also allow for explicitly setting the _time_now attribute for testing purposes. r)ÚTimeN)Z astropy.timerlZ _time_nowr@Znow)rCrlr)r)r*Útime_now€s  z IERS.time_now)NF)r<)r<F)r<F)r<F)N)r+r,r-r.r6Ú classmethodr:r;rErJrNrRrWrGr]rHrMrQÚpropertyrmr)r)r)r*rus  *     9c@sReZdZdZdZedd„ƒZedd„ƒZeddd„ƒZd d „Z d d „Z d d„Z dS)ramIERS Table class targeted to IERS A, provided by USNO. These include rapid turnaround and predicted times. See http://maia.usno.navy.mil/ Notes ----- The IERS A file is not part of astropy. It can be downloaded from ``iers.IERS_A_URL`` or ``iers.IERS_A_URL_MIRROR``. See ``iers.__doc__`` for instructions on use in ``Time``, etc. NcCsÎ||dj|dj@}| |¡}t |dj|dj|dj¡|d<|dj|d_t |dj|djd¡|d<t |dj|d j|dj¡|d <|d j|d _t |d j|d j|d j¡|d <|d j|d _t |dj|djd¡|d<t |dj|dj|dj¡|d<|dj|d_t |dj|dj|dj¡|d<|dj|d_t |dj|djd¡|d<|ddk|ddkB}t t |¡¡|jd<|d|jd|jd<|S)zZ Return a new table with appropriate combination of IERS_A and B columns. Z UT1_UTC_AZ PolPMFlag_AÚ UT1_UTC_BrFZ UT1Flag_AÚBÚUT1FlagÚPM_X_BZPM_x_ArOÚPM_Y_BZPM_y_ArPÚ PolPMFlagZ dX_2000A_BZ dX_2000A_ArKZ dY_2000A_BZ dY_2000A_ArLZ NutFlag_AÚNutFlagÚPÚpredictive_indexrXÚpredictive_mjd) ÚmaskÚ_substitute_iers_brAÚwherer%ZunitÚminZ flatnonzeroÚmeta)r9Úiers_aÚtableZ is_predictiver)r)r*Ú_combine_a_b_columnsœsH             zIERS_A._combine_a_b_columnscCs|S)Nr))r9r€r)r)r*r{ÛszIERS_A._substitute_iers_bcCsR|dkr t}|dkrt}tj|d|d}| |¡}||jd<||jd<|| ¡ƒS)a¾Read IERS-A table from a finals2000a.* file provided by USNO. Parameters ---------- file : str full path to ascii file holding IERS-A data. Defaults to ``iers.IERS_A_FILE``. readme : str full path to ascii file holding CDS-style readme. Defaults to package version, ``iers.IERS_A_README``. Returns ------- ``IERS_A`` class instance NÚcds)ÚformatÚreadmeZ data_pathZ readme_path)rrrr8rr~Úfilled)r9r5r„rr€r)r)r*r8às   z IERS_A.readcCs6|d|}t |¡t}t||dk<t||dk<|S)z1Set UT1-UTC source flag for entries in IERS tablerrÚIrw)rAÚ ones_likerrr)rCrfZut1flagrer)r)r*rHs    zIERS_A.ut1_utc_sourcecCs6|d|}t |¡t}t||dk<t||dk<|S)z8Set CIP correction source flag for entries in IERS tablervr†rw)rAr‡rrr)rCrfZnutflagrer)r)r*rM s    zIERS_A.dcip_sourcecCs6|d|}t |¡t}t||dk<t||dk<|S)z6Set polar motion source flag for entries in IERS tablerur†rw)rAr‡rrr)rCrfZpmflagrer)r)r*rQs    zIERS_A.pm_source)NN) r+r,r-r.r6rnrr{r8rHrMrQr)r)r)r*rs  ?  !c@s:eZdZdZdZed dd„ƒZdd„Zdd „Zd d „Z dS) raIERS 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``. NécCs6|dkr t}|dkrt}tj|d||d}|| ¡ƒS)a9Read 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 Nr‚)rƒr„Ú data_start)rrrr8r…)r9r5r„r‰Úiers_br)r)r*r8(s z IERS_B.readcCst |¡tS)z1Set UT1-UTC source flag for entries in IERS table)rAr‡r)rCrfr)r)r*rHFszIERS_B.ut1_utc_sourcecCst |¡tS)z8Set CIP correction source flag for entries in IERS table)rAr‡r)rCrfr)r)r*rMJszIERS_B.dcip_sourcecCst |¡tS)z,Set PM source flag for entries in IERS table)rAr‡r)rCrfr)r)r*rQNszIERS_B.pm_source)NNrˆ) r+r,r-r.r6rnr8rHrMrQr)r)r)r*rs  c@s<eZdZdZdZedd„ƒZdd„Zdd„Zed d „ƒZ dS) rzp Provide most-recent IERS data and automatically handle downloading of updated values as necessary. Nc Csútjst ¡|_|jStjtjf}|jdk rD|jj d¡|krD|jSd}g}xT|D]L}yt |dd}Wn0t k r–}z|  t |ƒ¡Wdd}~XYqRXd}PqRW|sÖt td d |¡d |¡¡ƒƒt ¡|_|jS|j|d |_t |ƒ|jjd<|jS) aÖ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 NÚdata_urlFT)r4z-failed to download {}, using local IERS-B: {}z and ú;)r5)r r0rr:r6r2r3r~Úgetr&r@rbÚstrrr rƒÚjoinr8)r9Úall_urlsÚ dl_successÚerr_listÚurlÚfilenameÚerrr)r)r*r:Zs2      zIERS_Auto.opencCsP|jd}tjdk rtjn t t¡j}||krL|jj||krLt t   |¡ƒ‚dS)aFCheck 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. ryN) r~r r1rAÚfinfoÚfloatrcrmrDÚ ValueErrorÚINTERPOLATE_ERRORrƒ)rCrTrUrVryr1r)r)r*rW™s   z$IERS_Auto._check_interpolate_indicesc Csôt |¡}|jj}|jd}|jd}tjdk r6tjn t t¡j}|dkrRt dƒ‚||krð|||krðtj tj f}d}g} x\|D]T} yt | ƒt | dd} Wn0tk rÐ} z|  t| ƒ¡Wdd} ~ XYq„Xd}Pq„W|sttd  d  |¡d  | ¡¡ƒƒdS|jj| d } t| ƒ| jd <| dd|ddkrÜtj| dj|dd}t|ƒ|}| |||…||d…<| d|||dddtjkr¦t dƒ‚x$| ||d…D]}| |¡q¸W|j | j¡nttd tj¡ƒƒdS)a—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. rxryNé zAIERS auto_max_age configuration value must be larger than 10 daysFT)r4afailed 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.z and rŒ)r5r‹rXr"rY)rZgð?z0unexpected gap in MJD when refreshing IERS tablezmIERS_Auto predictive values are older than {} days but downloading the latest table did not find newer values)rArcrmrDr~r r1r–r—r˜r2r3r r&r@rbrŽrr rƒrÚ __class__r8r^r_r`ÚuÚdZadd_rowr7r)rCrDrVZnow_mjdZfpiryr1rr‘r’r“r”r•Z new_tableZnew_fpiZ n_replaceÚrowr)r)r*r]©sN         (z"IERS_Auto._refresh_table_as_neededcCsØt ¡}|d|dj}tj|dj|ddd}tj|dj|ddd}|||…}t|ƒ}|dkrÔt |dd|…|dj¡s’td ƒ‚|d j|dd|…<|d j|d d|…<|d j|dd|…<|S)aMSubstitute 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. rXrprÚleft)rZr"rYNzAunexpected mismatch when copying IERS-B values into IERS-A table.rFrOrsrPrt) rr:rzrAr^r_r`Zallcloser˜)r9r€rŠZmjd_brhrgZn_iers_br)r)r*r{s  zIERS_Auto._substitute_iers_b) r+r,r-r.r6rnr:rWr]r{r)r)r)r*rSs  ?W)1r.ÚwarningsrrÚ ImportErrorZ urllib.parseZnumpyrAZastropyrr/rrœZ astropy.tablerrZastropy.utils.datarr r Zastropy.utils.exceptionsr Ú__all__rrrrrrrrrrrrrBr™r&rZConfigNamespacer r Ú IndexErrorrrrrrr8r)r)r)r*Ú s^      8M