Hot-keys on this page

r m x p   toggle line displays

j k   next/prev highlighted chunk

0   (zero) top of page

1   (one) first highlighted chunk

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

32

33

34

35

36

37

38

39

40

41

42

43

44

45

46

47

48

49

50

51

52

53

54

55

56

57

58

59

60

61

62

63

64

65

66

67

68

69

70

71

72

73

74

75

76

77

78

79

80

81

82

83

84

85

86

87

88

89

90

91

92

93

94

95

96

97

98

99

100

101

102

103

104

105

# -*- coding: utf-8 -*- 

""" 

obspy.neries - NERIES Web service client for ObsPy 

================================================== 

The obspy.neries package contains a client for the Seismic Data Portal 

(http://www.seismicportal.eu) which was developed under the European 

Commission-funded NERIES project. The Portal provides a single point of access 

to diverse, distributed European earthquake data provided in a unique joint 

initiative by observatories and research institutes in and around Europe. 

 

:copyright: 

    The ObsPy Development Team (devs@obspy.org) 

:license: 

    GNU Lesser General Public License, Version 3 

    (http://www.gnu.org/copyleft/lesser.html) 

 

Basic Usage 

----------- 

(1) :meth:`~obspy.neries.client.Client.getEvents()`: The following example 

    illustrates how to request all earthquakes of a magnitude of >=9 in the 

    year 2004. 

 

    >>> from obspy.neries import Client 

    >>> client = Client(user='test@obspy.org') 

    >>> events = client.getEvents(min_datetime="2004-01-01", 

    ...                           max_datetime="2005-01-01", 

    ...                           min_magnitude=9) 

    >>> len(events) 

    1 

    >>> events #doctest: +SKIP 

    [{'author': u'CSEM', 'event_id': u'20041226_0000148', 'origin_id': 127773, 

      'longitude': 95.724, 'datetime': UTCDateTime(2004, 12, 26, 0, 58, 50), 

      'depth': -10.0, 'magnitude': 9.3, 'magnitude_type': u'mw', 

      'latitude': 3.498, 'flynn_region': u'OFF W COAST OF NORTHERN SUMATRA'}] 

 

(2) :meth:`~obspy.neries.client.Client.getLatestEvents()`: Returns only the 

    latest earthquakes. 

 

    >>> from obspy.neries import Client 

    >>> client = Client(user='test@obspy.org') 

    >>> events = client.getLatestEvents(num=5, format='list') 

    >>> len(events) 

    5 

    >>> events[0]  #doctest: +SKIP 

    [{'author': u'CSEM', 'event_id': u'20041226_0000148', 'origin_id': 127773, 

      'longitude': 95.724, 'datetime': u'2004-12-26T00:58:50Z', 'depth': -10.0, 

      'magnitude': 9.3, 'magnitude_type': u'mw', 'latitude': 3.498, 

      'flynn_region': u'OFF W COAST OF NORTHERN SUMATRA'}] 

 

(3) :meth:`~obspy.neries.client.Client.getEventDetail()`: Returns additional 

    information for each event by a given event_id. 

 

    >>> from obspy.neries import Client 

    >>> client = Client(user='test@obspy.org') 

    >>> result = client.getEventDetail("20041226_0000148", 'list') 

    >>> len(result)  # Number of calculated origins 

    11 

    >>> result[0]  # Details about first calculated origin #doctest: +SKIP 

    {'author': u'CSEM', 'event_id': u'20041226_0000148', 'origin_id': 127773, 

     'longitude': 95.724, 'datetime': UTCDateTime(2004, 12, 26, 0, 58, 50), 

     'depth': -10.0, 'magnitude': 9.3, 'magnitude_type': u'mw', 

     'latitude': 3.498, 'flynn_region': u'OFF W COAST OF NORTHERN SUMATRA'} 

 

(4) :meth:`~obspy.neries.client.Client.getWaveform()`: Wraps a NERIES Web 

    service build on top of the ArcLink protocol. Here we give a small example 

    how to fetch and display waveforms. 

 

    >>> from obspy.neries import Client 

    >>> from obspy import UTCDateTime 

    >>> client = Client(user='test@obspy.org') 

    >>> dt = UTCDateTime("2009-08-20 04:03:12") 

    >>> st = client.getWaveform("BW", "RJOB", "", "EH*", dt - 3, dt + 15) 

    >>> st.plot()  #doctest: +SKIP 

 

    .. plot:: 

 

        from obspy.neries import Client 

        from obspy import UTCDateTime 

        client = Client(user='test@obspy.org') 

        dt = UTCDateTime("2009-08-20 04:03:12") 

        st = client.getWaveform("BW", "RJOB", "", "EH*", dt - 3, dt + 15) 

        st.plot() 

 

(5) :meth:`~obspy.neries.client.Client.getTravelTimes()`: Wraps a Taup Web 

    service, an utility to compute arrival times using a few default velocity 

    models such as ``'iasp91'``, ``'ak135'`` or ``'qdt'``. 

 

    >>> from obspy.neries import Client 

    >>> client = Client(user='test@obspy.org') 

    >>> locations = [(48.0, 12.0), (48.1, 12.0)] 

    >>> result = client.getTravelTimes(latitude=20.0, longitude=20.0, 

    ...                                depth=10.0, locations=locations, 

    ...                                model='iasp91') 

    >>> len(result) 

    2 

    >>> result[0] # doctest: +SKIP 

    {'P': 356981.13561726053, 'S': 646841.5619481194} 

""" 

 

from client import Client 

 

 

if __name__ == '__main__': 

    import doctest 

    doctest.testmod(exclude_empty=True)