import h_transport_materials as htm
from h_transport_materials.plotting import plot
import matplotlib.pyplot as plt

# filter only tungsten and H
diffusivities = htm.diffusivities.filter(material="tungsten").filter(isotope="h")

plot(diffusivities)

plt.title("Tungsten diffusivity")
plt.yscale("log")
plt.legend()
plt.show()