README.txt

Execute the files in the C_code folder in the following order:
0) Run erosion_per_PFT.py and part 3 of data_processing_part2.ksh in the E_code folder to derive the erosion rates per PFT.
   However, one needs first to obtain the PFT map from Peng et al. (2017) and remap it to 5arcmin resolution. 
   Be aware that erosion rates provided in this study are only for the Rhine catchment while the erosion_per_PFT.py file needs the erosion input file for the entire globe.
1) f2py -c -m mod_append mod_append.f90
2) f2py -c -m mod_sorting mod_sorting.f90
3) f2py -c -m routing_main_basin routing_main_basin.f90
4) f2py -c -m routing_main_all routing_main_all.f90
5) Enter your input/output directories and choose simulation type (S0,S1,S2) in the files Run_eq.py, func1.py, func2.py for the equilibrium simulation
6) python Run_eq.py
7) Enter your input/output directories and choose simulation type (S0,S1,S2) in the files Run_tr.py, func1_tr.py, func2_tr.py for the transient simulation
The derivation of the C pools and fluxes to construct the turnover rates (K's) of ORCHIDEE land surface model is described in detail by Naipal et al. (2018) published in Biogeosciences.

Data needed to run the model are available at the repository Zenodo:
1) Erosion rates under climate change and land use change for the period 1850-2005 for the Rhine catchment per land cover type
2) The turnover rates (K's) of ORCHIDEE land surface model needed for the Emulator extracted for the Rhine catchment for the period 1850-2005
3) Data on the coefficient 're' of the soil organic carbon profile (equation 13 in paper Naipal et al. (2019))
4) Scaled topographical slope for the Rhine
5) Aggregated flow-accumulation for the Rhine
6) Derived stream-length for the Rhine
7) Grid cell area at 5 arcmin resolution
8) Rhine subbasin division
9) Mask files 

Data used from other studies are not included in the repository but are needed to run the model:
1) Vegetation fractions from the PFT map of Peng et al. (2017)
2) Soil information (bulk density)
3) Hillslope and lowland fractions from Pelletier et al. (2016)
Please refer to Naipal et al. (2019) for information on how to obtain this data. When using this data please refer to the respective study.

When using the carbon erosion model (CE-DYNAM) or some of its methods or modeldata please cite the paper Naipal et al. (2019). 