sofars 0.6.0

Pure Rust implementation of the IAU SOFA library
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36

use crate::consts::D2PI;
use std::ops::Rem;

/// Fundamental argument, IERS Conventions (2003):
/// mean longitude of Saturn.
///
/// Status:  canonical model.
///
/// Given:
///    t     f64    TDB, Julian centuries since J2000.0 (Note 1)
///
/// Returned (function value):
///          f64    mean longitude of Saturn, radians (Note 2)
///
/// Notes:
///
/// 1) Though t is strictly TDB, it is usually more convenient to use
///    TT, which makes no significant difference.
///
/// 2) The expression used is as adopted in IERS Conventions (2003) and
///    comes from Souchay et al. (1999) after Simon et al. (1994).
///
/// References:
///
///    McCarthy, D. D., Petit, G. (eds.), IERS Conventions (2003),
///    IERS Technical Note No. 32, BKG (2004)
///
///    Simon, J.-L., Bretagnon, P., Chapront, J., Chapront-Touze, M.,
///    Francou, G., Laskar, J. 1994, Astron.Astrophys. 282, 663-683
///
///    Souchay, J., Loysel, B., Kinoshita, H., Folgueira, M. 1999,
///    Astron.Astrophys.Supp.Ser. 135, 111
pub fn fasa03(t: f64) -> f64 {
    /* Mean longitude of Saturn (IERS Conventions 2003). */
    (0.874016757 + 21.3299104960 * t).rem(D2PI)
}