haz_real.check() haz_synth.check() if plot_img: 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()
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"]) imp_fun.mdd[:] = 0.1 ifset_hail.append(imp_fun) if True: #plot_img: axis = None if True: num_plts = ifset_hail.size() num_row, num_col = u_plot._get_row_col_size(num_plts) _, axis = plt.subplots(num_row, num_col) for ax in axis: ax.set_ylim([9.5, 10.5]) # ifset_hail.plot(axis = axis) ifset_hail.plot(axis=axis) #%% Exposure exp_infr_meshs = fct.load_exp_infr(force_new_hdf5_generation, name_hdf5_file, input_folder, haz_real) exp_infr_dur = exp_infr_meshs.copy() exp_infr_dur["if_HL"] = 8 #change if_HL to match the corresponding imp_id 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()