settings.cfg

[Settings]
# Above line is REQUIRED, do not edit!

### Configuration file for REDPy ###

# Most important settings are at the top
# Settings not explicitly defined will be assigned a default value

### RUN PARAMETERS ###

# Describe the run
title=REDPy Catalog
groupDesc=Default Test Run

# groupName may not have spaces, corresponds to folder where outputs are stored
groupName=default

# filename must end in .h5
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
channel=EHZ,EHZ,EHZ,EHZ,EHZ,EHZ,EHZ,EHZ
network=UW,UW,UW,UW,UW,UW,UW,UW
location=--,--,--,--,--,--,--,--
# If so, how many stations need to be triggered to count?
nstaC=5

# Sampling rate to store all waveforms at (will resample if data do not match)
samprate=100.

# Bandpass filter
fmin=1.
fmax=10.

# Frequency Index windows (should be within bounds of filtering)
filomin=1.
filomax=2.5
fiupmin=5.
fiupmax=10.

# Accepts any FDSN web service that has waveforms (see obspy.clients.fdsn for full list)
# Server may also be a waveserver (e.g., mazama.ess.washington.edu), SAC, or miniSEED
# If using SAC, please be sure your headers match what you use for the SCNL above!!
server=IRIS
# Port is only necessary if using a waveserver
port=16017
# If using SAC files, define the path to the directory where they exist, ending in /
sacdir=./
# Same with miniSEED
mseeddir=./


# STALTA trigger settings
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.


### 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
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 'active' subplot (default is 1 day)
dybin=1
# Width (in hours) of bins for 'active' subplot on 'recent' timeline (default is 1 hour)
hrbin=1.
# Number of days prior to last repeater to show on 'recent' timeline (default is 2 weeks)
recplot=14. 


### 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!

# Use or not use this feature (defaults to not used)
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 (default 0.5 degrees)
locdeg = 0.5
# Distance in degrees to search for regional matches (default 2 degrees)
regdeg = 2.
# Minimum magnitude to allow for regional matches (default M2.5)
regmag = 2.5
# Minimum magnitude to allow for teleseismic matches (default M4.5)
telemag = 4.5


### SPIKE REMOVAL ###

# Kurtosis of waveform; ~80-100 is appropriate for 5 s window, ~130 for 15 s,
# and ~200 for 25 s, updated testing suggests these values may be too high?
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.06

### TELESEISM REMOVAL ###

# 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
# more than one station 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.
# Number of stations allowed to pass that exceed threshold
teleok = 1