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)