PhreeqcRM SetConcentrations

[strofp]

Hello,
If a component value is low (e.g. 10^-23) then GetConcentrations (immediately after SetConcentrations) returns zero. Is there any switch how to fix it? What is the MINIMUM VALUE LIMIT?

Thanks for the help

[dlparkhurst]

Looks like a value of 1e-18 is coded in the Update method of solution.cpp.

PHREEQC is somewhat inconsistent on censoring values. The value in the Update method should probably be set to MIN_TOTAL, which is used in most of the code (although I won't guarantee that it is used consistently either). MIN_TOTAL is set to 1e-25 in phreeqc.cpp.

[strofp]

Thanks David,
it was exactly as you wrote and it was enough to change solution.cpp at line 1114.
Best Regards

Pavel