import numpy as np
import matplotlib.pyplot as plt
fig, ax = plt.subplots()
y = np.linspace(700.0, 900.0, 401)
offsets = [94.0, 95.0, 96.0, 97.0]
times = [800.0, 790.0, 780.0, 770.0]
for offset, time in zip(offsets, times):
    x = offset + 10 * (np.sin(y / 2.0) * 1 / (10.0 * np.sqrt(2 * np.pi)) *
        np.exp(-(y - time) ** 2 / (2 * 10.0 ** 2)))
ax.plot(x, y, 'k-')
ax.fill_betweenx(y, offset, x, where=x > offset, color='k')
ax.set_xlim(93, 98)
plt.show()
plotting seismic wiggle traces using matplotlib