| 6 | | 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) |
| 7 | | == 1) refCN3238Palmi01 == |
| 8 | | * STOMATE_IMPOSE_CN = y |
| 9 | | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 10 | | * Resolution = 2 degree |
| 11 | | * !CyclicBegin = 1901 ; !CyclicEnd = 1910 |
| 12 | | * LAT = [-85] ; LON = ![33] |
| 13 | | * PERIOD LENGTH = 100 years |
| 14 | | * FREQUENCY OUTPUT = yearly |
| 15 | | * START TIME = 1901-01-01 |
| 16 | | * '''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. |
| 17 | | |
| 18 | | |
| 19 | | == 2) refCN3238Palmi02 == |
| 20 | | |
| 21 | | * STOMATE_IMPOSE_CN = y |
| 22 | | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 23 | | * Resolution = 2 degree |
| 24 | | * Imposing the PFT C3Ag (PFT10) for the whole grid. |
| 25 | | * !CyclicBegin = 1901 ; !CyclicEnd = 1910 |
| 26 | | * LAT = [-85] ; LON = ![33] |
| 27 | | * PERIOD LENGTH = 50 years |
| 28 | | * FREQUENCY OUTPUT = yearly |
| 29 | | * START TIME = 1901-01-01 |
| 30 | | * '''COMMENTS:''' The CN_ratio is not stable for the C3Ag PFT and stops exactly at the same time of the refCN3238Palmi01 simulation. |
| 31 | | |
| 32 | | |
| 33 | | == 3) refCN3238Palmi03 == |
| 34 | | |
| 35 | | * STOMATE_IMPOSE_CN = y |
| 36 | | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 37 | | * Resolution = 2 degree |
| 38 | | * Imposing the PFT C3Ag (PFT10) for the whole grid. |
| 39 | | * !CyclicBegin = 1901 ; !CyclicEnd = 1910 |
| 40 | | * LAT = [-85] ; LON = ![33] |
| 41 | | * PERIOD LENGTH = 50 years |
| 42 | | * FREQUENCY OUTPUT = daily |
| 43 | | * START TIME = 1901-01-01 |
| 44 | | * '''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. |
| 45 | | |
| 46 | | |
| 47 | | == 4) refCN3238Palmi04 == |
| 48 | | |
| 49 | | * STOMATE_IMPOSE_CN = y |
| 50 | | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 51 | | * Resolution = 2 degree |
| 52 | | * Imposing the PFT C3Ag (PFT10) for the whole grid. |
| 53 | | * !CyclicBegin = 1901 ; !CyclicEnd = 1910 |
| 54 | | * LAT = [-85] ; LON = ![33] |
| 55 | | * PERIOD LENTH = 20 years |
| 56 | | * FREQUENCY OUTPUT = daily |
| 57 | | * START TIME = 1901-01-01 |
| 58 | | * '''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. |
| 59 | | |
| 60 | | |
| 61 | | == 5) refCN3238Palmi05 == |
| 62 | | |
| 63 | | * STOMATE_IMPOSE_CN = y |
| 64 | | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 65 | | * Resolution = 2 degree |
| 66 | | * Imposing the PFT C3Ag (PFT10) for the whole grid. |
| 67 | | * !CyclicBegin = 1911 ; !CyclicEnd = 1911 |
| 68 | | * LAT = [-85] ; LON = ![33] |
| 69 | | * PERIOD LENTH = 5 years |
| 70 | | * FREQUENCY OUTPUT = daily |
| 71 | | * START TIME = 1911-01-01 |
| 72 | | * !RestartPath = /ccc/scratch/cont003/dsm/p24palm/IGCM_OUT/OL2/TEST/refCN3238Palmi04 |
| 73 | | * !RestartJobName = FG1low.CN.Palmi04 |
| 74 | | * '''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. |
| 75 | | |
| 76 | | |
| 77 | | == 6) refCN3238Palmi06 == |
| 78 | | |
| 79 | | * STOMATE_IMPOSE_CN = y |
| 80 | | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 81 | | * Resolution = 2 degree |
| 82 | | * Imposing the PFT C3Ag (PFT10) for the whole grid. |
| 83 | | * !CyclicBegin = 1911 ; !CyclicEnd = 1911 |
| 84 | | * LAT = [-85] ; LON = ![33] |
| 85 | | * PERIOD LENTH = 5 years |
| 86 | | * FREQUENCY OUTPUT = daily |
| 87 | | * START TIME = 1911-01-01 |
| 88 | | * !RestartPath = /ccc/scratch/cont003/dsm/p24palm/IGCM_OUT/OL2/TEST/refCN3238Palmi04 |
| 89 | | * !RestartJobName = FG1low.CN.Palmi04 |
| 90 | | * '''COMMENTS:''' If plant_uptake(:,m,initrate) is negative, we impose plant_uptake(:,m,initrate) = 0. If we impose it the CN_ratio is stable. |
| 91 | | |
| 92 | | |
| 93 | | == 7) refCN3238Palmi07 == |
| 94 | | |
| 95 | | * STOMATE_IMPOSE_CN = y |
| 96 | | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 97 | | * Resolution = 2 degree |
| 98 | | * Imposing the PFT C3Ag (PFT10) for the whole grid. |
| 99 | | * !CyclicBegin = 1911; !CyclicEnd = 1911 |
| 100 | | * LAT = [-85] ; LON = ![33] |
| 101 | | * PERIOD LENTH = 5 years |
| 102 | | * FREQUENCY OUTPUT = daily |
| 103 | | * START TIME = 1911-01-01 |
| 104 | | * !RestartPath = /ccc/scratch/cont003/dsm/p24palm/IGCM_OUT/OL2/TEST/refCN3238Palmi04 |
| 105 | | * !RestartJobName = FG1low.CN.Palmi04 |
| 106 | | * '''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. |
| 107 | | |
| 108 | | |
| 109 | | == 8) refCN3238Palmi08 == |
| 110 | | |
| 111 | | * STOMATE_IMPOSE_CN = y |
| 112 | | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 113 | | * Resolution = 2 degree |
| 114 | | * !CyclicBegin = 1901 ; !CyclicEnd = 1910 |
| 115 | | * LAT = [-179.0, 179.0] ; LON = [-89.0, 89.0] |
| 116 | | * PERIOD LENTH = 30 years |
| 117 | | * FREQUENCY OUTPUT = monthly |
| 118 | | * START TIME = 1901-01-01 |
| 119 | | * '''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). |
| 120 | | |
| 121 | | |
| 122 | | == 9) refCN3238Palmi09 == |
| 123 | | |
| 124 | | * STOMATE_IMPOSE_CN = y |
| 125 | | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 126 | | * Resolution = 2 degree |
| 127 | | * !CyclicBegin = 1901 ; !CyclicEnd = 1910 |
| 128 | | * LAT = [-179.0, 179.0] ; LON = [-89.0, 89.0] |
| 129 | | * PERIOD LENTH = 30 years |
| 130 | | * FREQUENCY OUTPUT = monthly |
| 131 | | * START TIME = 1901-01-01 |
| 132 | | * '''COMMENTS:''' Global simulation, with no code change compared to refCN3238Palmi07 simulation. It is the reference run to compare with refCN3238Palmi08 and refCN3238Palmi10 simulations. |
| 133 | | |
| 134 | | |
| 135 | | == 10) refCN3238Palmi10 == |
| 136 | | |
| 137 | | * STOMATE_IMPOSE_CN = y |
| 138 | | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 139 | | * Resolution = 2 degree |
| 140 | | * !CyclicBegin = 1901 ; !CyclicEnd = 1910 |
| 141 | | * LAT = [-179.0, 179.0] ; LON = [-89.0, 89.0] |
| 142 | | * PERIOD LENTH = 30 years |
| 143 | | * FREQUENCY OUTPUT = monthly |
| 144 | | * START TIME = 1901-01-01 |
| 145 | | * '''COMMENTS:''' Global simulation, leaching is calculated after n_adsorbed in stomate_soilcarbon.f90. The CN_ratio is more stable, but anyway starts to change. |
| 146 | | |
| 147 | | == 11) refCN3238Palmi11 == |
| 148 | | |
| 149 | | * STOMATE_IMPOSE_CN = y |
| 150 | | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 151 | | * Resolution = 2 degree |
| 152 | | * !CyclicBegin = 1901 ; !CyclicEnd = 1910 |
| 153 | | * LAT = [-179.0, 179.0] |
| 154 | | * LON = [-89.0, 89.0] |
| 155 | | * PERIOD LENTH = 200 years |
| 156 | | * FREQUENCY OUTPUT = yearly |
| 157 | | * START TIME = 1901-01-01 |
| 158 | | * '''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. |
| 159 | | |
| 160 | | == 12) refCN3238Palmi12 == |
| 161 | | |
| 162 | | * STOMATE_IMPOSE_CN = y |
| 163 | | * Set-up as Nicolas run FG1low.Choi.CN3238 |
| 164 | | * Resolution = 2 degree |
| 165 | | * !CyclicBegin = 1901 ; !CyclicEnd = 1910 |
| 166 | | * LAT = ![47] |
| 167 | | * LON = ![83] |
| 168 | | * Imposing the PFT !BoBrDe (PFT08) for the whole grid. |
| 169 | | * PERIOD LENTH = 200 years |
| 170 | | * FREQUENCY OUTPUT = daily |
| 171 | | * START TIME = 1901-01-01 |
| 172 | | * '''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. |
