Skip to content

Boundary interior module #668

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Draft
rdc49 wants to merge 6 commits into
base: main
Choose a base branch
from
Draft

Boundary interior module #668

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
6 commits
Select commit Hold shift + click to select a range
7305185
commit before pulling from main branch
rdc49 Jan 26, 2026
bd72d69
commit before pulling from main branch 2
rdc49 Jan 26, 2026
dd817c2
added git ignore to some .toml files
rdc49 Jan 26, 2026
5f202b4
working with dummy atm and no radio
rdc49 Apr 14, 2026
9c1c966
fixing merge conflicts
rdc49 Apr 14, 2026
e986c34
bug fixes
rdc49 Apr 14, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
71 changes: 30 additions & 41 deletions input/all_options.toml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -24,7 +24,6 @@
# [observe] synthetic observations

# ----------------------------------------------------

version = "2.0"

# Parameters
Expand Down Expand Up @@ -100,7 +99,6 @@ version = "2.0"
spin_enabled = true
offset_spin = 0 # correction to calculated Breakup period [s]


# ----------------------------------------------------
# Star
[star]
Expand Down Expand Up @@ -128,7 +126,6 @@ version = "2.0"
phoenix_radius = "none" # Stellar radius [R_sun] used for PHOENIX spectrum scaling
phoenix_log_g = "none" # Stellar surface gravity [dex]
phoenix_Teff = "none" # Stellar effective temperature [K]

[star.dummy]
radius = 1.0 # Constant stellar radius [R_sun]
calculate_radius = false # Calculate star radius using scaling from Teff?
Expand All @@ -137,11 +134,11 @@ version = "2.0"
# Orbital system
[orbit]
instellation_method = 'sma' # whether to define orbit using semi major axis ('sma') or instellation flux ('inst')
instellationflux = 1.0 # instellation flux received from the planet in [Earth units]
semimajoraxis = 1.0 # initial semi-major axis of planet's orbit [AU]
eccentricity = 0.0 # initial eccentricity of planet's orbit [dimensionless]
zenith_angle = 48.19 # characteristic zenith angle [degrees]
s0_factor = 0.375 # instellation scale factor [dimensionless]
instellationflux = 1.0 # instellation flux received from the planet in Earth units.
semimajoraxis = 1.0 # AU
eccentricity = 0.0 # dimensionless
zenith_angle = 48.19 # degrees
s0_factor = 0.375 # dimensionless

evolve = false # whether to evolve the SMaxis and eccentricity
module = "none" # module used to calculate tidal heating
Expand Down Expand Up @@ -173,8 +170,7 @@ version = "2.0"
update_dtmagma_frac = 0.03 # trigger on relative T_magma change (3%)
update_dphi_abs = 0.05 # trigger on absolute Phi_global change

[struct.zalmoxis]
# Per-layer EOS: "<source>:<material>". See Zalmoxis docs for valid options.
[struct.zalmoxis]# Per-layer EOS: "<source>:<material>". See Zalmoxis docs for valid options.
# Tabulated: Seager2007:iron, Seager2007:MgSiO3, WolfBower2018:MgSiO3, RTPress100TPa:MgSiO3, Seager2007:H2O
# Analytic: Analytic:iron, Analytic:MgSiO3, Analytic:MgFeSiO3, Analytic:H2O, Analytic:graphite, Analytic:SiC
core_eos = "Seager2007:iron" # core layer EOS
Expand Down Expand Up @@ -212,15 +208,14 @@ version = "2.0"
cloud_alpha = 0.0 # condensate retention fraction (1 -> fully retained)
surf_state = "fixed" # surface scheme: "mixed_layer" | "fixed" | "skin"
surf_greyalbedo = 0.1 # surface grey albedo
albedo_pl = 0.0 # Enforced input value of bond albedo; float (0 to 1) or str (path to CSV)
rayleigh = true # enable rayleigh scattering
albedo_pl = 0.0 # Bond albedo (scattering)
rayleigh = false # enable rayleigh scattering
tmp_minimum = 0.5 # temperature floor on solver
tmp_maximum = 5000.0 # temperature ceiling on solver

module = "agni" # Which atmosphere module to use

[atmos_clim.agni]
verbosity = 1 # output verbosity for agni (0:none, 1:info, 2:debug)
p_top = 1.0e-5 # bar, top of atmosphere grid pressure
p_obs = 0.02 # bar, level probed in transmission
spectral_group = "Honeyside" # which gas opacities to include
Expand All @@ -245,19 +240,15 @@ version = "2.0"
sens_heat = true # include sensible heat flux near surface, with TKE scheme
real_gas = false # use real-gas equations of state
psurf_thresh = 0.1 # bar, surface pressure where we switch to 'transparent' mode
dx_max_ini = 300.0 # initial maximum temperature step [kelvin] allowed by solver
dx_max = 35.0 # maximum temperature step [kelvin] allowed by solver
max_steps = 70 # max steps allowed by solver during each iteration
perturb_all = true # updated entire jacobian each step?
mlt_criterion = "s" # MLT convection stability criterion; (l)edoux or (s)chwarzschild
mlt_criterion = "l" # MLT convection stability criterion; (l)edoux or (s)chwarzschild
fastchem_floor = 150.0 # Minimum temperature allowed to be sent to FC
fastchem_maxiter_chem = 60000 # Maximum FC iterations (chemistry)
fastchem_maxiter_solv = 20000 # Maximum FC iterations (internal solver)
fastchem_xtol_chem = 1e-4 # FC solver tolerance (chemistry)
fastchem_xtol_elem = 1e-4 # FC solver tolerance (elemental)
ini_profile = 'isothermal' # Initial guess for temperature profile shape
ls_default = 2 # Default linesearch method (0:none, 1:gs, 2:bt)
fdo = 2 # finite-difference order (options: 2, 4)

[atmos_clim.janus]
p_top = 1.0e-6 # bar, top of atmosphere grid pressure
Expand All @@ -277,7 +268,7 @@ version = "2.0"
[escape]

module = "zephyrus" # Which escape module to use
reservoir = "outgas" # Reservoir that sets escaping composition when fractionation disabled
reservoir = "outgas" # Reservoir that sets escaping gas composition

[escape.zephyrus]
Pxuv = 5e-5 # Pressure at which XUV radiation become opaque in the planetary atmosphere [bar]
Expand All @@ -287,26 +278,6 @@ version = "2.0"
[escape.dummy]
rate = 0.0 # Bulk unfractionated escape rate [kg s-1]

[escape.boreas]
fractionate = true # Include fractionation in outflow?
efficiency = 0.1 # Escape efficiency factor
sigma_H = 1.89e-18 # H absorption cross-section in XUV [cm2]
sigma_O = 2.00e-18 # O absorption ^
sigma_C = 2.50e-18 # C absorption ^
sigma_N = 3.00e-18 # N absorption ^
sigma_S = 6.00e-18 # S absorption ^
kappa_H2 = 0.01 # H2 opacity in IR, grey [cm2 g-1]
kappa_H2O = 1.0 # H2O opacity ^
kappa_O2 = 1.0 # O2 opacity ^
kappa_CO2 = 1.0 # CO2 opacity ^
kappa_CO = 1.0 # CO opacity ^
kappa_CH4 = 1.0 # CH4 opacity ^
kappa_N2 = 1.0 # N2 opacity ^
kappa_NH3 = 1.0 # NH3 opacity ^
kappa_H2S = 1.0 # H2S opacity ^
kappa_SO2 = 1.0 # SO2 opacity ^
kappa_S2 = 1.0 # S2 opacity ^

# Interior - physics table
[interior]
grain_size = 0.1 # crystal settling grain size [m]
Expand All @@ -318,7 +289,6 @@ version = "2.0"
melting_dir = "Monteux-600" # Melting curve set used by all interior modules (Zalmoxis, Aragog, SPIDER)
eos_dir = "WolfBower2018_MgSiO3" # Dynamic EOS for SPIDER + Aragog (in FWL_DATA/interior_lookup_tables/EOS/dynamic/)


module = "spider" # Which interior module to use

[interior.spider]
Expand Down Expand Up @@ -366,6 +336,25 @@ version = "2.0"
mantle_cp = 1792.0 # Mantle heat capacity [J K-1 kg-1]
H_radio = 0.0 # Radiogenic heating rate [W/kg]

[interior.boundary]
rtol = 0.0001 # ODE solver relative tolerance.
atol = 1e-08 # ODE solver absolute tolerance.
T_surf_0 = 3900.0 # Initial surface temperature [K]
T_p_0 = 4000.0 # Initial mantle potential temperature [K]
T_solidus = 1420.0 # Mantle solidus temperature [K]
T_liquidus = 2020.0 # Mantle liquidus temperature [K]
critical_melt_fraction = 0.4 # Critical melt fraction for rheological transition.
critical_rayleigh_number = 1100.0 # Critical Rayleigh number for onset of convection
heat_fusion_silicate = 400000.0 # Latent heat of fusion for silicates [J/kg]
nusselt_exponent = 0.33 # Nusselt-Rayleigh scaling exponent
silicate_heat_capacity = 1200.0 # Silicate heat capacity [J/kg/K]
thermal_conductivity = 4.2 # Thermal conductivity [W/m/K]
thermal_diffusivity = 1e-06 # Thermal diffusivity [m^2/s]
thermal_expansivity = 2e-05 # Thermal expansivity [1/K]
transition_width = 0.15 # Width of viscosity transition in melt fraction space
eta_solid_const = 1e+21 # Constant solid viscosity for aggregate formulation [Pa s]
eta_melt_const = 100.0 # Constant melt viscosity for aggregate formulation [Pa s]

# Outgassing - physics table
[outgas]
fO2_shift_IW = 4 # atmosphere/interior boundary oxidation state [log10(ΔIW)]
Expand Down Expand Up @@ -461,7 +450,7 @@ version = "2.0"
ini_mix = "profile" # Initial mixing ratios (profile, outgas)
fix_surf = false # Fixed surface mixing ratios
network = "SNCHO" # Class of chemical network to use (CHO, NCHO, SNCHO)
save_frames = false # Plot frames during iterations
save_frames = true # Plot frames during iterations
yconv_cri = 0.05 # Convergence criterion, value of mixing ratios
slope_cri = 0.0001 # Convergence criterion, rate of change of mixing ratios

Expand Down
4 changes: 2 additions & 2 deletions input/chili/intercomp/_base.grid.toml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3,13 +3,13 @@ symlink = ""
ref_config = "input/chili/intercomp/_base.toml"
use_slurm = true

max_jobs = 9 # maximum number of concurrent tasks (e.g. 500 on Habrok)
max_jobs = 16 # maximum number of concurrent tasks (e.g. 500 on Habrok)
max_days = 1 # maximum number of days to run (e.g. 1)
max_mem = 3 # maximum memory per CPU in GB (e.g. 3)

["delivery.elements.H_kg"]
method = "direct"
values = [1.60e20, 7.80e20, 16e20]
values = [1.60e20, 4.70e20, 7.80e20, 16e20]

["delivery.elements.C_kg"]
method = "direct"
Expand Down
29 changes: 13 additions & 16 deletions input/chili/intercomp/_base.toml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -12,7 +12,7 @@ archive_mod = "none"
remove_sf = true

[params.dt]
starspec = 1e9
starspec = 10000000.0
starinst = 10.0
method = "adaptive"
minimum = 100.0
Expand Down Expand Up @@ -64,11 +64,10 @@ age_ini = 0.05

[star.mors]
age_now = 4.567
spec = "stellar_spectra/Named/sun.txt"
rot_pcntle = 50.0
rot_period = "none"
tracks = "spada"
star_name = "sun"
spectrum_source = "solar"

Comment on lines 65 to 71
[star.dummy]
Teff = 5772.0
Expand Down Expand Up @@ -99,7 +98,7 @@ radius_int = "none"

[atmos_clim]
module = "agni"
surf_state = "skin"
surf_state = "fixed"
prevent_warming = true
surface_d = 0.01
surface_k = 2.0
Expand All @@ -114,22 +113,20 @@ tmp_maximum = 5000.0
[atmos_clim.agni]
spectral_group = "Dayspring"
spectral_bands = "48"
p_top = 1e-04
p_top = 1e-05
p_obs = 0.02
surf_material = "surface_albedos/Hammond24/lunarmarebasalt.dat"
num_levels = 45
chemistry = "eq"
surf_material = "greybody"
num_levels = 42
chemistry = "none"
solve_energy = true
solution_atol = 0.08
solution_atol = 0.1
solution_rtol = 0.1
overlap_method = "ee"
rainout = false
condensation = false
latent_heat = false
real_gas = true
real_gas = false
psurf_thresh = 0.1
dx_max = 15.0
dx_max_ini = 200.0
ls_default = 1
dx_max = 10.0
max_steps = 70
perturb_all = true
mlt_criterion = "s"
Expand Down Expand Up @@ -159,9 +156,9 @@ rheo_phi_loc = 0.5
rheo_phi_wid = 0.2

[interior.spider]
ini_entropy = 3900.0
ini_entropy = 3000.0
ini_dsdr = -4.698e-06
num_levels = 150
num_levels = 140
mixing_length = 2
tolerance = 1e-08
tolerance_rel = 1e-08
Expand Down
Loading
Loading