Hi All;

I have an ascii file that looks like this:

strike dip rake depth VR Mag

130 74 74 19 8.437e+01 4.17

129 74 72 20 8.459e+01 4.17

130 75 73 21 8.458e+01 4.19

129 74 72 22 8.490e+01 4.19

129 74 70 23 8.505e+01 4.19

…

After plot the time series x,y, where x=depth and y =VR, I would like to plot the fault plane solution for each

x,y pair in the same plot, so that, I can see the change in fault geometry as a function of depth.

I was trying something like this:

import numpy as np

import matplotlib.pyplot as plt

from obspy.imaging.beachball import Beach

file = open(‘meca.dat’,‘r’)

lines = file.readlines()

file.close()

x1 =[]

x2=[]

x3=[]

x4=[]

x5=[]

x6=[]

strike=np.array(x1)

dip=np.array(x2)

rake=np.array(x3)

dep=np.array(x4)

vrd=np.array(x5)

mag=np.array(x6)

fig, ax = plt.subplots(figsize=(10,6))

plt.scatter(dep, vrd, s=mag*15, c=mag, alpha=0.75)

plt.colorbar()

ax.set_ylabel(‘Variance reduction (%)’)

ax.set_xlabel(‘Depth (km)’)

bball = Beach([strike, dip, rake], xy=([dep,vrd]), width=(0.80,0.45), linewidth=1, alpha=0.85)

bball.set_zorder(1)

ax.add_collection(bball)

plt.show()

But it doesn’t work: ValueError: The truth value of an array with more than one element is ambiguous. Use a.any() or a.all()

Is there an smart way to do this? Any help will be appreciated.

Cheers;

-Esteban