Added >0 constraints on gamma tleaf and par

noaa-oar-arl · Mar 18, 2024 · a8d1dea · a8d1dea
1 parent d23126f
commit a8d1dea
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/canopy_bioemi_mod.F90 b/src/canopy_bioemi_mod.F90
@@ -165,6 +165,8 @@ SUBROUTINE CANOPY_BIO( ZK, FCLAI, FCH, LAI, FSUN, &
 
         GammaTLEAF_AVE = (GammaTLEAF_SUN*FSUN) + (GammaTLEAF_SHADE*(1.0-FSUN)) ! average = sum sun and shade weighted by sunlit fraction
 
+        GammaTLEAF_AVE = MAX( GammaTLEAF_AVE, 0.0_rk )
+
 ! Calculate gamma (activity) values for average PPFD (Clifton et al., 2022; based on Guenther et al. 2012)
         ALPHA_P_SUN = 0.004 - 0.0005*log(PPFD240_SUN)
         ALPHA_P_SHADE = 0.004 - 0.0005*log(PPFD240_SHADE)
@@ -175,6 +177,8 @@ SUBROUTINE CANOPY_BIO( ZK, FCLAI, FCH, LAI, FSUN, &
 
         GammaPPFD_AVE = (GammaPPFD_SUN*FSUN) + (GammaPPFD_SHADE*(1.0-FSUN)) ! average = sum sun and shade weighted by sunlit fraction
 
+        GammaPPFD_AVE = MAX( GammaPPFD_AVE, 0.0_rk )
+
 ! Get CO2 inhibition factor for isoprene only
 
         if (EMI_IND .eq. 1) then  !Isoprene