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


Ignore:
Timestamp:
2016-05-10T12:01:58+02:00 (8 years ago)
Author:
pmessina
Comment:

--

Legend:

Unmodified
Added
Removed
Modified

  • DevelopmentActivities/ORCHIDEE-CN/Global

    v1 v2  
    11= Global evaluation of ORCHIDEE-CN = 
     2 
     3 
     4We 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) 
     5== 1) refCN3238Palmi01 ==  
     6* STOMATE_IMPOSE_CN = y 
     7* Set-up as Nicolas run FG1low.Choi.CN3238 
     8* !CyclicBegin = 1901 ; !CyclicEnd = 1910 
     9* LAT = [-85] ; LON = ![33] 
     10* PERIOD LENGTH = 100 years 
     11* FREQUENCY OUTPUT = yearly 
     12* START TIME = 1901-01-01 
     13* '''COMMENTS:''' CN_ratio is not stable for all PFTs in the grid cell. Depending on PFTs, the departure form the fixed value happens in different moment of the simulation. 
     14 
     15 
     16== 2) refCN3238Palmi02 == 
     17 
     18* STOMATE_IMPOSE_CN = y 
     19* Set-up as Nicolas run FG1low.Choi.CN3238 
     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* Imposing the PFT C3Ag (PFT10) for the whole grid. 
     34* !CyclicBegin = 1901 ; !CyclicEnd = 1910 
     35* LAT = [-85] ; LON = ![33] 
     36* PERIOD LENGTH = 50 years 
     37* FREQUENCY OUTPUT = daily 
     38* START TIME = 1901-01-01 
     39* '''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. 
     40 
     41 
     42== 4) refCN3238Palmi04 ==   
     43 
     44* STOMATE_IMPOSE_CN = y 
     45* Set-up as Nicolas run FG1low.Choi.CN3238 
     46* Imposing the PFT C3Ag (PFT10) for the whole grid. 
     47* !CyclicBegin = 1901 ; !CyclicEnd = 1910 
     48* LAT = [-85] ; LON = ![33]  
     49* PERIOD LENTH = 20 years  
     50* FREQUENCY OUTPUT = daily  
     51* START TIME = 1901-01-01 
     52* '''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.   
     53 
     54 
     55== 5) refCN3238Palmi05 == 
     56 
     57* STOMATE_IMPOSE_CN = y 
     58* Set-up as Nicolas run FG1low.Choi.CN3238 
     59* Imposing the PFT C3Ag (PFT10) for the whole grid. 
     60* !CyclicBegin = 1911 ; !CyclicEnd = 1911 
     61* LAT = [-85] ; LON = ![33] 
     62* PERIOD LENTH = 5 years 
     63* FREQUENCY OUTPUT = daily 
     64* START TIME = 1911-01-01 
     65* !RestartPath = /ccc/scratch/cont003/dsm/p24palm/IGCM_OUT/OL2/TEST/refCN3238Palmi04 
     66* !RestartJobName = FG1low.CN.Palmi04 
     67* '''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.   
     68 
     69 
     70== 6) refCN3238Palmi06 ==   
     71 
     72* STOMATE_IMPOSE_CN = y 
     73* Set-up as Nicolas run FG1low.Choi.CN3238 
     74* Imposing the PFT C3Ag (PFT10) for the whole grid. 
     75* !CyclicBegin = 1911 ; !CyclicEnd = 1911 
     76* LAT = [-85] ; LON = ![33] 
     77* PERIOD LENTH = 5 years 
     78* FREQUENCY OUTPUT = daily 
     79* START TIME = 1911-01-01 
     80* !RestartPath = /ccc/scratch/cont003/dsm/p24palm/IGCM_OUT/OL2/TEST/refCN3238Palmi04 
     81* !RestartJobName = FG1low.CN.Palmi04 
     82* '''COMMENTS:''' If plant_uptake(:,m,initrate) is negative, we impose plant_uptake(:,m,initrate) = 0. If we impose it the CN_ratio is stable.  
     83 
     84 
     85== 7) refCN3238Palmi07 == 
     86 
     87* STOMATE_IMPOSE_CN = y 
     88* Set-up as Nicolas run FG1low.Choi.CN3238 
     89* Imposing the PFT C3Ag (PFT10) for the whole grid. 
     90* !CyclicBegin = 1911; !CyclicEnd = 1911 
     91* LAT = [-85] ; LON = ![33] 
     92* PERIOD LENTH = 5 years 
     93* FREQUENCY OUTPUT = daily 
     94* START TIME = 1911-01-01 
     95* !RestartPath = /ccc/scratch/cont003/dsm/p24palm/IGCM_OUT/OL2/TEST/refCN3238Palmi04 
     96* !RestartJobName = FG1low.CN.Palmi04 
     97* '''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. 
     98 
     99 
     100== 8) refCN3238Palmi08 ==   
     101 
     102* STOMATE_IMPOSE_CN = y 
     103* Set-up as Nicolas run FG1low.Choi.CN3238 
     104* !CyclicBegin = 1901 ; !CyclicEnd = 1910 
     105* LAT = [-179.0, 179.0] ; LON = [-89.0, 89.0]  
     106* PERIOD LENTH = 30 years 
     107* FREQUENCY OUTPUT = monthly  
     108* START TIME = 1901-01-01 
     109* '''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).  
     110 
     111 
     112== 9) refCN3238Palmi09 ==  
     113 
     114* STOMATE_IMPOSE_CN = y 
     115* Set-up as Nicolas run FG1low.Choi.CN3238 
     116* !CyclicBegin = 1901 ; !CyclicEnd = 1910 
     117* LAT = [-179.0, 179.0] ; LON = [-89.0, 89.0] 
     118* PERIOD LENTH = 30 years 
     119* FREQUENCY OUTPUT = monthly 
     120* START TIME = 1901-01-01 
     121* '''COMMENTS:''' Global simulation, with no code change compared to refCN3238Palmi07 simulation. It is the reference run to compare with refCN3238Palmi08 and refCN3238Palmi10 simulations.  
     122 
     123 
     124== 10) refCN3238Palmi10 == 
     125 
     126* STOMATE_IMPOSE_CN = y 
     127* Set-up as Nicolas run FG1low.Choi.CN3238 
     128* !CyclicBegin = 1901 ; !CyclicEnd = 1910 
     129* LAT = [-179.0, 179.0] ; LON = [-89.0, 89.0] 
     130* PERIOD LENTH = 30 years 
     131* FREQUENCY OUTPUT = monthly 
     132* START TIME = 1901-01-01 
     133* '''COMMENTS:''' Global simulation, leaching is calculated after n_adsorbed in stomate_soilcarbon.f90. The CN_ratio is more stable, but anyway starts to change.  
     134 
     135== 11) refCN3238Palmi11 ==   
     136 
     137* STOMATE_IMPOSE_CN = y 
     138* Set-up as Nicolas run FG1low.Choi.CN3238 
     139* !CyclicBegin = 1901 ; !CyclicEnd = 1910 
     140* LAT = [-179.0, 179.0] 
     141* LON = [-89.0, 89.0]  
     142* PERIOD LENTH = 200 years 
     143* FREQUENCY OUTPUT = yearly  
     144* START TIME = 1901-01-01 
     145* '''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.