organic carbon model in R