earths 0.0.3

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
70

use crate::lighting::solar_position::SolarPosition;
#[derive(Debug, Clone, Copy, PartialEq)]
pub enum DaylightState {
    Day,
    CivilTwilight,
    NauticalTwilight,
    AstronomicalTwilight,
    Night,
}
pub struct DayNightCycle {
    pub juliandate: f64,
}
impl DayNightCycle {
    pub fn new(juliandate: f64) -> Self {
        Self { juliandate }
    }
    pub fn stateat(&self, latdeg: f64, londeg: f64) -> DaylightState {
        let sun = SolarPosition::compute(self.juliandate, latdeg, londeg);
        match sun.elevationdeg {
            e if e > 0.0 => DaylightState::Day,
            e if e > -6.0 => DaylightState::CivilTwilight,
            e if e > -12.0 => DaylightState::NauticalTwilight,
            e if e > -18.0 => DaylightState::AstronomicalTwilight,
            e if e <= -18.0 => DaylightState::Night,
            e => {
                debug_assert!(e.is_nan());
                DaylightState::Night
            }
        }
    }
    pub fn terminatorpoints(&self, numpoints: usize) -> Vec<[f64; 2]> {
        let mut points = Vec::with_capacity(numpoints);
        for i in 0..numpoints {
            let lon = -180.0 + 360.0 * i as f64 / numpoints as f64;
            let mut latlow = -90.0f64;
            let mut lathigh = 90.0f64;
            for iter in 0..30 {
                let latmid = (latlow + lathigh) / 2.0;
                let sun = SolarPosition::compute(self.juliandate, latmid, lon);
                if sun.elevationdeg > 0.0 {
                    latlow = latmid;
                } else {
                    lathigh = latmid;
                }
                if (lathigh - latlow).abs() < 1e-8 || iter == 29 {
                    break;
                }
            }
            points.push([(latlow + lathigh) / 2.0, lon]);
        }
        points
    }
    pub fn ambientlight(&self, latdeg: f64, londeg: f64) -> f64 {
        let sun = SolarPosition::compute(self.juliandate, latdeg, londeg);
        let e = sun.elevationdeg;
        if e > 10.0 {
            1.0
        } else if e > 0.0 {
            0.5 + 0.5 * (e / 10.0) * (e / 10.0)
        } else if e > -6.0 {
            0.5 * (10.0f64).powf(e / 6.0)
        } else if e > -12.0 {
            0.05 * (10.0f64).powf((e + 6.0) / 6.0)
        } else if e > -18.0 {
            0.005 * (10.0f64).powf((e + 12.0) / 6.0)
        } else {
            0.001
        }
    }
}