-
Notifications
You must be signed in to change notification settings - Fork 16
/
setup_parameters.m
90 lines (81 loc) · 3.93 KB
/
setup_parameters.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
% Script to setup parameters used for the whole project
%%%% %%%% %%%% %%%% %%%% %%%% %%%% %%%% %%%% %%%%
%%%% Global settings
parameters.proj_name = 'TA';
parameters.component = 'LHZ'; % determined by filenames
parameters.lalim=[25 50];
parameters.lolim=[-125 -65];
parameters.gridsize=0.3; % in degrees
parameters.periods = [20 25 32 40 50 60 80 100]; % in seconds
%%%% %%%% %%%% %%%% %%%% %%%% %%%% %%%% %%%% %%%%
% parameters for data downloading (if using IRIS DMC)
parameters.start_time = '2009-01-07 00:00:00';
parameters.end_time = '2009-06-08 00:00:00'; % put '' for using 4 days before current date
parameters.is_use_timestamp = 0;
parameters.network = '_US-ALL';
parameters.minMw = 6;
parameters.maxdepth = 50;
parameters.datalength = 7200; % in second
parameters.resample_delta = 1; % in second
%%%% %%%% %%%% %%%% %%%% %%%% %%%% %%%% %%%% %%%%
% parameters for the auto_win_select.m
parameters.largest_epidist_range = 3000;
parameters.cycle_before = 2;
parameters.cycle_after = 5;
parameters.min_dist_tol = deg2km(20);
parameters.max_dist_tol = deg2km(160);
parameters.min_groupv = 2;
parameters.max_groupv = 5;
parameters.cent_freq = 0.025;
%%%% %%%% %%%% %%%% %%%% %%%% %%%% %%%% %%%% %%%%
% parameters for the cross-correlation measurement
% (gsdfmain.m)
parameters.minstadist = 5;
parameters.maxstadist = 200; % station cross-correlation distance in km
parameters.is_rm_resp = 0;
parameters.periods = sort(parameters.periods); % make sure periods are ascending
parameters.refv = 4; % to select the correct cycle
parameters.refphv = ones(size(parameters.periods))*4;
parameters.min_width = 0.06; % to build up gaussian filters
parameters.max_width = 0.10;
parameters.wintaperlength = 30; % taper to build up the isolation filter
parameters.prefilter = [10,200];
parameters.xcor_win_halflength = 150; % window for the cross-correlation
parameters.xcor_win_iter = zeros(size(parameters.periods)); % re-apply the xcor window due to measured group delay, should be same length as periods, not used anymore
parameters.Nfit = 2;
parameters.Ncircle = 5;
parameters.cohere_tol = 0.5; % minimum coherenecy between two stations
parameters.tp_tol = 10; % seconds away from averaged phase velocity
%%%% %%%% %%%% %%%% %%%% %%%% %%%% %%%% %%%% %%%%
% parameters for the tomography
% (eikonal_eq.m helmholtz_eq.m)
parameters.smweight_array = 3*[0.4 0.3 0.2 0.2 0.2 0.5 1 2]; % smoothing weight for the deltaSx and delta Sy
parameters.raydensetol=deg2km(parameters.gridsize)*2;
parameters.Tdumpweight = 0; % dumping the ray to the girgle circle path
parameters.Rdumpweight = 0; % dumping the region to have the same phase velocity
parameters.fiterrtol = 3; % error allowed in the wavelet fitting
parameters.isRsmooth = 1; % smoothing due to Sx and Sy or Sr and S_theta
parameters.dterrtol = 2; % largest variance of the inversion error allowed
parameters.inverse_err_tol = 2; % count be number of standard devition
parameters.min_amp_tol = 0.1; % station with amplitude smaller than this ratio of average amplitude will not be used.
parameters.amp_var_tol = 2; % how much times change of amplitude of single station to the mean value of nearby stations should be considered as bad measurement
parameters.alpha_range = [1 1];
parameters.alpha_search_grid = 0.1;
%%%% %%%% %%%% %%%% %%%% %%%% %%%% %%%% %%%% %%%%
% parameter for stacking
% (stack_phv.m stack_helm.m)
parameters.min_csgoodratio=[3 3 3 3 5 10 15 15]; % minimum radio between good and bad measurements for a good event
parameters.min_phv_tol = 3;
parameters.max_phv_tol = 5;
parameters.is_raydense_weight = 1;
parameters.min_event_num = 10;
parameters.err_std_tol = 2;
parameters.issmoothmap = 1;
parameters.smooth_wavelength = 0.25;
parameters.event_bias_tol = 2;
%%%% %%%% %%%% %%%% %%%% %%%% %%%% %%%% %%%% %%%%
% parameters for azimuthal anisotropy inversion
% beta version
parameters.smsize = 1; % averaging nearby grid number
parameters.off_azi_tol = 30; % differ from great circle path
parameters.is_one_phi = 1;