1### Miscellaneous
2misc:
3 # Use direct- or coda-wave?
4 wavetype: 'direct-wave'
5 # wavetype: 'coda-wave'
6 # Number of cpus
7 ncpu: 20
8 # Logging level with increasing information:
9 # 'CRITICAL', 'ERROR', 'WARNING', 'INFO', 'DEBUG'
10 log_level: 'INFO'
11
12### Directories
13dir:
14 # Everything will be relative to this directory
15 project: '/home/szhang/Desktop/ThreeStation/example/data'
16 # Two-statoin ambient noise interferogram
17 I2: 'I2'
18 # Output
19 out: 'test_I3_hyp'
20 # Store paths to I2 and common source stations
21 meta: 'metadata'
22 # Write stacked I3
23 I3: 'stack_all'
24 # Write randomly stacked I3 for uncertainty estimation
25 I3_rand: 'stack_rand'
26
27### Station files
28fstation:
29 # Receiver-station names, [group number]
30 receiver: 'receiver.csv'
31 # Source-station names
32 source: 'source.csv' # [stnm]
33 # All stations: network, station, longitude, latitude
34 all: 'all.csv'
35
36### Signal-to-noise ratio
37# Signal window: [tmin - bfact*Tmax, tmax + efact*Tmax]
38# Noise window: [tmax + efact*Tmax + dsn, + nlen]
39snr:
40 # Minimum SNR of I2 for accepting source-stations
41 min: 10
42 vmax: 5
43 vmin: 1
44 Tmin: 5
45 Tmax: 150
46 bfact: 1
47 efact: 1
48 dsn: 500
49 nlen: 500
50
51### Cutting signals from two-station interferograms
52# Direct-wave window: [tmin - bfact_dw*Tmax, tmax + efact_dw*Tmax]
53# Coda window: [tmax + efact_cw*Tmax, te]
54cut:
55 # Desired sampling interval
56 delta: 1
57 # Minimum length of window (s)
58 min_len: 1500
59 # If mute outside surface wave window
60 mute: False
61 # Only mute precursory for calculating SNR
62 mute_prc_only: True
63 vmin: 1
64 vmax: 5
65 Tmax: 150
66 # Maximum lag time
67 te: 3000
68 bfact_dw: 1
69 efact_dw: 1
70 efact_cw: 1
71 # Ensure that trailing noise in SNR is non-zero
72 restrict_src_dist: True
73 cw4snr: 1000
74
75### Preprocess two-station interferograms
76preproc:
77 # One-bit time-domain normalization
78 onebit: False
79 # Spectral whitening
80 whiten:
81 val: False
82 plot: False
83 # Brute force, i.e., make spectral amplitudes to 1
84 brute: False
85 # Window width for running average (Hz)
86 freq_width: 0.001
87 Tmin: 5
88 Tmax: 80
89
90### Three-station interferometry
91interferometry:
92 # Interferometry operator: 'corrlation' or 'convolution'
93 # operator: 'convolution'
94 operator: 'correlation'
95 # Whether de-bias source-specific interferogram via phase shift
96 phase_shift: False
97 # Whether pick the appropriate lag (causal or acausal) for
98 # direct-wave correlations
99 pick_lag: True
100 # 1-D predicted phase velocity in two columns: [period, phase velocity]
101 pred_pv_1d:
102 # 2-D predicted phase velocity in a dictionary
103 # mapping names of station-pair to [period, phase velocity]
104 pred_pv_2d:
105 # Whether trim the convolutional interferogram at the end
106 # to make it the same as the input length
107 trim_conv: True
108 # Whether save source-station direction to SAC header
109 # 1 or 2 correspond to whether the source is close to station1 or station2
110 return_srcdir: True
111 # Limit source-stations to stationary phase zone
112 spz: True
113 # Define SPZ using a constant, a fraction of inter-receiver distance, or an azimuth
114 method: 'const'
115 # Maximum deviation from inter-receiver distance (km)
116 max_dr: 30
117 # Maximum deviation from inter-receiver distance (%)
118 max_drpct: 1
119 # Minimum source-receiver distance (km)
120 min_srdist: 30
121 # Whether use Welch's method in coda correlation
122 Welch: False
123 # Sub-window length for Welch's method (s)
124 subwin: 800
125 # Sub-window overlapping fraction for Welch's method
126 lap: .5
127 # Whether flip negative lag
128 flip_nlag: True
129 # No. of lags from two-station interferograms to use
130 # 2: positive-positive and negative-negative
131 # 4: also include positive-negative and negative-positive
132 nlag: 4
133 # Demean data before correlation
134 demean: False
135 # Normalize correlation
136 normalize: False
137 # Correlate in 'time' or 'freq' domain
138 domain: 'freq'
139 # Symmetrize correlation
140 symmetric: True
141
142### Stack source-specific interferograms
143stack:
144 # Minimum source-stations required for stack
145 min_src: 1
146 # How to weigh individual C3 in stack
147 # 'RMS': reciprocal of RMS of trailing noise
148 # 'unit': reciprocal of absolute maximum
149 # 'None': no weight
150 # 'SNR': signal-to-noise ratio
151 weight: 'rms_noise'
152 # Whether output symmetric component
153 symmetric: False
154 # Whether sort coda correlations by length before stacking
155 sort: False
156 # Whehter compute random stacks for uncertainty estimation
157 rand: False
158 # Number of random stacks
159 nsub: 4
160
161### Whether skip certain processing
162skip:
163 # Skip finding paths to correlations if files already exist
164 find_I2: False
165 # Skip cutting lags if already exists
166 cut: False
167 # Skip preprocessing lags if already exists
168 preprocess: False
169 # Skip finding source stations if files already exist
170 find_source: False
171 # Skip computing C3 if already exists
172 C3: False
173 # Skip stacking for I3 if already exists
174 I3: False
175
176### Whether write intermediate files
177write:
178 # Path to I2 & common source-stations
179 meta: True
180 # I2 lags
181 lag: True
182 # C3
183 C3: True
184 # I3
185 stack: True
186
187### Suffix and prefix for file names
188# Suffix of metadata
189sfx:
190 # Appended to files of paths
191 I2: '_I2'
192 # Appended to files of sources
193 source: '_source'
194# Prefix
195pfx:
196 plag: 'P'
197 nlag: 'N'
198 sym: 'S'