pub fn ice_tectonic_rate_mm_myr(tidal_stress_mpa: f64) -> f64 {
0.8 + 2.5 * tidal_stress_mpa.max(0.0)
}
pub fn sublimation_rate_mm_yr(temperature_k: f64) -> f64 {
(temperature_k - 50.0).max(0.0) * 0.01
}
pub fn space_weathering_index(exposure_myr: f64, solar_wind_factor: f64) -> f64 {
(1.0 - (-exposure_myr / 60.0).exp()) * solar_wind_factor.max(0.0)
}