6 Apr
2017
6 Apr
'17
12:42 p.m.
Hi all,
I have the following error trying to calculate the DTT with the command
(master version)
msnoise compute_dtt -i 999
for 2 stations with autocorrelations. It seems that the script stacks in a
some days (eg. 2011-02-27).
Any idea?
Cheers,
Christos
-----------------------------------------------
2017-04-06 10:37:45 [INFO] *** Starting: Compute DT/T ***
2017-04-06 10:37:45 [DEBUG] Found 2229 updated days
2017-04-06 10:37:45 [INFO] Loading mov=5 days for filter=02
2017-04-06 10:37:45 [DEBUG] Processing 2011-02-03 - 02 - 05 mov
/home/noise/anaconda2/lib/python2.7/site-packages/msnoise/s06compute_dtt.py:140:
RuntimeWarning: invalid value encountered in double_scalars
wavg = (d * w).sum() / w.sum()
/home/noise/anaconda2/lib/python2.7/site-packages/msnoise/s06compute_dtt.py:142:
RuntimeWarning: invalid value encountered in double_scalars
wstd = np.sqrt(np.sum(w * (d - wavg) ** 2) / ((N - 1) * np.sum(w) / N))
/home/noise/anaconda2/lib/python2.7/site-packages/msnoise/s06compute_dtt.py:340:
RuntimeWarning: invalid value encountered in less_equal
errindex = np.where(errArray[i] <= maxErr)[0]
/home/noise/anaconda2/lib/python2.7/site-packages/msnoise/s06compute_dtt.py:341:
RuntimeWarning: invalid value encountered in less_equal
dtindex = np.where(np.abs(dtArray[i]) <= maxDt)[0]
2017-04-06 10:37:45 [DEBUG] 2011-02-03: exporting: 4 pairs
2017-04-06 10:37:45 [DEBUG] Processing 2011-02-04 - 02 - 05 mov
2017-04-06 10:37:45 [DEBUG] 2011-02-04: exporting: 4 pairs
2017-04-06 10:37:45 [DEBUG] Processing 2011-02-05 - 02 - 05 mov
2017-04-06 10:37:45 [DEBUG] 2011-02-05: exporting: 4 pairs
2017-04-06 10:37:45 [DEBUG] Processing 2011-02-06 - 02 - 05 mov
2017-04-06 10:37:45 [DEBUG] 2011-02-06: exporting: 4 pairs
2017-04-06 10:37:45 [DEBUG] Processing 2011-02-07 - 02 - 05 mov
2017-04-06 10:37:45 [DEBUG] 2011-02-07: exporting: 4 pairs
2017-04-06 10:37:45 [DEBUG] Processing 2011-02-08 - 02 - 05 mov
2017-04-06 10:37:45 [DEBUG] 2011-02-08: exporting: 4 pairs
2017-04-06 10:37:45 [DEBUG] Processing 2011-02-09 - 02 - 05 mov
/home/noise/anaconda2/lib/python2.7/site-packages/msnoise/s06compute_dtt.py:142:
RuntimeWarning: divide by zero encountered in double_scalars
wstd = np.sqrt(np.sum(w * (d - wavg) ** 2) / ((N - 1) * np.sum(w) / N))
2017-04-06 10:37:46 [DEBUG] 2011-02-09: exporting: 4 pairs
2017-04-06 10:37:46 [DEBUG] Processing 2011-02-10 - 02 - 05 mov
2017-04-06 10:37:46 [DEBUG] 2011-02-10: exporting: 4 pairs
2017-04-06 10:37:46 [DEBUG] Processing 2011-02-11 - 02 - 05 mov
2017-04-06 10:37:46 [DEBUG] 2011-02-11: exporting: 4 pairs
2017-04-06 10:37:46 [DEBUG] Processing 2011-02-13 - 02 - 05 mov
2017-04-06 10:37:46 [DEBUG] 2011-02-13: exporting: 4 pairs
2017-04-06 10:37:46 [DEBUG] Processing 2011-02-14 - 02 - 05 mov
2017-04-06 10:37:46 [DEBUG] 2011-02-14: exporting: 4 pairs
2017-04-06 10:37:46 [DEBUG] Processing 2011-02-15 - 02 - 05 mov
2017-04-06 10:37:46 [DEBUG] 2011-02-15: exporting: 4 pairs
2017-04-06 10:37:46 [DEBUG] Processing 2011-02-16 - 02 - 05 mov
2017-04-06 10:37:46 [DEBUG] 2011-02-16: exporting: 4 pairs
2017-04-06 10:37:46 [DEBUG] Processing 2011-02-17 - 02 - 05 mov
2017-04-06 10:37:46 [DEBUG] 2011-02-17: exporting: 4 pairs
2017-04-06 10:37:46 [DEBUG] Processing 2011-02-18 - 02 - 05 mov
2017-04-06 10:37:46 [DEBUG] 2011-02-18: exporting: 4 pairs
2017-04-06 10:37:46 [DEBUG] Processing 2011-02-19 - 02 - 05 mov
2017-04-06 10:37:46 [DEBUG] 2011-02-19: exporting: 4 pairs
2017-04-06 10:37:46 [DEBUG] Processing 2011-02-20 - 02 - 05 mov
2017-04-06 10:37:46 [DEBUG] 2011-02-20: exporting: 4 pairs
2017-04-06 10:37:46 [DEBUG] Processing 2011-02-21 - 02 - 05 mov
2017-04-06 10:37:46 [DEBUG] 2011-02-21: exporting: 4 pairs
2017-04-06 10:37:46 [DEBUG] Processing 2011-02-22 - 02 - 05 mov
2017-04-06 10:37:46 [DEBUG] 2011-02-22: exporting: 4 pairs
2017-04-06 10:37:46 [DEBUG] Processing 2011-02-23 - 02 - 05 mov
2017-04-06 10:37:46 [DEBUG] 2011-02-23: exporting: 4 pairs
2017-04-06 10:37:46 [DEBUG] Processing 2011-02-24 - 02 - 05 mov
2017-04-06 10:37:46 [DEBUG] 2011-02-24: exporting: 4 pairs
2017-04-06 10:37:46 [DEBUG] Processing 2011-02-25 - 02 - 05 mov
2017-04-06 10:37:46 [DEBUG] 2011-02-25: exporting: 4 pairs
2017-04-06 10:37:46 [DEBUG] Processing 2011-02-26 - 02 - 05 mov
2017-04-06 10:37:46 [DEBUG] 2011-02-26: exporting: 4 pairs
2017-04-06 10:37:46 [DEBUG] Processing 2011-02-27 - 02 - 05 mov
Traceback (most recent call last):
File "/home/noise/anaconda2/bin/msnoise", line 11, in <module>
load_entry_point('msnoise==1.5a0', 'console_scripts',
'msnoise')()
File
"/home/noise/anaconda2/lib/python2.7/site-packages/msnoise/scripts/msnoise.py",
line 708, in run
cli(obj={})
File "/home/noise/anaconda2/lib/python2.7/site-packages/click/core.py",
line 722, in __call__
return self.main(*args, **kwargs)
File "/home/noise/anaconda2/lib/python2.7/site-packages/click/core.py",
line 697, in main
rv = self.invoke(ctx)
File "/home/noise/anaconda2/lib/python2.7/site-packages/click/core.py",
line 1066, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/home/noise/anaconda2/lib/python2.7/site-packages/click/core.py",
line 895, in invoke
return ctx.invoke(self.callback, **ctx.params)
File "/home/noise/anaconda2/lib/python2.7/site-packages/click/core.py",
line 535, in invoke
return callback(*args, **kwargs)
File
"/home/noise/anaconda2/lib/python2.7/site-packages/msnoise/scripts/msnoise.py",
line 421, in compute_dtt
main(interval)
File
"/home/noise/anaconda2/lib/python2.7/site-packages/msnoise/s06compute_dtt.py",
line 356, in main
intercept_origin=False)
File "/home/noise/anaconda2/lib/python2.7/site-packages/msnoise/api.py",
line 1508, in linear_regression
xtol=1e-20)
File
"/home/noise/anaconda2/lib/python2.7/site-packages/scipy/optimize/minpack.py",
line 740, in curve_fit
raise RuntimeError("Optimal parameters not found: " + errmsg)
RuntimeError: Optimal parameters not found: xtol=0.000000 is too small, no
further improvement in the approximate
solution is possible.