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

Active crown fire (Scott & Reinhardt 2001) #1

Closed
wants to merge 1,379 commits into from

Conversation

slevis-lmwg
Copy link
Owner

@slevis-lmwg slevis-lmwg commented Apr 19, 2022

Description:

This PR follows in the footsteps of
#573
#572
jkshuman#5
https://github.com/jkshuman/fates/tree/active_crown_Scott_2001
in getting the FATES fire module closer to having an active crown fire parameterization.

Collaborators:

@jkshuman @lmkueppers @pollybuotte @ckoven @glemieux

Expectation of Answer Changes:

At least when fates_fire_active_crown_fire = 1 and active crown fires do occur.

Checklist:

  • My change requires a change to the documentation.
  • I have updated the in-code documentation .AND. (the technical note .OR. the wiki) accordingly.
  • I have read the CONTRIBUTING document.
  • FATES PASS/FAIL regression tests were run
  • If answers were expected to change, evaluation was performed and provided

Test Results:

CTSM (or) E3SM (specify which) test hash-tag:

CTSM (or) E3SM (specify which) baseline hash-tag:

FATES baseline hash-tag:

Test Output:

JunyanDing and others added 30 commits July 21, 2021 13:30
These changes add new parameters to the fates parameter file. No existing parameters were changed or removed.  See issue 753 for details on new parameters.
This needs to be not equal since when in any mode other than no comp
or its derivatives, the label will be -999
…to print statements and new file for parameters at barro colorado island panama.
Refactored the do_patch_dynamics check added for SP mode as well
as some older duplicate logic checks that can be combined.
Updating the element id from a number to the element parameter name to clarify.

Co-authored-by: Charlie Koven <[email protected]>
glemieux and others added 23 commits March 10, 2022 12:06
Update FATES PFT optical to Majasalmi Bright 2019
We never added the logo that we had decided on a while ago to our GitHub repository, so here it is.
move frac_lai and frac_sai calculation outside do ib=1,hlm_numSWb loop
Resolved conflicts:
fire/SFMainMod.F90
main/EDPftvarcon.F90
main/EDTypesMod.F90
main/FatesInterfaceMod.F90
parameter_files/fates_params_default.cdl
This and next few commits resolve conflicts missed by git and caught by
github
Committing only to get ./manage_externals/checkout_externals to finish
Using existing test:
SMS_Lm12_D_Mmpi-serial.1x1_brazil.I2000Clm50FatesCruRsGs.izumi_intel.clm-FatesFireLightningPopDens
There were two copies of this parameter: a scalar and an array and the
model was getting confused when trying to read them
I found references to both currentCohort%active_crown_fire_flg and
currentPatch%active_crown_fire_flg and I think the latter is correct
and the former incorrect
Calculate required intensity to scorch crown and canopy, then re-calculate
rate of spread and fire intensity

Fixes: 573

User interface changes?: No

Code review:
Fix mistake in equation to solve for phi_wind as part of steps towards calculation of canopy fraction burnt.
To be used as on/off switch.
Rm the version that varies by pft.
@slevis-lmwg slevis-lmwg changed the base branch from read_lightning_data to active_crown_fire April 19, 2022 21:05
! TODO slevis: Are we missing a -phi_s after the -1?
ROS_torch = 1._r8 / (54.683_r8 * wind_reduce) * &
((60._r8 * passive_crown_FI * currentPatch%fuel_bulkd * eps * q_ig / &
(heat_per_area * ir * xi) - 1._r8) / (c * beta_ratio**-e))**(1._r8 / b)
Copy link
Owner Author

Choose a reason for hiding this comment

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

@jkshuman
commenting out lines 408 and 495 got us past the divide-by-zero in eq 18; however, are we missing a -phi_s term after the -1 or did you omit it intentionally?

Copy link
Contributor

Choose a reason for hiding this comment

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

that is for the slope. I meant to add it in as a placeholder with a value of zero as we do not have slope currently. So you can add a parameter "slope_factor" and set it to zero.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

10 participants