enceladuss 0.0.3

Enceladus celestial simulation crate for the MilkyWay SolarSystem workspace
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

use crate::{ENCELADUS_DAY_HOURS, MAX_SURFACE_TEMP_K, MIN_SURFACE_TEMP_K};

#[derive(Debug, Clone, Copy, PartialEq)]
pub struct ThermalEnvironment {
    pub local_time_hours: f64,
    pub latitude_deg: f64,
    pub shadowed: bool,
}

impl ThermalEnvironment {
    pub fn surface_temperature_k(&self) -> f64 {
        if self.shadowed {
            return MIN_SURFACE_TEMP_K;
        }

        let diurnal_phase = (self.local_time_hours / ENCELADUS_DAY_HOURS).fract();
        let solar_term = -((diurnal_phase * std::f64::consts::TAU).cos());
        let latitude_cooling = self.latitude_deg.abs() * 0.15;
        let tiger_stripe_heating = if self.latitude_deg < -55.0 { 40.0 } else { 0.0 };
        (55.0 + 20.0 * solar_term.max(-0.95) - latitude_cooling + tiger_stripe_heating)
            .clamp(MIN_SURFACE_TEMP_K, MAX_SURFACE_TEMP_K)
    }
}

pub fn thermal_cycle_span_k(latitude_deg: f64) -> f64 {
    let noon = ThermalEnvironment {
        local_time_hours: ENCELADUS_DAY_HOURS * 0.5,
        latitude_deg,
        shadowed: false,
    }
    .surface_temperature_k();
    let midnight = ThermalEnvironment {
        local_time_hours: 0.0,
        latitude_deg,
        shadowed: false,
    }
    .surface_temperature_k();
    noon - midnight
}