earths 0.0.4

High-fidelity Earth simulation engine — orbit, atmosphere, geology, hydrology, biosphere, terrain, lighting, rendering, satellites, and temporal systems with full scientific coupling
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
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69

use crate::temporal::calendar::SECONDSPERDAY;
use sciforge::hub::domain::astronomy::celestial::gravitational_force;
use sciforge::hub::domain::astronomy::orbits::{
    angular_momentum, apoapsis, escape_velocity, kepler_period, orbital_energy, periapsis, vis_viva,
};
use sciforge::hub::domain::common::constants::{AU, EARTH_MASS, G, SOLAR_MASS};
pub const SEMIMAJORAXIS: f64 = 1.000001018 * AU;
pub const ECCENTRICITY: f64 = 0.0167086;
pub const INCLINATIONDEG: f64 = 0.00005;
pub const LONGITUDEASCENDINGNODEDEG: f64 = -11.26064;
pub const ARGUMENTPERIHELIONDEG: f64 = 114.20783;
const MUSUN: f64 = G * SOLAR_MASS;
pub struct EarthOrbit {
    pub semimajoraxism: f64,
    pub eccentricity: f64,
    pub trueanomalyrad: f64,
}
impl Default for EarthOrbit {
    fn default() -> Self {
        Self::new()
    }
}
impl EarthOrbit {
    pub fn new() -> Self {
        Self {
            semimajoraxism: SEMIMAJORAXIS,
            eccentricity: ECCENTRICITY,
            trueanomalyrad: 0.0,
        }
    }
    pub fn orbitalperiods(&self) -> f64 {
        kepler_period(self.semimajoraxism, MUSUN)
    }
    pub fn orbitalperioddays(&self) -> f64 {
        self.orbitalperiods() / SECONDSPERDAY
    }
    pub fn velocityatdistance(&self, r: f64) -> f64 {
        vis_viva(MUSUN, r, self.semimajoraxism)
    }
    pub fn perihelionm(&self) -> f64 {
        periapsis(self.semimajoraxism, self.eccentricity)
    }
    pub fn aphelionm(&self) -> f64 {
        apoapsis(self.semimajoraxism, self.eccentricity)
    }
    pub fn specific_orbital_energy(&self) -> f64 {
        orbital_energy(MUSUN, self.semimajoraxism)
    }
    pub fn specific_angular_momentum(&self) -> f64 {
        angular_momentum(MUSUN, self.semimajoraxism, self.eccentricity)
    }
    pub fn escape_velocity_at_surface() -> f64 {
        escape_velocity(G * EARTH_MASS, EARTH_RADIUS)
    }
    pub fn gravitational_force_sun(&self) -> f64 {
        let r = self.current_radius();
        gravitational_force(SOLAR_MASS, EARTH_MASS, r)
    }
    pub fn current_radius(&self) -> f64 {
        let a = self.semimajoraxism;
        let e = self.eccentricity;
        a * (1.0 - e * e) / (1.0 + e * self.trueanomalyrad.cos())
    }
    pub fn mean_orbital_velocity(&self) -> f64 {
        let period = self.orbitalperiods();
        2.0 * std::f64::consts::PI * self.semimajoraxism / period
    }
}
use sciforge::hub::domain::common::constants::EARTH_RADIUS;