## Chapter 3 Scripts: Elementary Radiation Balance

• `Planck.py`` `: Computes and plots the Planck function, and the cumulative emission (integral of the Planck function with respect to wavenumber, over the wavenumber range from 0 to nu)
• `SpectralBalance.py :`Plots the spectrum of incoming vs. outgoing radiation for the Solar System planets
• `EnergyBalance.py`` :` Generates the plot illustrating how a planet's surface temperature is determined by the balance between absorbed solar radiation and outgoing longwave radiation
• `PlanetT.py`` :` Computes and plots the isothermal blackbody temperatures of the Solar system planets
• `ERBEplot.py`` :` Makes the plot comparing the Earth's observed OLR with the OLR the surface would emit if the atmosphere were transparent to infrared.
• `IceAlbedoZeroD.py :` This makes plots showing various ways of analyzing the solutions of the zero-dimensional energy balance model with ice-albedo feedback. Among other things, it gives an example of how to compute the ice-albedo bifurcation diagram. This version uses solar constant (rather than radiating pressure) as the control variable for most of the calculations, but at the end illustrates how to use the radiating pressure (or some other parameter characterizing the greenhouse effect) as a control parameter.
• `Stability.py :`This script does various things to illustrate stability and hysteresis in the ice-albedo feedback system. It implements a time-dependent version of the ice-albedo feedback system, and shows how equilibrium is approached from various different initial conditions. It also shows how a hysteresis loop manifests itself in the time-dependent system with a time-varying control parameter. It illustrates an alternate method of generating a plot of the hysteresis loop, using a Newton solver to get the full set of solutions. This script isn't specifically referenced in the text, or used to make any of the plots, but it provides some additional insight.