[MSNoise] Problems with s03compute_cc

Rueda Nuñez Juán José jjrueda at fomento.es
Thu Feb 19 16:17:32 UTC 2015


Dear Thomas.
In this moment the filter configuration I use :

low : 0.01 or 0.1
mwcs_low : 0.12
high : 0.98
mwcs_high : 1.0
rms_treshold : 0
mwcs_len : 10
mwcs_step : 5

And the result is:

2015-02-19 17:03:04,831 [DEBUG] ES.CGIN.Z Reading 1 Files
2015-02-19 17:03:06,029 [DEBUG] ES.CGIN.Z Merging Stream
2015-02-19 17:03:06,030 [DEBUG] ES.CGIN.Z Slicing Stream to 2015-01-05T00:00:00.000000Z:2015-01-05T23:59:59.990000Z
2015-02-19 17:03:06,109 [DEBUG] ES.CGIN.Z Lowpass at 8.00 Hz
2015-02-19 17:03:08,506 [DEBUG] ES.CGIN.Z Highpass at 0.01 Hz
2015-02-19 17:03:09,214 [DEBUG] ES.CGIN.Z Decimate by a factor of 5
2015-02-19 17:03:09,263 [DEBUG] ES.CGUI.Z Reading 1 Files
2015-02-19 17:03:10,589 [DEBUG] ES.CGUI.Z Merging Stream
2015-02-19 17:03:10,590 [DEBUG] ES.CGUI.Z Slicing Stream to 2015-01-05T00:00:00.000000Z:2015-01-05T23:59:59.990000Z
2015-02-19 17:03:10,590 [DEBUG] ES.CGUI.Z Lowpass at 8.00 Hz
2015-02-19 17:03:11,266 [DEBUG] ES.CGUI.Z Highpass at 0.01 Hz
2015-02-19 17:03:11,972 [DEBUG] ES.CGUI.Z Decimate by a factor of 5
2015-02-19 17:03:12,019 [DEBUG] ES.CHIE.Z Reading 1 Files
2015-02-19 17:03:13,379 [DEBUG] ES.CHIE.Z Merging Stream
2015-02-19 17:03:13,379 [DEBUG] ES.CHIE.Z Slicing Stream to 2015-01-05T00:00:00.000000Z:2015-01-05T23:59:59.980000Z
2015-02-19 17:03:13,380 [DEBUG] ES.CHIE.Z Lowpass at 8.00 Hz
2015-02-19 17:03:13,717 [DEBUG] ES.CHIE.Z Highpass at 0.01 Hz
2015-02-19 17:03:14,066 [DEBUG] ES.CHIE.Z Decimate by a factor of 5
Traceback (most recent call last):
  File "s03compute_cc.py", line 285, in <module>
    tramef_Z[istation] = trame
ValueError: could not broadcast input array from shape (864000) into shape (1728000)

Thanks.
Juan


-----Mensaje original-----
De: msnoise-bounces at mailman-as.oma.be [mailto:msnoise-bounces at mailman-as.oma.be] En nombre de Thomas Lecocq
Enviado el: jueves, 19 de febrero de 2015 9:19
Para: msnoise at mailman-as.oma.be
Asunto: Re: [MSNoise] Problems with s03compute_cc

Hi,

what are the frequency bounds for the configured filter(s) ?

Thomas

Le 18/02/2015 18:20, Rueda Nuñez Juán José a écrit :
> Dear all,
> I have a problem with the cross-correlation computation.
> The output for the program is:
>
>
> 2015-02-18 17:56:31,075 [INFO] Will compute ZZ
> 2015-02-18 17:56:31,306 [INFO] New CC Job: 2015-01-01 (4 pairs with 5 
> stations)
> 2015-02-18 17:56:31,408 [DEBUG] ES.CJUL.Z Reading 1 Files
> 2015-02-18 17:56:33,837 [DEBUG] ES.CJUL.Z Merging Stream
> 2015-02-18 17:56:33,839 [DEBUG] ES.CJUL.Z Slicing Stream to 
> 2015-01-01T00:00:00.000000Z:2015-01-01T23:59:59.990000Z
> 2015-02-18 17:56:33,937 [DEBUG] ES.CJUL.Z Lowpass at 8.00 Hz
> 2015-02-18 17:56:38,981 [DEBUG] ES.CJUL.Z Highpass at 0.01 Hz
> 2015-02-18 17:56:39,705 [DEBUG] ES.CJUL.Z Decimate by a factor of 5
> 2015-02-18 17:56:39,758 [DEBUG] ES.CORC.Z Reading 1 Files
> 2015-02-18 17:56:42,276 [DEBUG] ES.CORC.Z Merging Stream
> 2015-02-18 17:56:42,277 [DEBUG] ES.CORC.Z Slicing Stream to 
> 2015-01-01T00:00:00.000000Z:2015-01-01T23:59:59.990000Z
> 2015-02-18 17:56:42,277 [DEBUG] ES.CORC.Z Lowpass at 8.00 Hz
> 2015-02-18 17:56:42,924 [DEBUG] ES.CORC.Z Highpass at 0.01 Hz
> 2015-02-18 17:56:43,625 [DEBUG] ES.CORC.Z Decimate by a factor of 5
> 2015-02-18 17:56:43,670 [DEBUG] ES.CRST.Z Reading 1 Files
> 2015-02-18 17:56:46,168 [DEBUG] ES.CRST.Z Merging Stream
> 2015-02-18 17:56:46,168 [DEBUG] ES.CRST.Z Slicing Stream to 
> 2015-01-01T00:00:00.000000Z:2015-01-01T23:59:59.990000Z
> 2015-02-18 17:56:46,169 [DEBUG] ES.CRST.Z Lowpass at 8.00 Hz
> 2015-02-18 17:56:46,817 [DEBUG] ES.CRST.Z Highpass at 0.01 Hz
> 2015-02-18 17:56:47,478 [DEBUG] ES.CRST.Z Decimate by a factor of 5
> 2015-02-18 17:56:47,525 [DEBUG] ES.CTAB.Z Reading 1 Files
> 2015-02-18 17:56:50,316 [DEBUG] ES.CTAB.Z Merging Stream
> 2015-02-18 17:56:50,443 [DEBUG] ES.CTAB.Z Slicing Stream to 
> 2015-01-01T00:00:00.000000Z:2015-01-01T23:59:59.990000Z
> 2015-02-18 17:56:50,445 [DEBUG] ES.CTAB.Z Lowpass at 8.00 Hz
> 2015-02-18 17:56:51,115 [DEBUG] ES.CTAB.Z Highpass at 0.01 Hz
> 2015-02-18 17:56:51,775 [DEBUG] ES.CTAB.Z Decimate by a factor of 5
> 2015-02-18 17:56:51,815 [DEBUG] ES.CTIG.Z Reading 1 Files
> 2015-02-18 17:56:54,611 [DEBUG] ES.CTIG.Z Merging Stream
> 2015-02-18 17:56:55,737 [DEBUG] ES.CTIG.Z Slicing Stream to 
> 2015-01-01T00:00:00.000000Z:2015-01-01T23:59:59.990000Z
> 2015-02-18 17:56:55,738 [DEBUG] ES.CTIG.Z Lowpass at 8.00 Hz
> 2015-02-18 17:56:56,403 [DEBUG] ES.CTIG.Z Highpass at 0.01 Hz
> 2015-02-18 17:56:57,074 [DEBUG] ES.CTIG.Z Decimate by a factor of 5
> 2015-02-18 17:56:57,113 [DEBUG] Processing pair: ES.CJUL vs ES.CTIG 
> Traceback (most recent call last):
>    File "s03compute_cc.py", line 402, in <module>
>      trame2h[i]*cp, Nfft, dt, low, high, plot=False)
>    File "/home/juan/MSNoise-1.2.5/whiten.py", line 58, in whiten
>      low = J[0] - Napod
> IndexError: index 0 is out of bounds for axis 0 with size 0
>
>
> Any solution?
>
> Thank you very much for your time.
>
> Juan Rueda
> Instituto Geográfico Nacional
> Spain
> _______________________________________________
> MSNoise mailing list
> MSNoise at mailman-as.oma.be
> http://mailman-as.oma.be/mailman/listinfo/msnoise

_______________________________________________
MSNoise mailing list
MSNoise at mailman-as.oma.be
http://mailman-as.oma.be/mailman/listinfo/msnoise


More information about the MSNoise mailing list