#
# +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
# photoferrotrophy.R
# +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
#
1-D water column model for competition between oxygenic and anoxygenic phototrophs
#
This script models competition between oxygenic (cyanobacterial) and anoxygenic (photoferrotrophic) photosynthetic bacteria in an idealized 1-D water column.  The model parameters and equations are described in K. Ozaki, K.J. Tompson, R.L. Simister, S.A. Crowe, and C.T. Reinhard [2019] "Anoxygenic photosynthesis and the delayed oxygenation of Earth's atmosphere", Nature Communications, doi:10.1038/s41467-019-10872-z
#
As a courtesy, we request that others who use this code please cite Ozaki et al. [2019].  We also request that those who use/modify the code please send publications and/or modified code to the corresponding author [chris.reinhard@eas.gatech.edu]
#
REQUIREMENTS: R and/or R Studio, packages ReacTran, rootSolve, deSolve, shape, marelac, tictoc
#
TO RUN THE CODE
From the command line:
[1] navigate to the directory hosting photoferrotrophy.R
[2] execute: R script photoferrotrophy.R
[3] output will be stored as "photosynth.out.csv" in the directory specified in Line 200 of photoferrotrophy.R
#
Using R or R Studio:
[1] open new console window
[2] execute: source('~/{DIRECTORY}/photoferrotrophy.R')
[3] output will be stored as above
#
# +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++