From eaf440cc4bab13bb29d12f604fae54994468c4b5 Mon Sep 17 00:00:00 2001 From: Rich Fiorella Date: Sun, 31 Mar 2024 15:34:13 -0600 Subject: [PATCH] Update VegStructUpdateMod to get right array sizes --- components/elm/src/biogeochem/VegStructUpdateMod.F90 | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/components/elm/src/biogeochem/VegStructUpdateMod.F90 b/components/elm/src/biogeochem/VegStructUpdateMod.F90 index e2abeb4ef4d7..9b73344fd540 100644 --- a/components/elm/src/biogeochem/VegStructUpdateMod.F90 +++ b/components/elm/src/biogeochem/VegStructUpdateMod.F90 @@ -164,11 +164,11 @@ subroutine VegStructUpdate(num_soilp, filter_soilp, & ! taper and stocking density can be set as input variables now to ! change from default values set in pftvarcon.F90 if (spinup_state >= 1) then - htop(p) = ((3._r8 * deadstemc(p) * spinup_mortality_factor * taper * taper)/ & - (SHR_CONST_PI * stocking * dwood(ivt(p))))**(1._r8/3._r8) + htop(p) = ((3._r8 * deadstemc(p) * spinup_mortality_factor * taper(p) * taper(p))/ & + (SHR_CONST_PI * stocking(p) * dwood(ivt(p))))**(1._r8/3._r8) else - htop(p) = ((3._r8 * deadstemc(p) * taper * taper)/ & - (SHR_CONST_PI * stocking * dwood(ivt(p))))**(1._r8/3._r8) + htop(p) = ((3._r8 * deadstemc(p) * taper(p) * taper(p))/ & + (SHR_CONST_PI * stocking(p) * dwood(ivt(p))))**(1._r8/3._r8) end if ! Peter Thornton, 5/3/2004 @@ -240,7 +240,7 @@ subroutine VegStructUpdate(num_soilp, filter_soilp, & ! Liston and Hiemstra, 2011; Belke-Brea et al. 2020 if (ivt(p) > noveg .and. ivt(p) <= nbrdlf_dcd_brl_shrub ) then ol = min( max(snow_depth(c)-hbot(p), 0._r8), htop(p)-hbot(p)) - fb = 1._r8 - (ol / max(1.e-06_r8, bendresist(p)*(htop(p)-hbot(p))) ** vegshape + fb = 1._r8 - (ol / max(1.e-06_r8, bendresist(p) * (htop(p)-hbot(p))) ** vegshape(p) else fb = 1._r8 - max(min(snow_depth(c),0.2_r8),0._r8)/0.2_r8 ! 0.2m is assumed !depth of snow required for complete burial of grasses