From 3be78ec57603af2f81bda159c05a17cd61f538aa Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Ludger=20K=C3=BCperkoch?= Date: Tue, 19 Nov 2019 15:12:10 +0100 Subject: [PATCH] Some debugging. --- pylot/core/analysis/magnitude.py | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/pylot/core/analysis/magnitude.py b/pylot/core/analysis/magnitude.py index 99e7cc37..920ff0a0 100644 --- a/pylot/core/analysis/magnitude.py +++ b/pylot/core/analysis/magnitude.py @@ -625,7 +625,7 @@ def calcsourcespec(wfstream, onset, vp, delta, azimuth, incidence, "calcsourcespec: Determined corner frequency: %f Hz" % (w01, fc1)) # use of conventional fitting - [w02, fc2] = fitSourceModel(F, YYcor, Fcin, 2, verbosity) + [w02, fc2] = fitSourceModel(F, YYcor, Fcin, iplot, verbosity) # get w0 and fc as median of both # source spectrum fits @@ -652,7 +652,7 @@ def calcsourcespec(wfstream, onset, vp, delta, azimuth, incidence, plt.xlabel('Time since %s' % zdat[0].stats.starttime) plt.ylabel('Displacement [mm]') - if iplot >= 1: + if iplot > 1: plt.subplot(2, 1, 2) p1, = plt.loglog(f, Y.real, 'k') p2, = plt.loglog(F, YY.real) @@ -780,8 +780,8 @@ def fitSourceModel(f, S, fc0, iplot, verbosity=False): w0 = max(S) if verbosity: print( - "fitSourceModel: best fc: {0} Hz, best w0: {1} m/Hz".format(fc, w0)) - if iplot > 0: + "fitSourceModel: best fc: {0} Hz, best w0: {1} m/Hz".format(Fc, w0)) + if iplot >= 1: plt.figure() # iplot) plt.loglog(f, S, 'k') plt.loglog([f[0], Fc], [w0, w0], 'g')