def test_cb_one_meas_fut_pass(self):
"""Test _cost_ben_one with same future"""
meas_name = 'Mangroves'
meas_val = dict()
meas_val['cost'] = (1.3117683608515418e+09, 1)
meas_val['risk'] = 4.850407096284983e+09
meas_val['efc'] = None
meas_val['risk_transf'] = 0
cb = CostBenefit()
cb.present_year = 2018
cb.future_year = 2040
cb.imp_meas_future['no measure'] = dict()
cb.imp_meas_future['no measure']['risk'] = 6.51220115756442e+09
disc_rates = DiscRates()
disc_rates.years = np.arange(2000, 2051)
disc_rates.rates = np.ones(disc_rates.years.size) * 0.02
time_dep = cb._time_dependency_array()
cb._cost_ben_one(meas_name, meas_val, disc_rates, time_dep)
self.assertAlmostEqual(cb.benefit[meas_name],
3.100583368954022e+10,
places=3)
self.assertAlmostEqual(cb.cost_ben_ratio[meas_name],
0.04230714690616641)
def test_cb_one_meas_pres_pass(self):
"""Test _cost_ben_one with different future"""
meas_name = 'Mangroves'
meas_val = dict()
meas_val['cost'] = (1.3117683608515418e+09, 1)
meas_val['risk'] = 4.826231151473135e+10
meas_val['efc'] = None
meas_val['risk_transf'] = 0
cb = CostBenefit()
cb.present_year = 2018
cb.future_year = 2040
cb.imp_meas_present['no measure'] = dict()
cb.imp_meas_present['no measure']['risk'] = 6.51220115756442e+09
cb.imp_meas_present['Mangroves'] = dict()
cb.imp_meas_present['Mangroves']['risk'] = 4.850407096284983e+09
cb.imp_meas_present['Mangroves']['risk_transf'] = 0
cb.imp_meas_future['no measure'] = dict()
cb.imp_meas_future['no measure']['risk'] = 5.9506659786664024e+10
disc_rates = DiscRates()
disc_rates.years = np.arange(2016, 2051)
disc_rates.rates = np.ones(disc_rates.years.size) * 0.02
time_dep = cb._time_dependency_array(1)
cb._cost_ben_one(meas_name, meas_val, disc_rates, time_dep)
self.assertAlmostEqual(cb.benefit[meas_name],
113345027690.81276,
places=3)
self.assertAlmostEqual(cb.cost_ben_ratio[meas_name],
0.011573232523528404)
def test_calc_cb_no_change_pass(self):
"""Test _calc_cost_benefit without present value against reference value"""
hazard = Hazard('TC')
hazard.read_mat(HAZ_TEST_MAT)
entity = Entity()
entity.read_mat(ENT_TEST_MAT)
entity.measures._data['TC'] = entity.measures._data.pop('XX')
for meas in entity.measures.get_measure('TC'):
meas.haz_type = 'TC'
entity.check()
cost_ben = CostBenefit()
cost_ben._calc_impact_measures(hazard, entity.exposures, entity.measures,
entity.impact_funcs, when='future', risk_func=risk_aai_agg, save_imp=True)
cost_ben.present_year = 2018
cost_ben.future_year = 2040
cost_ben._calc_cost_benefit(entity.disc_rates)
self.assertEqual(cost_ben.imp_meas_present, dict())
self.assertEqual(len(cost_ben.imp_meas_future), 5)
self.assertEqual(cost_ben.present_year, 2018)
self.assertEqual(cost_ben.future_year, 2040)
self.assertEqual(cost_ben.cost_ben_ratio['Mangroves'], 0.04230714690616641)
self.assertEqual(cost_ben.cost_ben_ratio['Beach nourishment'], 0.06998836431681373)
self.assertEqual(cost_ben.cost_ben_ratio['Seawall'], 0.2679741183248266)
self.assertEqual(cost_ben.cost_ben_ratio['Building code'], 0.30286828677985717)
self.assertEqual(cost_ben.benefit['Mangroves'], 3.100583368954022e+10)
self.assertEqual(cost_ben.benefit['Beach nourishment'], 2.468981832719974e+10)
self.assertEqual(cost_ben.benefit['Seawall'], 3.3132973770502796e+10)
self.assertEqual(cost_ben.benefit['Building code'], 3.0376240767284798e+10)
self.assertEqual(cost_ben.tot_climate_risk, 1.2150496306913972e+11)
def test_calc_cb_change_pass(self):
"""Test _calc_cost_benefit with present value against reference value"""
hazard = Hazard('TC')
hazard.read_mat(HAZ_TEST_MAT)
entity = Entity()
entity.read_mat(ENT_TEST_MAT)
entity.measures._data['TC'] = entity.measures._data.pop('XX')
for meas in entity.measures.get_measure('TC'):
meas.haz_type = 'TC'
entity.check()
cost_ben = CostBenefit()
cost_ben._calc_impact_measures(hazard, entity.exposures, entity.measures,
entity.impact_funcs, when='present', risk_func=risk_aai_agg, save_imp=False)
ent_future = Entity()
ent_future.read_excel(ENT_DEMO_FUTURE)
ent_future.check()
haz_future = copy.deepcopy(hazard)
haz_future.intensity.data += 25
cost_ben._calc_impact_measures(haz_future, ent_future.exposures, ent_future.measures,
ent_future.impact_funcs, when='future', risk_func=risk_aai_agg, save_imp=False)
cost_ben.present_year = 2018
cost_ben.future_year = 2040
cost_ben._calc_cost_benefit(entity.disc_rates, imp_time_depen=1)
self.assertEqual(cost_ben.present_year, 2018)
self.assertEqual(cost_ben.future_year, 2040)
self.assertEqual(cost_ben.tot_climate_risk, 5.768659152882021e+11)
self.assertEqual(cost_ben.imp_meas_present['no measure']['risk'], 6.51220115756442e+09)
self.assertEqual(cost_ben.imp_meas_present['Mangroves']['risk'], 4.850407096284983e+09)
self.assertEqual(cost_ben.imp_meas_present['Beach nourishment']['risk'], 5.188921355413834e+09)
self.assertEqual(cost_ben.imp_meas_present['Seawall']['risk'], 4.736400526119911e+09)
self.assertEqual(cost_ben.imp_meas_present['Building code']['risk'], 4.884150868173321e+09)
self.assertEqual(cost_ben.imp_meas_future['no measure']['risk'], 5.9506659786664024e+10)
self.assertEqual(cost_ben.imp_meas_future['Mangroves']['risk'], 4.826231151473135e+10)
self.assertEqual(cost_ben.imp_meas_future['Beach nourishment']['risk'], 5.0647250923231674e+10)
self.assertEqual(cost_ben.imp_meas_future['Seawall']['risk'], 21089567135.7345)
self.assertEqual(cost_ben.imp_meas_future['Building code']['risk'], 4.462999483999791e+10)
self.assertAlmostEqual(cost_ben.benefit['Mangroves'], 113345027690.81276)
self.assertAlmostEqual(cost_ben.benefit['Beach nourishment'], 89444869971.53653)
self.assertAlmostEqual(cost_ben.benefit['Seawall'], 347977469896.1333)
self.assertAlmostEqual(cost_ben.benefit['Building code'], 144216478822.05154)
self.assertAlmostEqual(cost_ben.cost_ben_ratio['Mangroves'], 0.011573232523528404)
self.assertAlmostEqual(cost_ben.cost_ben_ratio['Beach nourishment'], 0.01931916274851638)
self.assertAlmostEqual(cost_ben.cost_ben_ratio['Seawall'], 0.025515385913577368)
self.assertAlmostEqual(cost_ben.cost_ben_ratio['Building code'], 0.06379298728650741)
self.assertEqual(cost_ben.tot_climate_risk, 576865915288.2021)
def test_time_array_no_pres_pass(self):
"""Test _time_dependency_array"""
cb = CostBenefit()
cb.present_year = 2018
cb.future_year = 2030
time_arr = cb._time_dependency_array()
n_years = cb.future_year - cb.present_year + 1
self.assertEqual(time_arr.size, n_years)
self.assertTrue(np.array_equal(time_arr, np.ones(n_years)))
def test_npv_unaverted_no_pres_pass(self):
""" Test _npv_unaverted_impact """
cb = CostBenefit()
cb.present_year = 2018
cb.future_year = 2030
risk_future = 1000
disc_rates = DiscRates()
disc_rates.years = np.arange(cb.present_year, cb.future_year+1)
disc_rates.rates = np.ones(disc_rates.years.size)*0.025
time_dep = np.linspace(0, 1, disc_rates.years.size)
res = cb._npv_unaverted_impact(risk_future, disc_rates, time_dep,
risk_present=None)
self.assertEqual(res, disc_rates.net_present_value(cb.present_year, \
cb.future_year, time_dep * risk_future))
def test_time_array_pres_pass(self):
""" Test _time_dependency_array """
cb = CostBenefit()
cb.present_year = 2018
cb.future_year = 2030
imp_time_depen = 1.0
time_arr = cb._time_dependency_array(imp_time_depen)
n_years = cb.future_year - cb.present_year + 1
self.assertEqual(time_arr.size, n_years)
self.assertTrue(np.allclose(time_arr[:-1], np.arange(0, 1, 1/(n_years-1))))
self.assertEqual(time_arr[-1], 1)
imp_time_depen = 0.5
time_arr = cb._time_dependency_array(imp_time_depen)
n_years = cb.future_year - cb.present_year + 1
self.assertEqual(time_arr.size, n_years)
self.assertTrue(np.allclose(time_arr, np.arange(n_years)**imp_time_depen / \
(n_years-1)**imp_time_depen))
