Skip to content

Commit

Permalink
simplify calculation of biomass derivatives using biomass
Browse files Browse the repository at this point in the history
  • Loading branch information
XiulinGao committed Jun 14, 2024
1 parent 812c7cc commit 9857e2f
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions biogeochem/FatesAllometryMod.F90
Original file line number Diff line number Diff line change
Expand Up @@ -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

Expand Down Expand Up @@ -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
Expand Down

0 comments on commit 9857e2f

Please sign in to comment.