GASTLI
Classes
Class defining objects for carrying out interior structure calculations |
Module Contents
- class GASTLI.int_planet(j_max=30, cnt_conv_max=3, conv_prec=1e-05, pow_law=0.32, chk_EOS=0, EOS_lim_P=[500000000000.0, 500000000000.0, 500000000000.0, 500000000000.0, 500000000000.0, 500000000000.0, 500000000000.0, 500000000000.0, 500000000000.0, 500000000000.0], corEOS=1)
Class defining objects for carrying out interior structure calculations for a given set of mass, composition and surface conditions
- Args:
- path_to_file:
path to input data directory
- j_max (Optional):
Maximum number of iterations - must be < 100.
- cnt_conv_max (Optional):
Maximum number of detected oscillations.
- conv_prec (Optional):
Precision of the convergence condition.
- pow_law (Optional):
Power exponent for planet radius estimation. Default is 0.32. Increase if planet is very massive
- chk_EOS (Optional):
Index to check the validity range of the used EOS. Do not change
- EOS_lim_P (Optional, array):
Upper limit in pressure. Do not change
- corEOS (Optional):
Type of correction on the thermodynamical parameters of Vinet EOS. (=0: no correction / =1: range [1:1.5] / =2: range [1:5] / =3: range [1:10]). Do not change
- path_to_file
- j_max = 30
- cnt_conv_max = 3
- conv_prec = 1e-05
- pow_law = 0.32
- chk_EOS = 0
- EOS_lim_P = [500000000000.0, 500000000000.0, 500000000000.0, 500000000000.0, 500000000000.0, 500000000000.0,...
- corEOS = 1
- n_lay
- n_EOS
- n_pts
- n_int
- n_mat_max
- f_alloy = 0.0
- MgD = 0.0
- MgSi = 0.0
- M_P = 0.0
- x_core = 0.0
- x_H2O = 0.0
- T_surf = 0.0
- P_surf = 0.0
- x_H2Oc = 0.0
- x_corec = 0.0
- M_Pc = 0.0
- R_P = 0.0
- FeSi = 0.0
- rho_p = 0.0
- intrf
- intrf_hist
- iter_num
- g
- T
- P
- rho
- cv
- r
- setup_parameters(f_alloy=cte.constants.f_alloy_e, MgD=cte.constants.mgd_e, MgSi=cte.constants.mgsi_e)
- Set up parameters for chemistry of mantle and core
Do not change any of the parameters below
- Args:
- f_alloy (Optional):
if no Fe in core, not applicable
- MgD (Optional):
if not Vinet Eq. for mantle, not applicable
- MgSi (Optional):
if not Vinet Eq. for mantle, not applicable
- Returns:
- x:
matrix with material properties
- calc_radius(M_P, x_core, x_H2O, T_surf, P_surf, Zenv)
Function that calculates planet radius and interior structure
- Args:
M_P: Planet mass in Earth masses x_core: Core mass fraction x_H2O: Envelope mass fraction. Defined as 1-x_core in this case T_surf: Outer surface (boundary) temperature in K P_surf: Outer surface (boundary) pressure in Pa
- Returns:
- x_H2Oc:
Re-computed envelope mass fraction (for checking)
- x_corec:
Re-computed core mass fraction (for checking)
- M_Pc:
Re-computed interior mass in Earth masses (for checking)
- R_P:
Planet radius in Earth radii
- FeSi:
Fe-to-Si ratio (not applicable if no Fe core)
- rho_p:
Planet density in g/cm3
- OtoH:
Envelope O:H ratio
- metal:
Envelope metallicity ([Fe/H] in x solar units)
- intrf:
Array with indexes of interfaces between layers
- g:
Gravity acceleration profile array, g(r), in cm/s2
- T:
Temperature profile array, T(r), in K
- P:
Pressure profile array, P(r), in Pa
- rho:
Density profile array, rho(r), in kg/m3
- cv:
Specific heat capacity at constant volume profile array, Cv(r), in J/kg/K
- entropy:
Entropy profile array, S(r), in J/kg/K
- r:
Radius profile array, r, in m