Source code for download_events_script

#!/usr/bin/python3
# ##!/usr/bin/python2.7
## #!/usr/bin/env python
import sys
from obspy import UTCDateTime
#requires python version2.7
#Requires fdsnwsscripts for fdsn=True (pip install fdsnwsscripts)
#arclink_fetch client no longer works
#Create a working directory
#Update paths and options (default is to download example events from example stations)
#Run script
#View readme for further instructions

#############Main parameters
#Working directory where data will be saved (requires the trailing "/")

[docs]
wd='/home/kmichall/Downloads/obs/sac/'


#Directory containing the functions document

[docs]
fd="/home/kmichall/Desktop/Codes/github/rfmpy/rfmpy/run_examples/download_data/get_wav_data/download_events_functions.py"


###########EVENTS###################
#Path to the events csv file
#eventcsv='./0events/example_events.csv'
#eventcsv='' 

[docs]
eventcsv='./0events/USGS-fjg_rainfall_niwa.csv'

##OR using an available catalog

[docs]
useclient=True


[docs]
cl="USGS"

#APS starttime=UTCDateTime("2009-01-01")
#APS endtime=UTCDateTime("2018-12-01")

[docs]
starttime=UTCDateTime("2017-05-01")


[docs]
endtime=UTCDateTime("2018-03-01")

###Create unique event names YYYY.DDD.HH.MM.SS.NW.STA.C.SAC

[docs]
cnames=True

####################################


###########STATIONS###################
#Client (routing clients are "iris-federator" and "eida-routing")

[docs]
client_name="eida-routing"

#Is this a routing client?

[docs]
rclient=True

#Path to the stations csv file. Note: A "*" entry means download all stations available for that network (_ALPARRAY if no network name is provided)
#stationcsv='./0stations/PA-test.txt'

[docs]
stationcsv='/home/kmichall/Desktop/Codes/github/rfmpy/rfmpy/run_examples/download_data/get_wav_data/Z3_obs.csv'

network="_ALPARRAY"
#Set c_inv equal to an obspy inventory (e.g. from read_inventory or get_stations) if you want to use your own station inventory rather than download a new one

[docs]
c_inv=[]

#Or get stations from client in lat/longbox (leave True if using c_inv)

[docs]
usestatclient=False

#APS usestatclient=True   #APS


[docs]
network="Z3"


[docs]
minlatitude=-90


[docs]
minlongitude=-180


[docs]
maxlatitude=90


[docs]
maxlongitude=180


[docs]
includeZS=False #Include the ZS network (defaults as part of _ALPARRAY)

#ZNE rotations (correct station azimuths with .csv file?)

[docs]
znepath='./rotations.csv'


[docs]
rotzne=False

####################################


#Phase (see obspy for detailed options)

[docs]
phase="P"


#minimum magnitude (events below this value will be completely ignored)

[docs]
minmag=5.5


#minimum epicentral distance

[docs]
minepi=30


#maximum epicentral distance

[docs]
maxepi=95


#Window start time (seconds relative to predicted P-wave arrival time)

[docs]
ws=-120  


#Window end time (seconds relative to predicted P-wave arrival time)

[docs]
we=120


#Sort by events by either "station" or "event"

[docs]
sortby="event" 


#High value (Hz) for 2-pole butterworth bandpass filter (zerophase) (No filter if set to None)

[docs]
fhi=None #


# Low value (Hz) for 2-pole butterworth bandpass filter (zerophase) (no filter if set to None)
#APS flo=0.03 

[docs]
flo=None


#Downsample to 20Hz? (using decimate unless the actual sampling rate isn't an integer multiple of 20)

[docs]
downsample=True


#Rotate coordinates? "NE->RT" or "ZNE->LQT" or None (Warning: does not correct ZNE misorientations listed in metadata)

[docs]
rotrt= None


#Earth model for predicting P-wave travel time (see obspy for options)

[docs]
model="iasp91" 


#Mode for running; options include "all" to download everything, "continue" to continue, or "retry" to retry failed events
#APS mode="continue"

[docs]
mode="all"



###########
##Source functions and start download
exec(open(fd).read())