def test_estimateMagnitude(self): """ Tests against PITSA. Note that PITSA displays microvolt, that is the amplitude values must be computed back into counts (for this stations .596microvolt/count was used). Pitsa internally calculates with the sensitivity 2800 of the WA. Using this we get for the following for event 2009-07-19 23:03:: RTSH PITSA 2.263 ObsPy 2.294 RTBE PITSA 1.325 ObsPy 1.363 RMOA PITSA 1.629 ObsPy 1.675 """ paz = {'poles': [-4.444 + 4.444j, -4.444 - 4.444j, -1.083 + 0j], \ 'zeros': [0 + 0j, 0 + 0j, 0 + 0j], \ 'gain': 1.0, \ 'sensitivity': 671140000.0} mag_RTSH = estimateMagnitude(paz, 3.34e6, 0.065, 0.255) self.assertAlmostEqual(mag_RTSH, 2.1328727151723488) mag_RTBE = estimateMagnitude(paz, 3.61e4, 0.08, 2.197) self.assertAlmostEqual(mag_RTBE, 1.1962687721890191) mag_RNON = estimateMagnitude(paz, 6.78e4, 0.125, 1.538) self.assertAlmostEqual(mag_RNON, 1.4995311686507182)
def test_estimateMagnitude(self): """ Tests against PITSA. Note that PITSA displays microvolt, that is the amplitude values must be computed back into counts (for this stations .596microvolt/count was used). Pitsa internally calculates with the sensitivity 2800 of the WA. Using this we get for the following for event 2009-07-19 23:03:: RTSH PITSA 2.263 ObsPy 2.294 RTBE PITSA 1.325 ObsPy 1.363 RMOA PITSA 1.629 ObsPy 1.675 """ paz = {'poles': [-4.444 + 4.444j, -4.444 - 4.444j, -1.083 + 0j], \ 'zeros': [0 + 0j, 0 + 0j, 0 + 0j], \ 'gain': 1.0, \ 'sensitivity': 671140000.0} mag_RTSH = estimateMagnitude(paz, 3.34e6, 0.065, 0.255) self.assertAlmostEqual(mag_RTSH, 2.1328727151723488) mag_RTBE = estimateMagnitude(paz, 3.61e4, 0.08, 2.197) self.assertAlmostEqual(mag_RTBE, 1.1962687721890191) mag_RNON = estimateMagnitude(paz, 6.78e4, 0.125, 1.538) self.assertAlmostEqual(mag_RNON, 1.4995311686507182)