forked from ahotovec/REDPy
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathsettings.cfg
executable file
·209 lines (187 loc) · 9.7 KB
/
settings.cfg
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
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
[Settings]
# Above line is REQUIRED, do not edit!
############################# Configuration file for REDPy #############################
# Settings not explicitly defined will be assigned a default value, so feel free to #
# delete all of the comment lines and any parameters you wish to keep at default for #
# a cleaner configuration file. #
########################################################################################
###### RUN PARAMETERS ######
# Describe the run - this will be used in the title of the webpages and master plot title
title=REDPy Catalog
# Path to the folder that will contain outputs, not including groupName (below); leave
# empty to use current directory, and ensure ends with / or \ as appropriate for your
# operating system
outputPath=
# groupName may not have spaces, corresponds to folder where outputs are stored
groupName=default
# filename must end in .h5, and can point to a relative or absolute path to the file
filename=redpytable.h5
# Amount of days to keep orphans in the queue based on size before permanent deletion
minorph=0.05
maxorph=7.
# nsec should be the maximum amount of time you want to process at once (in seconds)
# This depends on how much data your computer can hold in memory
nsec=3600
###### STATION PARAMETERS ######
# Number of stations
nsta=8
# List of stations in order of storage preference, NO SPACES
station=SEP,YEL,HSR,SHW,EDM,STD,JUN,SOS
# List of channels, etc., one per station even if they're all the same
channel=EHZ,EHZ,EHZ,EHZ,EHZ,EHZ,EHZ,EHZ
network=UW,UW,UW,UW,UW,UW,UW,UW
location=--,--,--,--,--,--,--,--
# Sampling rate to store all waveforms at (will resample if data do not match)
samprate=100.
# Bandpass filter (Hz)
fmin=1.
fmax=10.
# Frequency Index windows (should be within bounds of filtering above)
# Please note that these settings are ALSO used for teleseism filtering (below); if you
# change these from the default you will likely need to change 'telefi' or change
# 'teleok' to the same value as 'nsta' if you do not want to exclude teleseisms.
# FI lower bounds
filomin=1.
filomax=2.5
# FI upper bounds
fiupmin=5.
fiupmax=10.
###### DATA SOURCE ######
# Server accepts any FDSN web service that has waveforms (see obspy.clients.fdsn for full
# list or specify a server, e.g., server=fdsnws://http://service.iris.edu).
# It may also be a waveserver (e.g., server=waveserver://mazama.ess.washington.edu:16024),
# or a seedlink server (e.g., server=seedlink://rtserve.iris.washington.edu:18000) or
# server=file for local files in any format ObsPy can read. If using local files like SAC
# or miniSEED, please be sure your headers EXACLTY match what you use for the SCNL above,
# as this is what is used to associate station information and not the filename.
# Seedlink default port : 18000
# WaveServer default port : 16017
server=IRIS
# Port is only necessary if using a waveserver with legacy format
# Now deprecated, place port in server value, see above!
port=16017
# If using local files, define the path to the top directory where they exist,
# ending in / or \ as appropriate for your operating system. If there are files in any
# subdirectories within this directory, they will be found
searchdir=./
# You can specify a pattern for your files to reduce the files within the directory
# searched. For example, filepattern=2019.06.*.mseed if your files are miniSEED files
# named by date and you only want those from June 2019. Simple wildcarding is supported
# (i.e., * and ?, [] for ranges of values or lists) but not full regular expressions.
filepattern=*
###### TRIGGERING SETTINGS ######
# Choose the STALTA algorithm. Options include a selection of trigger algorithms
# supported by ObsPy's network coincidence trigger (suggested values for 'lwin' and 'swin'
# are in parentheses, followed by suggested values for 'trigon' and 'trigoff'):
# * 'classicstalta' : Classic STA/LTA (8., 0.7), (3., 2.)
# * 'recstalta' : Recursive STA/LTA (3., 8.), (2., 1.5)
# * 'delayedstalta' : Delayed STA/LTA (5., 10.), (5., 10,)
# For more details, please visit:
# https://docs.obspy.org/packages/autogen/obspy.core.trace.Trace.trigger.html#obspy.core.trace.Trace.trigger
# and https://docs.obspy.org/tutorial/code_snippets/trigger_tutorial.html
trigalg=classicstalta
# How many stations need to be triggered in order to be considered using a coincidence
# trigger? I like using just over half of nsta. Using less allows more triggers through,
# but that may not be real events or are of dubious quality. Requiring all stations to
# trigger is not advised unless you have 100% uptime on all stations.
nstaC=5
# STALTA trigger settings for long and short windows
lwin=8.
swin=0.7
# Higher trigon will reduce number of triggers, but may miss small events
trigon=3.
# trigoff should be slightly less than trigon
trigoff=2.
# Optional time offset
# Accepts a list (in same order as stations) of positive numbers corresponding to number
# of seconds to move the trace to the left. The way I think of this is as a list of
# arrival times of waves from an expected location. Can be used to quickly assess whether
# an event of interest is located near a known point (e.g., beneath the summit of a
# volcano or near the main shock in an aftershock sequence), or to bring arrivals from
# distant stations to be within a short correlation window. Be aware that the trigger time
# listed in the catalog will also be offset by the smallest value given. If used, should
# have values for each station. If not used, can just be 0.
# Example for default Mount St. Helens network might be:
# offset=0.00,0.17,0.49,0.60,0.50,0.82,1.09,1.05
offset=0
### CROSS-CORRELATION PARAMETERS ###
# Cross-correlation window length (in samples, 2^n is best)
winlen=1024
# Minimum correlation coefficient to be considered a repeater
cmin=0.7
# Number of stations this value must be exceeded on to be counted as a repeater. Like
# nstaC, I like to use about half.
ncor=4
### PLOTTING PARAMETERS ###
# Station number to use for plotting cores & amplitudes, beginning at 0, in station list
printsta=2
# Minimum number of members within a family needed to appear on the timeline plot
minplot=5
# Width (in days) of bins for histogram subplot (default is 1 day)
dybin=1
# Width (in hours) of bins for histogram subplot on 'recent' timeline (default is 1 hour)
hrbin=1.
# Width (in hours) of bars in occurrence plot (default is 1 hour) - if this is set for
# a day or longer, the color bar will be scaled automatically up to a maximum of 1000
# events per bin; recommend customizing this setting if your run extends >1 year
occurbin=1
# Width (in hours) of bars in recent occurrence plot (default is 1 hour)
recbin=1
# Number of days prior to last repeater to show on 'recent' timeline (default is 2 weeks)
recplot=14.
# Print a more verbose version of the catalog
verbosecatalog=False
# Annotation file (see included annotation.csv for format)
anotfile=annotation.csv
# Use 'global' or 'family' to define amplitude plot limits
amplims=global
### CHECK COMCAT FOR MATCHES ###
# If used, will search the ANSS Comprehensive Catalog for potential located seismicity
# that has arrival times consistent with the trigger time of repeaters, and if found,
# will write those matches to the end of the family HTML file with a map if at least one
# match is local to the study area. This can be time consuming for very large families!
# If the location of interest is in Northern California, it will also check NCEDC for
# located seismicity that may not have made it to ComCat.
checkComCat=False
# Station latitudes and longitudes (can also be a single value near the center of your
# network, here these match the stations above)
stalats=46.200210,46.209550,46.174280,46.193470,46.197170,46.237610,46.147060,46.243860
stalons=-122.190600,-122.188990,-122.180650,-122.236350,-122.151210,-122.223960,-122.152430,-122.137870
# Seconds of error allowable to be considered a match
serr = 5.
# Distance in degrees to search for local matches
locdeg = 0.5
# Distance in degrees to search for regional matches
regdeg = 2.
# Minimum magnitude to allow for regional matches
regmag = 2.5
# Minimum magnitude to allow for teleseismic matches
telemag = 4.5
# Maximum number of events to try to match, chosen in descending order by amplitude on
# 'printSta' station. Defaults to 0, which matches all. I like 50 as a tradeoff between
# matching everything in most clusters and limiting the amount of time spent trying to
# match everything in extremely large clusters with 1000+ members
matchMax = 0
### AUTOMATED SPIKE AND TELESEISM REMOVAL ###
# Spikes, harmonic noise, and teleseisms are automatically put into the 'junk' table and
# never considered as a repeater if they exceed the following thresholds. You can view the
# contents of this table with the plotJunk.py helper script and experiment with settings
# that might be more appropriate for your data. So far these settings have done well at
# not throwing away too many real events and not letting too many bad triggers through.
# Kurtosis of waveform in a window: ~80-100 is appropriate for 5 s window, ~130 for 15 s,
# and ~200 for 25 s
kurtwin=5.
kurtmax=80.
# Kurtosis of amplitude spectrum to eliminate calibration pulses with unnaturally harmonic
# signals--be careful not to set too low or you could eliminate real harmonic events
kurtfmax=150.
# Maximum ratio of outliers to total number of datapoints in trace
oratiomax=0.15
# Teleseisms are filtered based on frequency index (uses above settings)
# Some tests with the default windows suggest that using a cutoff of -1.0 on
# multiple stations should catch a decent number of teleseisms without throwing
# out too many local LF events
# Frequency index minimum threshold (above this is OK)
telefi = -1.0
# Number of stations allowed to pass that exceed threshold
teleok = 2