#[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 stickney_bias =
        if (-5.0..=10.0).contains(&latitude_deg) && (40.0..=60.0).contains(&longitude_deg) {
            -3_500.0
        } else {
            0.0
        };
    let limtoc_bias =
        if (15.0..=30.0).contains(&latitude_deg) && (150.0..=180.0).contains(&longitude_deg) {
            -1_200.0
        } else {
            0.0
        };
    let wave = latitude_deg.to_radians().sin() * 300.0 + longitude_deg.to_radians().cos() * 200.0;
    ElevationSample {
        latitude_deg,
        longitude_deg,
        elevation_m: stickney_bias + limtoc_bias + wave,
    }
}
