When I open the thermo.dat file, I saw there are some parameters for pressure after those for the principle temperatures. How are the pressure effect considered in GWB? I am using GWB to calculate the solution under a constant CO2 pressure of 102 atm, and at temperature of 95 deg C. Under this condition, CO2 is at supercritical phase. I am wondering how could I calculate the fugacity coefficient under this pressure?
Also, the equilibrium constant for other aqueous species will change under high pressure also, how would GWB consider this effect? thank you.
Duan, Z., and Sun, R. (2003) An improved model calculating CO2 solubility in pure water and aqueous NaCl solutions from 273 to 533 K and from 0 to 2000 bar, Chemical Geology, v193, p. 257-271.