stf_falloff_rate = 6 #Only affects Dreger STF, choose 4-8 are reasonable values
num_modes = 500
rake = 180
rise_time_depths = [5, 8] #Transition depths for rise time scaling
buffer_factor = 0.5
###############################################################################
#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)
# 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,
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,
ncpus,
mean_slip_name=mean_slip_name,
force_magnitude=force_magnitude,
force_area=force_area,
hypocenter=hypocenter,
force_hypocenter=force_hypocenter,
shear_wave_fraction=shear_wave_fraction,
max_slip_rule=max_slip_rule,
previous_runs=previous_runs)
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,
pmax,kmax,0,time_epi,0,ncpus,custom_stf,impulse=True)
# Synthesize the waveforms
if make_waveforms==1:
forward.waveforms_fakequakes(home,project_name,fault_name,rupture_list,GF_list,
model_name,run_name,dt,NFFT,G_from_file,G_name,source_time_function)