30 Mar
2017
30 Mar
'17
11:33 a.m.
Hi, Thomas
I'm Sun Jie and i am a python beginner.Following your suggestion, I have solved the
'pws compute_cc' problem ,many thanks!
However, i have another question when i run 'msnoise stack -r -m ' I
'config' the stack_method=pws ,pws power=2 ,
pws_timegate=10.0,autocorr=Y,ZZ=Y....,
There is following error when I 'stack'. How can i solve this problem? many
thanks!
$ msnoise stack -r -m
rusult in
Lets STACK !
2017-03-30 16:46:58 [DEBUG] Starting the ref stack
[]
2017-03-30 16:46:58 [DEBUG] Processing YP_NE49:YP_NE49-ZZ-1
2017-03-30 16:46:58 [DEBUG] Found 278 updated days
2017-03-30 16:46:58 [DEBUG] New Data for YP_NE49:YP_NE49-ZZ-1
/home/sj/anaconda2/lib/python2.7/site-packages/obspy/io/sac/util.py:370: UserWarning: Old
header has invalid reftime.
warnings.warn(msg)
2017-03-30 16:47:01 [DEBUG] Processing YP_NE49:YP_NE49-RR-1
2017-03-30 16:47:01 [DEBUG] Found 278 updated days
2017-03-30 16:47:01 [DEBUG] New Data for YP_NE49:YP_NE49-RR-1
Traceback (most recent call last):
File "/home/sj/anaconda2/bin/msnoise", line 11, in <module>
sys.exit(run())
File
"/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/scripts/msnoise.py",
line 614, in run
cli(obj={})
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 722, in
__call__
return self.main(*args, **kwargs)
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 697, in
main
rv = self.invoke(ctx)
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 1066,
in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 895, in
invoke
return ctx.invoke(self.callback, **ctx.params)
File "/home/sj/.local/lib/python2.7/site-packages/click/core.py", line 535, in
invoke
return callback(*args, **kwargs)
File
"/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/scripts/msnoise.py",
line 301, in stack
main('ref', interval)
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/s04stack.py",
line 209, in main
db, sta1, sta2, filterid, components, datelist, format=format)
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/api.py", line
1020, in get_results
corr = stack(session, stack_data)
File "/home/sj/anaconda2/lib/python2.7/site-packages/msnoise/api.py", line
978, in stack
coh = 1. / data.shape[0] * np.abs(phasestack)
ZeroDivisionError: float division by zero