def setup(self): self.catalogs = { 'Geofon': catalog.Geofon(), 'USGS/NEIC PDE': catalog.USGS('pde'), 'USGS/NEIC US': catalog.USGS('us'), 'Global-CMT': catalog.GlobalCMT(), 'Kinherd': catalog.Kinherd(), } catkeys = sorted(self.catalogs.keys()) self.set_name('Catalog Search') self.add_parameter(Choice('Catalog', 'catalog', catkeys[0], catkeys)) self.add_parameter(Param('Min Magnitude', 'magmin', 0, 0, 10)) self.set_live_update(False)
def testUSGS(self): def is_the_haiti_event(ev): assert near(ev.magnitude, 7.0, 0.1) assert near(ev.lat, 18.443, 0.01) assert near(ev.lon, -72.571, 0.01) assert near(ev.depth, 13000., 1.) cat = catalog.USGS() tmin = util.ctimegm('2010-01-12 21:50:00') tmax = util.ctimegm('2010-01-13 03:17:00') names = cat.get_event_names(time_range=(tmin, tmax), magmin=5.) assert len(names) == 13 for name in names: ev = cat.get_event(name) if ev.magnitude >= 7.: is_the_haiti_event(ev) ident = ev.name assert ident is not None cat.flush() ev = cat.get_event(ident) is_the_haiti_event(ev)
import numpy as num from pyrocko import trace, util, io, cake, catalog, automap, pile, model from pyrocko import orthodrome, weeding from pyrocko.client import fdsn from pyrocko.io import resp, enhanced_sacpz as epz, stationxml from beat import utility from beat import heart km = 1000. g_sites_available = sorted(fdsn.g_site_abbr.keys()) geofon = catalog.Geofon() usgs = catalog.USGS(catalog=None) gcmt = catalog.GlobalCMT() tfade_factor = 1.0 ffade_factors = 0.5, 1.5 fdsn.g_timeout = 60. class starfill(object): def __getitem__(self, k): return '*' def nice_seconds_floor(s): nice = [
def USGS(searchparameter): cat = catalog.USGS() getevents(searchparameter, cat)