From 7e1bf6056a3417f78edd1541ea2aeca4302e3e93 Mon Sep 17 00:00:00 2001 From: Sam Rabin Date: Thu, 3 Oct 2024 12:18:12 -0600 Subject: [PATCH] Bugfix: Add termination mortality to _SZAP outputs. --- main/FatesHistoryInterfaceMod.F90 | 13 ++++++++++--- 1 file changed, 10 insertions(+), 3 deletions(-) diff --git a/main/FatesHistoryInterfaceMod.F90 b/main/FatesHistoryInterfaceMod.F90 index 6da14a8de1..79c56c9334 100644 --- a/main/FatesHistoryInterfaceMod.F90 +++ b/main/FatesHistoryInterfaceMod.F90 @@ -4899,20 +4899,27 @@ subroutine update_history_dyn2_ageclass(this,nc,nsites,sites,bc_in) ! treat as such iscag = i_scls - ! TODO: Why not also *cpatch%area like other mortality terms? - weight = AREA_INV - ! add imort to other mortality terms. consider imort as understory mortality even if it happens in ! cohorts that may have been promoted as part of the patch creation, and use the pre-calculated site-level ! values to avoid biasing the results by the dramatically-reduced number densities in cohorts that are subject to imort + weight = AREA_INV hio_mortality_understory_si_scag(io_si,iscag) = hio_mortality_understory_si_scag(io_si,iscag) + & sites(s)%imort_rate(i_scls, ft) * weight ! add fire mortality to other mortality terms + weight = AREA_INV hio_mortality_canopy_si_scag(io_si,iscag) = hio_mortality_canopy_si_scag(io_si,iscag) + & sites(s)%fmort_rate_canopy(i_scls, ft) * weight hio_mortality_understory_si_scag(io_si,iscag) = hio_mortality_understory_si_scag(io_si,iscag) + & sites(s)%fmort_rate_ustory(i_scls, ft) * weight + + ! add termination mortality to other mortality terms + weight = AREA_INV * days_per_year + hio_mortality_canopy_si_scag(io_si,iscag) = hio_mortality_canopy_si_scag(io_si,iscag) + & + sum(sites(s)%term_nindivs_canopy(:,i_scls,ft)) * weight + hio_mortality_understory_si_scag(io_si,iscag) = hio_mortality_understory_si_scag(io_si,iscag) + & + sum(sites(s)%term_nindivs_ustory(:,i_scls,ft)) * weight + end do ! size class loop end do ! pft loop end do siteloop