R code related to SRB model Geobiology (2011), 9, 446–457