-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathfevlp_backup.m
More file actions
186 lines (156 loc) · 3.41 KB
/
fevlp_backup.m
File metadata and controls
186 lines (156 loc) · 3.41 KB
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
*sac macro file for automatically windowing and filtering all the stations
*in a directory for a particular event. Window parameters should be based
*on previous examination of filtered records at a near and far station of the
*array - then linearly interpolate windows based on distance
*
*To avoid useless files, should previously have labelled files that aren't
*to be used with prefix b for bad, as in b.D30.N07.HOT07.Z
*
*
echo on c m p
setbb ffile $1
rh %ffile
setbb dist &1,DIST
setbb phv20 3.0
evaluate to cute %dist / %phv20 + 1800
if %cute GT 3000
evaluate to cutb %cute - 3000
else
evaluate to cutb 100
endif
*do file list D05.N85.SUR.T
do file wild D*HZ*SAC
setbb fn1 $file
setbb fn2 %fn1%.w
r %fn1
setbb endtime &1,E
setbb npts &1,NPTS
if %cutb LT %endtime
cut %cutb %cute
r %fn1
w %fn2
cut off
xlim off
rmean
taper
do BANDP from 1 to 18 by 1
r %fn2
* p1
if $BANDP EQ 1
setbb SUFIX 'f050'
bp co .045 .055 n 4 p 2
w %fn1%.%SUFIX%
envelope
* p1
w %fn1%.%SUFIX%.evlp
elseif $BANDP EQ 2
setbb SUFIX 'f045'
bp co .04 .05 n 4 p 2
w %fn1%.%SUFIX%
envelope
w %fn1%.%SUFIX%.evlp
elseif $BANDP EQ 3
bp co .035 .045 n 4 p 2
setbb SUFIX 'f040'
w %fn1%.%SUFIX%
envelope
w %fn1%.%SUFIX%.evlp
elseif $BANDP EQ 4
bp co .032 .042 n 4 p 2
setbb SUFIX 'f037'
w %fn1%.%SUFIX%
envelope
w %fn1%.%SUFIX%.evlp
elseif $BANDP EQ 5
bp co .028 .038 n 4 p 2
setbb SUFIX 'f033'
w %fn1%.%SUFIX%
envelope
w %fn1%.%SUFIX%.evlp
elseif $BANDP EQ 6
bp co .024 .034 n 4 p 2
setbb SUFIX 'f029'
w %fn1%.%SUFIX%
envelope
w %fn1%.%SUFIX%.evlp
elseif $BANDP EQ 7
bp co .02 .03 n 4 p 2
setbb SUFIX 'f025'
w %fn1%.%SUFIX%
envelope
w %fn1%.%SUFIX%.evlp
elseif $BANDP EQ 8
bp co .017 .027 n 4 p 2
setbb SUFIX 'f022'
w %fn1%.%SUFIX%
envelope
w %fn1%.%SUFIX%.evlp
elseif $BANDP EQ 9
bp co .015 .025 n 4 p 2
setbb SUFIX 'f020'
w %fn1%.%SUFIX%
envelope
w %fn1%.%SUFIX%.evlp
elseif $BANDP EQ 10
bp co .012 .022 n 4 p 2
setbb SUFIX 'f017'
w %fn1%.%SUFIX%
envelope
w %fn1%.%SUFIX%.evlp
elseif $BANDP EQ 11
bp co .01 .02 n 4 p 2
setbb SUFIX 'f015'
w %fn1%.%SUFIX%
envelope
w %fn1%.%SUFIX%.evlp
elseif $BANDP EQ 12
bp co .008 .018 n 4 p 2
setbb SUFIX 'f013'
w %fn1%.%SUFIX%
envelope
w %fn1%.%SUFIX%.evlp
elseif $BANDP EQ 13
bp co .0065 .0165 n 4 p 2
setbb SUFIX 'f012'
w %fn1%.%SUFIX%
envelope
w %fn1%.%SUFIX%.evlp
elseif $BANDP EQ 14
bp co .005 .015 n 4 p 2
setbb SUFIX 'f010'
w %fn1%.%SUFIX%
envelope
w %fn1%.%SUFIX%.evlp
elseif $BANDP EQ 15
bp co .004 .014 n 4 p 2
setbb SUFIX 'f009'
w %fn1%.%SUFIX%
envelope
w %fn1%.%SUFIX%.evlp
elseif $BANDP EQ 16
bp co .003 .013 n 4 p 2
setbb SUFIX 'f008'
w %fn1%.%SUFIX%
envelope
w %fn1%.%SUFIX%.evlp
elseif $BANDP EQ 17
bp co .002 .012 n 4 p 2
setbb SUFIX 'f007'
w %fn1%.%SUFIX%
envelope
w %fn1%.%SUFIX%.evlp
elseif $BANDP EQ 18
bp co .001 .011 n 4 p 2
setbb SUFIX 'f006'
w %fn1%.%SUFIX%
envelope
w %fn1%.%SUFIX%.evlp
endif
enddo
sc rm %fn2
elseif %npts GT 5000
evaluate to cutb %cutb - 1000
else
message "failed"
enif
enddo