From 9857e2fdc05267724181698fb38f4fce52303bed Mon Sep 17 00:00:00 2001 From: Xiulin Gao Date: Thu, 13 Jun 2024 18:17:50 -0600 Subject: [PATCH] simplify calculation of biomass derivatives using biomass --- biogeochem/FatesAllometryMod.F90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/biogeochem/FatesAllometryMod.F90 b/biogeochem/FatesAllometryMod.F90 index 8b9fbebb39..d0957ef190 100644 --- a/biogeochem/FatesAllometryMod.F90 +++ b/biogeochem/FatesAllometryMod.F90 @@ -1732,7 +1732,7 @@ subroutine dh2blmax_3pwr_grass(d,h,dhdd,p1,p2,p3,dbh_maxh,c2b,blmax,dblmaxdd) if(d .ge. dbh_maxh)then dblmaxdd = 0._r8 else - dblmaxdd = p1 * (p2 * duse**(p2 - 1.0_r8) * h**p3 + p3 * h**(p3 - 1.0_r8) * dhdd * duse**p2) / c2b + dblmaxdd = blmax * (p2 / duse + p3 * dhdd / h) end if end if @@ -2177,7 +2177,7 @@ subroutine dh2bagw_3pwr_grass(d,h,dhdd,p1,p2,p3,c2b,bagw,dbagwdd) !----Compute the aboveground biomass derivative with basal diameter if needed if (present(dbagwdd)) then - dbagwdd = p1 * (p2 * d**(p2 - 1.0_r8) * h**p3 + p3 * h**(p3 - 1.0_r8) * dhdd * d**p2) / c2b + dbagwdd = p2 * bagw / d + p3 * bagw * dhdd / h end if return