Value Error when using compute_cc_rot

[JessieSchuler]

I received the error below when running compute_cc_rot for ZZ and TT components. When looking at the data from the stations and the day which threw the error, nothing was obviously wrong with the days waveforms. This error did not occur when running just the compute_cc for the ZZ component for the exact same data set. I am a little bit stumped as to what could possibly be causing this error, help would be extremely appreciated!

2022-07-15 01:01:54.693141 msnoise [pid 2890485][INFO]: New CC Job: 2019-12-03 (15 pairs with 6 stations)2022-07-15 01:03:49.805719 msnoise [pid 2890485][INFO]: Will do 15 pairs
2022-07-15 01:03:49.805913 msnoise [pid 2890485][INFO]: Processing pair: 2B.BEEZ vs 2B.JANZ
2022-07-15 01:03:49.810476 msnoise [pid 2890485][INFO]: Azimuth=61.5
2022-07-15 01:03:50.646469 msnoise [pid 2890485][INFO]: Finished processing this pair. It took 0.84 seconds
2022-07-15 01:03:50.650214 msnoise [pid 2890485][INFO]: Processing pair: 2B.BEEZ vs NZ.HATZ
2022-07-15 01:03:50.652911 msnoise [pid 2890485][INFO]: Azimuth=112.7
2022-07-15 01:03:51.468286 msnoise [pid 2890485][INFO]: Finished processing this pair. It took 0.82 seconds
2022-07-15 01:03:51.471671 msnoise [pid 2890485][INFO]: Processing pair: 2B.BEEZ vs NZ.RATZ
2022-07-15 01:03:51.473891 msnoise [pid 2890485][INFO]: Azimuth=192.1
2022-07-15 01:03:52.313519 msnoise [pid 2890485][INFO]: Finished processing this pair. It took 0.84 seconds
2022-07-15 01:03:52.317092 msnoise [pid 2890485][INFO]: Processing pair: 2B.BEEZ vs NZ.WATZ
2022-07-15 01:03:52.319722 msnoise [pid 2890485][INFO]: Azimuth=331.4
2022-07-15 01:03:53.293301 msnoise [pid 2890485][INFO]: Finished processing this pair. It took 0.98 seconds
2022-07-15 01:03:53.297333 msnoise [pid 2890485][INFO]: Processing pair: 2B.BEEZ vs NZ.WHTZ
2022-07-15 01:03:53.300828 msnoise [pid 2890485][INFO]: Azimuth=45.6
2022-07-15 01:03:54.149179 msnoise [pid 2890485][INFO]: Finished processing this pair. It took 0.85 seconds
2022-07-15 01:03:54.152763 msnoise [pid 2890485][INFO]: Processing pair: 2B.JANZ vs NZ.HATZ
2022-07-15 01:03:54.155113 msnoise [pid 2890485][INFO]: Azimuth=158.5
2022-07-15 01:03:54.995363 msnoise [pid 2890485][INFO]: Finished processing this pair. It took 0.84 seconds
2022-07-15 01:03:54.998658 msnoise [pid 2890485][INFO]: Processing pair: 2B.JANZ vs NZ.RATZ
2022-07-15 01:03:55.000942 msnoise [pid 2890485][INFO]: Azimuth=227.6
2022-07-15 01:03:55.900344 msnoise [pid 2890485][INFO]: Finished processing this pair. It took 0.90 seconds
2022-07-15 01:03:55.903964 msnoise [pid 2890485][INFO]: Processing pair: 2B.JANZ vs NZ.WATZ
2022-07-15 01:03:55.906422 msnoise [pid 2890485][INFO]: Azimuth=269.2
2022-07-15 01:03:56.805884 msnoise [pid 2890485][INFO]: Finished processing this pair. It took 0.90 seconds
2022-07-15 01:03:56.809332 msnoise [pid 2890485][INFO]: Processing pair: 2B.JANZ vs NZ.WHTZ
2022-07-15 01:03:56.811663 msnoise [pid 2890485][INFO]: Azimuth=320.7
2022-07-15 01:03:57.651677 msnoise [pid 2890485][INFO]: Finished processing this pair. It took 0.84 seconds
2022-07-15 01:03:57.655009 msnoise [pid 2890485][INFO]: Processing pair: NZ.HATZ vs NZ.RATZ
2022-07-15 01:03:57.657426 msnoise [pid 2890485][INFO]: Azimuth=275.4
2022-07-15 01:03:58.685092 msnoise [pid 2890485][INFO]: Finished processing this pair. It took 1.03 seconds
2022-07-15 01:03:58.688469 msnoise [pid 2890485][INFO]: Processing pair: NZ.HATZ vs NZ.WATZ
2022-07-15 01:03:58.691099 msnoise [pid 2890485][INFO]: Azimuth=303.0
2022-07-15 01:03:59.732047 msnoise [pid 2890485][INFO]: Finished processing this pair. It took 1.04 seconds
2022-07-15 01:03:59.735727 msnoise [pid 2890485][INFO]: Processing pair: NZ.HATZ vs NZ.WHTZ
2022-07-15 01:03:59.738349 msnoise [pid 2890485][INFO]: Azimuth=334.8
2022-07-15 01:04:00.791191 msnoise [pid 2890485][INFO]: Finished processing this pair. It took 1.06 seconds
2022-07-15 01:04:00.795035 msnoise [pid 2890485][INFO]: Processing pair: NZ.RATZ vs NZ.WATZ
2022-07-15 01:04:00.797669 msnoise [pid 2890485][INFO]: Azimuth=349.1
2022-07-15 01:04:01.835085 msnoise [pid 2890485][INFO]: Finished processing this pair. It took 1.04 seconds
2022-07-15 01:04:01.838538 msnoise [pid 2890485][INFO]: Processing pair: NZ.RATZ vs NZ.WHTZ
2022-07-15 01:04:01.841013 msnoise [pid 2890485][INFO]: Azimuth=35.9
2022-07-15 01:04:02.967900 msnoise [pid 2890485][INFO]: Finished processing this pair. It took 1.13 seconds
2022-07-15 01:04:02.971752 msnoise [pid 2890485][INFO]: Processing pair: NZ.WATZ vs NZ.WHTZ
2022-07-15 01:04:02.974436 msnoise [pid 2890485][INFO]: Azimuth=76.6
2022-07-15 01:04:04.027373 msnoise [pid 2890485][INFO]: Finished processing this pair. It took 1.06 seconds
2022-07-15 01:04:04.051348 msnoise [pid 2890485][INFO]: Job Finished. It took 129.36 seconds (preprocess: 115.11 s & process 14.25 s)
2022-07-15 01:04:04.108529 msnoise [pid 2890485][INFO]: New CC Job: 2019-12-04 (15 pairs with 6 stations)
2022-07-15 01:05:58.293444 msnoise [pid 2890485][INFO]: Will do 15 pairs
2022-07-15 01:05:58.293655 msnoise [pid 2890485][INFO]: Processing pair: 2B.BEEZ vs 2B.JANZ
2022-07-15 01:05:58.298394 msnoise [pid 2890485][INFO]: Azimuth=61.5
2022-07-15 01:05:59.507618 msnoise [pid 2890485][INFO]: Finished processing this pair. It took 1.21 seconds
2022-07-15 01:05:59.511583 msnoise [pid 2890485][INFO]: Processing pair: 2B.BEEZ vs NZ.HATZ
2022-07-15 01:05:59.514418 msnoise [pid 2890485][INFO]: Azimuth=112.7
Traceback (most recent call last):
  File "/opt/conda/envs/msnoise/bin/msnoise", line 10, in
    sys.exit(run())
  File "/opt/conda/envs/msnoise/lib/python3.9/site-packages/msnoise/scripts/msnoise.py", line 1202, in run
    cli(obj={})
  File "/opt/conda/envs/msnoise/lib/python3.9/site-packages/click/core.py", line 1128, in call
    return self.main(*args, **kwargs)
  File "/opt/conda/envs/msnoise/lib/python3.9/site-packages/click/core.py", line 1053, in main
    rv = self.invoke(ctx)
  File "/opt/conda/envs/msnoise/lib/python3.9/site-packages/click/core.py", line 1659, in invoke
    return _process_result(sub_ctx.command.invoke(sub_ctx))
  File "/opt/conda/envs/msnoise/lib/python3.9/site-packages/click/core.py", line 1395, in invoke
    return ctx.invoke(self.callback, **ctx.params)
  File "/opt/conda/envs/msnoise/lib/python3.9/site-packages/click/core.py", line 754, in invoke
    return __callback(*args, **kwargs)
  File "/opt/conda/envs/msnoise/lib/python3.9/site-packages/click/decorators.py", line 26, in new_func
    return f(get_current_context(), *args, **kwargs)
  File "/opt/conda/envs/msnoise/lib/python3.9/site-packages/msnoise/scripts/msnoise.py", line 697, in compute_cc_rot
    main(loglevel=loglevel)
  File "/opt/conda/envs/msnoise/lib/python3.9/site-packages/msnoise/s03compute_cc.py", line 305, in main
    t2_novert = make_same_length(t2_novert)
  File "/opt/conda/envs/msnoise/lib/python3.9/site-packages/msnoise/api.py", line 1754, in make_same_length
    mask = np.logical_or(st[0].data.mask, st[1].data.mask)
ValueError: operands could not be broadcast together with shapes (2159959,) (2159960,)

[ThomasLecocq]

It looks like the rotation/or another processing step produces a trace object with 1 less (or more) sample than the others, which shouldn't occur. You could try to edit (temporarily) /opt/conda/envs/msnoise/lib/python3.9/site-packages/msnoise/api.py with something to trim the data arrays to the same length, add this before line 1754:

minlen = np.min([tr.stats.npts for tr in st])
for tr in st:
    tr.data = tr.data[:minlen]