data = thermo.tdat verify
surface_data = Pb-Kd.sdat
# download the surface dataset and save it in your working directory for this script

interval start at 0 yr, fluid = contaminated
interval interval-2 at 2 yr, fluid = flush
time end at 10 yr

length = 1000 m
Nx = 400

discharge start = 30 m/yr
discharge interval-2 = 30 m/yr

porosity = .3
dispersivity = 1 m

scope = initial
   H2O          = 1 free kg
   Na+          = 10 mmolal
   balance on Cl-
   Pb++         = 1e-12 mmolal
   Br-          = 1e-12 mmolal
scope = contaminated
   H2O          = 1 free kg
   Na+          = 10 mmolal
   balance on Cl-
   Pb++         = .1 mmolal
   Br-          = .1 mmolal
scope = flush
   H2O          = 1 free kg
   Na+          = 10 mmolal
   balance on Cl-
   Pb++         = 1e-12 mmolal
   Br-          = 1e-12 mmolal
scope