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

Implemented read_metadata.py in calcsourcespec.py (was this intended?).

Browse Source
This commit is contained in:
Ludger Küperkoch 2016-09-20 15:26:55 +02:00
parent e4c1912896
commit a95a33abc9
1 changed files with 5 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -6,6 +6,7 @@ Created autumn/winter 2015.
:author: Ludger Küperkoch / MAGS2 EP3 working group :author: Ludger Küperkoch / MAGS2 EP3 working group
""" """
import pdb
import matplotlib.pyplot as plt import matplotlib.pyplot as plt
import numpy as np import numpy as np
from obspy.core import Stream, UTCDateTime from obspy.core import Stream, UTCDateTime
@ -14,7 +15,7 @@ from pylot.core.util.utils import getPatternLine
from scipy.optimize import curve_fit from scipy.optimize import curve_fit
from scipy import integrate, signal from scipy import integrate, signal
from pylot.core.io.data import Data from pylot.core.io.data import Data
from pylot.core.util.dataprocessing import restitute_data from pylot.core.util.dataprocessing import restitute_data, read_metadata
from pylot.core.util.utils import common_range from pylot.core.util.utils import common_range
@ -306,7 +307,7 @@ def calcMoMw(wfstream, w0, rho, vp, delta, inv):
return Mo, Mw return Mo, Mw
def calcsourcespec(wfstream, onset, metadata, vp, delta, azimuth, incidence, def calcsourcespec(wfstream, onset, invdir, vp, delta, azimuth, incidence,
qp, iplot): qp, iplot):
''' '''
Subfunction to calculate the source spectrum and to derive from that the plateau Subfunction to calculate the source spectrum and to derive from that the plateau
@ -321,7 +322,7 @@ def calcsourcespec(wfstream, onset, metadata, vp, delta, azimuth, incidence,
:param: onset, P-phase onset time :param: onset, P-phase onset time
:type: float :type: float
:param: inventory, path/name of inventory or dataless-SEED file :param: invdir, path to inventory or dataless-SEED file
:type: string :type: string
:param: vp, Vp-wave velocity :param: vp, Vp-wave velocity
@ -354,7 +355,7 @@ def calcsourcespec(wfstream, onset, metadata, vp, delta, azimuth, incidence,
data = Data() data = Data()
wf_copy = wfstream.copy() wf_copy = wfstream.copy()
invtype, inventory = metadata invtype, inventory = read_metadata(invdir)
[cordat, restflag] = restitute_data(wf_copy, invtype, inventory) [cordat, restflag] = restitute_data(wf_copy, invtype, inventory)
if restflag is True: if restflag is True: