Version 10 (modified by ajornet, 9 years ago) (diff) |
---|
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
modularisation of code used by nitrogen and phosphorus routines
There are many aspects which the nutrient cycles have in common, for example root uptake kinetics, stoichiometric considerations, etc. To avoid the risks of introducing inconsistencies and redundant code, I started to write subroutines which are listed in the following. All subroutines are in stomate_phosphorus.f90.
1. root_conductivity
This subroutine calculates the uptake capacity of N & P per mass root according to the uptake kinetics of ammonia and nitrate used in OCN (Zaehle & Friend, 2010).
2. f_XY_plant
This subroutine calculates the scaling functions based on the stoichiometric ratio of labile plant tissue (reserve, labile and leaf) which are used to scale root uptake, biological N2 fixation, and biochemical mineralization. It currently supports scaling functions based on the P-to-N ratio, N-to-C ratio and P-to-C ratio. All other combinations will cause a 'STOP', but could be easily implemented if needed.
SECTIONS
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. Also delete non necessary variables.
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
Attachments (11)
- missingPoints.png (20.5 KB) - added by dgoll 8 years ago.
-
USDA_map.png
(82.4 KB) -
added by dgoll 8 years ago.
The input file with the USDA soil orders in alphabetical order
-
USDA_map.2.png
(76.8 KB) -
added by dgoll 8 years ago.
input file with USDA soil orders in alphabetical order
-
USDA_map.3.png
(76.0 KB) -
added by dgoll 8 years ago.
The input file with the USDA soil orders in alphabetical order
-
lith_dom_map.png
(96.2 KB) -
added by dgoll 8 years ago.
dominant lithology
-
soilshield_map.png
(83.9 KB) -
added by dgoll 8 years ago.
soil shielding factor
-
ORC-CNP.png
(69.6 KB) -
added by dgoll 8 years ago.
Schematic representation of the key processes represented in ORCHIDEE-CNP
-
CNP_longterm.png
(46.6 KB) -
added by dgoll 8 years ago.
In red original formulation, in green the 3yr average. Shown is the C, the N and the CN of slow soil pool.
- BurnLabile.png (106.3 KB) - added by dgoll 8 years ago.
- Zoom.png (57.2 KB) - added by dgoll 8 years ago.
-
GPP_Hai.png
(48.2 KB) -
added by dgoll 8 years ago.
average annual cycle of DE-Hai site w & w/o N-limitation / w & w/o burning
Download all attachments as: .zip