class TimeSettings:
load_qr = 15
load_pg = 6
next_message = about(10)
next_contact = about(120)
start_time = 7 #starting at 07:00
stop_time = 0 #breaking at 00:00
def vol_inventory(mass, escape_vel, rms_vel, stellar_mass, zone,
greenhouse_effect, accreted_gas):
'''This implements Fogg's eq.17. The 'inventory' returned is unitless.'''
velocity_ratio = escape_vel / rms_vel
if velocity_ratio >= GAS_RETENTION_THRESHOLD:
if zone == Zone.ZONE_1:
proportion_ = 140000.0
'''100 . 140 JLB'''
elif zone == Zone.ZONE_2:
proportion_ = 75000.0
elif zone == Zone.ZONE_3:
proportion_ = 250.0
else:
raise NotImplementedError("orbital zone not initialized correctly")
earth_units = mass * SUN_MASS_IN_EARTH_MASSES
temp1 = (proportion_ * earth_units) / stellar_mass
temp2 = about(temp1, 0.2)
temp2 = temp1
if greenhouse_effect or accreted_gas:
return temp2
else:
return temp2 / 140.0 # 100 . 140 JLB
else:
return 0.0
def log_about_data():
logger.info('train:\n' + util.about(train, SINGLE_LINE=True))
logger.info('valid:\n' + util.about(valid, SINGLE_LINE=True))
logger.info('test:\n' + util.about(test, SINGLE_LINE=True))
logger.info('y_train:\n' + util.about(y_train, SINGLE_LINE=True))
logger.info('y_valid:\n' + util.about(y_valid, SINGLE_LINE=True))
logger.info('y_test:\n' + util.about(y_test, SINGLE_LINE=True))
def generate_planet(protoplanet, star, random_tilt=0, planet_id=None, do_gases=True, do_moons=True, is_moon=False):
planet = Planet(
sun=star,
orbit=protoplanet.orbit,
dust_mass=protoplanet.dust_mass,
gas_mass=protoplanet.gas_mass,
mass=protoplanet.mass,
axial_tilt=inclination(protoplanet.orbit.a) if random_tilt else 0,
atmosphere=None,
surf_temp=0,
high_temp=0,
low_temp=0,
max_temp=0,
min_temp=0,
greenhs_rise=0,
resonant_period=False,
orbit_zone=orb_zone(star.luminosity_ratio, protoplanet.orbit.a),
orb_period=period(protoplanet.orbit.a, protoplanet.mass, star.mass_ratio)
)
planet.exospheric_temp = EARTH_EXOSPHERE_TEMP / \
((planet.orbit.a / star.r_ecosphere) ** 2)
planet.rms_velocity = rms_vel(MOL_NITROGEN, planet.exospheric_temp)
planet.core_radius = kothari_radius(
planet.dust_mass, False, planet.orbit_zone)
# Calculate the radius as a gas giant, to verify it will retain gas.
# Then if mass > Earth, it's at least 5% gas and retains He, it's
# some flavor of gas giant.
planet.density = empirical_density(
planet.mass, planet.orbit.a, star.r_ecosphere, True)
planet.radius = volume_radius(planet.mass, planet.density)
planet.surf_accel = acceleration(planet.mass, planet.radius)
planet.surf_grav = gravity(planet.surf_accel)
planet.molec_weight = min_molec_weight(planet)
if (((planet.mass * SUN_MASS_IN_EARTH_MASSES) > 1.0)
and ((planet.gas_mass / planet.mass) > 0.05)
and (min_molec_weight(planet) <= 4.0)):
if ((planet.gas_mass / planet.mass) < 0.20):
planet.type = PlanetType.SUB_SUB_GAS_GIANT
elif ((planet.mass * SUN_MASS_IN_EARTH_MASSES) < 20.0):
planet.type = PlanetType.SUB_GAS_GIANT
else:
planet.type = PlanetType.GAS_GIANT
else: # If not, it's rocky.
planet.radius = kothari_radius(planet.mass, False, planet.orbit_zone)
planet.density = volume_density(planet.mass, planet.radius)
planet.surf_accel = acceleration(planet.mass, planet.radius)
planet.surf_grav = gravity(planet.surf_accel)
if ((planet.gas_mass / planet.mass) > 0.000001):
h2_mass = planet.gas_mass * 0.85
he_mass = (planet.gas_mass - h2_mass) * 0.999
h2_loss = 0.0
he_loss = 0.0
h2_life = gas_life(MOL_HYDROGEN, planet)
he_life = gas_life(HELIUM, planet)
if (h2_life < star.age):
#math.exp with a value above 709 results in a math range error
#this is a dumb fix. STH 2021-0131
if (star.age / h2_life)>709:
h2_loss = ((1.0 - (1.0 / exp(709.0))) * h2_mass)
else:
h2_loss = ((1.0 - (1.0 / exp(star.age / h2_life))) * h2_mass)
planet.gas_mass -= h2_loss
planet.mass -= h2_loss
planet.surf_accel = acceleration(planet.mass, planet.radius)
planet.surf_grav = gravity(planet.surf_accel)
if (he_life < star.age):
he_loss = ((1.0 - (1.0 / exp(star.age / he_life))) * he_mass)
planet.gas_mass -= he_loss
planet.mass -= he_loss
planet.surf_accel = acceleration(planet.mass, planet.radius)
planet.surf_grav = gravity(planet.surf_accel)
'''if (((h2_loss + he_loss) > .000001) and (flag_verbose & 0x0080)):
fprintf(stderr, "%s\tLosing gas: H2: %5.3Lf EM, He: %5.3Lf EM\n",
planet_id,
h2_loss * SUN_MASS_IN_EARTH_MASSES, he_loss * SUN_MASS_IN_EARTH_MASSES)'''
planet.day = day_length(planet) # Modifies planet.resonant_period
planet.esc_velocity = escape_vel(planet.mass, planet.radius)
if planet.type == PlanetType.GAS_GIANT or planet.type == PlanetType.SUB_GAS_GIANT or planet.type == PlanetType.SUB_SUB_GAS_GIANT:
planet.greenhouse_effect = False
planet.volatile_gas_inventory = INCREDIBLY_LARGE_NUMBER
planet.surf_pressure = INCREDIBLY_LARGE_NUMBER
planet.boil_point = INCREDIBLY_LARGE_NUMBER
planet.surf_temp = INCREDIBLY_LARGE_NUMBER
planet.greenhs_rise = 0
planet.albedo = about(GAS_GIANT_ALBEDO, 0.1)
planet.hydrosphere = 1.0
planet.cloud_cover = 1.0
planet.ice_cover = 0.0
planet.surf_grav = gravity(planet.surf_accel)
planet.molec_weight = min_molec_weight(planet)
planet.surf_grav = INCREDIBLY_LARGE_NUMBER
planet.estimated_temp = est_temp(
star.r_ecosphere, planet.orbit.a, planet.albedo)
planet.estimated_terr_temp = est_temp(
star.r_ecosphere, planet.orbit.a, EARTH_ALBEDO)
temp = planet.estimated_terr_temp
if (temp >= FREEZING_POINT_OF_WATER) and (temp <= EARTH_AVERAGE_KELVIN + 10.) and (star.age > 2.0E9):
pass
'''if (flag_verbose & 0x8000):
fprintf (stderr, "%s\t%s (%4.2LfEM %5.3Lf By)%s with earth-like temperature (%.1Lf C, %.1Lf F, %+.1Lf C Earth).\n",
planet_id,
planet.type == PlanetType.GAS_GIANT ? "Jovian" :
planet.type == PlanetType.SUB_GAS_GIANT ? "Sub-Jovian" :
planet.type == PlanetType.SUB_SUB_GAS_GIANT ? "Gas Dwarf" :
"Big",
planet.mass * SUN_MASS_IN_EARTH_MASSES,
star.age /1.0E9,
planet.first_moon == NULL ? "" : " WITH MOON",
temp - FREEZING_POINT_OF_WATER,
32 + ((temp - FREEZING_POINT_OF_WATER) * 1.8),
temp - EARTH_AVERAGE_KELVIN)'''
else:
planet.estimated_temp = est_temp(
star.r_ecosphere, planet.orbit.a, EARTH_ALBEDO)
planet.estimated_terr_temp = est_temp(
star.r_ecosphere, planet.orbit.a, EARTH_ALBEDO)
planet.surf_grav = gravity(planet.surf_accel)
planet.molec_weight = min_molec_weight(planet)
planet.greenhouse_effect = grnhouse(star.r_ecosphere, planet.orbit.a)
planet.volatile_gas_inventory = vol_inventory(planet.mass,
planet.esc_velocity,
planet.rms_velocity,
star.mass_ratio,
planet.orbit_zone,
planet.greenhouse_effect,
(planet.gas_mass
/ planet.mass) > 0.000001)
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)
# Sets:
# planet.surf_temp
# planet.greenhs_rise
# planet.albedo
# planet.hydrosphere
# planet.cloud_cover
# planet.ice_cover
iterate_surface_temp(planet)
if (do_gases and (planet.max_temp >= FREEZING_POINT_OF_WATER) and (planet.min_temp <= planet.boil_point)):
calculate_gases(star, planet, planet_id)
# Next we assign a type to the planet.
if (planet.surf_pressure < 1.0):
if (not is_moon) and ((planet.mass * SUN_MASS_IN_EARTH_MASSES) < ASTEROID_MASS_LIMIT):
planet.type = PlanetType.ASTEROIDS
else:
planet.type = PlanetType.ROCK
elif (planet.surf_pressure > 6000.0) and (planet.molec_weight <= 2.0): # Retains Hydrogen
planet.type = PlanetType.SUB_SUB_GAS_GIANT
planet.gases = 0
planet.atmosphere = None
else:
# Atmospheres:
if (int(planet.day) == int(planet.orb_period * 24.0)) or planet.resonant_period:
planet.type = PlanetType.ONE_FACE
elif (planet.hydrosphere >= 0.95):
planet.type = PlanetType.WATER # >95% water
elif (planet.ice_cover >= 0.95):
planet.type = PlanetType.ICE # >95% ice
elif (planet.hydrosphere > 0.05):
planet.type = PlanetType.TERRESTRIAL # Terrestrial
# else <5% water
elif (planet.max_temp > planet.boil_point):
planet.type = PlanetType.VENUSIAN # Hot = Venusian
elif ((planet.gas_mass / planet.mass) > 0.0001):
# Accreted gas
planet.type = PlanetType.ICE # But no Greenhouse
planet.ice_cover = 1.0 # or liquid water
# Make it an Ice World
elif (planet.surf_pressure <= 250.0): # Thin air = Martian
planet.type = PlanetType.MARTIAN
elif (planet.surf_temp < FREEZING_POINT_OF_WATER):
planet.type = PlanetType.ICE
else:
planet.type = PlanetType.UNKNOWN # TODO(woursler): Consider throwing an error here.
'''if (flag_verbose & 0x0001)
fprintf (stderr, "%12s\tp=%4.2Lf\tm=%4.2Lf\tg=%4.2Lf\tt=%+.1Lf\t%s\t Unknown %s\n",
type_string (planet.type),
planet.surf_pressure,
planet.mass * SUN_MASS_IN_EARTH_MASSES,
planet.surf_grav,
planet.surf_temp - EARTH_AVERAGE_KELVIN,
planet_id,
((int)planet.day == (int)(planet.orb_period * 24.0) or
(planet.resonant_period)) ? "(1-Face)" : ""
)'''
if do_moons and not is_moon:
for protomoon in protoplanet.moons:
if protomoon.mass * SUN_MASS_IN_EARTH_MASSES > .000001:
protomoon.orbit = planet.orbit
# Note: adjusts density, which is used in computing the roche limit.
moon = generate_planet(
protoplanet=protomoon,
star=star,
random_tilt=random_tilt,
do_gases=do_gases,
do_moons=do_moons,
is_moon=True
)
# TODO(woursler): these should be their own subroutines.
roche_limit_r = roche_limit(planet, moon)
hill_sphere_r = hill_sphere(planet, star)
if (roche_limit_r * 3.0) < hill_sphere_r:
moon_a = random_number(
roche_limit_r * 1.5, hill_sphere_r / 2.0) / KM_PER_AU
moon_e = random_eccentricity()
moon.orbit = Orbit(a=moon_a, e=moon_e)
else:
moon.orbit = Orbit(a=0, e=0)
planet.moons.append(moon)
'''
if (flag_verbose & 0x40000):
fprintf (stderr,
" Roche limit: R = %4.2Lg, rM = %4.2Lg, rm = %4.2Lg . %.0Lf km\n"
" Hill Sphere: a = %4.2Lg, m = %4.2Lg, M = %4.2Lg . %.0Lf km\n"
"%s Moon orbit: a = %.0Lf km, e = %.0Lg\n",
planet.radius, planet.density, ptr.density,
roche_limit,
planet.orbit.a * KM_PER_AU, planet.mass * SOLAR_MASS_IN_KILOGRAMS, star.mass_ratio * SOLAR_MASS_IN_KILOGRAMS,
hill_sphere,
moon_id,
ptr.moon_a * KM_PER_AU, ptr.moon_e
)
if (flag_verbose & 0x1000):
fprintf (stderr, " %s: (%7.2LfEM) %d %4.2LgEM\n",
planet_id,
planet.mass * SUN_MASS_IN_EARTH_MASSES,
n,
ptr.mass * SUN_MASS_IN_EARTH_MASSES)'''
return planet
def inclination(orb_radius):
'''The orbital radius is expected in units of Astronomical Units (AU).
Inclination is returned in units of degrees. '''
temp = int((orb_radius**0.2) * about(EARTH_AXIAL_TILT, 0.4))
return temp % 360
def generate_planet(protoplanet,
star,
random_tilt=0,
planet_id=None,
do_gases=True,
do_moons=True,
is_moon=False):
planet = Planet(
sun=star,
orbit=protoplanet.orbit,
dust_mass=protoplanet.dust_mass,
gas_mass=protoplanet.gas_mass,
mass=protoplanet.mass,
axial_tilt=inclination(protoplanet.orbit.a) if random_tilt else 0,
atmosphere=None,
surf_temp=0,
high_temp=0,
low_temp=0,
max_temp=0,
min_temp=0,
greenhs_rise=0,
resonant_period=False,
orbit_zone=orb_zone(star.luminosity_ratio, protoplanet.orbit.a),
orb_period=period(protoplanet.orbit.a, protoplanet.mass,
star.mass_ratio),
mapname="Unknown")
planet.exospheric_temp = EARTH_EXOSPHERE_TEMP / \
((planet.orbit.a / star.r_ecosphere) ** 2)
planet.rms_velocity = rms_vel(MOL_NITROGEN, planet.exospheric_temp)
planet.core_radius = kothari_radius(planet.dust_mass, False,
planet.orbit_zone)
# Calculate the radius as a gas giant, to verify it will retain gas.
# Then if mass > Earth, it's at least 5% gas and retains He, it's
# some flavor of gas giant.
planet.density = empirical_density(planet.mass, planet.orbit.a,
star.r_ecosphere, True)
planet.radius = volume_radius(planet.mass, planet.density)
planet.surf_accel = acceleration(planet.mass, planet.radius)
planet.surf_grav = gravity(planet.surf_accel)
planet.molec_weight = min_molec_weight(planet)
if (((planet.mass * SUN_MASS_IN_EARTH_MASSES) > 1.0)
and ((planet.gas_mass / planet.mass) > 0.05)
and (min_molec_weight(planet) <= 4.0)):
if ((planet.gas_mass / planet.mass) < 0.20):
planet.type = PlanetType.SUB_SUB_GAS_GIANT
elif ((planet.mass * SUN_MASS_IN_EARTH_MASSES) < 20.0):
planet.type = PlanetType.SUB_GAS_GIANT
else:
planet.type = PlanetType.GAS_GIANT
else: # If not, it's rocky.
planet.radius = kothari_radius(planet.mass, False, planet.orbit_zone)
planet.density = volume_density(planet.mass, planet.radius)
planet.surf_accel = acceleration(planet.mass, planet.radius)
planet.surf_grav = gravity(planet.surf_accel)
if ((planet.gas_mass / planet.mass) > 0.000001):
h2_mass = planet.gas_mass * 0.85
he_mass = (planet.gas_mass - h2_mass) * 0.999
h2_loss = 0.0
he_loss = 0.0
h2_life = gas_life(MOL_HYDROGEN, planet)
he_life = gas_life(HELIUM, planet)
if (h2_life < star.age):
h2_gamma = star.age / h2_life
if star.age / h2_life > 709: #This h2_gamma business is to avoid overflows.
h2_gamma = 709
h2_loss = ((1.0 - (1.0 / exp(h2_gamma))) * h2_mass)
planet.gas_mass -= h2_loss
planet.mass -= h2_loss
planet.surf_accel = acceleration(planet.mass, planet.radius)
planet.surf_grav = gravity(planet.surf_accel)
if (he_life < star.age):
he_gamma = star.age / he_life
if star.age / he_life > 709: #This he_gamma business is to avoid overflows.
he_gamma = 709
he_loss = ((1.0 - (1.0 / exp(he_gamma))) * he_mass)
planet.gas_mass -= he_loss
planet.mass -= he_loss
planet.surf_accel = acceleration(planet.mass, planet.radius)
planet.surf_grav = gravity(planet.surf_accel)
'''if (((h2_loss + he_loss) > .000001) and (flag_verbose & 0x0080)):
fprintf(stderr, "%s\tLosing gas: H2: %5.3Lf EM, He: %5.3Lf EM\n",
planet_id,
h2_loss * SUN_MASS_IN_EARTH_MASSES, he_loss * SUN_MASS_IN_EARTH_MASSES)'''
planet.day = day_length(planet) # Modifies planet.resonant_period
planet.esc_velocity = escape_vel(planet.mass, planet.radius)
if planet.type == PlanetType.GAS_GIANT or planet.type == PlanetType.SUB_GAS_GIANT or planet.type == PlanetType.SUB_SUB_GAS_GIANT:
planet.greenhouse_effect = False
planet.volatile_gas_inventory = INCREDIBLY_LARGE_NUMBER
planet.surf_pressure = INCREDIBLY_LARGE_NUMBER
planet.boil_point = INCREDIBLY_LARGE_NUMBER
planet.surf_temp = INCREDIBLY_LARGE_NUMBER
planet.greenhs_rise = 0
planet.albedo = about(GAS_GIANT_ALBEDO, 0.1)
planet.hydrosphere = 1.0
planet.cloud_cover = 1.0
planet.ice_cover = 0.0
planet.surf_grav = gravity(planet.surf_accel)
planet.molec_weight = min_molec_weight(planet)
planet.surf_grav = INCREDIBLY_LARGE_NUMBER
planet.estimated_temp = est_temp(star.r_ecosphere, planet.orbit.a,
planet.albedo)
planet.estimated_terr_temp = est_temp(star.r_ecosphere, planet.orbit.a,
EARTH_ALBEDO)
temp = planet.estimated_terr_temp
if (temp >= FREEZING_POINT_OF_WATER) and (
temp <= EARTH_AVERAGE_KELVIN + 10.) and (star.age > 2.0E9):
pass
'''if (flag_verbose & 0x8000):
fprintf (stderr, "%s\t%s (%4.2LfEM %5.3Lf By)%s with earth-like temperature (%.1Lf C, %.1Lf F, %+.1Lf C Earth).\n",
planet_id,
planet.type == PlanetType.GAS_GIANT ? "Jovian" :
planet.type == PlanetType.SUB_GAS_GIANT ? "Sub-Jovian" :
planet.type == PlanetType.SUB_SUB_GAS_GIANT ? "Gas Dwarf" :
"Big",
planet.mass * SUN_MASS_IN_EARTH_MASSES,
star.age /1.0E9,
planet.first_moon == NULL ? "" : " WITH MOON",
temp - FREEZING_POINT_OF_WATER,
32 + ((temp - FREEZING_POINT_OF_WATER) * 1.8),
temp - EARTH_AVERAGE_KELVIN)'''
else:
planet.estimated_temp = est_temp(star.r_ecosphere, planet.orbit.a,
EARTH_ALBEDO)
planet.estimated_terr_temp = est_temp(star.r_ecosphere, planet.orbit.a,
EARTH_ALBEDO)
planet.surf_grav = gravity(planet.surf_accel)
planet.molec_weight = min_molec_weight(planet)
planet.greenhouse_effect = grnhouse(star.r_ecosphere, planet.orbit.a)
planet.volatile_gas_inventory = vol_inventory(
planet.mass, planet.esc_velocity, planet.rms_velocity,
star.mass_ratio, planet.orbit_zone, planet.greenhouse_effect,
(planet.gas_mass / planet.mass) > 0.000001)
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)
# Sets:
# planet.surf_temp
# planet.greenhs_rise
# planet.albedo
# planet.hydrosphere
# planet.cloud_cover
# planet.ice_cover
iterate_surface_temp(planet)
if (do_gases and (planet.max_temp >= FREEZING_POINT_OF_WATER)
and (planet.min_temp <= planet.boil_point)):
logging.debug("Calculating gases...")
calculate_gases(star, planet, planet_id)
else:
logging.debug("No gases for this planet.")
logging.debug("Do gases?: %(dg)s", {'dg': do_gases})
logging.debug("Max temp high enough? %(mt)s >= %(fpw)s", {
'mt': planet.max_temp,
'fpw': FREEZING_POINT_OF_WATER
})
logging.debug("Min temp less than boiling? %(mt)s <= %(bp)s", {
'mt': planet.min_temp,
'bp': planet.boil_point
})
# Next we assign a type to the planet.
if (planet.surf_pressure < 1.0):
if (not is_moon) and ((planet.mass * SUN_MASS_IN_EARTH_MASSES) <
ASTEROID_MASS_LIMIT):
planet.type = PlanetType.ASTEROIDS
else:
planet.type = PlanetType.ROCK
elif (planet.surf_pressure > 6000.0) and (planet.molec_weight <=
2.0): # Retains Hydrogen
planet.type = PlanetType.SUB_SUB_GAS_GIANT
planet.gases = 0
planet.atmosphere = None
else:
# Atmospheres:
if (int(planet.day) == int(
planet.orb_period * 24.0)) or planet.resonant_period:
planet.type = PlanetType.ONE_FACE
elif (planet.hydrosphere >= 0.95):
planet.type = PlanetType.WATER # >95% water
elif (planet.ice_cover >= 0.95):
planet.type = PlanetType.ICE # >95% ice
elif (planet.hydrosphere > 0.05):
planet.type = PlanetType.TERRESTRIAL # Terrestrial
# else <5% water
elif (planet.max_temp > planet.boil_point):
planet.type = PlanetType.VENUSIAN # Hot = Venusian
elif ((planet.gas_mass / planet.mass) > 0.0001):
# Accreted gas
planet.type = PlanetType.ICE # But no Greenhouse
planet.ice_cover = 1.0 # or liquid water
# Make it an Ice World
elif (planet.surf_pressure <= 250.0): # Thin air = Martian
planet.type = PlanetType.MARTIAN
elif (planet.surf_temp < FREEZING_POINT_OF_WATER):
planet.type = PlanetType.ICE
else:
planet.type = PlanetType.UNKNOWN # TODO(woursler): Consider throwing an error here.
#Breathability index
planet.breath = breathability(planet)
'''if (flag_verbose & 0x0001)
fprintf (stderr, "%12s\tp=%4.2Lf\tm=%4.2Lf\tg=%4.2Lf\tt=%+.1Lf\t%s\t Unknown %s\n",
type_string (planet.type),
planet.surf_pressure,
planet.mass * SUN_MASS_IN_EARTH_MASSES,
planet.surf_grav,
planet.surf_temp - EARTH_AVERAGE_KELVIN,
planet_id,
((int)planet.day == (int)(planet.orb_period * 24.0) or
(planet.resonant_period)) ? "(1-Face)" : ""
)'''
if do_moons and not is_moon:
for protomoon in protoplanet.moons:
if protomoon.mass * SUN_MASS_IN_EARTH_MASSES > .000001:
protomoon.orbit = planet.orbit
# Note: adjusts density, which is used in computing the roche limit.
moon = generate_planet(protoplanet=protomoon,
star=star,
random_tilt=random_tilt,
do_gases=do_gases,
do_moons=do_moons,
is_moon=True)
# TODO(woursler): these should be their own subroutines.
roche_limit_r = roche_limit(planet, moon)
hill_sphere_r = hill_sphere(planet, star)
if (roche_limit_r * 3.0) < hill_sphere_r:
moon_a = random_number(roche_limit_r * 1.5,
hill_sphere_r / 2.0) / KM_PER_AU
moon_e = random_eccentricity()
moon.orbit = Orbit(a=moon_a, e=moon_e)
else:
moon.orbit = Orbit(a=0, e=0)
planet.moons.append(moon)
'''
if (flag_verbose & 0x40000):
fprintf (stderr,
" Roche limit: R = %4.2Lg, rM = %4.2Lg, rm = %4.2Lg . %.0Lf km\n"
" Hill Sphere: a = %4.2Lg, m = %4.2Lg, M = %4.2Lg . %.0Lf km\n"
"%s Moon orbit: a = %.0Lf km, e = %.0Lg\n",
planet.radius, planet.density, ptr.density,
roche_limit,
planet.orbit.a * KM_PER_AU, planet.mass * SOLAR_MASS_IN_KILOGRAMS, star.mass_ratio * SOLAR_MASS_IN_KILOGRAMS,
hill_sphere,
moon_id,
ptr.moon_a * KM_PER_AU, ptr.moon_e
)
if (flag_verbose & 0x1000):
fprintf (stderr, " %s: (%7.2LfEM) %d %4.2LgEM\n",
planet_id,
planet.mass * SUN_MASS_IN_EARTH_MASSES,
n,
ptr.mass * SUN_MASS_IN_EARTH_MASSES)'''
#Traveller map information added here
map_attributes = []
#Generate values.
map_attributes.append('########-') #map_location = '' #Not yet decided.
map_attributes.append(
''
) #map_starport = '' #Not yet decided. #Placeholder character here since starport needs population to calculate.
if (planet.radius * 2 * 6.369e-05) + (planet.surf_grav *
6.32775242) + 0.82893812 > 1000:
map_attributes.append("X")
else:
map_attributes.append(
int(
round((planet.radius * 2 * 6.369e-05) +
(planet.surf_grav * 6.32775242) +
0.82893812))) #map_size
if (planet.surf_pressure) > 100000:
map_attributes.append('G')
else:
map_attributes.append(
int(
round((planet.surf_pressure * 0.000986923) * 3.49523894 +
1.74382934))
) #map_atmosphere = '' 0.000986923 added to convert mb to atm
map_attributes.append(
int(round(planet.hydrosphere * 0.1 -
0.416666666666664))) #map_hydro = ''
if planet.type == PlanetType.GAS_GIANT or planet.type == PlanetType.SUB_GAS_GIANT or planet.type == PlanetType.SUB_SUB_GAS_GIANT:
map_attributes.append('X')
else:
map_attributes.append(int(round(random.gauss(7, 2.4145392935) -
2))) #map_population = ''
if map_attributes[5] < 0:
map_attributes[5] = 0
#Starport Now
if map_attributes[5] != 'X' and map_attributes[5] > 0:
temp = int(
round(
random.gauss(7, 2.4145392935) + (map_attributes[5] * 0.5 - 3)))
if temp <= 2:
temp = 'X'
elif 2 < temp <= 5:
temp = 'E'
elif 5 < temp <= 7:
temp = 'D'
elif 7 < temp <= 9:
temp = 'C'
elif 9 < temp <= 11:
temp = 'B'
elif 11 < temp:
temp = 'A'
map_attributes[1] = temp
temp = ''
else:
map_attributes[1] = 'X'
if map_attributes[5] != 'X' and map_attributes[5] > 0:
map_attributes.append(
int(round(random.gauss(7, 2.4145392935) - 7 + map_attributes[5])))
elif map_attributes[5] == 'X':
map_attributes.append('X')
else:
map_attributes.append(0) #map_government = ''
if map_attributes[5] != 'X' and map_attributes[5] > 0:
map_attributes.append(
int(round(random.gauss(7, 2.4145392935) - 7 + map_attributes[6])))
elif map_attributes[5] == 'X':
map_attributes.append('X')
else:
map_attributes.append(0) #map_law = ''
if map_attributes[5] != 'X' and map_attributes[5] > 0:
map_attributes.append(
int(round(random.gauss(3.5, 1.7058722109) + map_attributes[5])))
elif map_attributes[5] == 'X':
map_attributes.append('X')
else:
map_attributes.append(0) #map_tech = ''
#Now convert high numbers to single letter characters.
digitbank = ['A', 'B', 'C', 'D', 'E', 'F']
for e in map_attributes:
try:
if map_attributes[e] > 9:
if map_attributes[e] > 15:
map_attributes[e] = 'G'
map_attributes[e] = map_attributes[e] - 10
map_attributes[e] = digitbank[map_attributes[e]]
if map_attributes[e] < 0:
map_attributes[e] = 0
except:
pass
planet.mapname = ''.join(str(e) for e in map_attributes)
return planet