: Computes and plots the CO2 partial pressure in equilibrium with three different silicate minerals as a function of temperature.
WHAK.py :This script implements a variant of the Walker, Hayes And Kasting silicate weathering thermostat calculation. It calculations the equilibrium CO2 partial pressure and equilibrium surface temperature as a function of stellar constant and CO2 outgassing rate.
CarbonateEq.py :This is a collection of routines used in calculations involving ocean carbonate/bicarbonate equilibria, and partitioning of carbon between atmospheric and oceanic reservoirs. It carries out the calculations used in Section 8.4 and associated workbook problems.
:This is a simple utility script that uses an integral over the Planck function to compute the blackbody photon flux in various wavenumber ranges. It is primarily used in the discussion of the UV and EUV photon fluxes that drive photochemistry and atmospheric escape.
RayleighJeans.py :Computes atmospheric escape flux due to the Rayleigh-Jeans mechanism. It also computes exobase heights and properties. This was used to make the table of Rayleigh-Jeans escape rates in the text.
ColdTrap.py :This script estimates the homopause mixing ratio of a condensable substance by computing the mixing ratio at the cold trap (the point of minimum saturation mixing ratio encountered in an atmosphere in radiative-convective equilibrium). This is a simplified version of the calculation, which specifies a tropopause height and an isothermal stratosphere, without solving for the full stratospheric radiative equilibrium.
HydroEscape.py :This script calculates atmospheric structure and escape flux for steady, radially symmetric hydrodynamic escape. It can handle both the adiabatic case, and the case where escape is sustained by a specified radiative heating profile. Convergence of the iteration to satisfy the transonic boundary condition is quite delicate in the latter case. At some future point, I will supplement this code with a 1D time-stepping code, which provides more reliable convergence.
Impact.py :Carries out simple estimates of atmosphere loss by impact erosion. Note that this is not a detailed hydrodynamic calculation. It makes use of fairly crude scaling laws to estimate how much atmospheric mass is blown off as a function of impact energy, and integrates over a hypothetical spectrum of impactors. The scaling laws employed are not expected to work for giant impacts which significantly deplace the whole planet.