num_windows=20
reg_spatial=np.logspace(0.25,0.35,num=2) #Set to False if you don't want to use it
reg_temporal=np.logspace(0,0,num=1) #Set to False if don't want to use it
nstrike=21 ; ndip=9 ; nfaults=(nstrike,ndip)
beta=45 #Rotational offset (in degrees) applied to rake (0 for normal)
solver='nnls' # 'lstsq','nnls'
top='free' ; bottom='locked' ; left='locked' ; right='locked' #'locked' or 'free'
bounds=(top,bottom,left,right)
################################################################################
######## Run-time modifications to the time series ############
decimate=4 #Decimate by constant (=None for NO decimation)
lowpass=None #Low pass corner frequency in Hz =None if no filter is desired
################################################################################
#Initalize project folders
if init==1:
runslip.init(home,project_name)
# Run green functions
if make_green==1 or make_synthetics==1:
runslip.inversionGFs(home,project_name,GF_list,tgf_file,fault_name,model_name,
dt,tsun_dt,NFFT,tsunNFFT,coord_type,make_green,make_synthetics,dk,pmin,
pmax,kmax,beta,time_epi,hot_start)
#Run inversion
if invert==1:
runslip.run_inversion(home,project_name,run_name,fault_name,model_name,GF_list,G_from_file,
G_name,epicenter,rupture_speed,num_windows,coord_type,reg_spatial,reg_temporal,
nfaults,beta,decimate,lowpass,solver,bounds)
top = 'free'
bottom = 'locked'
left = 'locked'
right = 'locked' #'locked' or 'free'
bounds = (top, bottom, left, right)
################################################################################e=
######## Run-time modifications to the time series ############
weight = True
decimate = None #Decimate by constant (=None for NO decimation)
bandpass = None #np.array([0.5]) #Corner frequencies in Hz =None if no filter is desired
################################################################################
#Initalize project folders
if init == 1:
runslip.init(home, project_name)
# Run green functions
if make_green == 1 or make_synthetics == 1:
runslip.inversionGFs(home, project_name, GF_list, tgf_file, fault_name,
model_name, dt, tsun_dt, NFFT, tsunNFFT, make_green,
make_synthetics, dk, pmin, pmax, kmax, beta, time_epi,
hot_start, ncpus)
#Run inversion
if invert == 1:
runslip.run_inversion(home, project_name, run_name, fault_name, model_name,
GF_list, G_from_file, G_name, epicenter,
rupture_speed, num_windows, reg_spatial,
reg_temporal, nfaults, beta, decimate, bandpass,
solver, bounds, weight, Ltype)
reg_temporal=None #Set to None if don't want to use it
nstrike=20 ; ndip=15 ; nfaults=(nstrike,ndip) #set nstrike to total no. of faults and ndip to 1 if using Tikh
beta=45 #Rotational offset (in degrees) applied to rake (0 for normal)
Ltype=0 # 0 for Tikhonov and 2 for Laplacian
solver='nnls' # 'lstsq','nnls'
top='free' ; bottom='locked' ; left='locked' ; right='locked' #'locked' or 'free'
bounds=(top,bottom,left,right)
################################################################################e=
######## Run-time modifications to the time series ############
weight=True
decimate=None #Decimate by constant (=None for NO decimation)
bandpass=np.array([0.5]) #Corner frequencies in Hz =None if no filter is desired
################################################################################
#Initalize project folders
if init==1:
runslip.init(home,project_name)
# Run green functions
if make_green==1 or make_synthetics==1:
runslip.inversionGFs(home,project_name,GF_list,tgf_file,fault_name,model_name,
dt,tsun_dt,NFFT,tsunNFFT,make_green,make_synthetics,dk,pmin,
pmax,kmax,beta,time_hypo,hot_start,ncpus,custom_stf)
#Run inversion
if invert==1:
runslip.run_inversion(home,project_name,run_name,fault_name,model_name,GF_list,G_from_file,
G_name,hypocenter,rupture_speed,num_windows,reg_spatial,reg_temporal,
nfaults,beta,decimate,bandpass,solver,bounds,weight,Ltype)