marcel/pylot
3
0
Fork 0
Code Issues 7 Pull Requests 1 Releases Wiki Activity

Bugfix in calcsourcespec: take into account special case of components 1, 2, and Z.

Browse Source
This commit is contained in:
Ludger Küperkoch 2016-10-27 16:15:38 +02:00
parent 9e5b3f6916
commit 0318fe7c32
1 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
pylot/core/analysis/magnitude.py
View File

@ -459,12 +459,14 @@ def calcsourcespec(wfstream, onset, vp, delta, azimuth, incidence,
"no rotation of components possible!") "no rotation of components possible!")
# instead, use component 3 # instead, use component 3
ldat = LQT.select(component="3") ldat = LQT.select(component="3")
if len(ldat == 0:
# maybe component z available
ldat = LQT.select(component="Z")
# integrate to displacement # integrate to displacement
# unrotated vertical component (for comparison) # unrotated vertical component (for comparison)
inttrz = signal.detrend(integrate.cumtrapz(zdat[0].data, None, dt)) inttrz = signal.detrend(integrate.cumtrapz(zdat[0].data, None, dt))
# rotated component Z => L # rotated component Z => L
print ldat
Ldat = signal.detrend(integrate.cumtrapz(ldat[0].data, None, dt)) Ldat = signal.detrend(integrate.cumtrapz(ldat[0].data, None, dt))
# get window after P pulse for # get window after P pulse for