Commit 61ad6a4a authored by Sebastian Heimann's avatar Sebastian Heimann
Browse files

fix problem in 'grond check --waveforms'

parent c4011c92
...@@ -1526,7 +1526,6 @@ def solve(problem, ...@@ -1526,7 +1526,6 @@ def solve(problem,
xcandidates[icandidate, ipar] = v xcandidates[icandidate, ipar] = v
x = select_most_excentric( x = select_most_excentric(
xcandidates, xcandidates,
xhist[chains_i[jchoice, :], :], xhist[chains_i[jchoice, :], :],
...@@ -1908,16 +1907,20 @@ def check( ...@@ -1908,16 +1907,20 @@ def check(
freqlimits[3] = 0.5/deltat freqlimits[3] = 0.5/deltat
freqlimits = tuple(freqlimits) freqlimits = tuple(freqlimits)
trs_projected, trs_restituted, trs_raw = \ try:
ds.get_waveform( trs_projected, trs_restituted, trs_raw = \
target.codes, ds.get_waveform(
tmin=tmin+tobs_shift, target.codes,
tmax=tmax+tobs_shift, tmin=tmin+tobs_shift,
tfade=tfade, tmax=tmax+tobs_shift,
freqlimits=freqlimits, tfade=tfade,
deltat=deltat, freqlimits=freqlimits,
backazimuth=target.get_backazimuth_for_waveform(), deltat=deltat,
debug=True) backazimuth=target.get_backazimuth_for_waveform(),
debug=True)
except dataset.NotFound, e:
logger.warn(str(e))
continue
trs_projected = copy.deepcopy(trs_projected) trs_projected = copy.deepcopy(trs_projected)
trs_restituted = copy.deepcopy(trs_restituted) trs_restituted = copy.deepcopy(trs_restituted)
...@@ -2045,8 +2048,8 @@ def check( ...@@ -2045,8 +2048,8 @@ def check(
event.name or util.time_to_str(event.time), event.name or util.time_to_str(event.time),
str(e))) str(e)))
if show_waveforms: if show_waveforms:
trace.snuffle(trs_all, stations=ds.get_stations(), markers=markers) trace.snuffle(trs_all, stations=ds.get_stations(), markers=markers)
g_state = {} g_state = {}
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment