Changes between Version 32 and Version 33 of DevelopmentActivities/ORCHIDEE-CN/NH3


Ignore:
Timestamp:
2017-07-11T11:57:09+02:00 (7 years ago)
Author:
taudoux
Comment:

--

Legend:

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

    v32 v33  
    1 = Work of Thomas on NH3 = 
    2  
    3 == Evaluation of ammonia emissions results calculated by the model (Fertilizer =0) == 
    4 This has been done with rev 4164 
     1= Work of Thomas Audoux carried out during an internship on NH3 = 
     2This has been done with rev 4164. My work is divided in 3 parts, the first part corresponds to the work done for the internship report. Then, the second part is a brief bibliography about emission factors. Finally, the last part corresponds to the second part of my internship: analyzing ammonia emission factor. 
     3== I. Evaluation of ammonia emissions results calculated by the model ( Fertilizer =0 ) == 
     4 
     5[https://forge.ipsl.jussieu.fr/orchidee/attachment/wiki/DevelopmentActivities/ORCHIDEE-CN/NH3/Rapport%20de%20stage%20final.pdf Internship_report.pdf] 
     6 
    57=== Sites and PFT ===  
    68||= PFT =||= Sites =|| 
     
    2325 * Humidity 
    2426 * pH 
     27 * Nitrogen Pools (Soil_NH4, Soil_NO3, Soil_NOX and Soil_N2O) [[BR]] 
     28 * Emissions (NH3) and Depositions (NHX) [[BR]] 
    2529path : /home/users/taudoux/documents/Results/Parameters/ [[BR]] 
    26  
    27 Nitrogen Pools (Soil_NH4, Soil_NO3, Soil_NOX and Soil_N2O) [[BR]] 
    2830path : /home/users/taudoux/documents/Results/Pools/ [[BR]] 
    29  
    30 Emissions (NH3) and Depositions (NHX) [[BR]] 
    3131path : /home/users/taudoux/documents/Results/Emission-Deposition/  
    3232 
    33   * Files  
    34 . The script "script_calcul_parameters.ksh" used in order to create files of Temperature, Humidity and pH of each site in function of time listed in : [[BR]]/home/users/taudoux/documents/Results/Parameters/T2M/SITES/;[[BR]]/home/users/taudoux/documents/Results/Parameters/LITTERHUM/SITES/;[[BR]]/home/users/taudoux/documents/Results/Parameters/pH/SITES/.[[BR]] 
    35  
    36 . The script "script_calcul_pools.ksh" used in order to create files of soil_NH4+, Soil_NO3-, Soil_NOX and Soil_N2O for each site in function of time listed in : [[BR]] /home/users/taudoux/documents/Results/Pools/SOIL_NH4/SITES/ ; [[BR]] /home/users/taudoux/documents/Results/Pools/SOIL_NO3/SITES/ ; [[BR]] /home/users/taudoux/documents/Results/Pools/SOIL_NOX/SITES/ ;  [[BR]] /home/users/taudoux/documents/Results/Pools/SOIL_N2O/SITES/ . [[BR]] 
    37  
    38 . The Script "script_calcul_emission-deposition.ksh" used in order to create files of NH3 emission and NHX deposition for each site in function of time in : [[BR]] /home/users/taudoux/documents/Results/Emission-Deposition/NH3_EMISSION/Sites/ ; [[BR]] /home/users/taudoux/documents/Results/Emission-Deposition/NHX_DEPOSITION/Sites/.[[BR]]  
    39  
    40 For each parameter are listed :  
    41  
    42 . Figure of correlation with NH3 emissions (except pH) obtained using "script_calcul_r2.ksh" to obtain the R2 and then bar_plot.py in order to create the figure ; [[BR]] . Average emission of ammonia on each site in function of the parameter studied using "scatterplot.py" ; [[BR]] . Text documents in which are listed the different values. 
    43  
    44    * Comparison with data used by LMDz-INCA model 
     33==== Files ====  
     34* The script "script_calcul_parameters.ksh" used in order to create files of Temperature, Humidity and pH of each site in function of time listed in :  
     35 * /home/users/taudoux/documents/Results/Parameters/T2M/SITES/ 
     36 * /home/users/taudoux/documents/Results/Parameters/LITTERHUM/SITES/ 
     37 * /home/users/taudoux/documents/Results/Parameters/pH/SITES/ 
     38 
     39* The script "script_calcul_pools.ksh" used in order to create files of soil_NH4+, Soil_NO3-, Soil_NOX and Soil_N2O for each site in function of time listed in :  
     40 * /home/users/taudoux/documents/Results/Pools/SOIL_NH4/SITES/   
     41 * /home/users/taudoux/documents/Results/Pools/SOIL_NO3/SITES/   
     42 * /home/users/taudoux/documents/Results/Pools/SOIL_NOX/SITES/   
     43 * /home/users/taudoux/documents/Results/Pools/SOIL_N2O/SITES/  
     44 
     45* The Script "script_calcul_emission-deposition.ksh" used in order to create files of NH3 emission and NHX deposition for each site in function of time in :  
     46 * /home/users/taudoux/documents/Results/Emission-Deposition/NH3_EMISSION/Sites/  
     47 * /home/users/taudoux/documents/Results/Emission-Deposition/NHX_DEPOSITION/Sites/ 
     48 
     49* For each parameter are listed :  
     50 * Figure of correlation with NH3 emissions (except pH) obtained using "script_calcul_r2.ksh" to obtain the R2 and then bar_plot.py in order to create the figure  
     51 * Average emission of ammonia on each site in function of the parameter studied using "scatterplot.py"  
     52 * Text documents in which are listed the different values. 
     53 
     54   === Comparison with data used by LMDz-INCA model === 
    4555path : /home/users/taudoux/documents/Results/Comparison/  [[BR]] 
    4656 
     
    7585---- 
    7686 
    77 === Bibliography on ammonia emission factor 
     87== II. Bibliography on ammonia emission factor == 
    7888||||=  '''Bouwman et al., 1997 :''' A global high-resolution emission inventory for ammonia =|| 
    7989||= [[Image(EF_bouwman_1.png, 50%)]] =||= [[Image(EF_bouwman_2.png, 50%)]] =|| 
     
    8696---- 
    8797 
    88 === Evaluation of ammonia emissions results calculated by the model (Fertilizer added) 
    89  
    90  
    91 . Scripts used to generate the files.txt are in the folder : [[BR]]  
    92 /home/users/taudoux/documents/Results/B/Script [[BR]] 
    93 For the script : script_delta_X.ksh  [[BR]] 
    94 X depends on the simulation we are analyzing.[[BR]] 
    95 arg1= Name of the site [[BR]] 
    96 arg2= PFT [[BR]] 
    97 arg3= Last Time Step [[BR]] 
    98  
    99 Others scripts don't need arguments and generate files.txt according to their appellation : [[BR]] 
    100 . 0.5 --> Simuls_ratio0.5_v1/ ; 0.8 --> Simuls_ratio0.8_v1/ [[BR]]  
    101 . e --> EFe/ ; f --> EFf ; g --> EFg etc. 
    102  
    103 .Scripts python used to generate graphics are in the folder : [[BR]]  
    104 /home/users/taudoux/documents/Results/B/barplot_X.py [[BR]] 
    105  
    106 X depends on the simulation we are analyzing.[[BR]] 
    107  
    108 . files.txt are in the folders : [[BR]] 
    109 /home/users/taudoux/documents/Results/B/Simuls_ratio0.5_v1/[[BR]] 
    110 /home/users/taudoux/documents/Results/B/Simuls_ratio0.5_v2/ [[BR]] 
    111 /home/users/taudoux/documents/Results/B/Simuls_ratio0.8_v1/ [[BR]] 
    112 /home/users/taudoux/documents/Results/B/Simuls_ratio0.8_v2/ [[BR]] 
    113 /home/users/taudoux/documents/Results/B/Simuls_ratio_init/ [[BR]] 
    114 /home/users/taudoux/documents/Results/B/RundBis/ 
    115  
    116  
    117  
    118  
    119 Code was changed : separation of nitrogen input. Initially fertilizer, bnf and atmospheric deposition were in the same loop.  
     98== III. Evaluation of ammonia emissions results calculated by the model (Fertilizer added) == 
     99 
     100=== Where to find the files === 
     101 
     102* Scripts used to generate the files.txt are in the folder :  
     103 * /home/users/taudoux/documents/Results/B/Script  
     104 * For the script : script_delta_X.ksh  ( X depends on the simulation we are analyzing ) 
     105  * arg1= Name of the site  
     106  * arg2= PFT  
     107  * arg3= Last Time Step  
     108 
     109* Others scripts don't need arguments and generate files.txt according to their appellation :  
     110 * 0.5 --> Simuls_ratio0.5_v1/ ; 0.8 --> Simuls_ratio0.8_v1/   
     111 *  e --> EFe/ ; f --> EFf ; g --> EFg etc. 
     112 
     113* Scripts python used to generate graphics are in the folder :  
     114 * /home/users/taudoux/documents/Results/B/barplot_X.py ( X depends on the simulation we are analyzing) 
     115 
     116* files.txt are in the folders :  
     117 * /home/users/taudoux/documents/Results/B/Simuls_ratio0.5_v1/ 
     118 * /home/users/taudoux/documents/Results/B/Simuls_ratio0.5_v2/  
     119 * /home/users/taudoux/documents/Results/B/Simuls_ratio0.8_v1/  
     120 * /home/users/taudoux/documents/Results/B/Simuls_ratio0.8_v2/  
     121 * /home/users/taudoux/documents/Results/B/Simuls_ratio_init/  
     122 * /home/users/taudoux/documents/Results/B/RundBis/ 
     123 
     124 
     125 
     126 
     127=== Code was changed ===  
     128'''Separation of nitrogen input. Initially fertilizer, bnf and atmospheric deposition were in the same loop ''' 
    120129 
    121130 
     
    130139|| 13 || NL-Lan ; US-Bo1 ; US-Ne1 ; US-Ne2 ; US-Ne3 ||  
    131140 
     141/home/users/taudoux/documents/Results/B/Simul_ratio0.5_v1/ [[BR]] 
     142/home/users/taudoux/documents/Results/B/Simul_ratio0.8_v1/ 
    132143 
    133144* Parameters studied  
    134 General path : path : /home/users/taudoux/documents/Results/B/Simul_ratio0.5_v1/ or /home/users/taudoux/documents/Results/B/Simul_ratio0.8_v1/ [[BR]] 
    135 Temperature, Humidity in the folder Parameters/ [[BR]] 
    136  
    137 Nitrogen Pools (Soil_NH4, Soil_NO3, Soil_NOX and Soil_N2O) and Leaching (NH4 and NO3) in the folder Pools_fluxes/ [[BR]] 
    138  
    139 Emissions (NH3, N2, N2O, NOx) in the folder NH3_EM/ [[BR]] 
    140  
    141 GPP in the folder GPP/ [[BR]]  
    142  
    143 Nitrogen input summed (fertiliser, NHx and NOx depositions) in the folder DEP_FERT/ [[BR]] 
     145 * Temperature 
     146 * Humidity  
     147 * Nitrogen Pools (Soil_NH4, Soil_NO3, Soil_NOX and Soil_N2O) and Leaching (NH4 and NO3)  
     148 * Emissions (NH3, N2, N2O, NOx)  
     149 * GPP  
     150 * Nitrogen input summed (fertilizer, NHx and NOx depositions)  
     151 
     152Respectively in the folders Parameters (T° and Hum), Pools_fluxes, NH3_EM, DEP_FERT and GPP. 
    144153 
    145154We looked at how those parameters evolved in each site with the amount of fertilizer with both of ammonium's percentage in the fertilizer (0.5 and 0.8). [[BR]] 
    146 Script used : scatterplot_X.py in the folders [[BR]] /home/users/taudoux/documents/Results/B/Simuls_ratio0.5_v1/Graph/ [[BR]] 
    147 /home/users/taudoux/documents/Results/B/Simuls_ratio0.8_v1/Graph/ [[BR]] 
     155 
     156* Script used : scatterplot_X.py in the folders  
     157 * /home/users/taudoux/documents/Results/B/Simuls_ratio0.5_v1/Graph/  
     158 * /home/users/taudoux/documents/Results/B/Simuls_ratio0.8_v1/Graph/  
    148159X = em_n2 ; em_n2o ; em_nh3 ; em_nox ; GPP ; leaching_nh4 ; leaching_no3 ; pools_n2 ; pools_n2o ; pools_nh4 ; pools_no3 ; pools_nox. 
     160 
    149161||||= N2 emission =||||= N2O emission =||||= NH3 emission =||||= NOx emission =|| 
    150162||= 0.5 =||= 0.8 =||= 0.5 =||= 0.8 =||= 0.5 =||= 0.8 =||= 0.5 =||= 0.8 =|| 
     
    158170 
    159171 
    160 '''Legend :'''  [[BR]]  
    161 . Red = DK-Ris [[BR]] 
    162 . Yellow = FR-Gri [[BR]] 
    163 . Green = NL-Lan [[BR]] 
    164 . Blue= NL-Lut [[BR]] 
    165 . Black= US-Bo1 [[BR]] 
    166 . Gray = US-Ne1 [[BR]] 
    167 . Orange = US-Ne2 [[BR]] 
    168 . Purple = US-Ne3 [[BR]] 
     172* '''Legend :'''  
     173 * Red = DK-Ris  
     174 * Yellow = FR-Gri  
     175 * Green = NL-Lan  
     176 * Blue= NL-Lut  
     177 * Black= US-Bo1  
     178 * Gray = US-Ne1  
     179 * Orange = US-Ne2  
     180 * Purple = US-Ne3  
    169181 
    170182 
    171183We calculated emission factor (EF) of different fluxes of nitrogen (NH3, N2, N2O and NOx emission, but also NH4 and NO3 leaching) with different methods : [[BR]] 
    172 . The first (A.) consisted to normalized the fluxes minus emissions when nfert=0 by the amount of fertiliser added ; [[BR]]  
    173 . The second (B.) consisted to normalized the fluxes by the amount of fertiliser added ; [[BR]] 
    174 . The third (C.) consisted to normalized the fluxes by the amount of all nitrogen inputs. [[BR]] 
    175 . The third (C.) consisted to normalized the fluxes by all nitrogen inputs. 
    176  
    177 A. EF of each parameters on each site for NH4%=0.5 and 0.8 [[BR]] 
     184* The first (A.) consisted to normalized the fluxes minus emissions when nfert=0 by the amount of fertilizer added ; [[BR]]  
     185* The second (B.) consisted to normalized the fluxes by the amount of fertilizer added  
     186* The third (C.) consisted to normalized the fluxes by the amount of all nitrogen inputs 
     187* The third (C.) consisted to normalized the fluxes by all nitrogen inputs. 
     188 
     189 
     190=== A. EF of each parameters on each site for NH4%=0.5 and 0.8 === 
    178191path : /home/users/taudoux/documents/Results/B/EF/ 
    179192Python : 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]] 
    180 The first (A.) consisted to normalized the fluxes minus emissions when nfert=0 by the amount of fertiliser added. 
     193 
    181194 
    182195||||= DK-Ris =|| 
     
    216229 
    217230 
    218 B. EF of each parameters on each site for NH4%=0.5 and 0.8 
     231=== B. EF of each parameters on each site for NH4%=0.5 and 0.8 === 
    219232 
    220233path : /home/users/taudoux/documents/Results/B/EFb/ 
    221234 
    222 The second (B.) consisted to normalized the fluxes by the amount of fertiliser added.  
    223  
    224235||||= DK-Ris =|| 
    225236||= 0.5 =||= 0.8 =|| 
     
    258269 
    259270 
    260 C. EF of each parameters on each site for NH4%=0.5 and 0.8 
    261  
     271=== C. EF of each parameters on each site for NH4%=0.5 and 0.8 === 
    262272path : /home/users/taudoux/documents/Results/B/EFc/ [[BR]] 
    263273Python : barplot_c.py (arg1= Name of the site) 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).  
     
    303313 
    304314path : /home/users/taudoux/documents/Results/B/EFd/ [[BR]] 
    305 python : barplot_d.py arg1=name of the site 
    306 Difference between nitrogen pools normalized by the amount of all nitrogen inputs (fertilizer, NHx and NOx deposition). [[BR]] 
     315 
     316python : barplot_d.py arg1=name of the site [[BR]] 
     317 
     318'''Difference between nitrogen pools normalized by the amount of all nitrogen inputs (fertilizer, NHx and NOx deposition).''' [[BR]] 
    307319 
    308320