Changes between Version 1 and Version 2 of DevelopmentActivities/ORCHIDEE-CN/Global/Firstest


Ignore:
Timestamp:
2016-09-02T12:41:18+02:00 (8 years ago)
Author:
pmessina
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • DevelopmentActivities/ORCHIDEE-CN/Global/Firstest

    v1 v2  
    1 pp 
     1 
     2We 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.