Merge commit '429062528e015a0988b9f9aff4cdd715e9e73405' into sunt05/i…

…ssue234
UMEP-dev · Jun 16, 2024 · 334902d · 334902d
2 parents b8a5f39 + 4290625
commit 334902d
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/src/suews/src/suews_ctrl_driver.f95 b/src/suews/src/suews_ctrl_driver.f95
@@ -7634,7 +7634,10 @@ SUBROUTINE SUEWS_cal_surf_DTS( &
          ! get individual building fractions of each layer
          ! NB.: sum(sfr_roof) = building_frac(1)
          sfr_roof = 0.
-         IF (nlayer > 1) sfr_roof(1:nlayer - 1) = building_frac(1:nlayer - 1) - building_frac(2:nlayer)
+         IF (nlayer > 1) sfr_roof(1:nlayer - 1) = &
+            MAX( &
+            building_frac(1:nlayer - 1) - building_frac(2:nlayer), &
+            0.01) ! minimum value for sfr_roof to avoid zero fractions when adjacent layers have the same building fraction
          sfr_roof(nlayer) = building_frac(nlayer)
 
          ! get individual net building height of each layer