Added a condition in class Nature(HSIProcedure) so that hsi is not calculated...

Added a condition in class Nature(HSIProcedure) so that hsi is not calculated when environmental_factor_array is None and self._suitabilities[land_type} is below hsi_threshold.

if environmental_factor_array is None:
   if self._suitabilities[land_type] >= self._hsi_threshold:

mnp/species_models/habitat_suitability.py

@@ -60,9 +60,10 @@ class Nature(HSIProcedure):
                 if land_type in land_types.keys():
                     # Just add land type suitabilities if not running with environmentals
                     if environmental_factor_array is None:
-                        hsi += land_types[land_type].multiply(
-                            self._suitabilities[land_type]
-                        )
+                        if self._suitabilities[land_type] >= self._hsi_threshold:
+                            hsi += land_types[land_type].multiply(
+                                self._suitabilities[land_type]
+                            )
 
                     # Otherwise, multiply hsi with environmental array
                     else:

[Roelofsen, Hans]

First commit, welcome Jan!