TITLE Do buffer problem using reactive transport
        instead of titration
SOLUTION        0
        temp        25
        units       mol/kgw
        pH          2        charge
        S(6)        .01

SOLUTION        1-40
        temp        25
        units       mol/kgw
        ph          7
        Na          .1
        Cl          .1        charge
        S(6)         1e-7
        Si           1e-7
EQUILIBRIUM_PHASES        1-40
        CO2(g)      -3.5
        O2(g)       -0.7
        Fe(OH)3(a)   0.0        0.01
        Calcite      0.0        0.08
        Kaolinite    0.0        0.03
        Quartz       0.6        0.10
        Gypsum       0.0        0.0
END
TRANSPORT
# 20 m/shift, 15 m/y is 20/15 y/shift
# 31556925 s/y *20/15 y/shift = 42075900 sec/shift
# 200 shifts in 40 cells = 5 pore volumes
# 200 shifts * 1.333 y/shift = 267 y total flow time
        -cells          40
        -lengths        40*20
        -time           42075900
        -shifts         200
        -dispersivity   40*10
        -correct_disp   true
        -punch_cells    1-40
        -punch_frequency        10
PRINT
        reset                false
SELECTED_OUTPUT
        -file            buffering
        -simulation      false
        -solution        false
        -state           false
        -pe              false
        -si              Calcite Fe(OH)3(a) Kaolinite Gypsum Quartz
        -equilibrium_phases Calcite Fe(OH)3(a) Kaolinite Gypsum Quartz
END