Re:Re: problems with uncertainty term in inverse model

["David L Parkhurst" <dlpark@xxxxxxxx>]

What happen if I don´t specify "dissolve or
precip"? Will the program show me all the
possible models although some of them may
not be right?

Yes.

Here is a file that makes a model. You can start from there trying to
determine what you like or don't like about the model. You do have chloride
decreasing, which is not possible except by dilution with another water, or
you are not on the same flow path. You probably need phases like kaolinite
and chalcedony if you are reacting primary silicates. Alternatively, you
may have other clays that are forming. This also opens up the possibility
of cation exchange, which I have not included.

The -mineral_water false statement eliminates some models that generate
water by moles of reaction of minerals. I prefer not to use this option,
but it simplifies the result for you. If you take it out, you will find
models that have mixing fractions greater than 1.0 for solution 1 and mole
transfers of tens of moles for some of the minerals.

Hope this gets you started.

David

INVERSE_MODELING
        -solution 1 2
             -uncertainty 0.05
             -range
             -phases

                         Albite     dis
                         Anorthite     dis
                         Adularia        dis
                         Calcite
                         Dolomite
                         Fe(OH)3(a)
                         Halite
#                        H2O(g)
                         N2(g)
                         Rhodochrosite
                         Siderite
                         GYPSUM

                        CO2(g)
#                       O2(g)
                        Chalcedony
                        Kaolinite
      -mineral_water false
END

David Parkhurst (dlpark@xxxxxxxx)
U.S. Geological Survey
Box 25046, MS 413
Denver Federal Center
Denver, CO 80225