diff --git a/main/EDMainMod.F90 b/main/EDMainMod.F90 index 2b29157512..fea464c08d 100644 --- a/main/EDMainMod.F90 +++ b/main/EDMainMod.F90 @@ -195,16 +195,13 @@ subroutine ed_ecosystem_dynamics(currentSite, bc_in, bc_out) ! Zero mass balance call TotalBalanceCheck(currentSite, 0) - ! We do not allow phenology while in ST3 mode either, it is hypothetically - ! possible to allow this, but we have not plugged in the litter fluxes - ! of flushing or turning over leaves for non-dynamics runs - if (hlm_use_ed_st3.eq.ifalse)then - if(hlm_use_sp.eq.ifalse) then - call phenology(currentSite, bc_in ) - else - call satellite_phenology(currentSite, bc_in ) - end if ! SP phenology - end if + ! Phenology. For now ST3 mode also represents phenology, though it is + ! an incomplete representation as it lacks transfers to the litter layer + if (hlm_use_sp == ifalse) then + call phenology(currentSite, bc_in ) + else + call satellite_phenology(currentSite, bc_in ) + end if ! SP phenology if (hlm_use_ed_st3.eq.ifalse.and.hlm_use_sp.eq.ifalse) then ! Bypass if ST3 diff --git a/parteh/PRTLossFluxesMod.F90 b/parteh/PRTLossFluxesMod.F90 index 833bbc7a40..8d71011498 100644 --- a/parteh/PRTLossFluxesMod.F90 +++ b/parteh/PRTLossFluxesMod.F90 @@ -24,6 +24,7 @@ module PRTLossFluxesMod use FatesConstantsMod, only : nearzero use FatesConstantsMod, only : calloc_abs_error use FatesConstantsMod, only : itrue + use FatesInterfaceTypesMod, only : hlm_use_ed_st3 use FatesGlobals , only : endrun => fates_endrun use FatesGlobals , only : fates_log use shr_log_mod , only : errMsg => shr_log_errMsg @@ -567,6 +568,13 @@ subroutine DeciduousTurnoverSimpleRetranslocation(prt,ipft,organ_id,mass_fractio if( prt_params%organ_param_id(organ_id) < 1 ) then retrans = 0._r8 + else if (hlm_use_ed_st3 == itrue) then + ! For now we retranslocate all biomass to storage to prevent carbon balance errors. + ! This is not an ideal solution because tissue abscission may contribute to a + ! seasonal sign in litter and consequently heterotrophic respiration (though the + ! lack of phenology is also problematic for seasonal cycles of GPP, ET, sensible + ! heat, etc.). + retrans = 1._r8 else if ( element_id == carbon12_element ) then retrans = 0._r8