forked from NoisyLeon/NoisePy
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathquake_Alaska.py
52 lines (47 loc) · 2.02 KB
/
quake_Alaska.py
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
import quakedbase
import numpy as np
import timeit
import matplotlib.pyplot as plt
import pyaftan
#------------------------------
# Initialize ASDF dataset
#------------------------------
dset = quakedbase.quakeASDF('/work1/leon/ALASKA_work/ASDF_data/surf_Alaska.h5')
# # # dset.add_quakeml('/scratch/summit/life9360/ALASKA_work/quakeml/alaska_2017_aug.ml')
# # # print dset.events[0]
#------------------------------
# Retrieving earthquake catalog
#------------------------------
# ISC catalog
# dset.get_events(startdate='1991-01-01', enddate='2015-02-01', Mmin=5.5, magnitudetype='mb', gcmt=True)
# gcmt catalog
# dset.get_events(startdate='1991-01-01', enddate='2017-08-31', Mmin=5.5, magnitudetype='mb', gcmt=True)
# gcmt catalog, updated on 20180831
# dset.get_events(startdate='2017-09-01', enddate='2018-04-30', Mmin=5.5, magnitudetype='mb', gcmt=True)
# gcmt catalog, updated on 20190218
# dset.get_events(startdate='2018-05-01', enddate='2018-10-30', Mmin=5.5, magnitudetype='mb', gcmt=True)
#------------------------------
# Getting station information
#------------------------------
# # # dset.get_stations(channel='LHZ', minlatitude=52., maxlatitude=72.5, minlongitude=-172., maxlongitude=-122.)
#------------------------------
# Download/Read data
#------------------------------
# # # t1=timeit.default_timer()
# #
# # # dset.read_surf_waveforms_DMT(datadir='/scratch/summit/life9360/ALASKA_work/surf_19950101_20170831', verbose=False)
# dset.get_surf_waveforms(startdate='2018-05-01', verbose=False)
#------------------------------
# aftan analysis
#------------------------------
# dset.quake_prephp(outdir='/work1/leon/ALASKA_work/quake_working_dir/pre_disp')
# inftan = pyaftan.InputFtanParam()
# inftan.tmax = 100.
# inftan.tmin = 5.
# # #
# # inftan.tmax = 120.
# # inftan.tmin = 20.
# dset.quake_aftan(prephdir='/work1/leon/ALASKA_work/quake_working_dir/pre_disp_R', inftan=inftan)
# pers = np.append( np.arange(11.)*2.+20., np.arange(10.)*5.+45.)
# # dset.interp_disp(verbose=True, pers=pers)
# dset.quake_get_field(pers=pers)