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,