def plottern(plt, alpha, vp):
gran = np.array([6.208, poisson(0.25, True)])
mafc = np.array([6.942, poisson(0.283, True)])
gray = np.array([6.302, poisson(0.257, True)])
vp = np.array([vp, alpha])
a = gran
b = mafc
c = gray
d = vp
endmembers = ("Granite Gneiss", "Mafic Granulite", "Gray Gneiss")
datalabels = ["Vp", r"$\sigma$"]
colors = ["b", "g", "c", "m", "r"]
lines = terntransform(a, b, c, d)
step = 0.2
for line, color, label in zip(lines, colors, datalabels):
tern.plot(line, color=color, linewidth=2.0, label=label)
ip = baryIntersect(lines)
tern.plotIntersect(ip)
tern.gridlines(step, "0.7")
tern.draw_boundary()
tern.addlabels(endmembers)
plt.legend()
# pyplot.box(on='off')
h = plt.gca()
# plt.axis('off')
xmin = -0.1
xmax = 1.1
ymin = -0.1
ymax = 1
v = [xmin, xmax, ymin, ymax]
plt.axis(v)
h.axes.get_xaxis().set_visible(False)
h.axes.get_yaxis().set_visible(False)
def plottern(plt, alpha, vp):
gran = np.array([6.208, 0.25])
mafc = np.array([6.942, 0.283])
gray = np.array([6.302, 0.257])
vp = np.array([vp, alpha])
a = gran
b = mafc
c = gray
d = vp
endmembers = ('Granite Gneiss', 'Mafic Granulite', 'Gray Gneiss')
datalabels = ['Vp',r'$\sigma$']
colors = ['b','g','c','m','r']
lines = terntransform(a, b, c, d)
step = 0.2
for line,color,label in zip(lines,colors,datalabels):
tern.plot(line, color = color, linewidth=2.0, label = label)
ip = baryIntersect(lines)
tern.plotIntersect(ip)
tern.gridlines(step, '0.7')
tern.draw_boundary()
tern.addlabels(endmembers)
plt.legend()
#pyplot.box(on='off')
h = plt.gca()
#plt.axis('off')
xmin = -0.1
xmax = 1.1
ymin = -0.1
ymax = 1
v = [xmin, xmax, ymin, ymax]
plt.axis(v)
h.axes.get_xaxis().set_visible(False)
h.axes.get_yaxis().set_visible(False)
def plottern(plt, alpha, vp):
gran = np.array([6.208, poisson(0.25, True)])
mafc = np.array([6.942, poisson(0.283, True)])
gray = np.array([6.302, poisson(0.257, True)])
vp = np.array([vp, alpha])
a = gran
b = mafc
c = gray
d = vp
endmembers = ('Granite Gneiss', 'Mafic Granulite', 'Gray Gneiss')
datalabels = ['Vp',r'$\sigma$']
colors = ['b','g','c','m','r']
lines = terntransform(a, b, c, d)
step = 0.2
for line,color,label in zip(lines,colors,datalabels):
tern.plot(line, color = color, linewidth=2.0, label = label)
ip = baryIntersect(lines)
tern.plotIntersect(ip)
tern.gridlines(step, '0.7')
tern.draw_boundary()
tern.addlabels(endmembers)
plt.legend()
#pyplot.box(on='off')
h = plt.gca()
#plt.axis('off')
xmin = -0.1
xmax = 1.1
ymin = -0.1
ymax = 1
v = [xmin, xmax, ymin, ymax]
plt.axis(v)
h.axes.get_xaxis().set_visible(False)
h.axes.get_yaxis().set_visible(False)
endmembers = ('Granite Gneiss', 'Mafic Granulite', 'Gray Gneiss')
datalabels = ['Vp',r'$\sigma$']
colors = ['b','g','c','m','r']
lines = terntransform(a, b, c, d)
step = 0.2
plt.figure()
for line,color,label in zip(lines,colors,datalabels):
tern.plot(line, color = color, linewidth=2.0, label = label)
ip = baryIntersect(lines)
tern.plotIntersect(ip)
tern.gridlines(step, '0.7')
tern.draw_boundary()
tern.addlabels(endmembers)
plt.legend()
#pyplot.box(on='off')
h = plt.gca()
#plt.axis('off')
xmin = -0.1
xmax = 1.1
ymin = -0.1
ymax = 1
v = [xmin, xmax, ymin, ymax]
data = lithdata[3:6]
data.append(region[1])
datalabels = ['Vp', r'$\sigma$']
colors = ['b', 'g', 'c', 'm', 'r', 'y']
lines = terntransform(*data)
for line, color, label in zip(lines, colors, datalabels):
tern.plot(line,
color=color,
linewidth=2.0,
label=label if (ind == 0) else None)
ip = baryIntersect(lines)
tern.plotIntersect(ip, marker='o' + colors[ind], label=region[0])
step = 0.2
tern.gridlines(step, '0.7')
tern.draw_boundary()
tern.addlabels(endmembers)
plt.legend()
#pyplot.box(on='off')
h = plt.gca()
#plt.axis('off')
xmin = -0.1
xmax = 1.1
ymin = -0.1
ymax = 1
