The code in model.demo.Rmd demonstrates the use of loadflex for model fitting, predictions, and biogeochemical analyses as described in the text, figures, and tables 1 and 2 of the main article. The code in model.assessment.Rmd demonstrates the computation of model performance metrics as described in the text and table 3 of the main article.