Re: [MSNoise] Compute MWCS

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 _______________________________________________ MSNoise mailing list MSNoise(a)mailman-as.oma.be http://mailman-as.oma.be/mailman/listinfo/msnoise