ioss 0.0.3

Io 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

use crate::{IO_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 / IO_DAY_HOURS).fract();
        let solar_term = -((diurnal_phase * std::f64::consts::TAU).cos());
        let latitude_cooling = self.latitude_deg.abs() * 0.50;
        (130.0 + 50.0 * solar_term.max(-0.95) - latitude_cooling)
            .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: IO_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
}