ORCHIDEE-CNP

This page describes the work done to include the phosphorus cycle into ORCHIDEE. It is based on MERGE-OCN, which was extended and corrected as described here: ​https://forge.ipsl.jussieu.fr/orchidee/wiki/Branches/MergeOCN/Goll as well as a recruitment routine by S. Peng (ask him if you need infos) The phosphorus cycle is an adaptation of the the model described by ​http://www.biogeosciences.net/9/3547/2012/bg-9-3547-2012.html

Soil type problem

The phosphorus model needs soil type specific parameters. Therefore, a new input routine is going to be written which reads in USDA soil orders. According to the soil order the parameter values are set for each grid box.

The problem is that ORCHIDEE already reads in fields of soil texture classes which were originally derived from soil orders. According to the soil texture classes, parameter related to soil hydrology are set. As the original data cannot be found anymore, ORCHIDEE will run with two independent sets of soil order specific parameters. Thereby inconsistencies between the phosphorus related parameters and the the hydrology related parameters can be introduced. The practice of reading in soil texture classes instead of soil orders is also problematic in respect of future developments and I suggest to change that.

Work for Albert

We need a routine which reads the soil orders (n=12) from this file: /home/users/dgoll/ORC_data/USDA_soilorders/USDA_SoilSuborder.nc .The file itself does not respect the netCDF standards, you would need to fix this.

Currently, all input reading routines are located in src_stomate/stomate_io.f90; for example the routines to read ref_temp, deposition and the lithological fractions. There is quite some doubling of code which could be minimized by re-writting the routines, for example outsourcing the remap procedure in a own subroutine/function.

The soil orders as read in from the file should be assigned to a new variable of your choice (following the texture classes or lithological fractions). I handle the rest.

You find the deposition input file here: /home/scratch01/dgoll/ORCHIDEE/ORC-CNP/simulations/simple_forcing_DSG/deposition.nc