Changes between Version 28 and Version 29 of DevelopmentActivities/ORCHIDEE-CN/NH3
- Timestamp:
- 2017-07-10T13:30:41+02:00 (7 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
DevelopmentActivities/ORCHIDEE-CN/NH3
v28 v29 84 84 85 85 === 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]] 90 For the script : script_delta_X.ksh [[BR]] 91 X depends on the simulation we are analyzing.[[BR]] 92 arg1= Name of the site [[BR]] 93 arg2= PFT [[BR]] 94 arg3= Last Time Step [[BR]] 95 96 Others 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 103 X 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 86 115 87 116 Code was changed : separation of nitrogen input. Initially fertilizer, bnf and atmospheric deposition were in the same loop. … … 141 170 . The second (B.) consisted to normalized the fluxes by the amount of fertiliser added ; [[BR]] 142 171 . 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. 153 173 154 174 A. EF of each parameters on each site for NH4%=0.5 and 0.8 [[BR]] 155 175 path : /home/users/taudoux/documents/Results/B/EF/ 156 176 Python : 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]] 157 177 The first (A.) consisted to normalized the fluxes minus emissions when nfert=0 by the amount of fertiliser added. 158 178 … … 238 258 239 259 path : /home/users/taudoux/documents/Results/B/EFc/ 240 241 The third (C.) consisted to normalized the fluxes by all nitrogen inputs. 260 Python : 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 242 262 243 263 ||||= DK-Ris =|| … … 279 299 The 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. 280 300 281 path : /home/users/taudoux/documents/Results/B/EFd/ 282 301 path : /home/users/taudoux/documents/Results/B/EFd/ [[BR]] 302 python : barplot_d.py arg1=name of the site 283 303 Difference between nitrogen pools normalized by the amount of all nitrogen inputs (fertilizer, NHx and NOx deposition). [[BR]] 284 304 … … 323 343 Loop 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]] 324 344 path : /home/users/taudoux/documents/Results/B/EFe/ [[BR]] 345 python : barplot_e.py arg1=name of the site [[BR]] 325 346 326 347 Values 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. … … 362 383 363 384 Loop 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]] 385 path : /home/users/taudoux/documents/Results/B/EFebis/ [[BR]] 386 barplot_ebis.py arg1=name of the site 364 387 365 388 ||||= DK-Ris =|| … … 400 423 We 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]] 401 424 Here are the results : [[BR]] 402 path : /home/users/taudoux/documents/Results/B/EFi/ 425 path : /home/users/taudoux/documents/Results/B/EFi/ [[BR]] 426 python : barplot_init.py arg1=name of the site 403 427 404 428 ||= FR-Gri =|| … … 421 445 422 446 423 Code was whanged : run_off, the variable "harvest_above_n" were added and ''the denitrifiers activity was changed''. [[BR]] 447 Code was whanged : run_off, the variable "harvest_above_n" were added and ''the denitrifiers activity was changed'' (UPPER BOUND). [[BR]] 448 449 path of files.txt : /home/users/taudoux/documents/Results/B/Simul 424 450 We looked again at the differents loops. [[BR]] 425 451 … … 427 453 428 454 path : /home/users/taudoux/documents/Results/B/EFf/ 455 python : barplot_f.py arg1=name of the site 429 456 430 457 Values 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. … … 465 492 466 493 Loop 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/ 494 path : /home/users/taudoux/documents/Results/B/EFg/[[BR]] 495 python : barplot_g.py arg1=name of the site 468 496 ||= DK-Ris =|| 469 497 ||= 0.5 =|| … … 501 529 502 530 503 531 Code was changed : ''the denitrifiers activity was changed'' (LOWER BOUND) 532 path of files.txt : /home/users/taudoux/documents/Results/B/RundBis [[BR]] 533 534 535 We looked again at the differents loops. [[BR]] 536 537 Loop 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 539 path : /home/users/taudoux/documents/Results/B/EFfbis/ [[BR]] 540 python : barplot_fbis.py arg1=name of the site 541 542 Values 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 578 Loop 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 580 path : /home/users/taudoux/documents/Results/B/EFgbis/ [[BR]] 581 python : 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 616 Those 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 618 path : /home/users/taudoux/documents/Results/B/EFz/ [[BR]] 619 python : 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%)]] =||