if plot_img:
exp_infr_meshs.plot_basemap()
exp_infr_meshs.plot_hexbin()
exp_infr_meshs.plot_scatter()
exp_infr_meshs.plot_raster()
#%% Impact
imp_infr_meshs = Impact()
imp_infr_meshs.calc(exp_infr_meshs, ifset_hail, haz_real, save_mat=True)
# imp_infr.plot_raster_eai_exposure()
freq_curve_infr_meshs = imp_infr_meshs.calc_freq_curve() #[1, 2, 5, 10, 20])
if plot_img:
freq_curve_infr_meshs.plot()
imp_infr_meshs.plot_basemap_eai_exposure()
imp_infr_meshs.plot_hexbin_eai_exposure()
imp_infr_meshs.plot_scatter_eai_exposure()
imp_infr_meshs.plot_raster_eai_exposure()
imp_infr_dur = Impact()
imp_infr_dur.calc(exp_infr_dur, ifset_hail, haz_dur, save_mat=True)
# imp_infr.plot_raster_eai_exposure()
freq_curve_infr_dur = imp_infr_dur.calc_freq_curve()
if plot_img:
freq_curve_infr_dur.plot()
imp_infr_dur.plot_basemap_eai_exposure()
imp_infr_dur.plot_hexbin_eai_exposure()
imp_infr_dur.plot_scatter_eai_exposure()
imp_infr_dur.plot_raster_eai_exposure()
print("dmg infr meshs {} Mio CHF, dmg infr dur {} Mio CHF".format(
imp_infr_meshs.aai_agg / 1e6, imp_infr_dur.aai_agg / 1e6))
exp_synth_agr_meshs.tag = Tag(file_name = "exp_agr", description="Exposure_description")
exp_synth_agr_meshs.plot_basemap()
exp_synth_agr_meshs.plot_hexbin()
exp_synth_agr_meshs.plot_scatter()
exp_synth_agr_meshs.plot_raster(raster_res = 0.001)
#%% Impact
imp_synth_infr_meshs = Impact()
imp_synth_infr_meshs.calc(exp_synth_infr_meshs, ifset_hail, haz_synth,save_mat=True)
freq_curve_synth_infr_meshs = imp_synth_infr_meshs.calc_freq_curve()
if plot_img:
freq_curve_synth_infr_meshs.plot()
imp_synth_infr_meshs.plot_basemap_eai_exposure()
imp_synth_infr_meshs.plot_hexbin_eai_exposure()
imp_synth_infr_meshs.plot_scatter_eai_exposure()
imp_synth_infr_meshs.plot_raster_eai_exposure(raster_res = 0.001)
imp_synth_agr_meshs = Impact()
imp_synth_agr_meshs.calc(exp_synth_agr_meshs, ifset_hail, haz_synth, save_mat=True)
freq_curve_synth_agr_meshs = imp_synth_agr_meshs.calc_freq_curve()
if plot_img:
freq_curve_synth_agr_meshs.plot()
imp_synth_agr_meshs.plot_basemap_eai_exposure()
imp_synth_agr_meshs.plot_hexbin_eai_exposure()
imp_synth_agr_meshs.plot_scatter_eai_exposure()
imp_synth_agr_meshs.plot_raster_eai_exposure(raster_res = 0.001)
print("dmg synth_infr_meshs {} Mio CHF, dmg synth_agr_meshs {} Mio CHF"
.format(imp_synth_infr_meshs.aai_agg/1e6, imp_synth_agr_meshs.aai_agg/1e6))
