void generate_stellar_system()
{
planet_pointer planet;
radians_per_rotation = 2.0 * PI;
stellar_mass_ratio = random_number(0.6,1.3);
stellar_luminosity_ratio = luminosity(stellar_mass_ratio);
planet = distribute_planetary_masses(stellar_mass_ratio,stellar_luminosity_ratio,0.0,stellar_dust_limit(stellar_mass_ratio));
main_seq_life = 1.0E10 * (stellar_mass_ratio / stellar_luminosity_ratio);
if ((main_seq_life >= 6.0E9))
age = random_number(1.0E9,6.0E9);
else
age = random_number(1.0E9,main_seq_life);
r_ecosphere = sqrt(stellar_luminosity_ratio);
r_greenhouse = r_ecosphere * GREENHOUSE_EFFECT_CONST;
while (planet != NULL)
{
planet->orbit_zone = orbital_zone(planet->a);
if (planet->gas_giant)
{
planet->density = empirical_density(planet->mass,planet->a,planet->gas_giant);
planet->radius = volume_radius(planet->mass,planet->density);
}
else
{
planet->radius = kothari_radius(planet->mass,planet->a,planet->gas_giant,planet->orbit_zone);
planet->density = volume_density(planet->mass,planet->radius);
}
planet->orbital_period = period(planet->a,planet->mass,stellar_mass_ratio);
planet->day = day_length(planet->mass,planet->radius,planet->orbital_period,planet->e,planet->gas_giant);
planet->resonant_period = spin_resonance;
planet->axial_tilt = inclination(planet->a);
planet->escape_velocity = escape_vel(planet->mass,planet->radius);
planet->surface_accel = acceleration(planet->mass,planet->radius);
planet->rms_velocity = rms_vel(MOLECULAR_NITROGEN,planet->a);
planet->molecule_weight = molecule_limit(planet->a,planet->mass,planet->radius);
if ((planet->gas_giant))
{
planet->surface_grav = INCREDIBLY_LARGE_NUMBER;
planet->greenhouse_effect = FALSE;
planet->volatile_gas_inventory = INCREDIBLY_LARGE_NUMBER;
planet->surface_pressure = INCREDIBLY_LARGE_NUMBER;
planet->boil_point = INCREDIBLY_LARGE_NUMBER;
planet->hydrosphere = INCREDIBLY_LARGE_NUMBER;
planet->albedo = about(GAS_GIANT_ALBEDO,0.1);
planet->surface_temp = INCREDIBLY_LARGE_NUMBER;
}
else
{
planet->surface_grav = gravity(planet->surface_accel);
planet->greenhouse_effect = greenhouse(planet->orbit_zone,planet->a,r_greenhouse);
planet->volatile_gas_inventory = vol_inventory(planet->mass,planet->escape_velocity,planet->rms_velocity,stellar_mass_ratio,planet->orbit_zone,planet->greenhouse_effect);
planet->surface_pressure = pressure(planet->volatile_gas_inventory,planet->radius,planet->surface_grav);
if ((planet->surface_pressure == 0.0))
planet->boil_point = 0.0;
else
planet->boil_point = boiling_point(planet->surface_pressure);
iterate_surface_temp(&(planet));
}
planet = planet->next_planet;
}
display_system( );
}
stellar_system* generate_stellar_system(unsigned long random_seed)
{
planet* planet;
double outer_dust_limit;
stellar_system *system = malloc(sizeof(stellar_system));
system->first_planet = NULL;
setup_seed(system, random_seed);
system->star_mass_r = random_number(0.6, 1.3); /* was 0.6, 1.3 */
system->star_radius_r = about(pow(system->star_mass_r, 1.0 / 3.0), 0.05);
/* for some unknown reason, only 3 digits wanted... */
system->star_radius_r = floor(system->star_radius_r * 1000.0) / 1000.0;
system->star_lum_r = luminosity(system->star_mass_r);
/* luminosity is proportional to T^4 and to area of star */
/* so temp is Tsol * 4th-root ( Lum / r^2 ) */
system->star_temp = 5650 * sqrt(sqrt(system->star_lum_r) / system->star_radius_r);
/* ignore fractional degrees */
system->star_temp = floor(system->star_temp);
sprintf(system->star_class, "%.16s", find_star_class(system->star_temp));
outer_dust_limit = stellar_dust_limit(system->star_mass_r);
system->first_planet = distribute_planetary_masses(system, 0.0, outer_dust_limit);
system->main_seq_life = 1.0E10 * (system->star_mass_r / system->star_lum_r);
if (system->main_seq_life > 6.0E9)
system->star_age = random_number(1.0E9, 6.0E9);
else if (system->main_seq_life > 1.0E9)
system->star_age = random_number(1.0E9, system->main_seq_life);
else
system->star_age = random_number(system->main_seq_life/10, system->main_seq_life);
system->r_ecosphere = sqrt(system->star_lum_r);
system->r_greenhouse = system->r_ecosphere * GREENHOUSE_EFFECT_CONST;
for (planet = system->first_planet; planet != NULL; planet = planet->next_planet)
{
planet->orbit_zone = orbital_zone(system, planet->a);
if (planet->gas_giant)
{
planet->density = empirical_density(system, planet->mass, planet->a,
planet->gas_giant);
planet->radius = volume_radius(planet->mass, planet->density);
}
else
{
planet->radius = kothari_radius(planet->mass, planet->gas_giant,
planet->orbit_zone);
planet->density = volume_density(planet->mass, planet->radius);
}
planet->orb_period = period(planet->a, planet->mass, system->star_mass_r);
planet->day = day_length(system, planet->mass, planet->radius, planet->e,
planet->density, planet->a,
planet->orb_period, planet->gas_giant,
system->star_mass_r);
planet->resonant_period = system->resonance;
planet->axial_tilt = inclination(planet->a);
planet->esc_velocity = escape_velocity(planet->mass, planet->radius);
planet->surf_accel = acceleration(planet->mass, planet->radius);
planet->rms_velocity = rms_velocity(system, MOL_NITROGEN, planet->a);
planet->molec_weight = molecule_limit(planet->mass, planet->radius);
if ((planet->gas_giant))
{
planet->surf_grav = INCREDIBLY_LARGE_NUMBER;
planet->greenhouse_effect = false;
planet->volatile_gas_inventory = INCREDIBLY_LARGE_NUMBER;
planet->surf_pressure = INCREDIBLY_LARGE_NUMBER;
planet->boil_point = INCREDIBLY_LARGE_NUMBER;
planet->hydrosphere = INCREDIBLY_LARGE_NUMBER;
planet->albedo = about(GAS_GIANT_ALBEDO, 0.1);
planet->surf_temp = INCREDIBLY_LARGE_NUMBER;
}
else
{
planet->surf_grav = gravity(planet->surf_accel);
planet->greenhouse_effect = greenhouse(planet->orbit_zone, planet->a,
system->r_greenhouse);
planet->volatile_gas_inventory = vol_inventory(planet->mass,
planet->esc_velocity,
planet->rms_velocity,
system->star_mass_r,
planet->orbit_zone,
planet->greenhouse_effect);
planet->surf_pressure = pressure(planet->volatile_gas_inventory,
planet->radius, planet->surf_grav);
if (planet->surf_pressure == 0.0)
planet->boil_point = 0.0;
else
planet->boil_point = boiling_point(planet->surf_pressure);
iterate_surface_temp(system, &(planet));
}
#ifdef MOON
if (args.make_moon)
{
#ifdef PROPER_MOON
planet->first_moon = dist_moon_masses(planet->mass,
star_lum_r, planet->e,
0.0, planet_dust_limit(planet->mass));
#else
planet->first_moon = do_dist_moon_masses(planet->mass, planet->radius);
{
planet* moon = planet->first_moon;
while (moon)
{
moon->radius = kothari_radius(moon->mass, 0, planet->orbit_zone);
moon->density = volume_density(moon->mass, moon->radius);
moon->density = random_number(1.5, moon->density * 1.1);
if (moon->density < 1.5)
moon->density = 1.5;
moon->radius = volume_radius(moon->mass, moon->density);
moon->orb_period = period(moon->a, moon->mass, planet->mass);
moon->day = day_length(system, moon->mass, moon->radius, moon->e,
moon->density, moon->a,
moon->orb_period, moon->gas_giant,
planet->mass);
moon->resonant_period = system->resonance;
moon->axial_tilt = inclination(moon->a);
moon->esc_velocity = escape_vel(moon->mass, moon->radius);
moon->surf_accel = acceleration(moon->mass, moon->radius);
moon->rms_velocity = rms_vel(system, MOL_NITROGEN, planet->a);
moon->molec_weight = molecule_limit(moon->mass, moon->radius);
moon->surf_grav = gravity(moon->surf_accel);
moon->greenhouse_effect = grnhouse(planet->orbit_zone,
planet->a,
system->r_greenhouse);
moon->volatile_gas_inventory = vol_inventory(moon->mass,
moon->esc_velocity,
moon->rms_velocity,
system->star_mass_r,
planet->orbit_zone,
moon->greenhouse_effect);
moon->surf_pressure = pressure(moon->volatile_gas_inventory,
moon->radius, moon->surf_grav);
if ((moon->surf_pressure == 0.0))
moon->boil_point = 0.0;
else
moon->boil_point = boiling_point(moon->surf_pressure);
iterate_surface_temp_moon(system, &planet, &moon);
moon = moon->next_planet;
}
}
#endif /* CC_MOON */
}
#endif /* MOON */
}
return system;
}
