Skip to content

Add 95th height history variable #40

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.

Sign up for GitHub

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

Open
wants to merge 2 commits into
base: noresm
Choose a base branch
from
Open

Add 95th height history variable #40

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
d43f1c6
add 95th pc height history variable
Oct 2, 2025
f9d41e1
fix typo in 95th height history var
Oct 9, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
29 changes: 29 additions & 0 deletions main/FatesHistoryInterfaceMod.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -306,6 +306,7 @@ module FatesHistoryInterfaceMod
integer :: ih_seedling_pool_si ! carbon only
integer :: ih_ba_weighted_height_si
integer :: ih_ca_weighted_height_si
integer :: ih_patch_weighted_95thpctile_height_si
integer :: ih_seeds_in_local_elem
integer :: ih_seeds_in_extern_elem
integer :: ih_seed_decay_elem
Expand Down Expand Up @@ -2413,6 +2414,7 @@ subroutine update_history_dyn_sitelevel(this,nc,nsites,sites,bc_in)
real(r8) :: leaf_herbivory ! mass of leaves eaten by herbivores [kg/yr]
real(r8) :: n_perm2 ! abundance per m2
real(r8) :: patch_fracarea ! Fraction of area for this patch
real(r8) :: crown_area_covered ! accumulator variable for patch crown area

associate( hio_npatches_si => this%hvars(ih_npatches_si)%r81d, &
hio_ncohorts_si => this%hvars(ih_ncohorts_si)%r81d, &
Expand All @@ -2424,6 +2426,7 @@ subroutine update_history_dyn_sitelevel(this,nc,nsites,sites,bc_in)
hio_fates_fraction_si => this%hvars(ih_fates_fraction_si)%r81d, &
hio_ba_weighted_height_si => this%hvars(ih_ba_weighted_height_si)%r81d, &
hio_ca_weighted_height_si => this%hvars(ih_ca_weighted_height_si)%r81d, &
hio_patch_weighted_95thpctile_height_si => this%hvars(ih_patch_weighted_95thpctile_height_si)%r81d, &
hio_canopy_spread_si => this%hvars(ih_canopy_spread_si)%r81d, &
hio_nesterov_fire_danger_si => this%hvars(ih_nesterov_fire_danger_si)%r81d, &
hio_fire_nignitions_si => this%hvars(ih_fire_nignitions_si)%r81d, &
Expand Down Expand Up @@ -2954,6 +2957,26 @@ subroutine update_history_dyn_sitelevel(this,nc,nsites,sites,bc_in)
ccohort => ccohort%taller
enddo cohortloop ! cohort loop


! mean 95th percentile height. loop through cohorts on patch again, this time from tallest to shortest
crown_area_covered = 0._r8
ccohort => cpatch%tallest
do while(associated(ccohort))
if (ccohort%canopy_layer .eq. 1) then ! ignore anything not in the canopy
if (((ccohort%c_area + crown_area_covered)/cpatch%area) .ge. 0.05_r8 ) then
hio_patch_weighted_95thpctile_height_si(io_si) = &
hio_patch_weighted_95thpctile_height_si(io_si) + ccohort%height * cpatch%area/AREA
ccohort => null() ! exit the cohort loop
else
crown_area_covered = crown_area_covered + ccohort%c_area
ccohort => ccohort%shorter
endif
else
ccohort => ccohort%shorter
endif
enddo ! cohort loop


cpatch => cpatch%younger
end do patchloop !patch loop

Expand Down Expand Up @@ -6304,6 +6327,12 @@ subroutine define_history_vars(this, initialize_variables)
upfreq=group_dyna_simple, ivar=ivar, initialize=initialize_variables, &
index=ih_ca_weighted_height_si)

call this%set_history_var(vname='FATES_MEAN_95PCTILE_HEIGHT', units='m', &
long='The patch-area-weighted mean of 95th percentile height of canopy plants within a given patch', use_default='active', &
avgflag='A', vtype=site_r8, hlms='CLM:ALM', &
upfreq=group_dyna_simple, ivar=ivar, initialize=initialize_variables, &
index=ih_patch_weighted_95thpctile_height_si)

call this%set_history_var(vname='FATES_COLD_STATUS', units='', &
long='site-level cold status, 0=not cold-dec, 1=too cold for leaves, 2=not too cold', &
use_default='active', avgflag='A', vtype=site_r8, hlms='CLM:ALM', &
Expand Down