SOLUTION_MASTER_SPECIES; Tetra Tetra 0.0 Tetra 1.0
SOLUTION_SPECIES; Tetra = Tetra; log_k 0.0
SURFACE_MASTER_SPECIES; Sor Sor
SURFACE_SPECIES
 Sor = Sor; log_k 0.0
 Sor + Tetra = SorTetra; log_k -99.52              # log(Kf = 3) - log(m_Sor)

SOLUTION 1-201                                     # The sand and clay aquifer

SURFACE  102-201                             # Sorption in stagnant cells (clay layer)
 Sor 1e100 1 1
END

SOLUTION 0; Cl 1; Tetra 1                          # The polluted water
PRINT; -reset false
TRANSPORT; -cells 100; -shifts 25 1; -timest 0.73e7      # 5 years pollution
 -length 10
# include dispersion...
# -disp 4.999
#  and stagnant zones...
# -stagnant 1 1.55e-9 0.24 0.12     # 1 stagnant layer, alpha, mobile por, immobile por
END

SOLUTION 0                                         # Clean water enters
TRANSPORT; -shifts 50 1                            # 10 years cleaning
 -punch_fr 50                                      # punch/graph only 50th shift
 -punch  1-100                                     # punchout mobile flow
# -punch 102-201                                   # punchout stagnant area
USER_GRAPH
 -initial_solutions false
 -headings  Dist CCl4 Cl
 -axis_titles "Distance / m" "mmol / L"
 -start
 10 graph_x Dist
 20 graph_y tot("Tetra") * 1e3, tot("Cl") * 1e3
 -end
END