haz_real.plot_intensity(event=0)
haz_real.plot_fraction(event=0)
#%% Impact_function
# Set impact function (see tutorial climada_entity_ImpactFuncSet)
ifset_hail = ImpactFuncSet()
for imp_fun_dict in imp_fun_parameter:
imp_fun = fct.create_impact_func(haz_type, imp_fun_dict["imp_id"],
imp_fun_dict["L"], imp_fun_dict["x_0"],
imp_fun_dict["k"])
ifset_hail.append(imp_fun)
ifset_hail.plot()
#%% Exposure
exp_infr = fct.load_exp_infr(force_new_hdf5_generation, name_hdf5_file,
input_folder, haz_real)
exp_meshs = fct.load_exp_agr(force_new_hdf5_generation, name_hdf5_file,
input_folder, haz_real)
exp_dur = exp_meshs.copy()
exp_dur["if_HL"] = exp_dur[
"if_HL"] + 3 #change if_HL to match the corresponding imp_id
if plot_img:
exp_infr.plot_basemap()
#This takes to long. Do over night!!!
#exp_agr.plot_basemap()
#%% Impact
imp_infr = Impact()
imp_infr.calc(exp_infr, ifset_hail, haz_real, save_mat=True)
# imp_infr.plot_raster_eai_exposure()
freq_curve_infr = imp_infr.calc_freq_curve()
#%% Impact_function
# Set impact function (see tutorial climada_entity_ImpactFuncSet)
ifset_hail = ImpactFuncSet()
for imp_fun_dict in imp_fun_parameter:
imp_fun = fct.create_impact_func(haz_type,
imp_fun_dict["imp_id"],
imp_fun_dict["L"],
imp_fun_dict["x_0"],
imp_fun_dict["k"])
ifset_hail.append(imp_fun)
if plot_img:
ifset_hail.plot()
#%% Exposure
exp_synth_infr_meshs = fct.load_exp_infr(force_new_hdf5_generation, name_hdf5_file, input_folder, haz_synth)
exp_synth_agr_meshs = fct.load_exp_agr(force_new_hdf5_generation, name_hdf5_file, input_folder, haz_synth)
if plot_img:
exp_synth_infr_meshs.plot_basemap()
exp_synth_infr_meshs.plot_hexbin()
exp_synth_infr_meshs.plot_scatter()
exp_synth_infr_meshs.plot_raster()
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()