PhreeqcUsers Discussion Forum
Processes => Mixing => Topic started by: Seven on April 24, 2017, 05:11:27 AM

Dear users,
I try to mix the ratio of stable isotopes, δ18O and δD by two solutions. It's perhaps a simple question.
But in the output, there are not the data of δ18O and δD directly.
However, take the example 18, there are the data of δ34S in the output.
How could I show the result of δ18O and δD in the mixing model?
Thanks for you help.
Seven
#Input as follow:
Solution 1 rain
isotope 18O 6.27 0.1
isotope 2H 40.42 1
save solution 1
end
Solution 2 river
isotope 18O 10.57 0.1
isotope 2H 72.66 1
save solution 2
end
Mix 1
1 40.49
2 59.57
save solution 3
end

The isotope option is used only in INVERSE_MODELING and does not result in a consideration of the isotopes in the initial solution calculation or in reaction calculations. Your script does not use INVERSE_MODELING, so the isotope definitions have no effect.
If you want to consider the isotopes as separate chemical components, you must use iso.dat database, and define the isotopes as [18O] and D.
Solution 1 rain
[18O] 6.27 0.1
D 40.42 1
save solution 1
end
Solution 2 river
[18O] 10.57 0.1
D 72.66 1
save solution 2
end
Mix 1
1 40.49
2 59.57
save solution 3
end

Thanks a lot! It's really a great help to me.