data = thermo.dat verify swap CO2(g) for HCO3- f CO2(g) = 10^-3 balance on H+ kinetic CO2(g) f_ext = 10^-3.5, rate_con = 2e-7, contact = 8 time end 24 hr suffix _Gastransfer