Belinda E. Medlyn, Martin G. De Kauwe, Yan-Shih Lin, Jürgen Knauer, Remko A. Duursma, Christopher A. Williams, Almut Arneth, Rob Clement, Peter Isaac, Jean-Marc Limousin, Maj-Lena Linderson, Patrick Meir, Nicolas Martin-StPaul and Lisa Wingate.

Repository containing all the code associated with our paper.

• Leaf gas exchange data is from Lin et al. (2015) Optimal stomatal behaviour around the world. Nature Climate Change, 5, 459–464.
• Stable isotope dataset: Dryad Digital Repository
• Eddy covariance dataset: FLUXNET
• Elevation data: GTOPO30
• NACP Synergetic (SYNMAP) C₄ relative fraction: Jung, M., Henkel, K., Herold, M., and Churkina, G.: Exploiting synergies of global land cover products for carbon cycle modeling. Remote Sens. Environ., 101, 534–553, doi:10.1016/j.rse.2006.01.020, 2006..
• Global CO₂ data: Ed Dlugokencky and Pieter Tans, NOAA/ESRL

To generate all g1 fits from leaf, isotope and fluxnet data and manuscript figures, simple type:

$ make

You will need a few python packages, namely numpy, pandas, matplotlib scipy and lmfit. The fitting code also exploits the python MPI library multiprocessing to speed things up. The R code depends on the doby, multcomp, lme4 and LMERConvenienceFunctions packages.

• Belinda Medlyn: B.Medlyn at westernsydney.edu.au
• Martin De Kauwe: mdekauwe at gmail.com