diff --git a/input/default.cfg b/input/default.cfg index fe0150dc..205f8941 100644 --- a/input/default.cfg +++ b/input/default.cfg @@ -10,7 +10,6 @@ star_luminosity_modern = 1 # L_sun, luminosity of star (TOD star_age_modern = 4.567e9 # yr, age of star (TODAY) star_rot_pctle = 50.0 star_spectrum = stellar_spectra/Named/sun.txt # Stellar spectrum file at 1 AU (TODAY) -star_omega = 1.0 # Star rotation rate [Omega_sun, rad s-1] ##### Planet configuration mean_distance = 1 # AU, star-planet distance @@ -22,8 +21,6 @@ asf_scalefactor = 0.375 # absorbed solar flux scale factor = 3/8 albedo_s = 0.1 # surface albedo albedo_pl = 0.1 # Bond albedo (scattering) -eccentricity = 0.017 # Earth eccentricity [dimensionless] - P_top = 1.0e-5 # bar, Pressure at TOA ##### PROTEUS settings @@ -78,12 +75,9 @@ emit_stop = 1 # Enable this break condition F_crit = 0.2 # Critical flux, below which the model will terminate # Atmospheric escape -escape_model = 1 # Escape model to be used, 0: None | 1: ZEPHYRUS | 2: Dummy +escape_model = 0 # Escape model to be used, 0: None | 1: ZEPHYRUS | 2: Dummy escape_stop = 3e-4 # Terminate when atm mass drops below this fraction of its initial mass escape_dummy_rate = 0.0 # Bulk escape rate for dummy escape model [kg s-1] -escape_el_rate = 0.15 # Efficiency factor for energy-limited escape (varies typically 0.1 < epsilon < 0.6) [dimensionless] -escape_el_tidal_correction = 0 # Tidal correction factor, 0: None | 1: yes [dimensionless] - # Method for solving for T(p) profile atmosphere_model = 0 # Atmosphere model to be used, 0: JANUS | 1: AGNI @@ -155,13 +149,13 @@ H2O_included = 1 H2O_initial_bar = 30.0 CO2_included = 1 -CO2_initial_bar = 30.0 +CO2_initial_bar = 0.0 N2_included = 1 -N2_initial_bar = 30.0 +N2_initial_bar = 0.0 S2_included = 1 -S2_initial_bar = 30.0 +S2_initial_bar = 0.0 SO2_included = 0 SO2_initial_bar = 0.0 diff --git a/input/escape.cfg b/input/escape.cfg new file mode 100644 index 00000000..fe0150dc --- /dev/null +++ b/input/escape.cfg @@ -0,0 +1,177 @@ +########## PROTEUS configuration file + +##### Star configuration + +star_model = 0 # Evolution model to use for star (0: Spada, 1: Baraffe) +star_mass = 1 # M_sun, mass of star +star_radius_modern = 1 # R_sun, radius of star (TODAY) +star_temperature_modern = 5772.0 # K, temperature of star (TODAY) +star_luminosity_modern = 1 # L_sun, luminosity of star (TODAY) used for Baraffe model +star_age_modern = 4.567e9 # yr, age of star (TODAY) +star_rot_pctle = 50.0 +star_spectrum = stellar_spectra/Named/sun.txt # Stellar spectrum file at 1 AU (TODAY) +star_omega = 1.0 # Star rotation rate [Omega_sun, rad s-1] + +##### Planet configuration +mean_distance = 1 # AU, star-planet distance +mass = 1 # M_earth +radius = 1 # R_earth + +zenith_angle = 48.19 # Cronin+14 Earth-like case with clear sky [arccos(2/3)] +asf_scalefactor = 0.375 # absorbed solar flux scale factor = 3/8 +albedo_s = 0.1 # surface albedo +albedo_pl = 0.1 # Bond albedo (scattering) + +eccentricity = 0.017 # Earth eccentricity [dimensionless] + +P_top = 1.0e-5 # bar, Pressure at TOA + +##### PROTEUS settings + +iter_max = 9000 +log_level = INFO + +# Output subdirectory name (relative to output folder) +dir_output = earth_demo + +# Choose times +time_star = 100.0e6 # yr, time since star formation +time_target = 4.567e+9 # yr, target time for MO evolution + +# SOCRATES spectral file to use (relative to FWL_DATA folder) +spectral_file = spectral_files/Frostflow/256/Frostflow.sf + +# Stellar heating toggle, 0: disabled | 1: enabled +stellar_heating = 1 + +# Time-stepping stuff +plot_iterfreq = 1 # Plotting frequency, 0: wait until completion | n: every n iterations +plot_format = pdf # Plotting image file format +sspec_dt_update = 1e9 # Time intervals at which to re-calculate the stellar spectrum +sinst_dt_update = 1e2 # Time intervals at which to re-calculate the instellation +dt_maximum = 3e7 # Maximum time-step +dt_minimum = 3e2 # Minimum time-step +dt_method = 1 # Time-stepping method, 0: proportional | 1: adaptive | 2: maximum +dt_propconst = 52.0 # Proportionality constant for dt_method=0 +dt_atol = 0.02 # Step size atol +dt_rtol = 0.07 # Step size rtol +dt_initial = 1e4 # Inital step size + +# Flux convergence scheme and tolerances for surface equilibration +shallow_ocean_layer = 0 # 0: off | 1: on +F_atm_bc = 0 # Boundary condition choice for F_atm, 0: TOA | 1: Surface +skin_d = 0.01 # m +skin_k = 2.0 # W m-1 K-1 +prevent_warming = 1 # Require that the planet only cool down over time, 0: disabled | 1: enabled + +# Break at solidification? +solid_stop = 1 # Enable this break condition +phi_crit = 0.005 # Model terminates with Phi_global < phi_crit + +# Break at steady state? +steady_stop = 1 # Enable this break condition +steady_flux = 0.8 # Maximum absolute value of F_atm allowed for convergence +steady_dprel = 1.0e-9 # Percentage change in melt fraction over time (dp/p)/dt*100 + +# Break at small flux? +emit_stop = 1 # Enable this break condition +F_crit = 0.2 # Critical flux, below which the model will terminate + +# Atmospheric escape +escape_model = 1 # Escape model to be used, 0: None | 1: ZEPHYRUS | 2: Dummy +escape_stop = 3e-4 # Terminate when atm mass drops below this fraction of its initial mass +escape_dummy_rate = 0.0 # Bulk escape rate for dummy escape model [kg s-1] +escape_el_rate = 0.15 # Efficiency factor for energy-limited escape (varies typically 0.1 < epsilon < 0.6) [dimensionless] +escape_el_tidal_correction = 0 # Tidal correction factor, 0: None | 1: yes [dimensionless] + + +# Method for solving for T(p) profile +atmosphere_model = 0 # Atmosphere model to be used, 0: JANUS | 1: AGNI +atmosphere_surf_state = 2 # Atmosphere bottom edge boundary condition, 0: free | 1: fixed at T_surf | 2: conductive skin + +# Number of levels +atmosphere_nlev = 250 + +# Temperature limits throughout atmosphere [K] +min_temperature = 0.5 +max_temperature = 5000.0 + +# Tropopause type, 0: none | 1: skin temperature | 2: dynamic +tropopause = 1 + +# Clouds +water_cloud = 0 # enable water cloud radiative effects? (1: yes, 0: no) +alpha_cloud = 0.0 # condensate retention fraction (1 -> fully retained) + +# Rayleigh scattering, 0: Disabled | 1: Enabled +rayleigh = 1 + +# Atmospheric chemistry, 0: Disabled | 1: Equilibrium | 2: Kinetics +atmosphere_chemistry = 0 + +##### SPIDER-specific settings + +interior_nlev = 220 +grain_size = 0.1 + +# Mixing length parameterization: 1: variable | 2: constant +mixing_length = 2 + +# Standard tolerance for solvers +solver_tolerance = 1.0e-10 + +# Maximum absolute surface temperature change [K] +tsurf_poststep_change = 20.0 + +# Maximum fractional surface temperature change [fraction] +tsurf_poststep_change_frac = 0.01 + +# Fractional core radius +planet_coresize = 0.55 + +# Initial entropy conditions for interior +ic_adiabat_entropy = 2600.0 +ic_dsdr = -4.698e-06 + +# Atmosphere heat flux start guess, adjusted during runtime [W/m^2] +F_atm = 8.0E4 + +# Oxygen fugacity offset relative to the IW buffer (log10 units) +fO2_shift_IW = 4 + +# Enable solving for initial partial pressures (0: off | 1: on) +solvevol_use_params = 0 + +# Parameters used to solve for initial partial pressures (when solvepp_use_params = 1) +Phi_global = 1.0 # Mantle melt fraction initial guess +CH_ratio = 1.0 # C/H ratio +hydrogen_earth_oceans = 6.0 # Hydrogen inventory in units of equivalent Earth oceans +nitrogen_ppmw = 2.0 # Nitrogen inventory in ppmw relative to mass of melt +sulfur_ppmw = 200.0 # Sulfur inventory in ppmw relative to mass of melt + +# Prescribed injected partial pressures [bar] +# Summed with solvepp results when solvepp_enabled = 1 +H2O_included = 1 +H2O_initial_bar = 30.0 + +CO2_included = 1 +CO2_initial_bar = 30.0 + +N2_included = 1 +N2_initial_bar = 30.0 + +S2_included = 1 +S2_initial_bar = 30.0 + +SO2_included = 0 +SO2_initial_bar = 0.0 + +H2_included = 0 +H2_initial_bar = 0.0 + +CH4_included = 0 +CH4_initial_bar = 0.0 + +CO_included = 0 +CO_initial_bar = 0.0 +