diff --git a/biogeochem/FatesAllometryMod.F90 b/biogeochem/FatesAllometryMod.F90 index 41980b001c..432e1efb1b 100644 --- a/biogeochem/FatesAllometryMod.F90 +++ b/biogeochem/FatesAllometryMod.F90 @@ -671,7 +671,7 @@ subroutine bsap_allom(d,ipft,canopy_trim,bsap,dbsapdd) bsap_cap = max_frac*(bagw+bbgw) if(bsap>bsap_cap) then - bsap = min( bsap_cap,bsap) + bsap = bsap_cap if(present(dbsapdd))then dbsapdd = max_frac*(dbagwdd+dbbgwdd) end if