fn main() {
let species = enceladuss::exosphere::species::reference_composition();
for s in &species {
println!(
"{}: {:.1}% abundance, scale height {} km",
s.name,
s.abundance_fraction * 100.0,
s.scale_height_km
);
}
let dominant = enceladuss::exosphere::species::dominant_species(&species).unwrap();
println!("Dominant species: {}", dominant.name);
let day = enceladuss::exosphere::charging::SurfaceCharging::dayside();
let night = enceladuss::exosphere::charging::SurfaceCharging::nightside();
println!("Dayside potential: {:.1} V", day.surface_potential_v());
println!("Nightside potential: {:.1} V", night.surface_potential_v());
let jeans = enceladuss::exosphere::escape::jeans_parameter(80.0, 0.018);
println!("Jeans parameter (H2O, 80K): {:.2}", jeans);
}