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

# http://pyrocko.org - GPLv3 

# 

# The Pyrocko Developers, 21st Century 

# ---|P------/S----------~Lg---------- 

from __future__ import absolute_import, division 

 

from pyrocko import util, model 

from .base_catalog import EarthquakeCatalog 

from .kinherd import ws_request 

 

import logging 

 

logger = logging.getLogger('pyrocko.client.saxony') 

 

km = 1000. 

 

 

class Saxony(EarthquakeCatalog): 

 

def __init__(self): 

self._events = None 

 

def retrieve(self): 

url = 'http://home.uni-leipzig.de/collm/auswertung_temp.html' 

 

f = ws_request(url) 

text = f.read() 

sec = 0 

events = {} 

for line in text.splitlines(): 

line = line.strip() 

if line == '<PRE>': 

sec += 1 

continue 

 

if sec == 1 and not line: 

sec += 1 

continue 

 

if sec == 1: 

t = line.split(' ', 1) 

name = t[0] 

sdate = t[1][0:11] 

stime = t[1][12:22] 

sloc = t[1][23:36] 

sdepth = t[1][37:42] 

smag = t[1][51:55] 

region = t[1][60:] 

 

sday, smon, syear = sdate.split('-') 

smon = { 

'Jan': '01', 'Feb': '02', 'Mar': '03', 'Apr': '04', 

'Mai': '05', 'Jun': '06', 'Jul': '07', 'Aug': '08', 

'Sep': '09', 'Okt': '10', 'Nov': '11', 'Dez': '12'}[smon] 

 

time = util.str_to_time( 

'%s-%s-%s %s' % (syear, smon, sday, stime)) 

 

slat, slon = sloc.split(';') 

 

ev = model.event.Event( 

time=time, 

lat=float(slat), 

lon=float(slon), 

depth=float(sdepth) * 1000., 

magnitude=float(smag), 

magnitude_type='Ml', 

name=name, 

region=region, 

catalog='Saxony') 

 

events[name] = ev 

 

self._events = events 

 

def iter_event_names(self, time_range=None, **kwargs): 

if self._events is None: 

self.retrieve() 

 

for name in sorted(self._events.keys()): 

if time_range is not None: 

ev = self._events[name] 

if time_range[0] <= ev.time and ev.time <= time_range[1]: 

yield name 

else: 

yield name 

 

def get_event(self, name): 

if self._events is None: 

self.retrieve() 

 

return self._events[name]