ex13a

TITLE Example 13A.--1 mmol/l NaCl/NO3 enters column with stagnant zones.
                    Implicit definition of first-order exchange model.
SOLUTION 0    # 1 mmol/l NaCl
        units   mmol/l
        pH       7.0
        pe      13.0    O2(g)   -0.7
        Na       1.0    # Na has Retardation = 2
        Cl       1.0    # Cl has Retardation = 1, stagnant exchange
        N(5)     1.0    # NO3 is conservative
#       charge imbalance is no problem ...
END
SOLUTION 1-41  # Column with KNO3
        units   mmol/l
        pH       7.0
        pe      13.0   O2(g)    -0.7
        K        1.0
        N(5)     1.0
EXCHANGE 1-41
        -equil  1
        X       1.e-3
EXCHANGE_SPECIES # For linear exchange, make KX exch. coeff. equal to NaX
        K+ + X- = KX
        log_k   0.0
        -gamma  3.5     0.015
END
TRANSPORT
        -cells  20
        -shifts 5
        -flow_d   forward
        -timest 3600
        -bcon   flux  flux
        -diffc  0.0
        -length 0.1
        -disp   0.015
        -stag       1  6.8e-6  0.3        0.1
#   1 stagnant layer^, ^alpha, ^theta(m), ^theta(im)
PRINT
        -reset false
	-echo true
END
SOLUTION 0  # Original solution reenters
        units   mmol/l
        pH       7.0
        pe      13.0   O2(g)    -0.7
        K        1.0
        N(5)     1.0
TRANSPORT
        -shifts 10
        -punch_frequency    10
        -punch_cells        1-20
SELECTED_OUTPUT
        -file   ex13a.sel
        -reset  false
        -solution
        -distance       true
USER_PUNCH
        -head Cl_mmol Na_mmol
10 PUNCH TOT("Cl")*1000, TOT("Na")*1000
END