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")

plot(diffusivities, alpha=0.5)
plot(diffusivities.mean(), color="black", linewidth=3)

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