1 | | pp |
| 1 | |
| 2 | We list here the simulations done for testing ORCHIDEE-CN at global scale, maintaining the ratio LEAF_M_c/LEAF_M_n fix (hereafter CN_ratio) |
| 3 | == 1) refCN3238Palmi01 == |
| 4 | * STOMATE_IMPOSE_CN = y |
| 5 | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 6 | * Resolution = 2 degree |
| 7 | * !CyclicBegin = 1901 ; !CyclicEnd = 1910 |
| 8 | * LAT = [-85] ; LON = ![33] |
| 9 | * PERIOD LENGTH = 100 years |
| 10 | * FREQUENCY OUTPUT = yearly |
| 11 | * START TIME = 1901-01-01 |
| 12 | * '''COMMENTS:''' CN_ratio is not stable for all PFTs in the grid cell. Depending on PFTs, the departure from the fixed value happens in different moment of the simulation. |
| 13 | |
| 14 | |
| 15 | == 2) refCN3238Palmi02 == |
| 16 | |
| 17 | * STOMATE_IMPOSE_CN = y |
| 18 | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 19 | * Resolution = 2 degree |
| 20 | * Imposing the PFT C3Ag (PFT10) for the whole grid. |
| 21 | * !CyclicBegin = 1901 ; !CyclicEnd = 1910 |
| 22 | * LAT = [-85] ; LON = ![33] |
| 23 | * PERIOD LENGTH = 50 years |
| 24 | * FREQUENCY OUTPUT = yearly |
| 25 | * START TIME = 1901-01-01 |
| 26 | * '''COMMENTS:''' The CN_ratio is not stable for the C3Ag PFT and stops exactly at the same time of the refCN3238Palmi01 simulation. |
| 27 | |
| 28 | |
| 29 | == 3) refCN3238Palmi03 == |
| 30 | |
| 31 | * STOMATE_IMPOSE_CN = y |
| 32 | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 33 | * Resolution = 2 degree |
| 34 | * Imposing the PFT C3Ag (PFT10) for the whole grid. |
| 35 | * !CyclicBegin = 1901 ; !CyclicEnd = 1910 |
| 36 | * LAT = [-85] ; LON = ![33] |
| 37 | * PERIOD LENGTH = 50 years |
| 38 | * FREQUENCY OUTPUT = daily |
| 39 | * START TIME = 1901-01-01 |
| 40 | * '''COMMENTS:''' The CN_ratio is not stable, we better investigate why it diverges. Probably some problems in stomate_growth_fun_all.f90 where the "impose_cn" condition is treated. |
| 41 | |
| 42 | |
| 43 | == 4) refCN3238Palmi04 == |
| 44 | |
| 45 | * STOMATE_IMPOSE_CN = y |
| 46 | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 47 | * Resolution = 2 degree |
| 48 | * Imposing the PFT C3Ag (PFT10) for the whole grid. |
| 49 | * !CyclicBegin = 1901 ; !CyclicEnd = 1910 |
| 50 | * LAT = [-85] ; LON = ![33] |
| 51 | * PERIOD LENTH = 20 years |
| 52 | * FREQUENCY OUTPUT = daily |
| 53 | * START TIME = 1901-01-01 |
| 54 | * '''COMMENTS:''' Some changes in stomate_growth_fun_all.f90 to optimize the part of the code linked to "impose_cn" condition. Things are not really improved. The origin of problems is more likely that N_UPTAKE_NO3 becomes negative and N_UPTAKE_NH4 becomes zero. |
| 55 | |
| 56 | |
| 57 | == 5) refCN3238Palmi05 == |
| 58 | |
| 59 | * STOMATE_IMPOSE_CN = y |
| 60 | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 61 | * Resolution = 2 degree |
| 62 | * Imposing the PFT C3Ag (PFT10) for the whole grid. |
| 63 | * !CyclicBegin = 1911 ; !CyclicEnd = 1911 |
| 64 | * LAT = [-85] ; LON = ![33] |
| 65 | * PERIOD LENTH = 5 years |
| 66 | * FREQUENCY OUTPUT = daily |
| 67 | * START TIME = 1911-01-01 |
| 68 | * !RestartPath = /ccc/scratch/cont003/dsm/p24palm/IGCM_OUT/OL2/TEST/refCN3238Palmi04 |
| 69 | * !RestartJobName = FG1low.CN.Palmi04 |
| 70 | * '''COMMENTS:''' We investigate when CN_ratio starts to change. Printed many variables in stomate_soilcarbon.f90 into the stardand output at each time step (every 30 minutes) to see where the problem comes from. We see that plant_uptake(:,m,initrate) is negative. |
| 71 | |
| 72 | |
| 73 | == 6) refCN3238Palmi06 == |
| 74 | |
| 75 | * STOMATE_IMPOSE_CN = y |
| 76 | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 77 | * Resolution = 2 degree |
| 78 | * Imposing the PFT C3Ag (PFT10) for the whole grid. |
| 79 | * !CyclicBegin = 1911 ; !CyclicEnd = 1911 |
| 80 | * LAT = [-85] ; LON = ![33] |
| 81 | * PERIOD LENTH = 5 years |
| 82 | * FREQUENCY OUTPUT = daily |
| 83 | * START TIME = 1911-01-01 |
| 84 | * !RestartPath = /ccc/scratch/cont003/dsm/p24palm/IGCM_OUT/OL2/TEST/refCN3238Palmi04 |
| 85 | * !RestartJobName = FG1low.CN.Palmi04 |
| 86 | * '''COMMENTS:''' If plant_uptake(:,m,initrate) is negative, we impose plant_uptake(:,m,initrate) = 0. If we impose it the CN_ratio is stable. |
| 87 | |
| 88 | |
| 89 | == 7) refCN3238Palmi07 == |
| 90 | |
| 91 | * STOMATE_IMPOSE_CN = y |
| 92 | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 93 | * Resolution = 2 degree |
| 94 | * Imposing the PFT C3Ag (PFT10) for the whole grid. |
| 95 | * !CyclicBegin = 1911; !CyclicEnd = 1911 |
| 96 | * LAT = [-85] ; LON = ![33] |
| 97 | * PERIOD LENTH = 5 years |
| 98 | * FREQUENCY OUTPUT = daily |
| 99 | * START TIME = 1911-01-01 |
| 100 | * !RestartPath = /ccc/scratch/cont003/dsm/p24palm/IGCM_OUT/OL2/TEST/refCN3238Palmi04 |
| 101 | * !RestartJobName = FG1low.CN.Palmi04 |
| 102 | * '''COMMENTS:''' Printed many more variables in stomate_soilcarbon.f90 in the stardand output at each time step (every 30 minutes), to see why plant_uptake(:,m,initrate) becomes negative. We find that when immob(:,:) is bigger than the nitrogen in the two soil pools (soil_n_min pools (nitrate and ammonium) and soil_n_min(:,:,initrate)), plant_uptake(:,m,initrate) becomes negative. |
| 103 | |
| 104 | |
| 105 | == 8) refCN3238Palmi08 == |
| 106 | |
| 107 | * STOMATE_IMPOSE_CN = y |
| 108 | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 109 | * Resolution = 2 degree |
| 110 | * !CyclicBegin = 1901 ; !CyclicEnd = 1910 |
| 111 | * LAT = [-179.0, 179.0] ; LON = [-89.0, 89.0] |
| 112 | * PERIOD LENTH = 30 years |
| 113 | * FREQUENCY OUTPUT = monthly |
| 114 | * START TIME = 1901-01-01 |
| 115 | * '''COMMENTS:''' Global simulation imposing leaching = 0 in stomate_soilcarbon.f90. We test the impact of leaching on CN_ratio. This simulation gives weird results. The CN ratio change very little in comparison to reference simulations (refCN3238Palmi09) and performs worse that refCN3238Palmi10 (see below). |
| 116 | |
| 117 | |
| 118 | == 9) refCN3238Palmi09 == |
| 119 | |
| 120 | * STOMATE_IMPOSE_CN = y |
| 121 | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 122 | * Resolution = 2 degree |
| 123 | * !CyclicBegin = 1901 ; !CyclicEnd = 1910 |
| 124 | * LAT = [-179.0, 179.0] ; LON = [-89.0, 89.0] |
| 125 | * PERIOD LENTH = 30 years |
| 126 | * FREQUENCY OUTPUT = monthly |
| 127 | * START TIME = 1901-01-01 |
| 128 | * '''COMMENTS:''' Global simulation, with no code change compared to refCN3238Palmi07 simulation. It is the reference run to compare with refCN3238Palmi08 and refCN3238Palmi10 simulations. |
| 129 | |
| 130 | |
| 131 | == 10) refCN3238Palmi10 == |
| 132 | |
| 133 | * STOMATE_IMPOSE_CN = y |
| 134 | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 135 | * Resolution = 2 degree |
| 136 | * !CyclicBegin = 1901 ; !CyclicEnd = 1910 |
| 137 | * LAT = [-179.0, 179.0] ; LON = [-89.0, 89.0] |
| 138 | * PERIOD LENTH = 30 years |
| 139 | * FREQUENCY OUTPUT = monthly |
| 140 | * START TIME = 1901-01-01 |
| 141 | * '''COMMENTS:''' Global simulation, leaching is calculated after n_adsorbed in stomate_soilcarbon.f90. The CN_ratio is more stable, but anyway starts to change. |
| 142 | |
| 143 | == 11) refCN3238Palmi11 == |
| 144 | |
| 145 | * STOMATE_IMPOSE_CN = y |
| 146 | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 147 | * Resolution = 2 degree |
| 148 | * !CyclicBegin = 1901 ; !CyclicEnd = 1910 |
| 149 | * LAT = [-179.0, 179.0] |
| 150 | * LON = [-89.0, 89.0] |
| 151 | * PERIOD LENTH = 200 years |
| 152 | * FREQUENCY OUTPUT = yearly |
| 153 | * START TIME = 1901-01-01 |
| 154 | * '''COMMENTS:''' in stomate_soilcarbon.f90 we add nitrogen when immob(:,:) is bigger than the nitrogen that can be given by the soil pools (nitrate and ammonium) and so soil_n_min(:,:,initrate) become negative. We take memory of the added nitrogen in N_support. We calculate leaching after n_adsorbed calculation (as in Zaehle version), in other words the ammonium available to be leached is decreased. This makes sense as the ammonium in less soluble than nitrate. |
| 155 | |
| 156 | == 12) refCN3238Palmi12 == |
| 157 | |
| 158 | * STOMATE_IMPOSE_CN = y |
| 159 | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 160 | * Resolution = 2 degree |
| 161 | * !CyclicBegin = 1901 ; !CyclicEnd = 1910 |
| 162 | * LAT = ![47] |
| 163 | * LON = ![83] |
| 164 | * Imposing the PFT !BoBrDe (PFT08) for the whole grid. |
| 165 | * PERIOD LENTH = 200 years |
| 166 | * FREQUENCY OUTPUT = daily |
| 167 | * START TIME = 1901-01-01 |
| 168 | * '''COMMENTS:''' Corrected a bug in stomate_phenology.f90 related to the calculation of labile biomass calculation. There are still some points that do not maintain a fixed CN_ratio in refCN3238Palmi11. We investigate in one of these points (China) to see if now it is fixed. |