data = thermo.dat verify balance off swap CO2(g) for HCO3- CO2(g) = -3.5 log fugacity pH = 2 slide pH to 10 fix f CO2(g) suffix _Gassoly