#[derive(Debug, Clone, Copy, PartialEq)]
pub struct ElevationSample {
pub latitude_deg: f64,
pub longitude_deg: f64,
pub elevation_m: f64,
}
pub fn sample_dem(latitude_deg: f64, longitude_deg: f64) -> ElevationSample {
let basin_bias = if latitude_deg < -70.0 { -4_500.0 } else { 0.0 };
let maria_bias = if longitude_deg.abs() < 60.0 {
-900.0
} else {
700.0
};
let wave = latitude_deg.to_radians().sin() * 650.0 + longitude_deg.to_radians().cos() * 420.0;
ElevationSample {
latitude_deg,
longitude_deg,
elevation_m: basin_bias + maria_bias + wave,
}
}