10, 15
] #Transition depths for rise time scaling (if slip shallower than first index, rise times are twice as long as calculated)
buffer_factor = 0.5 # I don't think this does anything anymore-- remove?
mean_slip_name = home + project_name + '/forward_models/mentawai_fine.rupt'
#mean_slip_name=None
shear_wave_fraction = 0.8
kappa = None
force_area = True
force_magnitude = False
force_hypocenter = True
###############################################################################
#Initalize project folders
if init == 1:
fakequakes.init(home, project_name)
#Generate rupture models
if make_ruptures == 1:
fakequakes.generate_ruptures(home,
project_name,
run_name,
fault_name,
slab_name,
mesh_name,
load_distances,
distances_name,
UTM_zone,
target_Mw,
model_name,
hurst,
G_name='GFs'
# Displacement and velocity waveform parameters
NFFT=1024 ; dt=1.0
#fk-parameters
dk=0.1 ; pmin=0 ; pmax=1 ; kmax=20
custom_stf=None
###############################################################################
#Initalize project folders
if init==1:
fakequakes.init(home,project_name)
#Generate rupture models
if make_ruptures==1:
fakequakes.generate_ruptures(home,project_name,run_name,fault_name,slab_name,
mesh_name,load_distances,distances_name,UTM_zone,target_Mw,model_name,
hurst,Ldip,Lstrike,num_modes,Nrealizations,rake,buffer_factor,
rise_time_depths,time_epi,max_slip,source_time_function,lognormal,
slip_standard_deviation,scaling_law,force_magnitude=force_magnitude,
hypocenter=hypocenter,force_hypocenter=force_hypocenter)
# Prepare waveforms and synthetics
if make_GFs==1 or make_synthetics==1:
runslip.inversionGFs(home,project_name,GF_list,None,fault_name,model_name,
dt,None,NFFT,None,make_GFs,make_synthetics,dk,pmin,