Hi there,

First of all, thank you so much for this fantastic tool! It's really helping us in our research.

Second, I wanted to see if you had a solution to this issue I'm seeing:

I would like to find the conductivity of high molality solution (1 - 3 molal) for potassium carbonate solution. When I use PHREEQC Interactive (the USGS software) directly I get results that match chemistry handbooks. So, for example, entering this input and running the Pitzer model:

SOLUTION 0 units mmol/kgw pH 10 charge temp 20 K 1000 C 500 END 

Gives 73.3 mS/cm and pH = 11.79

BUT, when I use your phreeqpy function

pp = PhreeqPython(database='pitzer.dat') solution = pp.add_solution({'units':'mol/kgw', #set the units (moles per kg of water) 'pH': '10 charge', 'temp': 20, 'K':2, 'C': 1 }) pH = solution.pH cond = solution.sc print('pH:',pH) print('Cond:',cond) 

I get 50.5 mS/cm and the same pH as above (11.79)

There is a substantial difference in conductivity and I'd really like to be able to reconcile that. The phreeqpy function is outputting a lower conductivity than physical measurements and than PHREEQC. Do you have any suggestions for why this may be happening? Note that at molality below ~0.5, I don't see the discrepancy anymore.

Thanks so much, and again, really appreciate this tool!