TITLE Leuchs 2D tritium_helium dating # Note... PHAST 1.1 calculates results with less numerical dispersion than version 1, used for Fig. 11.13 UNITS -time year -horizontal_grid meters -vertical_grid meters -head meters -hydraulic_conductivity m/day -flux m/year -specific_storage 1/m -dispersivity m SOLUTION_METHOD -direct_solver true -space_differencing 0. -time_differencing 1.0 SOLUTE_TRANSPORT true -diffusivity 0.0 #m^2/s STEADY_FLOW true GRID -uniform x 0.0 2500. 26 -uniform y 0.0 1.0 2 -uniform z 0.0 35. 141 -print_orientation XZ -chemistry_dimensions XZ MEDIA -zone 0. 0. 0. 2500. 1. 35. -porosity 0.3 -longitudinal_dispersivity 3.3333 -horizontal_dispersivity 0.03333 -vertical_dispersivity 0.003333 -Kx 40 -Ky 40 -Kz 40 # m/d -storage 0 # for homogeneous aquifer, cross out the following lines... -zone 1600. 0. 11. 2250. 1. 16. -vertical_dispersivity 0.003333e2 -Kx 240 FLUX_BC -zone 0.0 0. 0. 2500. 1. 0. -flux 0 0.3 -associated_solution 0 1950 1 1951 2 1952 3 1953 4 1954 5 1955 6 1956 7 1957 8 1958 9 1959 10 1960 11 1961 12 1962 13 1963 14 1964 15 1965 16 1966 17 1967 18 1968 19 1969 20 1970 21 1971 22 1972 23 1973 24 1974 25 1975 26 1976 27 1977 28 1978 29 1979 30 1980 31 1981 32 1982 33 1983 34 1984 35 1985 36 1986 SPECIFIED_HEAD_BC -zone 2500. 0. 0.2501 2500. 1. 35. -head 0 0.0 -associated_solution 0 1 FREE_SURFACE_BC false HEAD_IC -zone 0. 0. 0. 2500. 1. 35.0 -head X 0.1 0. 2.5 2500.0 CHEMISTRY_IC -zone 0. 0. 0. 2500. 1. 35.0 -solution 1 -kinetics 1 # for C_He is = 0 at the groundwater table ... # -zone 0.0 0. 0. 2500. 1. 0.1 # -equilibrium_phases 1 TIME_CONTROL -time_step; 0 0.25 year -time_end; 37. year PRINT_INITIAL -conduct true -hdf_steady true PRINT_LOCATIONS # Turn off xyz_print -zone 0 0 0 2500. 1. 35.0 -print 0 # Turn on xyz print for DFG1 -xyz_chemistry -zone 2295. 0 0 2305. 0 35.0 -print 1 -zone 1895. 0 0 1905. 0 35.0 -print 1 -zone 1295. 0 0 1305. 0 35.0 -print 1 -zone 295. 0 0 305. 0 35.0 -print 1 PRINT_FREQUENCY; 0 -hdf_chemistry 1 year -xyz_chemistry 37 year END