PaleoBudyko
PaleoBudyko is a MATLAB program allowing the user to investigate the last five million years of climate data and Earth orbital parameters and compare them to the output of the Budyko model. The program consists of a MATLAB GUI (PaleoBudyko.m) and a data file (data.mat).
You can download the files here:
PaleoBudyko.m
data.mat
To run the program, type 'PaleoBudyko' at the MATLAB command line. Make sure that both files are in the MATHLAB path, preferably in the current directory.
Credits
Computation of the Milankovitch Cycles are due to Laskar et al, "A longterm numerical solution for the insolation quantities of the Earth," Astronomy & Astrophysics 428 (2004), 261–285.
The computations are posted on the Web at
http://www.imcce.fr/Equipes/ASD/insola/earth/La2004/index.html.
Climate data are due to Lisiecki and Raymo, "A PliocenePleistocene stack of 57 globally distributed benthic δ^{18}O records," Paleoceanography 20 (2005), PA1003.
The data are posted on the Web at
http://lorrainelisiecki.com/LR04stack.txt.
A description of Budyko's model and the computation of its stable equilibria can be found in McGehee & Lehman, "A Paleoclimate Model of IceAlbedo Feedback Forced by Variations in Earth's Orbit," SIAM J. Applied Dynamical Systems 11 (2012), 684707.
DOI 10.1137/10079879X
