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

[Fang summer '24 PR 2.1] More improvements to fire model #2711

Merged
merged 13 commits into from
Aug 22, 2024

Conversation

samsrabin
Copy link
Collaborator

@samsrabin samsrabin commented Aug 16, 2024

Description of changes

(Supersedes #2615, which was the same as this but on the 2021 module instead of the new 2024 module.)

  • Improves peat burning dependence on climate
  • Improves dependence of deforestation burning on climate.
  • Allows lightning ignitions in tropical closed forest. (Human ignitions, other than deforestation, still not allowed.)
  • Adds effect of landscape fragmentation on ignitions and fire duration.
  • Recalibrates some parameters.
  • Adds fire_method li2024crujra.

Specific notes

Contributors other than yourself, if any: @lifang0209

CTSM Issues Fixed (include github issue #): None

Are answers expected to change (and if so in what way)? Yes, but only for fire_method li2024gswpfrc.

Any User Interface Changes (namelist or namelist defaults changes)?

  • Several namelist parameters get new default values for fire_method li2024gswpfrc.
  • Adds fire_method li2024crujra.

Does this create a need to change or add documentation? Did you do so? No

Testing performed, if any:

…ppression.

- Previously, neither lightning nor human ignitions (other than deforestation fire) were allowed there.
- Landscape fragmentation suppression effect has to do with lfwt.
Same as li2024gswpfrc except rh_hgh now 85 (was 80).
@samsrabin samsrabin added enhancement new capability or improved behavior of existing capability science Enhancement to or bug impacting science labels Aug 16, 2024
@samsrabin samsrabin self-assigned this Aug 16, 2024
Copy link
Collaborator

@ekluzek ekluzek left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This is great, thanks for working on this and getting it in quickly.

I have a question and a couple suggestions on changes to the namelist items at top. Those are the things I think you should do now. I'm marking as approve just so you don't have to wait on me to reapprove later.

There's a couple things looking at the code got me thinking about for longer term. There's two things that maybe you could do now.

The simple thing is redo code indentation because of an if statement that we removed. That's easy enough in VS code.

The only thing that might be good to do now would be to add some comments to explain what's going on in the code bits that change. If you could add something on that now -- that would be good to explain. But, if you need to get with Fang that could be done later.

bld/namelist_files/namelist_definition_ctsm.xml Outdated Show resolved Hide resolved
bld/namelist_files/namelist_definition_ctsm.xml Outdated Show resolved Hide resolved
bld/namelist_files/namelist_defaults_ctsm.xml Show resolved Hide resolved
src/biogeochem/CNFireFactoryMod.F90 Show resolved Hide resolved
src/biogeophys/Wateratm2lndBulkType.F90 Show resolved Hide resolved
src/biogeochem/CNFireLi2024Mod.F90 Outdated Show resolved Hide resolved
src/biogeochem/CNFireLi2024Mod.F90 Show resolved Hide resolved
src/biogeochem/CNFireLi2024Mod.F90 Show resolved Hide resolved
@samsrabin samsrabin changed the title [Fang summer '24 PR 2] More improvements to fire model [Fang summer '24 PR 2.1] More improvements to fire model Aug 19, 2024
Copy link
Collaborator Author

@samsrabin samsrabin left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Thanks for the review, Erik! I've made issue #2714 with my ideas for the fire refactor, and I've included the cleanup things you mentioned. Shortly I will make the changes you requested, then start aux_clm.

bld/namelist_files/namelist_defaults_ctsm.xml Show resolved Hide resolved
src/biogeochem/CNFireLi2024Mod.F90 Show resolved Hide resolved
src/biogeochem/CNFireLi2024Mod.F90 Outdated Show resolved Hide resolved
src/biogeophys/Wateratm2lndBulkType.F90 Show resolved Hide resolved
Value above which 30-day running relative humidity has no effect on fuel combustibility.
"PFT-dependent thresholds of P60d and P10d" in Li et al. (2013, doi:10.5194/bg-10-2293-2013):
- defo_fire_precip_thresh_bet: Value (mm/d) above which running mean daily precipitation (10 or 60 days) does not allow deforestation fire for a column with broadleaf evergreen tropical trees but no broadleaf deciduous tropical trees.
- defo_fire_precip_thresh_bdt: Vice versa
- nonborpeat_fire_precip_denom: Denominator of precipitation in equation relating that to non-boreal peat fire (unitless). Eq. 9 in Li et al. (2013, doi:10.5194/bg-10-2293-2013).
- borpeat_fire_soilmoist_denom: Denominator of exponential in soil moisture term of equation relating that and temperature to boreal peat fire (unitless). Eq. 10 in Li et al. (2013, doi:10.5194/bg-10-2293-2013).
@samsrabin samsrabin requested a review from ekluzek August 21, 2024 16:49
Copy link
Collaborator

@ekluzek ekluzek left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

These additions are really good. Thanks for bringing it in, and it looks like it was required for the new Li2024 work. Love to see more fire magic numbers turned into variables.

The one change I have is to add d00 to the end of the namelist parameters so that the constants are all double precision. You might not be able to replicate previous results without it. And it's an important practice to have.

But, nice work here. Thanks.

bld/namelist_files/namelist_defaults_ctsm.xml Outdated Show resolved Hide resolved
bld/namelist_files/namelist_defaults_ctsm.xml Show resolved Hide resolved
@samsrabin
Copy link
Collaborator Author

Made those changes @ekluzek. Could you re-approve? Tests are running now.

Copy link
Collaborator

@ekluzek ekluzek left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I saw one missing d00. But, approving so you don't have to wait for me.

@samsrabin samsrabin merged commit 50edd8e into ESCOMP:master Aug 22, 2024
2 checks passed
@lifang0209
Copy link

lifang0209 commented Aug 23, 2024

@samsrabin
Hi Sam,
(1)
if(( kmo == 1 .and. kda == 1 .and. mcsec == 0) .or. &
dtrotr_col(c) <=0._r8 )then
fbac1(c) = 0._r8
farea_burned(c) = baf_crop(c)+baf_peatf(c)
the farea_burned(c) = baf_crop(c)+baf_peatf(c) should be removed, so lightning ignition can occur in the tropical closed forests.

(2)
farea_burned(c) = farea_burned(c)+fb * cli * (0.01_r8/secspday) , 0.01 should be changed to cli_scale

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
enhancement new capability or improved behavior of existing capability science Enhancement to or bug impacting science
Projects
Status: Ready to eat (Done!)
Status: Done (non release/external)
Status: Done
Development

Successfully merging this pull request may close these issues.

3 participants