use if_ impact functions if no if_hazard_type present in exposures

Author: gabrielaznar

climada/engine/impact.py

@@ -161,10 +161,13 @@ def calc(self, exposures, impact_funcs, hazard, save_mat=False):
         # Get damage functions for this hazard
         if_haz = INDICATOR_IF + hazard.tag.haz_type
         haz_imp = impact_funcs.get_func(hazard.tag.haz_type)
-        if if_haz not in exposures:
-            LOGGER.error('Missing exposures column %s. No exposures with impact'\
-                         +' functions for peril %s.', if_haz, hazard.tag.haz_type)
+        if if_haz not in exposures and INDICATOR_IF not in exposures:
+            LOGGER.error('Missing exposures impact functions %s.', INDICATOR_IF)
             raise ValueError
+        elif if_haz not in exposures:
+            LOGGER.info('Missing exposures impact functions for hazard %s. ' +\
+                        'Using impact functions in %s.', if_haz, INDICATOR_IF)
+            if_haz = INDICATOR_IF
 
         # Check if deductible and cover should be applied
         insure_flag = False
@@ -203,7 +206,7 @@ def calc(self, exposures, impact_funcs, hazard, save_mat=False):
             self.imp_mat = self.imp_mat.tocsr()
 
     def plot_eai_exposure(self, mask=None, ignore_zero=True,
-                          pop_name=True, buffer_deg=0.0, extend='neither',
+                          pop_name=True, buffer=0.0, extend='neither',
                           var_name=None, **kwargs):
         """Plot expected annual impact of each exposure.
 
@@ -212,7 +215,7 @@ def plot_eai_exposure(self, mask=None, ignore_zero=True,
             ignore_zero (bool, optional): flag to indicate if zero and negative
                 values are ignored in plot. Default: False
             pop_name (bool, optional): add names of the populated places
-            buffer_deg (float, optional): border to add to coordinates.
+            buffer (float, optional): border to add to coordinates.
                 Default: 1.0.
             extend (str, optional): extend border colorbar with arrows.
                 [ 'neither' | 'both' | 'min' | 'max' ]
@@ -235,7 +238,7 @@ def plot_eai_exposure(self, mask=None, ignore_zero=True,
             kwargs['reduce_C_function'] = np.sum
         return u_plot.geo_bin_from_array(self.eai_exp[mask][pos_vals], \
             self.coord_exp[mask][pos_vals], var_name, title, pop_name, \
-            buffer_deg, extend, **kwargs)
+            buffer, extend, **kwargs)
 
     def _exp_impact(self, exp_iimp, exposures, hazard, imp_fun, insure_flag):
         """Compute impact for inpute exposure indexes and impact function.

climada/engine/test/test_impact.py

@@ -239,6 +239,48 @@ def test_calc_imp_mat_pass(self):
         self.assertTrue(np.allclose(np.array(np.sum(np.multiply(impact.imp_mat.todense(),
             impact.frequency.reshape(-1, 1)), axis=0)).reshape(-1), impact.eai_exp))
 
+    def test_calc_if_pass(self):
+        """ Execute when no if_HAZ present, but only if_ """
+        ent = Entity()
+        ent.read_excel(ENT_DEMO_TODAY)
+        ent.exposures.rename(columns={'if_TC':'if_'}, inplace=True)
+        ent.check()
+        
+        # Read default hazard file
+        hazard = Hazard('TC')
+        hazard.read_mat(HAZ_TEST_MAT)
+        # Create impact object
+        impact = Impact()
+        impact.calc(ent.exposures, ent.impact_funcs, hazard)
+
+        # Check result
+        num_events = len(hazard.event_id)
+        num_exp = ent.exposures.shape[0]
+        # Check relative errors as well when absolute value gt 1.0e-7
+        # impact.at_event == EDS.damage in MATLAB
+        self.assertEqual(num_events, len(impact.at_event))
+        self.assertEqual(0, impact.at_event[0])
+        self.assertEqual(0, impact.at_event[int(num_events/2)])
+        self.assertAlmostEqual(1.472482938320243e+08, impact.at_event[13809])
+        self.assertEqual(7.076504723057619e+10, impact.at_event[12147])
+        self.assertEqual(0, impact.at_event[num_events-1])
+        # impact.eai_exp == EDS.ED_at_centroid in MATLAB
+        self.assertEqual(num_exp, len(impact.eai_exp))
+        self.assertAlmostEqual(1.518553670803242e+08, impact.eai_exp[0])
+        self.assertAlmostEqual(1.373490457046383e+08, \
+                               impact.eai_exp[int(num_exp/2)], 6)
+        self.assertTrue(np.isclose(1.373490457046383e+08, \
+                                          impact.eai_exp[int(num_exp/2)]))
+        self.assertAlmostEqual(1.066837260150042e+08, \
+                               impact.eai_exp[num_exp-1], 6)
+        self.assertTrue(np.isclose(1.066837260150042e+08, \
+                                          impact.eai_exp[int(num_exp-1)]))
+        # impact.tot_value == EDS.Value in MATLAB
+        # impact.aai_agg == EDS.ED in MATLAB
+        self.assertAlmostEqual(6.570532945599105e+11, impact.tot_value)
+        self.assertAlmostEqual(6.512201157564421e+09, impact.aai_agg, 5)
+        self.assertTrue(np.isclose(6.512201157564421e+09, impact.aai_agg))
+
 class TestIO(unittest.TestCase):
     ''' Test impact calc method.'''