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=htm.TUNGSTEN)

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

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