From bf99a0b6ba0901411875ce566c8a69d2b507c180 Mon Sep 17 00:00:00 2001 From: Ryan Knox Date: Tue, 13 Sep 2022 10:07:37 -0600 Subject: [PATCH] use lai instead of treelai for b4b --- biogeochem/EDCanopyStructureMod.F90 | 9 ++++++--- 1 file changed, 6 insertions(+), 3 deletions(-) diff --git a/biogeochem/EDCanopyStructureMod.F90 b/biogeochem/EDCanopyStructureMod.F90 index b3175202ba..3305d2d857 100644 --- a/biogeochem/EDCanopyStructureMod.F90 +++ b/biogeochem/EDCanopyStructureMod.F90 @@ -1499,7 +1499,8 @@ subroutine leaf_area_profile( currentSite ) real(r8) :: dh ! vertical detph of height class (m) real(r8) :: min_chite ! bottom of cohort canopy (m) real(r8) :: max_chite ! top of cohort canopy (m) - real(r8) :: lai ! summed lai for checking m2 m-2 + real(r8) :: lai ! we use this to preserve b4b right now (m2/m2) + real(r8) :: sai ! "" !---------------------------------------------------------------------- @@ -1552,9 +1553,11 @@ subroutine leaf_area_profile( currentSite ) ! How much of each tree is stem area index? Assuming that there is ! This may indeed be zero if there is a sensecent grass ! ---------------------------------------------------------------- - + lai = currentCohort%treelai * currentCohort%c_area/currentPatch%total_canopy_area + sai = currentCohort%treesai * currentCohort%c_area/currentPatch%total_canopy_area if( (currentCohort%treelai+currentCohort%treesai) > nearzero)then - fleaf = currentCohort%treelai / (currentCohort%treelai + currentCohort%treesai) + !fleaf = currentCohort%treelai / (currentCohort%treelai + currentCohort%treesai) + fleaf = lai / (lai+sai) else fleaf = 0._r8 endif