diff --git a/parteh/PRTAllometricCarbonMod.F90 b/parteh/PRTAllometricCarbonMod.F90 index 35d39a07e3..266ae8ccde 100644 --- a/parteh/PRTAllometricCarbonMod.F90 +++ b/parteh/PRTAllometricCarbonMod.F90 @@ -485,6 +485,16 @@ subroutine DailyPRTAllometricCarbon(this) call bstore_allom(dbh,ipft,canopy_trim,target_store_c) + + ! ----------------------------------------------------------------------------------- + ! II 1/2. Update target biomass based on the leaf status. + ! ----------------------------------------------------------------------------------- + select case (leaf_status) + case (leaves_off) + target_leaf_c = 0.0_r8 + end select + + ! ----------------------------------------------------------------------------------- ! III. Prioritize some amount of carbon to replace leaf/root turnover ! Make sure it isnt a negative payment, and either pay what is available @@ -1120,6 +1130,16 @@ subroutine DailyPRTAllometricCarbonSimpler(this) call bstore_allom(dbh,ipft,canopy_trim,target_store_c) + + ! ----------------------------------------------------------------------------------- + ! II 1/2. Update target biomass based on the leaf status. + ! ----------------------------------------------------------------------------------- + select case (leaf_status) + case (leaves_off) + target_leaf_c = 0.0_r8 + end select + + ! ----------------------------------------------------------------------------------- ! III. If carbon is available, bring all the pools as close to the allometry ! as possible. This also includes the storage pool, even though carbon may