[MSNoise] error in pws ' stack'

[刘志坤]

> -----Original Messages-----
> From: "孙杰" <1010131334 at cugb.edu.cn>
> Sent Time: Thursday, March 30, 2017
> To: MSNoise at mailman-as.oma.be
> Cc: 
> Subject: [MSNoise] error in pws  ' stack'
> 
> 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
> 
> _______________________________________________
> MSNoise mailing list
> MSNoise at mailman-as.oma.be
> http://mailman-as.oma.be/mailman/listinfo/msnoise