

import numpy as np
import matplotlib.pyplot as plt

# ====== Parameters ======
temps = np.array([0, 20, 40, 60, 80, 100])

solubility_data = {
    "Ammonia": [1176, 702, np.nan, np.nan, np.nan, np.nan],
    "Hydrogen chloride": [507, 442, 386, 339, np.nan, np.nan],
    "Chlorine": [4.61, 2.30, 1.44, 1.02, 0.68, 0.00],
    "Sulfur dioxide": [80, 39, 19, np.nan, np.nan, np.nan],
    "Carbon dioxide": [1.71, 0.88, 0.53, 0.36, np.nan, np.nan],
    "Oxygen": [0.049, 0.031, 0.023, 0.019, 0.018, 0.017],
    "Air": [0.029, 0.019, 0.014, 0.012, 0.011, 0.011],
    "Nitrogen": [0.024, 0.016, 0.012, 0.010, 0.0096, 0.0095],
    "Hydrogen": [0.022, 0.018, 0.016, 0.016, 0.016, 0.016]
}

# ====== Plot ======
plt.figure(figsize=(10,6))
for gas, values in solubility_data.items():
    plt.plot(temps, values, marker='o', label=gas)
    for x, y in zip(temps, values):
        if not np.isnan(y):
            # プロットに座標表示
            plt.text(x, y, f"({x}, {y})", fontsize=8, ha='center', va='bottom')
            # コンソールに出力
            print(f"{gas}: Temperature={x}°C, Solubility={y}")

plt.title("Gas Solubility vs Temperature", fontsize=14)
plt.xlabel("Temperature [°C]", fontsize=12)
plt.ylabel("Solubility (cm³ gas / g water)", fontsize=12)
plt.grid(True)
plt.legend()
plt.show()