Add if statement to prevent negative leafc_slamax

NGEET · Jun 15, 2018 · a389cf8 · a389cf8
1 parent f50bb38
commit a389cf8
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/biogeochem/FatesAllometryMod.F90 b/biogeochem/FatesAllometryMod.F90
@@ -595,6 +595,9 @@ real(r8) function tree_lai( bl, status_coh, pft, c_area, n, cl, canopy_layer_tai
        ! Leafc_per_unitarea at which sla_max is reached due to exponential sla profile in canopy:
        leafc_slamax = (slat - sla_max * exp(-1.0_r8 * kn * laican)) / &
             (-1.0_r8 * kn * tai_to_lai * slat * sla_max)
+       if(leafc_slamax < 0.0_r8)then
+            leafc_slamax = 0.0_r8
+       endif
 
        ! Calculate tree_lai (m2 leaf area /m2 ground) = unitless LAI
        !----------------------------------------------------------------------