Changes between Version 40 and Version 41 of DevelopmentActivities/ORCHIDEE-CN/NH3


Ignore:
Timestamp:
2017-07-12T10:32:52+02:00 (7 years ago)
Author:
taudoux
Comment:

--

Legend:

Unmodified
Added
Removed
Modified

  • DevelopmentActivities/ORCHIDEE-CN/NH3

    v40 v41  
    44 
    55 
    6 This 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.[[BR]] 
     6This has been done with rev 4164. My work is divided in 4 parts, the first part corresponds to the work done for the internship report. Then, the second part is a brief bibliography about emission factors. Then, the third part corresponds to the second part of my internship: analyzing ammonia emission factor. Finally, the last part corresponds to a brief look at a global simulation of the model made by P.Messina.[[BR]] 
    77 
    88 
    … …  
    3434 * Emissions (NH3) and Depositions (NHX) [[BR]] 
    3535 
    36 >''' Detaille ce qu'il y a dans ces 3 différents répertoires''' 
    37 path : /home/users/taudoux/documents/Results/Parameters/ [[BR]] 
    38 path : /home/users/taudoux/documents/Results/Pools/ [[BR]] 
    39 path : /home/users/taudoux/documents/Results/Emission-Deposition/  
     36path : /home/users/taudoux/documents/Results/Parameters/ (Temperature, Humidity and pH)[[BR]] 
     37path : /home/users/taudoux/documents/Results/Pools/ (Nitrogen Pools)[[BR]] 
     38path : /home/users/taudoux/documents/Results/Emission-Deposition/ (NH3 emission and NHX deposition) 
    4039 
    4140==== Files ==== 
    42 >''' Où se trouvent les scripts que tu mentionnes'''  
     41All the script are in the folders /home/users/taudoux/documents/Results/Scripts 
    4342* 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 :  
    4443 * /home/users/taudoux/documents/Results/Parameters/T2M/SITES/ 
    … …  
    6968||  [[Image(nh3_em_nh4_soil.png, 50%)]]  ||  [[Image(nh3_em_n2o_soil.png, 50%)]]  ||  [[Image(nh3_em_nox_soil.png, 50%)]]  ||  [[Image(nh3_em_no3_soil.png, 50%)]]  || 
    7069 
    71 > '''Where is bar_plot_donneeS.py''' 
    7270Figure, created using "bar_plot_donnees.py", in which for each site are listed the average emission of NH3 for both model. 
    7371 
    … …  
    105103=== Where to find the files === 
    106104 
    107 * Scripts used to generate the files.txt are in the folder :  
     105 
     106* All the scripts used to generate the files.txt and figures are in the folder :  
    108107 * /home/users/taudoux/documents/Results/B/Script  
    109108 * For the script : script_delta_X.ksh  ( X depends on the simulation we are analyzing ) 
    … …  
    111110  * arg2= PFT  
    112111  * arg3= Last Time Step  
    113 > ''' Donne un exemple de la manière dont tu appelles ce script surtout pour le last time step''' 
     112Example : ./script_delta_0.5e.ksh FR-Gri 12 730  
     113 
     114* Scripts python used to generate graphics are named :  
     115 * barplot_X.py ( X depends on the simulation we are analyzing) 
    114116 
    115117* Others scripts don't need arguments and generate files.txt according to their appellation :  
    116118 * 0.5 --> Simuls_ratio0.5_v1/ ; 0.8 --> Simuls_ratio0.8_v1/   
    117  *  e --> EFe/ ; f --> EFf ; g --> EFg etc. 
    118 > ''' Ca correspond à quoi EFe , EFf EFg ??? Mentionne que ratio0.5 et ratio0.8 correspondent aux apports d'engrais et au rapport NH4+/NO3- c'est ça ?''' 
    119 * Scripts python used to generate graphics are in the folder :  
    120  * /home/users/taudoux/documents/Results/B/barplot_X.py ( X depends on the simulation we are analyzing) 
     119 *  x --> EFx/ ;  (x = b to g ; those letters correspond to the simulation I called "x") 
     120 
     121 
     122 
     123* Runs were launched in the folders /home/orchidee03/taudoux/ 
    121124 
    122125* files.txt are in the folders :  
    123126 * /home/users/taudoux/documents/Results/B/Simuls_ratio0.5_v1/ 
    124127 * /home/users/taudoux/documents/Results/B/Simuls_ratio0.5_v2/  
    125  * /home/users/taudoux/documents/Results/B/Simuls_ratio0.8_v1/  
    126  * /home/users/taudoux/documents/Results/B/Simuls_ratio0.8_v2/  
     128 * /home/users/taudoux/documents/Results/B/Simuls_ratio0.8_v1/ 
    127129 * /home/users/taudoux/documents/Results/B/Simuls_ratio_init/  
    128130 * /home/users/taudoux/documents/Results/B/RundBis/ 
    129  
    130 > '''Annonce brièvement ici ce que sont les deux versions modifiées avec lesquelles tu as travaillé, V1 et V2.''' 
    131 > ''' Plutot que tes deux sections 'code was changed', renomme les de manière plus explicite.''' 
    132 > Code was changed (V1) -> Sensitivity of NH3 emissions to the intensity of fertilization 
    133 > Code was changed (V1) -> Evaluation of a new model version for denitrification 
    134  
    135 === Code was changed (v1) ===  
     131V1/ -> Sensitivity of NH3 emissions to the intensity of fertilization [[BR]]  
     132V2/ and RundBis/ -> Evaluation of a new model version for denitrification v2 corresponds to the upper bound and RundBis to the lower bound. [[BR]] 
     133Init -> 
     134 
     135=== Sensitivity of NH3 emissions to the intensity of fertilization  ===  
    136136'''Separation of nitrogen input. Initially fertilizer, bnf and atmospheric deposition were in the same loop ''' 
    137137 
    … …  
    190190 
    191191 
    192 We 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]] 
    193 * The first (A.) consisted to normalized the fluxes minus emissions when nfert=0 by the amount of fertilizer added ; [[BR]]  
    194 * The second (B.) consisted to normalized the fluxes by the amount of fertilizer added  
    195 * The third (C.) consisted to normalized the fluxes by the amount of all nitrogen inputs 
     192We calculated emission factor (EF) of different fluxes of nitrogen (NH3, N2, N2O and NOx emission, but also NH4 and NO3 leaching) with different methods :  
     193* A. consisted to normalized the fluxes minus emissions when nfert=0 by the amount of fertilizer added   
     194* B. consisted to normalized the fluxes by the amount of fertilizer added  
     195* C. consisted to normalized the fluxes by the amount of all nitrogen inputs 
     196* D. Evaluation of the model v1  
     197* E. Evaluation of a new model version for denitrification v2 corresponds to the upper bound  
     198* F. Evaluation of a new model version for denitrification v2 corresponds to the lower bound  
    196199 
    197200 
    … …  
    325328python : barplot_d.py arg1=name of the site [[BR]] 
    326329 
    327 '''Difference between nitrogen pools normalized by the amount of all nitrogen inputs (fertilizer, NHx and NOx deposition).''' [[BR]] 
    328  
     330==== D. Difference between nitrogen pools normalized by the amount of all nitrogen inputs (fertilizer, NHx and NOx deposition) ==== 
    329331 
    330332||||= DK-Ris =|| 
    331333||= 0.5 =||= 0.8 =|| 
    332 ||= [[Image(DK-Ris_0.5d.png, 50%)]] =||= [[Image(DK-Ris_0.8d.png, 50%)]] =|| 
    333  
     334||= [[Image(DK-Ris_0.5d.png, 10%)]] =||= [[Image(DK-Ris_0.8d.png, 10%)]] =|| 
    334335 
    335336||||= FR-Gri =|| 
    336337||= 0.5 =||= 0.8 =|| 
    337 ||= [[Image(FR-Gri_0.5d.png, 50%)]] =||= [[Image(FR-Gri_0.8d.png, 50%)]] =|| 
     338||= [[Image(FR-Gri_0.5d.png, 10%)]] =||= [[Image(FR-Gri_0.8d.png, 10%)]] =|| 
    338339 
    339340||||= NL-Lan =|| 
    340341||= 0.5 =||= 0.8 =|| 
    341 ||= [[Image(NL-Lan_0.5d.png, 50%)]] =||= [[Image(NL-Lan_0.8d.png, 50%)]] =|| 
     342||= [[Image(NL-Lan_0.5d.png, 10%)]] =||= [[Image(NL-Lan_0.8d.png, 10%)]] =|| 
    342343 
    343344||||= NL-Lut =|| 
    344345||= 0.5 =||= 0.8 =|| 
    345 ||= [[Image(NL-Lut_0.5d.png, 50%)]] =||= [[Image(NL-Lut_0.8d.png, 50%)]] =|| 
     346||= [[Image(NL-Lut_0.5d.png, 10%)]] =||= [[Image(NL-Lut_0.8d.png, 10%)]] =|| 
    346347 
    347348||||= US-Bo1 =|| 
    348349||= 0.5 =||= 0.8 =|| 
    349 ||= [[Image(US-Bo1_0.5d.png, 50%)]] =||= [[Image(US-Bo1_0.8d.png, 50%)]] =|| 
     350||= [[Image(US-Bo1_0.5d.png, 10%)]] =||= [[Image(US-Bo1_0.8d.png, 10%)]] =|| 
    350351 
    351352||||= US-Ne1 =|| 
    352353||= 0.5 =||= 0.8 =|| 
    353 ||= [[Image(US-Ne1_0.5d.png, 50%)]] =||= [[Image(US-Ne1_0.8d.png, 50%)]] =|| 
    354  
     354||= [[Image(US-Ne1_0.5d.png, 10%)]] =||= [[Image(US-Ne1_0.8d.png, 10%)]] =|| 
    355355 
    356356||||= US-Ne2 =|| 
    357357||= 0.5 =||= 0.8 =|| 
    358 ||= [[Image(US-Ne2_0.5d.png, 50%)]] =||= [[Image(US-Ne2_0.8d.png, 50%)]] =|| 
    359  
     358||= [[Image(US-Ne2_0.5d.png, 10%)]] =||= [[Image(US-Ne2_0.8d.png, 10%)]] =|| 
    360359 
    361360||||= US-Ne3 =|| 
    362361||= 0.5 =||= 0.8 =|| 
    363 ||= [[Image(US-Ne3_0.5d.png, 50%)]] =||= [[Image(US-Ne3_0.8d.png, 50%)]] =|| 
    364  
    365  
     362||= [[Image(US-Ne3_0.5d.png, 10%)]] =||= [[Image(US-Ne3_0.8d.png, 10%)]] =|| 
     363 
     364This is just for the record on how we arrived at taking this decision for the analysis, this is included in the next part.  
    366365The values we found seem to be a little odd... So, in order to know if they are exploitable, we looked at the closure of different model components : 
    367366 
    368 ==== Closure of the inorganic nitrogen pools ==== 
     367===== Closure of the inorganic nitrogen pools ===== 
    369368 
    370369We 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]] 
    … …  
    410409 
    411410 
    412 ==== Closure of the N mass in ORCHIDEE model ====  
     411===== Closure of the N mass in ORCHIDEE model =====  
     412 
    413413We 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]] 
    414414 
    … …  
    476476 
    477477 
    478 === Code was changed (v2) ===  
     478==== E. Evaluation of a new model version for denitrification v2 corresponds to the upper bound  ====  
    479479Run_off, the variable "harvest_above_n" were added and ''the denitrifiers activity was changed'' (UPPER BOUND). [[BR]]  
    480480 
    … …  
    482482We looked again at the differents loops. [[BR]]  
    483483 
    484 * Loop at inorganic nitrogen scale  
     484===== Closure of the inorganic nitrogen pools ===== 
    485485 
    486486path : /home/users/taudoux/documents/Results/B/EFf/  
    … …  
    489489||= DK-Ris =|| 
    490490||= 0.5 =|| 
    491 ||= [[Image(DK-Ris_0.5f.png, 50%)]] =|| 
     491||= [[Image(DK-Ris_0.5f.png, 25%)]] =|| 
    492492 
    493493||= FR-Gri =|| 
    494494||= 0.5 =|| 
    495 ||= [[Image(FR-Gri_0.5f.png, 50%)]] =|| 
     495||= [[Image(FR-Gri_0.5f.png, 25%)]] =|| 
    496496 
    497497||= NL-Lan =|| 
    498498||= 0.5 =|| 
    499 ||= [[Image(NL-Lan_0.5f.png, 50%)]] =|| 
     499||= [[Image(NL-Lan_0.5f.png, 25%)]] =|| 
    500500 
    501501||= NL-Lut =|| 
    502502||= 0.5 =|| 
    503 ||= [[Image(NL-Lut_0.5f.png, 50%)]] =|| 
     503||= [[Image(NL-Lut_0.5f.png, 25%)]] =|| 
    504504 
    505505||= US-Bo1 =|| 
    506506||= 0.5 =|| 
    507 ||= [[Image(US-Bo1_0.5f.png, 50%)]] =|| 
     507||= [[Image(US-Bo1_0.5f.png, 25%)]] =|| 
    508508 
    509509||= US-Ne1 =|| 
    510510||= 0.5 =|| 
    511 ||= [[Image(US-Ne1_0.5f.png, 50%)]] =|| 
     511||= [[Image(US-Ne1_0.5f.png, 25%)]] =|| 
    512512 
    513513 
    514514||= US-Ne2 =|| 
    515515||= 0.5 =|| 
    516 ||= [[Image(US-Ne2_0.5f.png, 50%)]] =|| 
     516||= [[Image(US-Ne2_0.5f.png, 25%)]] =|| 
    517517 
    518518 
    … …  
    521521||= [[Image(US-Ne3_0.5f.png, 50%)]] =|| 
    522522 
    523 * Loop at ORCHIDEE scale  
     523===== Closure of the N mass in ORCHIDEE model ===== 
    524524 
    525525path : /home/users/taudoux/documents/Results/B/EFg/[[BR]] 
    … …  
    528528||= DK-Ris =|| 
    529529||= 0.5 =|| 
    530 ||= [[Image(DK-Ris_0.5g.png, 50%)]] =|| 
     530||= [[Image(DK-Ris_0.5g.png, 25%)]] =|| 
    531531 
    532532||= FR-Gri =|| 
    533533||= 0.5 =|| 
    534 ||= [[Image(FR-Gri_0.5g.png, 50%)]] =|| 
     534||= [[Image(FR-Gri_0.5g.png, 25%)]] =|| 
    535535 
    536536||= NL-Lan =|| 
    537537||= 0.5 =|| 
    538 ||= [[Image(NL-Lan_0.5g.png, 50%)]] =|| 
     538||= [[Image(NL-Lan_0.5g.png, 25%)]] =|| 
    539539 
    540540||= NL-Lut =|| 
    541541||= 0.5 =|| 
    542 ||= [[Image(NL-Lut_0.5g.png, 50%)]] =|| 
     542||= [[Image(NL-Lut_0.5g.png, 25%)]] =|| 
    543543 
    544544||= US-Bo1 =|| 
    545545||= 0.5 =|| 
    546 ||= [[Image(US-Bo1_0.5g.png, 50%)]] =|| 
     546||= [[Image(US-Bo1_0.5g.png, 25%)]] =|| 
    547547 
    548548||= US-Ne1 =|| 
    549549||= 0.5 =|| 
    550 ||= [[Image(US-Ne1_0.5g.png, 50%)]] =|| 
    551  
     550||= [[Image(US-Ne1_0.5g.png, 25%)]] =|| 
    552551 
    553552||= US-Ne2 =|| 
    554553||= 0.5 =|| 
    555 ||= [[Image(US-Ne2_0.5g.png, 50%)]] =|| 
    556  
     554||= [[Image(US-Ne2_0.5g.png, 25%)]] =|| 
    557555 
    558556||= US-Ne3 =|| 
    559557||= 0.5 =|| 
    560 ||= [[Image(US-Ne3_0.5g.png, 50%)]] =|| 
    561  
    562  
    563 === Code was changed  ===  
     558||= [[Image(US-Ne3_0.5g.png, 25%)]] =|| 
     559 
     560This is just for the record, part F is more important than part E. 
     561==== F. Evaluation of a new model version for denitrification v2 corresponds to the lower bound  ====  
    564562 
    565563''the denitrifiers activity was changed'' (LOWER BOUND) 
    … …  
    567565 
    568566 
     567 
    569568We looked again at the differents loops. [[BR]]  
    570569 
    571 * Loop at inorganic nitrogen scale  
     570===== Closure of the inorganic nitrogen pools =====  
     571 
    572572path : /home/users/taudoux/documents/Results/B/EFfbis/ [[BR]] 
    573573python : barplot_fbis.py arg1=name of the site 
    … …  
    607607||= [[Image(US-Ne3_0.5fbis.png, 50%)]] =|| 
    608608 
    609 * Loop at ORCHIDEE scale  
     609===== Closure of the N mass in ORCHIDEE model ===== 
    610610 
    611611path : /home/users/taudoux/documents/Results/B/EFgbis/ [[BR]] 
    … …  
    689689---- 
    690690 
    691 === Brief look at the global simulation done by Palmira MESSINA === 
     691== Brief look at the global simulation done by Palmira MESSINA == 
    692692This has been done with rev ?.  [[BR]] 
    693693