|
# C balance calculation
|
|
# C balance calculation
|
|
C balance is checked annually in the routine `src/lpj/check_fluxes.c`. The threshold is by default set to one (on absolute value) : `if(year>startyear && fabs(balance.carbon)>1)`. By running the exact same branch (nitrogen_cotton_tillage_newoutput) and with the same configurations and input, we obtained different C balance depending on which cluster we run the model: @PIK or @LAM
|
|
C balance is checked annually in the routine `src/lpj/check_fluxes.c`. The threshold is by default set to one (on absolute value) : `if(year>startyear && fabs(balance.carbon)>1)`. By running the exact same branch (nitrogen_cotton_tillage_newoutput) and with the same configurations and input, we obtained different C balance depending on which cluster we run the model: @PIK or @LAM
|
|
|
|
|
|
- add here plot of c balance
|
|
![image](uploads/7085b61b18ab23d08513e58a59514d0b/image.png)
|
|
|
|
Plot of carbon balance in two pixels : 24240 (south Morocco) and (north Morocco) (see [Regions and grid of Morocco ](Regions and grid of Morocco) for more information about pixel and region). C balance is higher in north pixel with value up to 20 which is a lot. And in comparaison, lpjml run @pik is under 0 anytime.
|
|
The threshold was desactivated to enable the run to continue and write the output without error.
|
|
These abnormal values seem to stop at the year 2000.
|
|
|
|
|
|
# Output comparison between run @PIK and @LAM
|
|
# Output comparison between run @PIK and @LAM
|
|
## Configuration
|
|
## Configuration
|
... | @@ -40,11 +40,20 @@ all configurations were checked for both several time |
... | @@ -40,11 +40,20 @@ all configurations were checked for both several time |
|
|
|
|
|
This error was avoided by disabled the threshold in check_fluxes routine.
|
|
This error was avoided by disabled the threshold in check_fluxes routine.
|
|
|
|
|
|
|
|
## output comparison
|
|
|
|
### output from @lam and @pik pixel 24240, south Maroc, Grand Oued Noun
|
|
|
|
[lpjml5_pixel24240_output_comparison_cluster.pdf](uploads/9bbc1ccc354d42d6100e819495940717/lpjml5_pixel24240_output_comparison_cluster.pdf)th Maroc, Grand oued noun
|
|
|
|
### difference (@lam-@pik) pixel [lpjml5_pixel24240_difference_output_comparison_cluster.pdf](uploads/a2b5886db7afa4a417a949c2214c9468/lpjml5_pixel24240_difference_output_comparison_cluster.pdf)
|
|
|
|
|
|
|
|
The major difference for this pixel is the soil carbon calculation which is bigger @lam and the difference is growing over the year.
|
|
|
|
|
|
|
|
### output from @lam and @pik pixel 24240, south Maroc, Grand Oued Noun
|
|
|
|
[lpjml5_pixel24897_output_comparison_cluster.pdf](uploads/2d07dbaa9262c3f271000a4c01510782/lpjml5_pixel24897_output_comparison_cluster.pdf)
|
|
|
|
### difference (@lam-@pik) pixel 24240, south Maroc, Grand oued noun
|
|
|
|
[lpjml5_pixel24897_difference_output_comparison_cluster.pdf](uploads/038abfb8cd87d20d4b721bf72ab44aef/lpjml5_pixel24897_difference_output_comparison_cluster.pdf)
|
|
|
|
|
|
# to do
|
|
# test to do
|
|
- [ ] plot c balance
|
|
- [ ] longer landuse spinup
|
|
- [ ] what about nitrogen ?
|
|
- [ ] 2 others pixels
|
|
- [ ] what about spinup
|
|
|
|
- [ ] add script in repository
|
|
|
|
|
|
|
|
|
|
|