/rust/registry/src/index.crates.io-1949cf8c6b5b557f/sofars-0.6.1/src/erst/ee00.rs
Line | Count | Source |
1 | | use super::eect00; |
2 | | |
3 | | /// Equation of the equinoxes, IAU 2000 |
4 | | /// |
5 | | /// The equation of the equinoxes, compatible with IAU 2000 resolutions, |
6 | | /// given the nutation in longitude and the mean obliquity. |
7 | | /// |
8 | | /// This function is part of the International Astronomical Union's |
9 | | /// SOFA (Standards of Fundamental Astronomy) software collection. |
10 | | /// |
11 | | /// Status: canonical model. |
12 | | /// |
13 | | /// Given: |
14 | | /// ``` |
15 | | /// date1,date2 double TT as a 2-part Julian Date (Note 1) |
16 | | /// epsa double mean obliquity (Note 2) |
17 | | /// dpsi double nutation in longitude (Note 3) |
18 | | /// ``` |
19 | | /// Returned (function value): |
20 | | /// ``` |
21 | | /// |
22 | | /// double equation of the equinoxes (Note 4) |
23 | | /// ``` |
24 | | /// Notes: |
25 | | /// |
26 | | /// 1) The TT date date1+date2 is a Julian Date, apportioned in any |
27 | | /// convenient way between the two arguments. For example, |
28 | | /// JD(TT)=2450123.7 could be expressed in any of these ways, |
29 | | /// among others: |
30 | | /// ``` |
31 | | /// date1 date2 |
32 | | /// |
33 | | /// 2450123.7 0.0 (JD method) |
34 | | /// 2451545.0 -1421.3 (J2000 method) |
35 | | /// 2400000.5 50123.2 (MJD method) |
36 | | /// 2450123.5 0.2 (date & time method) |
37 | | /// ``` |
38 | | /// The JD method is the most natural and convenient to use in |
39 | | /// cases where the loss of several decimal digits of resolution |
40 | | /// is acceptable. The J2000 method is best matched to the way |
41 | | /// the argument is handled internally and will deliver the |
42 | | /// optimum resolution. The MJD method and the date & time methods |
43 | | /// are both good compromises between resolution and convenience. |
44 | | /// |
45 | | /// 2) The obliquity, in radians, is mean of date. |
46 | | /// |
47 | | /// 3) The result, which is in radians, operates in the following sense: |
48 | | /// |
49 | | /// Greenwich apparent ST = GMST + equation of the equinoxes |
50 | | /// |
51 | | /// 4) The result is compatible with the IAU 2000 resolutions. For |
52 | | /// further details, see IERS Conventions 2003 and Capitaine et al. |
53 | | /// (2002). |
54 | | /// |
55 | | /// Called: |
56 | | /// ``` |
57 | | /// |
58 | | /// iauEect00 equation of the equinoxes complementary terms |
59 | | /// ``` |
60 | | /// References: |
61 | | /// |
62 | | /// Capitaine, N., Wallace, P.T. and McCarthy, D.D., "Expressions to |
63 | | /// implement the IAU 2000 definition of UT1", Astronomy & |
64 | | /// Astrophysics, 406, 1135-1149 (2003) |
65 | | /// |
66 | | /// McCarthy, D. D., Petit, G. (eds.), IERS Conventions (2003), |
67 | | /// IERS Technical Note No. 32, BKG (2004) |
68 | 0 | pub fn ee00(date1: f64, date2: f64, epsa: f64, dpsi: f64) -> f64 { |
69 | | /* Equation of the equinoxes. */ |
70 | 0 | dpsi * epsa.cos() + eect00(date1, date2) |
71 | 0 | } |