====== Transport Modeling with Mixing ======

<code>
TITLE Redox Reaction with transport

SOLUTION 0 # Recharge
        units   mmol/L
        pH      7.0
        temp    25.0
        Ca	 0.1
        Mg	 0.1
	  Na	 0.1 
        K	 0.11
        Alkalinity 0.4 as HCO3
        S(6) 0.05 as SO4
        Cl	 0.1 
        Zn   0.01

REACTION 0
        CH20 1
        0.000001 moles
SAVE SOLUTION 0

SOLUTION 1-10 # Aquifer
        units   ppm
        pH      8.22
        temp    25.0                  
        S(6) 3.3 as SO4
	  N(5) 0.58 as NO3
        P(5) 0.1
        Cl	 0.22 
	  Na	 0.38 
        K	 0.36
        Ca	 0.20
        Mg	 0.30
        Alkalinity 0.68 as HCO3

MIX 1 
    0    0.01
    1    0.99
MIX 2
    0    0.01
    2    0.99                   
MIX 3
    0    0.01
    3    0.99
MIX 4
    0    0.01
    4    0.99
MIX 5
    0    0.01
    5    0.99
MIX 6
    0    0.01
    6    0.99
MIX 7
    0    0.01
    7    0.99
MIX 8
    0    0.01
    8    0.99
MIX 9
    0    0.01
    9    0.99
MIX 10
    0    0.01
    10   0.99

TRANSPORT 
   -cells  10
   -shifts 10
   -flow_direction  forward 
   -time_step       86400*365
   -boundary_conditions   flux  flux 
   -lengths         10*100 
   -dispersivities  10*0.1
   -punch_cells	  5
   -punch_frequency 1

SELECTED_OUTPUT 
        -file aquifer.txt 
   	  -step true
        -pH true
        -pe true
        -reaction true
	  -temperature true
        -alkalinity true
        -si Calcite 
        -totals Ca Na Cl S(6) N(5) Na K Zn
        -molalities Ca+2 Na+ Cl- N2
END
</code>