fn main() {
    let composition = deimoss::exosphere::species::reference_composition();
    println!("Exosphere species count: {}", composition.len());

    let trace = deimoss::exosphere::species::trace_species();
    for s in &trace {
        println!(
            "Trace: {} abundance={:.2e} scale_height={:.3} km",
            s.name, s.abundance_fraction, s.scale_height_km
        );
    }

    let charging = deimoss::exosphere::charging::SurfaceCharging::dayside();
    println!("Dayside potential: {:.1} V", charging.surface_potential_v());
    println!("Dust lofting risk: {:.3}", charging.dust_lofting_risk());

    let hop = deimoss::exosphere::escape::ballistic_hop_range_m(4.0, 45.0);
    println!("Ballistic hop at 4 m/s 45 deg: {:.0} m", hop);

    let jeans = deimoss::exosphere::escape::jeans_parameter(200.0, 0.002);
    println!("Jeans parameter (H2 200K): {:.4}", jeans);
}
