'../green_function/G_He50km_Vis2.8E19_Rake85.h5'
],
nlin_par_names = ['log10(visM)','log10(He)','rake'],
file_d = '../obs/post_obs.h5',
file_sd = 'sd/sd_hor_1_ver_5_sea_inf.h5',
file_incr_slip0 = '../slip0/incr_slip0.h5',
file_co_slip = '../slip0/incr_slip0.h5',
filter_sites_file = 'sites_without_seafloor',
epochs = epochs,
regularization = None,
basis = basis,
)
inv.set_data_except_L()
reg_temp = 0
reg_edge = 0.06
for nrough, reg_rough in enumerate(reg_roughes):
outfname = 'outs/nrough_%02d.h5'%(nrough)
if exists(outfname):
print("Skip %s !"%outfname)
continue
print(outfname)
inv.regularization = \
reg = vj.create_temporal_edge_roughening(
fault_file, epochs, reg_temp, reg_edge, reg_rough)
inv.set_data_L()
inv.run()
inv.save(outfname, overwrite=True)
file_d = '../obs/cumu_post_with_seafloor.h5',
file_sd = 'sd/sd_uniform_seafloor_inf.h5',
file_incr_slip0 = 'slip0/incr_slip0.h5',
filter_sites_file = 'sites_2EXPs',
epochs = epochs,
regularization = None,
basis = basis,
)
inv.set_data_except(excepts=['L'])
reg_temp = 0
reg_edge = 0.06
for nrough, reg_rough in enumerate(reg_roughes):
outfname = 'outs/nrough_%02d.h5'%(nrough)
if exists(outfname):
print("Skip %s !"%outfname)
continue
print(outfname)
inv.regularization = \
reg = vj.create_temporal_edge_roughening(
fault_file, epochs, reg_temp, reg_edge, reg_rough)
inv.set_data_L()
inv.run()
inv.save(outfname, overwrite=True)