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

S = htm.Solubility(
   data_T=[673, 773, 873, 973, 1073] * htm.ureg.K,
   data_y=[3e+21, 9e+20, 5e+20, 3e+20, 1e+20]
   * htm.ureg.particle
   * htm.ureg.m**-3
   * htm.ureg.Pa**-0.5,
)

plot(S)

plt.yscale("log")
plt.show()