Running PHASTINPUT.

Processing flow and transport data file.

	TITLE
	 Central Oklahoma aquifer,
	 demonstration of PHAST
	UNITS
	        time                   years
	        horizontal_grid        meters
	        vertical_grid          meters
	        head                   meters
	        hydraulic_conductivity m/s
	        specific_storage       1/m
	        dispersivity           m
	        leaky_k                m/sec
	        leaky_thickness        m
	        flux                   m/day
	        river_bed_thickness    m
	        river_bed_hydraulic_conductivity       m/s     
	        well_diameter          in
	        well_flow_rate         L/day
	GRID
	        uniform x
	                0       90000 16
	        uniform y
	                0       48000 9
	        uniform z
	                0       400 5
	        print_orientation XY 
	        chemistry_dimensions XYZ
	FLOW_ONLY false
	STEADY_FLOW true
	        head_tolerance 1e-3
	        flow_balance_tolerance 1e-3
	FLUID_PROPERTIES
	        compressibility        0.0             #1/Pa   
	        density                1000.           #kg/m^3
	        viscosity              0.001           #Pa-s
	        diffusivity            0.0             #m^2/s
	MEDIA
	        zone   0. 0. 0. 90000.  48000. 400.
	                porosity               0.22
	                long_dispersivity      4000.
	                trans_dispersivity     50.
	                kx                     1.373e-5
	                ky                     1.373e-5
	                kz                     1.373e-7
	                storage                0
	        zone   48000. 0. 0. 90000. 48000. 100.
	                active 0
	RIVER 1 Canadian River
	        default_width  200.
	        default_depth  1.
	        default_bed_thickness  1.
	        default_bed_hydraulic_conductivity     1
	        default_solution 1
	        30000. 0
	                head   320.
	        90000. 0
	                head 275.
	RIVER 2 North Fork River
	        default_width  200.
	        default_depth  1.
	        default_bed_thickness  1.
	        default_bed_hydraulic_conductivity     1
	        default_solution 1
	        30000. 36000.
	                head   335.
	        30000. 48000.
	                head   320.
	        90000. 48000.
	                head 245.
	RIVER 3 North Canadian River
	        default_width  200.
	        default_depth  1.
	        default_bed_thickness  1.
	        default_bed_hydraulic_conductivity     1
	        default_solution 1
	        60000. 30000.
	                head   335.
	        90000. 20000.
	                head 290.
	FLUX_BC
	        zone   30000. 3000. 400.  90000. 45000. 400.
	                flux                   -10e-5
	                associated_solution         1
	SPECIFIED_VALUE_BC
	        zone   30000. 14000 300. 32000.  20000.  400.
	                head   348.
	                associated_solution_composition         1
	LEAKY_BC
	        zone   0. 48000. 0. 29000. 48000. 400.
	                hydraulic_conductivity      1.618e-5
	                thickness      30000.
	                head           305.0
	                associated_solution                  1
	FREE_SURFACE_BC true
	WELL 1 Observation well 1 in arsenic zone
	        12000 36000
	                diameter       2       # inches
	                pumping        1       # L/day
	                elevation      90 110  # m
	HEAD_IC
	        zone   0. 0. 0. 90000. 48000. 400.
	                head   380.
	CHEMISTRY_IC
	        zone   0. 0. 0. 90000. 48000. 400.
	                solution               2
	                equilibrium_phases     2
	                exchange               2
	                surface                2       
	SOLUTION_METHOD
	        direct_solver          false
	        tolerance              1e-10
	        space_differencing     0.
	        time_differencing      1.
	TIME_CONTROL
	        delta_time     2000    yr
	        end_time       100000  yr
	PRINT_INITIAL
		velocities		true
		xyz_head		true
		xyz_ss_velocities	true
	PRINT_FREQUENCY
		xyz_chemistry	        50000   yr
		hdf_chemistry	        2000    yr
	        xyz_well	        2000    yr
		save_final_heads	true
	PRINT_LOCATIONS
	       zone   0 0 0  90000 48000 0
	       print 0
	END
WARNING: Fluid compressibility has been set to zero for free surface calculation.
WARNING: Aquifer compressibility has been set to zero for free surface calculation.

PHASTINPUT done.

Running PHAST.

Running PHREEQC for initial conditions.

Processing database file.
	SOLUTION_MASTER_SPECIES
	SOLUTION_SPECIES
	PHASES
	EXCHANGE_MASTER_SPECIES
	EXCHANGE_SPECIES
	SURFACE_MASTER_SPECIES
	SURFACE_SPECIES
	END

Processing chemical data file.
	SURFACE_MASTER_SPECIES
	        Surf    SurfOH
	SURFACE_SPECIES
	SOLUTION_MASTER_SPECIES
	        As       H3AsO4        -1.0     74.9216         74.9216
	SOLUTION_SPECIES
	SOLUTION 1 20 x precipitation
	END
	USE solution 1
	EQUILIBRIUM_PHASES 1
	        Dolomite        0.0     2.7
	        Calcite         0.0     1.0
	        CO2(g)          -1.5    100.
	SAVE solution 1
	END
	SOLUTION 2 Brine
	END
	USE solution 2
	EQUILIBRIUM_PHASES 2
	        Dolomite        0.0     2.7
	        Calcite         0.0     1.0
	SAVE solution 2
	END
	EXCHANGE 2
	SURFACE 2
	END
	SELECTED_OUTPUT
	USER_PUNCH
	END
PHREEQC done.
. Central Oklahoma aquifer,
. demonstration of PHAST
reading static data for flow and transport simulation
WARNING: Multiple free surfaces in column of cells; Check initial condition on head field.
A4 array dimension requirement (D4ZGCG iterative solver)......... NRAL .        0 elements
Primary storage requirement (D4ZGCG iterative solver)............ NPRIST     6840 elements
Overhead storage requirement (D4ZGCG iterative solver)........... NOHST .   34264 elements
                    *** Calculation Information ***
          Direction index for d4 zig-zag renumbering ............ IDIR..    1
          Incomplete LU [f] or modified ILU [t] factorization ... MILU..     T
          Number of search directions before restart ............ NSDR..    5
          Tolerance on iterative solution ....................... EPSSLV 1.0E-10

Beginning flow calculation for steady-state i.c.
Reading transient data for simulation: Set  1
WARNING: 
Automatic Time Step Control Parameters for Steady Flow
Maximum pressure change allowed per time step.................... DPTAS    1.18E+06  (Pa )
Minimum time step required ...................................... DTIMMN    2.0      (yr )
Maximum time step allowed ....................................... DTIMMX   2.00E+06(yr )

Beginning steady state time step no.      1
     Maximum change in potentiometric head ..   -1.3312E+02     (m) at location ( 9.000E+04, 4.800E+04,  400.    )(m)
     Fractional flow residual ...............   -1.9439E+00

Beginning steady state time step no.      2
     Maximum change in potentiometric head ..   -2.3095E+01     (m) at location ( 9.000E+04, 4.800E+04,  200.    )(m)
     Fractional flow residual ...............   -1.9028E+00

Beginning steady state time step no.      3
     Maximum change in potentiometric head ..   -3.8663E+00     (m) at location ( 7.800E+04, 1.800E+04,  400.    )(m)
     Fractional flow residual ...............   -1.4759E+00

Beginning steady state time step no.      4
     Maximum change in potentiometric head ..   -6.3383E+00     (m) at location ( 7.800E+04, 1.800E+04,  400.    )(m)
     Fractional flow residual ...............   -1.0087E+00

Beginning steady state time step no.      5
     Maximum change in potentiometric head ..   -9.4795E+00     (m) at location ( 9.000E+04, 6.000E+03,  400.    )(m)
     Fractional flow residual ...............   -9.5331E-01

Beginning steady state time step no.      6
     Maximum change in potentiometric head ..   -1.3931E+01     (m) at location ( 9.000E+04, 6.000E+03,  400.    )(m)
     Fractional flow residual ...............   -8.6088E-01

Beginning steady state time step no.      7
     Maximum change in potentiometric head ..   -1.6528E+01     (m) at location ( 9.000E+04, 4.200E+04,  300.    )(m)
     Fractional flow residual ...............   -7.1811E-01

Beginning steady state time step no.      8
     Maximum change in potentiometric head ..   -1.7030E+01     (m) at location ( 9.000E+04, 3.600E+04,  300.    )(m)
     Fractional flow residual ...............   -5.2168E-01

Beginning steady state time step no.      9
     Maximum change in potentiometric head ..   -1.1431E+01     (m) at location ( 9.000E+04, 3.600E+04,  300.    )(m)
     Fractional flow residual ...............   -3.0728E-01

Beginning steady state time step no.     10
     Maximum change in potentiometric head ..   -1.1179E+01     (m) at location ( 0.000    , 6.000E+03,  400.    )(m)
     Fractional flow residual ...............   -1.4836E-01

Beginning steady state time step no.     11
     Maximum change in potentiometric head ..   -7.3995E+00     (m) at location ( 0.000    , 1.800E+04,  400.    )(m)
     Fractional flow residual ...............   -7.0581E-02

Beginning steady state time step no.     12
     Maximum change in potentiometric head ..   -5.4899E+00     (m) at location ( 0.000    , 0.000    ,  300.    )(m)
     Fractional flow residual ...............   -2.5324E-02

Beginning steady state time step no.     13
     Maximum change in potentiometric head ..   -1.4528E+00     (m) at location ( 0.000    , 1.200E+04,  300.    )(m)
     Fractional flow residual ...............   -4.0607E-03

Beginning steady state time step no.     14
     Maximum change in potentiometric head ..   -1.3801E-01     (m) at location ( 0.000    , 1.800E+04,  300.    )(m)
     Fractional flow residual ...............   -5.7376E-04

Beginning steady state time step no.     15
     Maximum change in potentiometric head ..   -1.1866E-02     (m) at location ( 0.000    , 3.000E+04,  300.    )(m)
     Fractional flow residual ...............   -4.1143E-05

Beginning steady state time step no.     16
     Maximum change in potentiometric head ..    9.0390E-04     (m) at location ( 2.400E+04, 1.200E+04,  300.    )(m)
     Fractional flow residual ...............   -5.3174E-06

Beginning steady state time step no.     17
     Maximum change in potentiometric head ..    2.7542E-05     (m) at location ( 0.000    , 1.200E+04,  300.    )(m)
     Fractional flow residual ...............   -2.7135E-09

Beginning steady state time step no.     18
     Current time step length ...................................   1.311E+05 (yr)
          No. of solver iterations, Relative residual:    3  3.4161738E-11
     Maximum change in potentiometric head ..    1.3214E-06     (m) at location ( 4.800E+04, 1.800E+04,  400.    )(m)
     Fractional flow residual ...............    2.8517E-09

Beginning steady state time step no.     19
     Current time step length ...................................   2.621E+05 (yr)
          No. of solver iterations, Relative residual:    3  5.1382054E-11
     Maximum change in potentiometric head ..   -3.5555E-08     (m) at location ( 3.600E+04, 3.000E+04,  400.    )(m)
     Fractional flow residual ...............    9.8500E-11

Done with steady-state flow.
Equilibration of cells for initial conditions.
          Total cells: 657
          Number of cells with same aqueous model: 529
          Total iterations all cells: 0
          Number of cells with zero iterations: 530
          Maximum iterations for one cell: 0
     Time at which next set of transient parameters will be read .......... TIMCHG  1.000E+05  (yr)

Beginning transient simulation.

Beginning time step no.      1
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.     1
     Time ..................................................   2.000E+03 (yr)

Beginning time step no.      2
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.     2
     Time ..................................................   4.000E+03 (yr)

Beginning time step no.      3
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.     3
     Time ..................................................   6.000E+03 (yr)

Beginning time step no.      4
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.     4
     Time ..................................................   8.000E+03 (yr)

Beginning time step no.      5
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.     5
     Time ..................................................   1.000E+04 (yr)

Beginning time step no.      6
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.     6
     Time ..................................................   1.200E+04 (yr)

Beginning time step no.      7
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.     7
     Time ..................................................   1.400E+04 (yr)

Beginning time step no.      8
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.     8
     Time ..................................................   1.600E+04 (yr)

Beginning time step no.      9
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.     9
     Time ..................................................   1.800E+04 (yr)

Beginning time step no.     10
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    10
     Time ..................................................   2.000E+04 (yr)

Beginning time step no.     11
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    11
     Time ..................................................   2.200E+04 (yr)

Beginning time step no.     12
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    12
     Time ..................................................   2.400E+04 (yr)

Beginning time step no.     13
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    13
     Time ..................................................   2.600E+04 (yr)

Beginning time step no.     14
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    14
     Time ..................................................   2.800E+04 (yr)

Beginning time step no.     15
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    15
     Time ..................................................   3.000E+04 (yr)

Beginning time step no.     16
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    16
     Time ..................................................   3.200E+04 (yr)

Beginning time step no.     17
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    17
     Time ..................................................   3.400E+04 (yr)

Beginning time step no.     18
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    18
     Time ..................................................   3.600E+04 (yr)

Beginning time step no.     19
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    19
     Time ..................................................   3.800E+04 (yr)

Beginning time step no.     20
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    20
     Time ..................................................   4.000E+04 (yr)

Beginning time step no.     21
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    21
     Time ..................................................   4.200E+04 (yr)

Beginning time step no.     22
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    22
     Time ..................................................   4.400E+04 (yr)

Beginning time step no.     23
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    23
     Time ..................................................   4.600E+04 (yr)

Beginning time step no.     24
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    24
     Time ..................................................   4.800E+04 (yr)

Beginning time step no.     25
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    25
     Time ..................................................   5.000E+04 (yr)

Beginning time step no.     26
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    26
     Time ..................................................   5.200E+04 (yr)

Beginning time step no.     27
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    27
     Time ..................................................   5.400E+04 (yr)

Beginning time step no.     28
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    28
     Time ..................................................   5.600E+04 (yr)

Beginning time step no.     29
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    29
     Time ..................................................   5.800E+04 (yr)

Beginning time step no.     30
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    30
     Time ..................................................   6.000E+04 (yr)

Beginning time step no.     31
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    31
     Time ..................................................   6.200E+04 (yr)

Beginning time step no.     32
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    32
     Time ..................................................   6.400E+04 (yr)

Beginning time step no.     33
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    33
     Time ..................................................   6.600E+04 (yr)

Beginning time step no.     34
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    34
     Time ..................................................   6.800E+04 (yr)

Beginning time step no.     35
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    35
     Time ..................................................   7.000E+04 (yr)

Beginning time step no.     36
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    36
     Time ..................................................   7.200E+04 (yr)

Beginning time step no.     37
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    37
     Time ..................................................   7.400E+04 (yr)

Beginning time step no.     38
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    38
     Time ..................................................   7.600E+04 (yr)

Beginning time step no.     39
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    39
     Time ..................................................   7.800E+04 (yr)

Beginning time step no.     40
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    40
     Time ..................................................   8.000E+04 (yr)

Beginning time step no.     41
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    41
     Time ..................................................   8.200E+04 (yr)

Beginning time step no.     42
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    42
     Time ..................................................   8.400E+04 (yr)

Beginning time step no.     43
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    43
     Time ..................................................   8.600E+04 (yr)

Beginning time step no.     44
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    44
     Time ..................................................   8.800E+04 (yr)

Beginning time step no.     45
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    45
     Time ..................................................   9.000E+04 (yr)

Beginning time step no.     46
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    46
     Time ..................................................   9.200E+04 (yr)

Beginning time step no.     47
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    47
     Time ..................................................   9.400E+04 (yr)

Beginning time step no.     48
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    48
     Time ..................................................   9.600E+04 (yr)

Beginning time step no.     49
     Beginning solute-transport calculation.
          H
          O
          As
          C
          Ca
          Cl
          Mg
          Na
          S
          Charge
     Beginning chemistry calculation.
     End of Time Step No.    49
     Time ..................................................   9.800E+04 (yr)

Beginning time step no.     50
     Current time step length ...................................   2.000E+03 (yr)
     Beginning solute-transport calculation.
          H
          No. of solver iterations, Relative residual:    7  9.9006280E-11
          O
          No. of solver iterations, Relative residual:    7  7.1245815E-11
          As
          No. of solver iterations, Relative residual:    7  9.8553809E-11
          C
          No. of solver iterations, Relative residual:    7  4.7869301E-11
          Ca
          No. of solver iterations, Relative residual:    7  4.7338936E-11
          Cl
          No. of solver iterations, Relative residual:    6  2.9764171E-11
          Mg
          No. of solver iterations, Relative residual:    7  5.6332877E-11
          Na
          No. of solver iterations, Relative residual:    7  1.8168590E-11
          S
          No. of solver iterations, Relative residual:    6  3.2690907E-11
          Charge
          No. of solver iterations, Relative residual:    7  1.8597661E-11
     Beginning chemistry calculation.
          Total cells: 657
          Number of cells with same aqueous model: 529
          Total iterations all cells: 3040
          Number of cells with zero iterations: 0
          Maximum iterations for one cell: 15
     End of Time Step No.    50
     Time ..................................................   1.000E+05 (yr)
     Maximum change in potentiometric head  -3.5555E-08      (m ) at location ( 3.600E+04, 3.000E+04,  400.    )(m)
     Maximum change in H .................. -6.9092E-04 (mol/kgw) at location ( 4.200E+04, 2.400E+04, 0.000    )(m)
     Maximum change in O .................. -3.0076E-03 (mol/kgw) at location ( 6.000E+03, 4.800E+04,  100.    )(m)
     Maximum change in As ................. -4.5911E-07 (mol/kgw) at location ( 2.400E+04, 4.200E+04, 0.000    )(m)
     Maximum change in C .................. -1.2292E-03 (mol/kgw) at location ( 6.000E+03, 4.800E+04, 0.000    )(m)
     Maximum change in Ca ................. -3.6729E-03 (mol/kgw) at location ( 0.000    , 6.000E+03,  300.    )(m)
     Maximum change in Cl ................. -5.6199E-02 (mol/kgw) at location ( 6.000E+03, 6.000E+03, 0.000    )(m)
     Maximum change in Mg ................. -1.8152E-03 (mol/kgw) at location ( 0.000    , 1.200E+04,  300.    )(m)
     Maximum change in Na ................. -4.5587E-02 (mol/kgw) at location ( 6.000E+03, 6.000E+03, 0.000    )(m)
     Maximum change in S .................. -3.7233E-05 (mol/kgw) at location ( 6.000E+03, 6.000E+03, 0.000    )(m)
     Maximum change in Charge ............. -1.9578E-03 (mol/kgw) at location ( 6.000E+03, 4.800E+04, 0.000    )(m)
Done with transient simulation.

                             ***** Simulation Completed *****
     Last time value calculated ...................................... 1.0000E+05 (yr)
     Last time step index ............................................      50