TITLE   Transport: Bear Creek example in Phreeqc 2.0
SOLUTION 0  incoming solution
        units   ppm
        pH      3.8
        pe      15
        temp    25.0
        Ca      310.00
        Mg      1000.
        Cl      550
        C(4)    6.3
        Al      1020.
        S(6)    16500.
        Fe(3)   1950.0
        Na      360.00   charge
        K       60.
        Si      40.5
END
SOLUTION 1-25 Initial solution for column
units   ppm
        pH      3.8
        pe      15
        temp    25.0
        Ca      310.00
        Mg      1000.
        Cl      550
        C(4)    6.3
        Al      1020.
        S(6)    16500.
        Fe(3)   1950.0
        Na      360.00   charge
        K       60.
        Si      40.5
EQUILIBRIUM_PHASES 1-25
        Al(OH)3(a)  -1.36  0.0
        Fe(OH)3(a)   1.69  0.2
        Calcite      0.0   0.0
        Gypsum       0.0   0.2
        SiO2(a)      0.0   0.0
        Illite       0.0   0.0
SAVE EQUILIBRIUM_PHASES 1-25
SAVE Solution 1-25
END
SOLUTION 26-60 Initial solution for column
        units   ppm
        pH      4.5
        pe      15
        temp    25.0
        Ca      420.00
        Mg      700.
        Cl      400
        C(4)    6.3
        Al      230.0
        S(6)    8100.00
        Fe(3)   926.00
        Na      278.00 Charge
        K       42.
        Si      10.
EQUILIBRIUM_PHASES 26-60
        Al(OH)3(a)  -1.36  0.1
        Fe(OH)3(a)   1.69  0.05
        Calcite      0.0   0.0
        Gypsum       0.0   0.2
        SiO2(a)      0.0   0.0
        Illite       0.0   0.0
SAVE EQUILIBRIUM_PHASES 26-60
SAVE Solution 26-60
END
SOLUTION 61-135 Initial solution for column
        units   ppm
        pH      6.5
        pe      15
        temp    25.0
        Ca      650.00
        Mg      250.
        Cl      375
        Alkalinity    1450
        Al      1.33
        S(6)    1650.00
        Fe(3)   2.18
        Na      212
        K       18.0
        Si      9.7
EQUILIBRIUM_PHASES 61-135
        Al(OH)3(a)  -1.36  0.01
        Fe(OH)3(a)   1.69  0.05
        Calcite      0.0   0.2
        Gypsum       0.0   0.0
        SiO2(a)      0.0   0.0
        Illite       0.0   0.0
SAVE EQUILIBRIUM_PHASES 61-135
SAVE Solution 61-135
END
SOLUTION 136-200 Initial solution for column
        units   ppm
        pH      6.7
        pe      15
        temp    25.0
        Ca      550.00
        Mg      150.0
        Cl      275
        Alkalinity    878.
        Al      1.150
        S(6)    1500 charge
        Fe(3)   0.69
        Na      265.0
        K       14.
        Si      8.40
EQUILIBRIUM_PHASES 136-200
        Al(OH)3(a)  -1.36  0.01
        Fe(OH)3(a)   1.69  0.05
        Calcite      0.0   0.2
        Gypsum       0.0   0.0
        SiO2(a)      0.0   0.0
        Illite       0.0   0.0
SAVE EQUILIBRIUM_PHASES 136-200
SAVE Solution 136-200
END
PRINT
        -reset  false
SELECTED_OUTPUT
        -file bearodu3.prn
        -totals H Cl Fe Al S Na Mg C Ca K Si
        -equilibrium_phases Al(OH)3(a) Fe(OH)3(a) Calcite Gypsum SiO2(a) Illite
        -pe false
        -time true
        -step false
        -state false
TRANSPORT
        -cells                200
        -length               4
        -shifts               62
        -time_step            2522880
        -flow_direction       forward
        -boundary_conditions  flux  flux
        -correct_disp         true
        -dispersivity         10.0
        -diffusion_coef       0.0e-0 #No diffusion
        -tempr  1.0    # No heat retardation
END
SOLUTION 0  Flush with background groundwater
        units   ppm
        pH      7.4
        temp    25.0
        Ca      158.00
        Mg      21.00
        Mn(3)   0.11
        Cl      25.    charge
        C(4)    153
        Al      .0100
        S(6)    425.00
        Fe(3)   0.1
        Na      61.0
        K       7.
        Si      5.60
END
SELECTED_OUTPUT
        -file bearout.prn
        -totals Cl Fe Al S Na Mg C Ca K Si
        -equilibrium_phases Al(OH)3(a) Fe(OH)3(a) Calcite Gypsum SiO2(a) Illite
        -pe false
        -time true
        -step false
        -state false
TRANSPORT
      -shifts 2500
END