Hi GWB Forum,
I am trying to mix two seawater samples with the same chemistry, but different temperatures and calculate a flash diagram for aragonite saturation state using React following the "Fluid Mixing and Scaling" instructions. I have also calculated the endmember fluid Log(Q/K) values for aragonite using GSS.
I am getting a very different answer for the end member compositions depending on which approach I take. I checked that I am using the same thermo database (phreeqpitz) for both, and am having trouble understanding why they are so different.
I've attached my React file, and my gss file for your reference.
Thanks for any insights,
Amanda
JFEndmembers.gss
JanFeb-Mixing-AssumeSWinboth.rea