@ -7,7 +7,7 @@ steps:
image: pyrocko-nest
pull: never
commands:
- flake8 src test setup.py
- flake8 src test examples setup.py
---
@ -1,5 +1,6 @@
import numpy as num
from pyrocko.dataset.gshhg import GSHHG
from matplotlib import pyplot as plt
gshhg = GSHHG.intermediate()
# gshhg = GSHHG.full()
@ -17,6 +18,5 @@ coordinates = num.array([lat_grid.ravel(), lon_grid.ravel()]).T
land_mask = gshhg.get_land_mask(coordinates).reshape(*lat_grid.shape)
plt.pcolormesh(lons, lats, land_mask, cmap='Greys')
plt.show()
@ -16,7 +16,7 @@ for network in station_xml.network_list:
frequency=1.0,
input_units=fdsn.Units('M'),
output_units=fdsn.Units('COUNTS')))
station_xml.validate()
# print(station_xml.dump_xml())
station_xml.dump_xml(filename='stations_flat_displacement.xml')
@ -6,12 +6,12 @@ exp = mtm.magnitude_to_moment(magnitude) # convert the mag to moment in [Nm]
# init pyrocko moment tensor
m = mtm.MomentTensor(
mnn = 2.34*exp,
mee = -2.64*exp,
mdd = 0.295*exp,
mne = 1.49*exp,
mnd = 0.182*exp,
med = -0.975*exp)
mnn=2.34*exp,
mee=-2.64*exp,
mdd=0.295*exp,
mne=1.49*exp,
mnd=0.182*exp,
med=-0.975*exp)
print(m) # print moment tensor
@ -8,4 +8,4 @@ get_example_data('responses.xml')
sx = stationxml.load_xml(filename='responses.xml')
# Extract Station objects from FDSNStationXML object
pyrocko_stations = sx.get_pyrocko_stations()
@ -20,8 +20,8 @@ easts = easts.reshape((lats.size, lons.size))
plt = plot.mpl_init(fontsize=10.)
fig = plt.figure(figsize=plot.mpl_papersize('a5', 'landscape'))
axes = fig.add_subplot(1, 1, 1, aspect=1.0)
cbar = axes.pcolormesh(easts, norths, tile.data,
cmap='gray', shading='gouraud')
fig.colorbar(cbar, label='Altitude [m]')
axes.set_title(dem_name)
axes.set_xlim(easts.min(), easts.max())