data = thermo.tdat verify length = 100 m Nx = 40 time end = 100 years discharge = 100 pore_volumes dispersivity = 10 cm scope = initial swap Quartz for SiO2(aq) Quartz = 70 volume% scope = inlet SiO2(aq) = 1 mg/kg scope kinetic Quartz rate_con = 4.2e-18 surface = 1000