[MSNoise] Compute MWCS

[weiyu]

Hi MSNoise group,
I recently ran into a problem when compute MWCS. I have 9 filters in my project, when I ran “msnoise compute-mwcs”,it worked successfully, but when calculated filter 8 and filter 9 something happened, the parameters were “0.04 0.05 0.0667 0.1” and “0.025 0.033 0.05 0.1” respectively. The code reported some Warnings: 

2019-06-10 15:39:13 [INFO] We will recompute all MWCS based on the new REF for SX.YDH:SX.YSG
2019-06-10 15:39:13 [DEBUG] Ready to compute
2019-06-10 15:39:14 [INFO] There are MWCS jobs for some days to recompute for CQ.CHK:CQ.CHS
2019-06-10 15:39:14 [DEBUG] Processing MWCS for: CQ_CHK_CQ_CHS.ZZ.08 - 2016-01-01 - 05 days
/data1/dell/Software/anaconda3/envs/msnoise/lib/python3.6/site-packages/scipy/optimize/minpack.py:779: OptimizeWarning: Covariance of the parameters could not be estimated
  category=OptimizeWarning)
/data1/dell/Software/anaconda3/envs/msnoise/lib/python3.6/site-packages/msnoise/move2obspy.py:349: RuntimeWarning: invalid value encountered in double_scalars
  e = np.sum((phi - m * v) ** 2) / (np.size(v) - 1)
/data1/dell/Software/anaconda3/envs/msnoise/lib/python3.6/site-packages/msnoise/move2obspy.py:349: RuntimeWarning: divide by zero encountered in double_scalars
  e = np.sum((phi - m * v) ** 2) / (np.size(v) - 1)
2019-06-10 15:39:16 [DEBUG] Processing MWCS for: CQ_CHK_CQ_CHS.ZZ.08 - 2016-01-01 - 15 days
2019-06-10 15:39:16 [DEBUG] Processing MWCS for: CQ_CHK_CQ_CHS.ZZ.08 - 2016-01-01 - 30 days
2019-06-10 15:39:17 [DEBUG] Processing MWCS for: CQ_CHK_CQ_CHS.ZZ.08 - 2016-01-01 - 45 days

With these warnings, the code still worked, and then new warnings occured:

2019-06-10 15:42:33 [DEBUG] Processing MWCS for: CQ_CHK_CQ_CHS.ZZ.08 - 2016-04-25 - 15 days
2019-06-10 15:42:33 [DEBUG] Processing MWCS for: CQ_CHK_CQ_CHS.ZZ.08 - 2016-04-25 - 30 days
/data1/dell/Software/anaconda3/envs/msnoise/lib/python3.6/site-packages/msnoise/move2obspy.py:327: RuntimeWarning: divide by zero encountered in true_divide
  w = 1.0 / (1.0 / (coh[index_range] ** 2) - 1.0)
/data1/dell/Software/anaconda3/envs/msnoise/lib/python3.6/site-packages/msnoise/move2obspy.py:327: RuntimeWarning: invalid value encountered in true_divide
  w = 1.0 / (1.0 / (coh[index_range] ** 2) - 1.0)
2019-06-10 15:42:34 [DEBUG] Processing MWCS for: CQ_CHK_CQ_CHS.ZZ.08 - 2016-04-25 - 45 days
2019-06-10 15:42:34 [DEBUG] Processing MWCS for: CQ_CHK_CQ_CHS.ZZ.08 - 2016-04-25 - 60 days

And the code was still running, the output was as follows: 

-3.300000000000000000e+02 -7.905952284653005258e+00 nan 9.937050757172795068e-01
-3.250000000000000000e+02 -8.138657988098142226e+00 nan 9.976592601521057091e-01
-3.200000000000000000e+02 -7.632117540664275168e+00 nan 9.894723420074960085e-01
-3.150000000000000000e+02 -4.285274200620775353e+00 nan 8.972202676215104145e-01
-3.100000000000000000e+02 -1.489405237146531613e+00 nan 9.666564114380370443e-01
-3.050000000000000000e+02 -1.821892139101924268e+00 nan 9.604215787607996990e-01
-3.000000000000000000e+02 8.124339913942842273e+00 nan 9.809556426874297408e-01
-2.950000000000000000e+02 1.436912500259247771e+00 nan 9.710551233524897796e-01
-2.900000000000000000e+02 6.558446922684406655e-01 nan 9.932296151167433473e-01
-2.850000000000000000e+02 5.593336746118490188e-01 inf 9.946785904215598162e-01
-2.800000000000000000e+02 1.623208176830735067e-02 inf 9.997050009388863678e-01
-2.750000000000000000e+02 -8.591600654862702413e+00 inf 9.924982060141236406e-01
-2.700000000000000000e+02 5.866839658532229329e+00 nan 9.421224113061066330e-01
-2.650000000000000000e+02 4.197498248311506863e+00 nan 9.094520333511715293e-01
-2.600000000000000000e+02 6.883655670067920340e+00 nan 9.634636129273534966e-01
-2.550000000000000000e+02 6.902610666221385749e+00 nan 9.689572985482811651e-01
-2.500000000000000000e+02 -2.224539350609665789e+00 nan 9.149593016370045628e-01
-2.450000000000000000e+02 -2.569921624096909163e+00 nan 9.193301735266846686e-01
-2.400000000000000000e+02 -1.787990836057278932e+00 inf 9.549244033726577152e-01
-2.350000000000000000e+02 -7.489977845697655123e-01 nan 9.923064807292454326e-01

I don’t know why this problem only occured at filter 8 and 9.
looking forward to know what’s going on

very very thanks!
weiyu