Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Adds a lower limit to stem allocation to avoid negative states #6589

Merged
merged 1 commit into from
Nov 1, 2024

Conversation

thorntonpe
Copy link
Contributor

@thorntonpe thorntonpe commented Sep 5, 2024

The stem:leaf allocation variable can go negative without this limiter, which
can cause negative stem carbon pools.

Fixes #6591
[non-BFB]

The stem:leaf allocation variable can go negative without this
limiter, which can cause negative stem carbon pools.
@thorntonpe thorntonpe added the non-BFB PR makes roundoff changes to answers. label Sep 5, 2024
@rljacob
Copy link
Member

rljacob commented Sep 5, 2024

Can you say more about why this is non-BFB ?

@thorntonpe
Copy link
Contributor Author

A condition can occur where the climate in a given gridcell is not supporting the growth of a woody vegetation type, and the carbon stock for that type as well as its productivity is declining over time. In some very rare instances this can lead to a negative value for a calculated parameter controlling how much new photosynthesis is put toward production of new wood. In even more rare circumstances this can lead to a calculation of unrealistically tall vegetation height, which can cause unrealistically high surface temperatures, which can trigger a model stop if the temperature exceeds a critical threshold. With this bug fix in place the calculated vegetation parameter for allocation of new wood is not allowed to go negative, and so there will be non-BFB results in a small number of gridcells for a small number of timesteps. The bug-fix will also prevent the occurrence of the even more rare canopy height and temperature threshold problems, and so runs that previously stopped will now continue.

@rljacob rljacob added this to the v3.1beta milestone Sep 5, 2024
@rljacob
Copy link
Member

rljacob commented Sep 5, 2024

Can you open an issue for the bug that this fixes?

@thorntonpe
Copy link
Contributor Author

Issue created (#6591 ) and PR description updated.

@thorntonpe
Copy link
Contributor Author

I see that this is connected to the v3.1beta milestone. Can it also be added to the v3.0.1 milestone, for integration with the V3 HES simulations?

@rljacob
Copy link
Member

rljacob commented Sep 23, 2024

v3 HES is going to use 3.1. It changes answers compared to 3.0.0 for watercycle cases so can't go in a 3.0.x tag.

@rljacob
Copy link
Member

rljacob commented Oct 31, 2024

@bishtgautam is this ready?

@bishtgautam
Copy link
Contributor

Yes, this is ready to be merged.

@rljacob
Copy link
Member

rljacob commented Oct 31, 2024

Go ahead and merge this to next today.

bishtgautam added a commit that referenced this pull request Oct 31, 2024
The stem:leaf allocation variable can go negative without this limiter, which
can cause negative stem carbon pools.

Fixes #6591
[non-BFB]
@bishtgautam bishtgautam merged commit 67ae4d5 into master Nov 1, 2024
20 checks passed
@bishtgautam bishtgautam deleted the thorntonpe/lnd/negc_debug branch November 1, 2024 19:16
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug fix PR Land non-BFB PR makes roundoff changes to answers.
Projects
None yet
Development

Successfully merging this pull request may close these issues.

Model stops with land temperature > 400K due to negative carbon stock
5 participants