modelling code used in paper of Marchingo et al., Science, 2014