Merge branch 'impactCalc_return_impactForecast' of https://github.com/CLIMADA-project/climada_python into impactCalc_return_impactForecast

Author: peanutfun

climada/engine/impact_forecast.py

@@ -98,3 +98,24 @@ def at_event(self, value):
         """Set the total exposure value close to a hazard"""
         LOGGER.warning("at_event for forecasts is not yet implemented.")
         self._at_event = value
+
+    def local_exceedance_impact(
+        self,
+        return_periods=(25, 50, 100, 250),
+        method="interpolate",
+        min_impact=0,
+        log_frequency=True,
+        log_impact=True,
+        bin_decimals=None,
+    ):
+        """Compution of local exceedance impact for given return periods is not
+        implemented for ImpactForecast. See climada.engine.impact.Impact for details.
+        Returns
+        -------
+        NotImplementedError
+        """
+
+        LOGGER.error("local_exceedance_impact is not defined for ImpactForecast")
+        raise NotImplementedError(
+            "local_exceedance_impact is not defined for ImpactForecast"
+        )

climada/engine/test/test_impact_forecast.py

@@ -103,3 +103,9 @@ def test_impact_forecast_concat(impact_forecast, member):
         [impact_forecast, impact_forecast], reset_event_ids=True
     )
     npt.assert_array_equal(impact_fc.member, np.concatenate([member, member]))
+
+
+def test_impact_forecast_exceedance_freq_curve_error(impact_forecast):
+    """Check if ImpactForecast.exceedance_freq_curve raises NotImplementedError"""
+    with pytest.raises(NotImplementedError):
+        impact_forecast.local_exceedance_impact(np.array([10, 50, 100]))