Changes between Version 28 and Version 29 of DevelopmentActivities/ORCHIDEE-CN/NH3


Ignore:
Timestamp:
2017-07-10T13:30:41+02:00 (7 years ago)
Author:
taudoux
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • DevelopmentActivities/ORCHIDEE-CN/NH3

    v28 v29  
    8484 
    8585=== Evaluation of ammonia emissions results calculated by the model (Fertilizer added) 
     86 
     87 
     88. Scripts used to generate the files.txt are in the folder : [[BR]]  
     89/home/users/taudoux/documents/Results/B/Script [[BR]] 
     90For the script : script_delta_X.ksh  [[BR]] 
     91X depends on the simulation we are analyzing.[[BR]] 
     92arg1= Name of the site [[BR]] 
     93arg2= PFT [[BR]] 
     94arg3= Last Time Step [[BR]] 
     95 
     96Others scripts don't need arguments and generate files.txt according to their appellation : [[BR]] 
     97. 0.5 --> Simuls_ratio0.5_v1/ ; 0.8 --> Simuls_ratio0.8_v1/ [[BR]]  
     98. e --> EFe/ ; f --> EFf ; g --> EFg etc. 
     99 
     100.Scripts python used to generate graphics are in the folder : [[BR]]  
     101/home/users/taudoux/documents/Results/B/barplot_X.py [[BR]] 
     102 
     103X depends on the simulation we are analyzing.[[BR]] 
     104 
     105. files.txt are in the folders : [[BR]] 
     106/home/users/taudoux/documents/Results/B/Simuls_ratio0.5_v1/[[BR]] 
     107/home/users/taudoux/documents/Results/B/Simuls_ratio0.5_v2/ [[BR]] 
     108/home/users/taudoux/documents/Results/B/Simuls_ratio0.8_v1/ [[BR]] 
     109/home/users/taudoux/documents/Results/B/Simuls_ratio0.8_v2/ [[BR]] 
     110/home/users/taudoux/documents/Results/B/Simuls_ratio_init/ [[BR]] 
     111/home/users/taudoux/documents/Results/B/RundBis/ 
     112 
     113 
     114 
    86115 
    87116Code was changed : separation of nitrogen input. Initially fertilizer, bnf and atmospheric deposition were in the same loop.  
     
    141170. The second (B.) consisted to normalized the fluxes by the amount of fertiliser added ; [[BR]] 
    142171. The third (C.) consisted to normalized the fluxes by the amount of all nitrogen inputs. [[BR]] 
    143  
    144  
    145 Script python used : [[BR]] 
    146 . barplot.py using files.txt in which values of nitrogen emissions and leaching are listed in order to plot EF in function of fertilizer added. [[BR]] 
    147 . barplot_c.py using files.txt in which values of nitrogen emissions and leaching are listed in order to plot EF in function of nitrogen inputs (fertilizer, NHx and NOx depositions). [[BR]] 
    148 . files.txt are in the folders : [[BR]] 
    149 /home/users/taudoux/documents/Results/B/Simuls_ratio0.5_v1/EF [[BR]] 
    150 /home/users/taudoux/documents/Results/B/Simuls_ratio0.5_v1/EFb [[BR]] 
    151 /home/users/taudoux/documents/Results/B/Simuls_ratio0.8_v1/EF [[BR]] 
    152 /home/users/taudoux/documents/Results/B/Simuls_ratio0.8_v1/EFb 
     172. The third (C.) consisted to normalized the fluxes by all nitrogen inputs. 
    153173 
    154174A. EF of each parameters on each site for NH4%=0.5 and 0.8 [[BR]] 
    155175path : /home/users/taudoux/documents/Results/B/EF/ 
    156  
     176Python : barplot.py using files.txt in which values of nitrogen emissions and leaching are listed in order to plot EF in function of fertilizer added. arg1= Name of the site [[BR]] 
    157177The first (A.) consisted to normalized the fluxes minus emissions when nfert=0 by the amount of fertiliser added. 
    158178 
     
    238258 
    239259path : /home/users/taudoux/documents/Results/B/EFc/ 
    240  
    241 The third (C.) consisted to normalized the fluxes by all nitrogen inputs. 
     260Python : barplot_c.py using files.txt in which values of nitrogen emissions and leaching are listed in order to plot EF in function of nitrogen inputs (fertilizer, NHx and NOx depositions). arg1= Name of the site[[BR]] 
     261 
    242262 
    243263||||= DK-Ris =|| 
     
    279299The sum of those EF is not equal to 1 for each site. So we looked at the difference between nitrogen pools (soil_nh4, soil_n2, soil_n2o, soil_no3 and soil_nox) at day1 and nitrogen pools on the last day.  
    280300 
    281 path : /home/users/taudoux/documents/Results/B/EFd/ 
    282  
     301path : /home/users/taudoux/documents/Results/B/EFd/ [[BR]] 
     302python : barplot_d.py arg1=name of the site 
    283303Difference between nitrogen pools normalized by the amount of all nitrogen inputs (fertilizer, NHx and NOx deposition). [[BR]] 
    284304 
     
    323343Loop at inorganic nitrogen scale : we stacked in barplot all the nitrogen outputs (N2, N2O, NH3, NOX emissions ; NH4 and NO3 Leaching; NH4 and NO3 Plant Uptake ; difference between nitrogen pools shown previously) divided by all nitrogen inputs (NOX and NHX deposition; Mineralisation but also Fertilizer) [[BR]] 
    324344path : /home/users/taudoux/documents/Results/B/EFe/ [[BR]]  
     345python : barplot_e.py arg1=name of the site [[BR]] 
    325346 
    326347Values in abscissa are all the nitrogen inputs but the fertilization is for each point : 0 , 54.8 , 82.2 , 109.6 , 140 , 164.4 , 191.8 , 219.2 , 365.25 kg/ha/y. 
     
    362383 
    363384Loop at ORCHIDEE scale : we stacked in barplot all the nitrogen outputs (N2, N2O, NH3, NOX emissions ; NH4 and NO3 Leaching; ; difference between nitrogen pools shown previously + difference between SOIL_ACTIVE_N , SOIL_SURF_N , SOIL_SLOW_N , SOIL_PASSIVE_N , TOTAL_M_N LITTER_STR_AB_N , LITTER_MET_AB_N , LITTER_STR_BE_N , LITTER_MET_BE_N , LITTER_WOD_AB_N , LITTER_WOD_BE_N); divided by all nitrogen inputs (NOX and NHX deposition and Fertilizer) [[BR]] 
     385path : /home/users/taudoux/documents/Results/B/EFebis/ [[BR]] 
     386barplot_ebis.py arg1=name of the site 
    364387 
    365388||||= DK-Ris =|| 
     
    400423We can see that there is a problem with NO3 stocks that are continually increasing or leached. It must be a problem with the denitrification. So we looked at the parameterization and saw that a variable called "anvf" takes an important role about denitrifiers. In order to know its effect, we ran a simulation in which this variable was deleted. [[BR]] 
    401424Here are the results : [[BR]] 
    402 path : /home/users/taudoux/documents/Results/B/EFi/ 
     425path : /home/users/taudoux/documents/Results/B/EFi/ [[BR]] 
     426python : barplot_init.py arg1=name of the site 
    403427 
    404428||= FR-Gri =|| 
     
    421445 
    422446 
    423 Code was whanged : run_off, the variable "harvest_above_n" were added and ''the denitrifiers activity was changed''. [[BR]]  
     447Code was whanged : run_off, the variable "harvest_above_n" were added and ''the denitrifiers activity was changed'' (UPPER BOUND). [[BR]]  
     448 
     449path of files.txt : /home/users/taudoux/documents/Results/B/Simul 
    424450We looked again at the differents loops. [[BR]]  
    425451 
     
    427453 
    428454path : /home/users/taudoux/documents/Results/B/EFf/  
     455python : barplot_f.py arg1=name of the site 
    429456 
    430457Values in abscissa are all the nitrogen inputs but the fertilization is for each point : 0 , 54.8 , 82.2 , 109.6 , 140 , 164.4 , 191.8 , 219.2 , 365.25 kg/ha/y. 
     
    465492 
    466493Loop at ORCHIDEE scale : we stacked in barplot all the nitrogen outputs (N2, N2O, NH3, NOX emissions ; NH4 and NO3 Leaching; ; difference between nitrogen pools shown previously + difference between SOIL_ACTIVE_N , SOIL_SURF_N , SOIL_SLOW_N , SOIL_PASSIVE_N , TOTAL_M_N LITTER_STR_AB_N , LITTER_MET_AB_N , LITTER_STR_BE_N , LITTER_MET_BE_N , LITTER_WOD_AB_N , LITTER_WOD_BE_N) and the variable harvest_above_n; divided by all nitrogen inputs (NOX and NHX deposition and Fertilizer) [[BR]] 
    467 path : /home/users/taudoux/documents/Results/B/EFg/ 
     494path : /home/users/taudoux/documents/Results/B/EFg/[[BR]] 
     495python : barplot_g.py arg1=name of the site 
    468496||= DK-Ris =|| 
    469497||= 0.5 =|| 
     
    501529 
    502530 
    503  
     531Code was changed : ''the denitrifiers activity was changed'' (LOWER BOUND) 
     532path of files.txt : /home/users/taudoux/documents/Results/B/RundBis [[BR]] 
     533 
     534 
     535We looked again at the differents loops. [[BR]]  
     536 
     537Loop at inorganic nitrogen scale : we stacked in barplot all the nitrogen outputs (N2, N2O, NH3, NOX emissions ; NH4 and NO3 Leaching; NH4 and NO3 Plant Uptake ; difference between nitrogen pools shown previously) divided by all nitrogen inputs (NOX and NHX deposition; Mineralisation but also Fertilizer) [[BR]] 
     538 
     539path : /home/users/taudoux/documents/Results/B/EFfbis/ [[BR]] 
     540python : barplot_fbis.py arg1=name of the site 
     541 
     542Values in abscissa are all the nitrogen inputs but the fertilization is for each point : 0 , 54.8 , 82.2 , 109.6 , 140 , 164.4 , 191.8 , 219.2 , 365.25 kg/ha/y. 
     543 
     544||= DK-Ris =|| 
     545||= 0.5 =|| 
     546||= [[Image(DK-Ris_0.5fbis.png, 50%)]] =|| 
     547 
     548||= FR-Gri =|| 
     549||= 0.5 =|| 
     550||= [[Image(FR-Gri_0.5fbis.png, 50%)]] =|| 
     551 
     552||= NL-Lan =|| 
     553||= 0.5 =|| 
     554||= [[Image(NL-Lan_0.5fbis.png, 50%)]] =|| 
     555 
     556||= NL-Lut =|| 
     557||= 0.5 =|| 
     558||= [[Image(NL-Lut_0.5fbis.png, 50%)]] =|| 
     559 
     560||= US-Bo1 =|| 
     561||= 0.5 =|| 
     562||= [[Image(US-Bo1_0.5fbis.png, 50%)]] =|| 
     563 
     564||= US-Ne1 =|| 
     565||= 0.5 =|| 
     566||= [[Image(US-Ne1_0.5fbis.png, 50%)]] =|| 
     567 
     568 
     569||= US-Ne2 =|| 
     570||= 0.5 =|| 
     571||= [[Image(US-Ne2_0.5fbis.png, 50%)]] =|| 
     572 
     573 
     574||= US-Ne3 =|| 
     575||= 0.5 =|| 
     576||= [[Image(US-Ne3_0.5fbis.png, 50%)]] =|| 
     577 
     578Loop at ORCHIDEE scale : we stacked in barplot all the nitrogen outputs (N2, N2O, NH3, NOX emissions ; NH4 and NO3 Leaching; ; difference between nitrogen pools shown previously + difference between SOIL_ACTIVE_N , SOIL_SURF_N , SOIL_SLOW_N , SOIL_PASSIVE_N , TOTAL_M_N LITTER_STR_AB_N , LITTER_MET_AB_N , LITTER_STR_BE_N , LITTER_MET_BE_N , LITTER_WOD_AB_N , LITTER_WOD_BE_N) and the variable harvest_above_n; divided by all nitrogen inputs (NOX and NHX deposition and Fertilizer) [[BR]] 
     579 
     580path : /home/users/taudoux/documents/Results/B/EFgbis/ [[BR]] 
     581python : barplot_gbis.py arg1=name of the site 
     582||= DK-Ris =|| 
     583||= 0.5 =|| 
     584||= [[Image(DK-Ris_0.5gbis.png, 50%)]] =| 
     585 
     586||= FR-Gri =|| 
     587||= 0.5 =|| 
     588||= [[Image(FR-Gri_0.5gbis.png, 50%)]] =|| 
     589 
     590||= NL-Lan =|| 
     591||= 0.5 =|| 
     592||= [[Image(NL-Lan_0.5gbis.png, 50%)]] =|| 
     593 
     594||= NL-Lut =|| 
     595||= 0.5 =|| 
     596||= [[Image(NL-Lut_0.5gbis.png, 50%)]] =|| 
     597 
     598||= US-Bo1 =|| 
     599||= 0.5 =|| 
     600||= [[Image(US-Bo1_0.5gbis.png, 50%)]] =|| 
     601 
     602||= US-Ne1 =|| 
     603||= 0.5 =|| 
     604||= [[Image(US-Ne1_0.5gbis.png, 50%)]] =|| 
     605 
     606 
     607||= US-Ne2 =|| 
     608||= 0.5 =|| 
     609||= [[Image(US-Ne2_0.5gbis.png, 50%)]] =|| 
     610 
     611 
     612||= US-Ne3 =|| 
     613||= 0.5 =|| 
     614||= [[Image(US-Ne3_0.5gbis.png, 50%)]] =|| 
     615 
     616Those results seem to suit good compared to others. So, in order to compare values of the emission with data we can find in the literature we calculate the emission factors of different fluxes. 
     617 
     618path : /home/users/taudoux/documents/Results/B/EFz/ [[BR]] 
     619python : barplot_EF.py arg1=name of the site 
     620 
     621||= DK-Ris =|| 
     622||= 0.5 =|| 
     623||= [[Image(DK-Ris_0.5_EF.png, 50%)]] =|| 
     624 
     625||= FR-Gri =|| 
     626||= 0.5 =|| 
     627||= [[Image(FR-Gri_0.5_EF.png, 50%)]] =|| 
     628 
     629||= NL-Lan =|| 
     630||= 0.5 =|| 
     631||= [[Image(NL-Lan_0.5_EF.png, 50%)]] =|| 
     632 
     633||= NL-Lut =|| 
     634||= 0.5 =|| 
     635||= [[Image(NL-Lut_0.5_EF.png, 50%)]] =|| 
     636 
     637||= US-Bo1 =|| 
     638||= 0.5 =|| 
     639||= [[Image(US-Bo1_0.5_EF.png, 50%)]] =|| 
     640 
     641||= US-Ne1 =|| 
     642||= 0.5 =|| 
     643||= [[Image(US-Ne1_0.5_EF.png, 50%)]] =|| 
     644 
     645 
     646||= US-Ne2 =|| 
     647||= 0.5 =|| 
     648||= [[Image(US-Ne2_0.5_EF.png, 50%)]] =|| 
     649 
     650 
     651||= US-Ne3 =|| 
     652||= 0.5 =|| 
     653||= [[Image(US-Ne3_0.5_EF.png, 50%)]] =||