-
Notifications
You must be signed in to change notification settings - Fork 1
/
neon_ternary.py
56 lines (43 loc) · 1.57 KB
/
neon_ternary.py
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
import matplotlib
matplotlib.use('Qt5Agg')
matplotlib.rcParams['backend.qt5'] = 'PySide2'
import ternary
import numpy as np
def plot_tern(data_point):
# Scatter Plot
scale = 100
figure, tax = ternary.figure(scale=scale)
figure.set_size_inches(10, 10)
tax.set_title("Neon Isotopes", fontsize=20)
tax.boundary(linewidth=2.0)
tax.gridlines(multiple=10, color="black")
tax.left_axis_label("Neon-22", offset=0.1)
tax.right_axis_label("Neon-21", offset=0.1)
tax.bottom_axis_label("Neon-20", offset=-0.05)
# Air
points = []
points.append([90.48, 00.27, 09.25])
tax.scatter(points, marker='s', color='red', label="Air")
# Cosmogenic (Quartz; Schaefer et al. 1999)
points = []
points.append([35.59, 28.47, 35.94])
tax.scatter(points, marker='s', color='green', label="Cosmogenic (Quartz)")
# Cosmogenic (Pyroxene; Schaefer et al. 1999)
points = []
points.append([29.67, 32.64, 37.69])
tax.scatter(points, marker='s', color='blue', label="Cosmogenic (Pyroxene)")
# Mantle
points = []
points.append([92.80, 00.47, 06.73])
tax.scatter(points, marker='s', color='purple', label="Mantle")
# Nucleogenic (Bulk Crustal)
points = []
points.append([08.39, 89.73, 01.88])
tax.scatter(points, marker='s', color='grey', label="Nucleogenic (Bulk Crust)")
points = []
points.append(data_point)
tax.scatter(points, marker='s', color='black', label="Data")
tax.legend()
tax.ticks(axis='lbr', multiple=10, linewidth=1, offset=0.012)
tax.clear_matplotlib_ticks()
tax.show()