Integration of CalcDeltaImpact class of unsequa into develop

CHANGELOG.md

@@ -14,6 +14,8 @@ Code freeze date: YYYY-MM-DD
 
 - Convenience method `api_client.Client.get_dataset_file`, combining `get_dataset_info` and `download_dataset`, returning a single file objet. [#821](https://github.com/CLIMADA-project/climada_python/pull/821)
 - Read and Write methods to and from csv files for the `DiscRates` class. [#818](ttps://github.com/CLIMADA-project/climada_python/pull/818)
+- Add `CalcDeltaClimate` to unsequa module to allow uncertainty and sensitivity analysis of impact change calculations. [#844](https://github.com/CLIMADA-project/climada_python/pull/844)
+- Addition of small function called `safe_divide` in util which handles division by zero and NaN values in the numerator or denominator. [#844](https://github.com/CLIMADA-project/climada_python/pull/844)
 
 ### Changed
 

climada/engine/unsequa/__init__.py

@@ -22,3 +22,4 @@
 from .calc_base import *
 from .calc_impact import *
 from .calc_cost_benefit import *
+from .calc_delta_climate import *

climada/engine/unsequa/calc_base.py

@@ -179,7 +179,6 @@ def est_comp_time(self, n_samples, time_one_run, processes=None):
                 "been assigned to exp before defining input_var, ..."
                 "\n If computation cannot be reduced, consider using"
                 " a surrogate model https://www.uqlab.com/", time_one_run)
-
         total_time = n_samples * time_one_run / processes
         LOGGER.info("\n\nEstimated computaion time: %s\n",
                     dt.timedelta(seconds=total_time))
@@ -323,15 +322,17 @@ def sensitivity(self, unc_output, sensitivity_method = 'sobol',
 
         Parameters
         ----------
-        unc_output : climada.engine.uncertainty.unc_output.UncOutput
+        unc_output : climada.engine.uncertainty.unc_output.UncOutput()
             Uncertainty data object in which to store the sensitivity indices
-        sensitivity_method : str, optional
-            Sensitivity analysis method from SALib.analyse. Possible choices: 'fast', 'rbd_fact',
-            'morris', 'sobol', 'delta', 'ff'. Note that in Salib, sampling methods and sensitivity
-            analysis methods should be used in specific pairs:
+        sensitivity_method : str
+            sensitivity analysis method from SALib.analyse
+            Possible choices:
+                'fast', 'rbd_fact', 'morris', 'sobol', 'delta', 'ff'
+            The default is 'sobol'.
+            Note that in Salib, sampling methods and sensitivity analysis
+            methods should be used in specific pairs.
             https://salib.readthedocs.io/en/latest/api.html
-            Default: 'sobol'
-        sensitivity_kwargs: dict, optional
+        sensitivity_kwargs: dict(), optional
             Keyword arguments of the chosen SALib analyse method.
             The default is to use SALib's default arguments.
 

climada/engine/unsequa/calc_cost_benefit.py

@@ -68,7 +68,9 @@ class CalcCostBenefit(Calc):
         ('haz_input_var', 'ent_input_var', 'haz_fut_input_var', 'ent_fut_input_var')
     _metric_names : tuple(str)
         Names of the cost benefit output metrics
-        ('tot_climate_risk', 'benefit', 'cost_ben_ratio', 'imp_meas_present', 'imp_meas_future')
+        ('tot_climate_risk', 'benefit', 'cost_ben_ratio',
+         'imp_meas_present', 'imp_meas_future')
+
     """
 
     _input_var_names = (