[MSNoise] error in pws ' stack'

[Thomas Lecocq]

Hi,

It's most probably because you don't have any RR CCF computed. I'll 
check, but I think this is corrected in the next version. In the 
meantime, you could disable stacking RR to test if all pass.

Thomas


Le 30/03/2017 à 11:33, 孙杰 a écrit :
> 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