Skip to content

Commit 0c775dc

Browse files
nts345045nts345045
authored
Merge branch 'cascadiaquakes:main' into main
2 parents 9677fa1 + 68467b4 commit 0c775dc

File tree

9 files changed

+46293
-1
lines changed

9 files changed

+46293
-1
lines changed

‎_toc.yml‎

Lines changed: 4 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -26,4 +26,7 @@ parts:
2626
- caption: Waveform Data QC and Catalog Wrangling
2727
chapters:
2828
- file: notebooks/Nate/waveform_qc.ipynb
29-
- file: notebooks/Nate/catalog_management.ipynb
29+
- file: notebooks/Nate/catalog_management.ipynb
30+
- caption: Event Relocation
31+
chapters:
32+
- file: notebooks/Felix/run_hypoDD_plot.ipynb

‎notebooks/Felix/fernd.sta‎

Lines changed: 92 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,92 @@
1+
GHO.NC 39.0445 -123.5411 687.0000
2+
GBL.NC 39.1171 -123.7088 144.0000
3+
GHGB.NC 39.1275 -122.8235 858.0000
4+
GCWB.NC 39.1300 -123.0765 1077.0000
5+
GMR.NC 39.1332 -123.4679 370.0000
6+
GSR.NC 39.1790 -122.6362 975.0000
7+
79666.CE 39.1853 -123.2109 193.0000
8+
GNR.NC 39.2007 -123.6351 327.0000
9+
GWR.NC 39.2057 -123.3010 654.0000
10+
BARR.BK 39.3048 -123.1975 264.7000
11+
GTSB.BK 39.3117 -122.6035 1102.4000
12+
MNDO.BK 39.3142 -123.7451 179.4000
13+
GVA.NC 39.3736 -123.1280 586.0000
14+
79668.CE 39.3907 -123.3433 432.0000
15+
GTC.NC 39.3994 -123.5553 369.0000
16+
79298.CE 39.4240 -123.8070 34.0000
17+
DCMP.BK 39.4690 -123.3375 450.0000
18+
GHM.NC 39.4955 -122.9310 1456.0000
19+
KRK.NC 39.5632 -123.1837 1259.0000
20+
KFP.NC 39.6389 -123.4251 768.0000
21+
KCPB.NC 39.6863 -123.5824 1261.0000
22+
79046.CE 39.7379 -123.8149 8.0000
23+
79036.CE 39.8070 -123.2464 425.0000
24+
KIP.NC 39.8085 -123.4813 1359.0000
25+
B934.PB 39.8445 -123.6914 331.0000
26+
KBN.NC 39.8924 -123.1951 1365.0000
27+
RVIT.BK 39.9063 -123.3786 781.8000
28+
GROB.NC 39.9168 -122.6708 1301.0000
29+
79435.CE 39.9238 -123.7616 245.0000
30+
89688.CE 40.0331 -124.0769 28.0000
31+
LVR.NC 40.0394 -122.6725 943.0000
32+
B933.PB 40.0600 -123.9690 268.1000
33+
KLR.NC 40.0629 -123.3308 1025.0000
34+
1591.NP 40.1054 -123.7903 178.0000
35+
PRDS.BK 40.1056 -124.0446 786.5000
36+
BRIC.BK 40.1122 -123.8965 180.5000
37+
LRB.NC 40.1432 -122.5577 308.0000
38+
LPG.NC 40.1451 -122.6879 574.0000
39+
PETY.BK 40.1766 -122.7884 962.0000
40+
HUNT.BK 40.1894 -123.1456 1540.0000
41+
KMR.NC 40.2022 -123.7086 915.0000
42+
KPR.NC 40.2196 -124.0573 450.0000
43+
B049.PB 40.2375 -123.8197 117.0000
44+
KBU.NC 40.2455 -123.8532 577.0000
45+
B047.PB 40.2499 -124.1698 80.0000
46+
KCO.NC 40.2567 -124.2660 892.0000
47+
B932.PB 40.2750 -124.2212 52.0000
48+
PETL.BK 40.2884 -124.3426 49.6000
49+
GHOP.BK 40.3064 -123.9780 1028.9000
50+
LBP.NC 40.3167 -122.8819 1019.0000
51+
B046.PB 40.3220 -124.2843 44.0000
52+
BLAS.BK 40.3227 -123.5325 1679.3000
53+
89101.CE 40.3250 -124.2877 37.0000
54+
RBOW.BK 40.3347 -124.1114 1100.0000
55+
KPPB.NC 40.3459 -123.3632 1738.0000
56+
KMPB.NC 40.4172 -124.1208 938.0000
57+
KCR.NC 40.4264 -123.8206 873.0000
58+
B045.PB 40.4360 -123.9965 29.4000
59+
DLIK.BK 40.4554 -122.8942 1293.0000
60+
KCT.NC 40.4756 -124.3375 378.0000
61+
B935.PB 40.4787 -123.5732 696.7000
62+
1584.NP 40.4976 -124.2945 753.0000
63+
1584B.NP 40.4986 -124.2923 720.0000
64+
89462.CE 40.5036 -124.1010 43.0000
65+
HAYF.BK 40.5227 -123.2320 776.6000
66+
KCSB.NC 40.5384 -123.5140 1658.5000
67+
BJES.BK 40.5436 -123.8800 1050.7000
68+
1023.NP 40.5761 -124.2637 16.0000
69+
89255.CE 40.5922 -124.1463 16.0000
70+
1586.NP 40.6442 -124.2207 38.0000
71+
LSF.NC 40.6582 -122.5237 1042.0000
72+
KHBB.NC 40.6599 -123.2197 1864.0000
73+
1582.NP 40.6990 -124.2010 19.0000
74+
DMOR.BK 40.7097 -123.7033 1305.0000
75+
KNEE.BK 40.7152 -123.9011 892.2000
76+
1581.NP 40.7350 -124.2076 8.0000
77+
WEAV.BK 40.7458 -122.9840 1066.8000
78+
89781.CE 40.7834 -124.1435 35.0000
79+
JCC.BK 40.8175 -124.0295 27.2000
80+
KHMB.NC 40.8747 -123.7326 1501.0000
81+
1580.NP 40.8777 -124.0774 52.0000
82+
LGP.NC 40.9123 -122.8295 1265.0000
83+
89464.CE 40.9567 -124.1164 25.0000
84+
TRIN.BK 40.9741 -122.7050 935.0000
85+
KTD.NC 41.0539 -124.1506 101.0000
86+
LBK.NC 41.0838 -122.6673 1304.0000
87+
KRP.NC 41.1577 -124.0233 813.0000
88+
GUMB.BK 41.2130 -122.5094 1849.8000
89+
KOM.NC 41.2787 -123.4531 1844.0000
90+
99700.CE 41.2906 -124.0569 10.0000
91+
SCOT.BK 41.3920 -122.8538 960.0000
92+
LPK.NC 41.4553 -122.5437 1880.0000

‎notebooks/Felix/hypoDD.inp‎

Lines changed: 108 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,108 @@
1+
hypoDD_2
2+
* RELOC.INP:
3+
*--- input file selection
4+
5+
*
6+
*catalog P diff times:
7+
dt.ct
8+
*
9+
* event file:
10+
event.sel
11+
*
12+
* station file:
13+
station.sel
14+
*
15+
*--- output file selection
16+
* original locations:
17+
hypoDD.loc
18+
* relocations:
19+
hypoDD.reloc
20+
* station information:
21+
hypoDD.sta
22+
* residual information:
23+
hypoDD.res
24+
* source paramater information:
25+
hypoDD.src
26+
*
27+
*
28+
*--- data type selection:
29+
* IDAT: 1= cross corr; 2= catalog; 3= cross & cat
30+
* IPHA: 1= P; 2= S; 3= P&S
31+
* DIST:max dist [km] between cluster centroid and station
32+
* IDAT IPHA MAXDIST
33+
2 3 150
34+
*
35+
*--- event clustering:
36+
* OBSCC: min # of obs/pair for crosstime data (0= no clustering)
37+
* OBSCT: min # of obs/pair for network data (0= no clustering)
38+
* MINDIST: min pair-station distance (-9= not active)
39+
* MAXDIST: max pair-station distance (-9= not active)
40+
* MAXGAP: max azimuthal gap (-9= not active)
41+
* OBSCC OBSCT MINDIST MAXDIST MAXGAP
42+
0 8 -999 -999 -999
43+
*
44+
*--- solution control:
45+
* ISTART: 1 = from single source; 2 = from network sources
46+
* ISOLV: 1 = SVD, 2=lsqr
47+
* IAQ: remove airquakes: 0=NO, 1=YES
48+
* NSET: number of sets of iteration with specifications following
49+
* ISTART ISOLV IAQ NSET
50+
2 2 0 3
51+
*
52+
*--- data weighting and re-weighting:
53+
* NITER: last iteration to used the following weights
54+
* WTCCP, WTCCS: weight cross P, S
55+
* WTCTP, WTCTS: weight catalog P, S
56+
* WRCC, WRCT: Cutoff threshold for correlation/catalog data outliers.
57+
* 0<wr<1: In seconds.
58+
* wr≥1: Factor to multiply standard deviation
59+
* (dynamic cutoff).
60+
* WDCC, WDCT: max dist [km] between cross, catalog linked pairs
61+
* DAMP: damping (for lsqr only)
62+
* --- CROSS DATA ----- ----CATALOG DATA ----
63+
* NITER WTCCP WTCCS WRCC WDCC WTCTP WTCTS WRCT WDCT DAMP
64+
5 -9 -9 -9 -9 1.0 0.5 10 10 100
65+
5 -9 -9 -9 -9 1.0 0.5 8 8 100
66+
5 -9 -9 -9 -9 1.0 0.5 6 6 100
67+
*
68+
*--- FORWARD MODEL:
69+
* IMOD 0= 1D local, layered, P velocity w/ fixed vp/vs ratio (hypoDD1.0);
70+
* 1= 1D local, layered, P velocity w/ variable vp/vs ratio
71+
* 9= 3D local model (simulps raytracer)
72+
1
73+
*
74+
*--- 1D model, fixed vp/vs ratio:
75+
* NLAY: number of model layers
76+
* RATIO: vp/vs ratio
77+
* TOP: depths of top of layer (km)
78+
* VEL: layer velocities (km/s)
79+
* IMOD
80+
* NLAY RATIO
81+
* MOD_TOP , MOD_V
82+
* 21 1.73
83+
*
84+
*--- 1D model, variable vp/vs ratio:
85+
* TOP: depths of top of layer (km)
86+
* VEL: layer velocities (km/s), end w/ -9
87+
* RATIO: layer ratios, end w/ -9
88+
* IMOD
89+
**17 28 MEN MENDOCINO
90+
0.0000 0.6000 1.2000 1.8000 2.4000 3.0000 3.6000 4.2000 4.8000 5.4000 6.0000 8.0000 10.0000 12.0000 14.0000 16.0000 18.0000 20.0000 22.0000 24.0000 26.0000
91+
3.5000 3.7214 3.9429 4.1643 4.3857 4.6071 4.8286 5.0500 5.1124 5.1747 5.2371 5.4449 5.6528 5.8607 6.0685 6.2764 6.4843 6.6921 6.9000 7.9000 7.9000
92+
1.7300 1.7300 1.7300 1.7300 1.7300 1.7300 1.7300 1.7300 1.7300 1.7300 1.7300 1.7300 1.7300 1.7300 1.7300 1.7300 1.7300 1.7300 1.7300 1.7300 1.7300
93+
*
94+
*--- Local 3D model:
95+
* IMOD
96+
* 3DMOD name of 3D velocity model
97+
* LAT_3D LON_3D ROT_3D origin of 3D velocity model and rotation of
98+
* cart coord. system (pos= anti clockwise)
99+
* raytracing parameters
100+
* IPHA NDIP ISKIP SCALE1 SCALE2 XFAX TLIM NITPB simulps ray tracing parameters
101+
*
102+
*
103+
*--- event selection:
104+
* CID: cluster to be relocated (0 = all)
105+
* ID: cuspids of event to be relocated (8 per line)
106+
* CID
107+
1
108+
* ID

0 commit comments

Comments
 (0)