fn main() {
    let species = titanss::exosphere::species::reference_composition();
    for s in &species {
        println!(
            "{}: {:.0}%, {} km",
            s.name,
            s.abundance_fraction * 100.0,
            s.scale_height_km
        );
    }
    let dom = titanss::exosphere::species::dominant_species(&species).unwrap();
    println!("Dominant: {}", dom.name);
    let day = titanss::exosphere::charging::SurfaceCharging::dayside();
    println!("Dayside potential: {:.1} V", day.surface_potential_v());
    println!(
        "Jeans parameter: {:.2}",
        titanss::exosphere::escape::jeans_parameter(100.0, 0.032)
    );
}
