<amanzi_input version="2.0.0" type="unstructured">

	<model_description name="SR Test Problem">
		<comments>Getting started on a simple 2D steady-state saturated test</comments>
		<model_id>Single-Phase Saturated</model_id>
		<author>Elizabeth Keating</author>
		<units>
			<length_unit>m</length_unit>
			<time_unit>s</time_unit>
			<mass_unit>kg</mass_unit>
			<conc_unit>molar</conc_unit>
		</units>
	</model_description>

	<definitions>
		<macros>
			<time_macro name="Steady State">
				<time>0.0</time>
			</time_macro>
		</macros>	
	</definitions>

	<process_kernels>
		<comments>Single phase flow only</comments>
      		<flow state = "on" model = "saturated"/>
      		<transport state = "off" algorithm="none"/>
      		<chemistry state = "off" process_model="none"/>
	</process_kernels>

	<phases>
		<liquid_phase name="water">
			<eos>false</eos>
			<viscosity>1.002E-03</viscosity>
			<density>998.2</density>
		</liquid_phase>
	</phases>

	<execution_controls>
		<verbosity level="high"/>
		<execution_control_defaults mode="steady" method="bdf1"/>
		<execution_control  start="0.0" end="0.0" mode = "steady"  method="bdf1"/>
	</execution_controls>

	<numerical_controls>
		<comments>Many of these were not in Elizabeth's original file</comments>
		<steady-state_controls>
			<comments>Steady-state controls comments here</comments>
          		<min_iterations>10</min_iterations>
		      	<max_iterations>20</max_iterations>
          		<max_preconditioner_lag_iterations>30</max_preconditioner_lag_iterations>
          		<nonlinear_tolerance>1.0e-5</nonlinear_tolerance>
          		<pseudo_time_integrator>
               			<method>picard</method>
          			<preconditioner>hypre_amg</preconditioner>
          			<linear_solver>aztec00</linear_solver>
          			<control_options>pressure</control_options>
          			<max_iterations>100</max_iterations>
          			<clipping_saturation>0.9</clipping_saturation>
          			<convergence_tolerance>1.0e-10</convergence_tolerance>
          			<initialize_with_darcy>true</initialize_with_darcy>
          		</pseudo_time_integrator>
		</steady-state_controls>

		<linear_solver>
		        <comments>linear solver method not specified in F-Area file; what is default? gmres?
			couldn't find "method" in spec so deleted it -vlf</comments>
		        <!-- max iterations for gmres? -->
			<max_iterations>10</max_iterations>
			<tolerance>1.0e-18</tolerance>
			<preconditioner name = "hypre_amg">
				<hypre_cycle_applications>10</hypre_cycle_applications>
			        <hypre_smoother_sweeps>3</hypre_smoother_sweeps>
			        <hypre_tolerance>0.1</hypre_tolerance>
			        <hypre_strong_threshold>0.4</hypre_strong_threshold>
			</preconditioner>
		</linear_solver>

	</numerical_controls>

        <mesh framework="mstk">
		<comments>Three-dimensional box 100m x 50m x 50m</comments>
		<dimension>3</dimension>
		<generate>
			<number_of_cells nx = "100"  ny = "5"  nz = "50"/>
			<box  low_coordinates = "0.0,0.0,0.0" high_coordinates = "100.0,50.0,50.0"/>
		</generate>
	</mesh>

	<regions>
  		<comments></comments>
		<region name ="Entire Domain">
			<comments>One region comprising the entire domain</comments>
			<box low_coordinates= "0.0,0.0,0.0"     high_coordinates="100.0,50.0,50.0"/>
		</region>
		<region name ="Upstream">
			<comments>Upstream boundary</comments>
			<box low_coordinates= "0.0,0.0,0.0"     high_coordinates="0.0,50.0,50.0"/>
		</region>
		<region name ="Downstream">
			<comments>Downstream boundary</comments>
			<box low_coordinates= "100.0,0.0,0.0"     high_coordinates="100.0,50.0,50.0"/>
		</region>
		<point name="Well 1" coordinate="(0.5,25.0,49.5)"/>
		<point name="Well 2" coordinate="(50.5,25.0,25.5)"/>
		<point name="Well 2t" coordinate="(50.5,25.0,49.5)"/>
		<point name="Well 2b" coordinate="(50.5,25.0,0.5)"/>
		<point name="Well 3" coordinate="(99.5,25.0,49.5)"/>
	</regions>
	
	<materials>
		<material name="Soil">
		       	<comments>Homogeneous material properties throughout the entire domain</comments>
               		<mechanical_properties>
                         	<porosity value="0.25"/>
 	 	      	</mechanical_properties>
 	 	        <permeability x="1.0E-12" y="1.0E-12" z="1.0E-12"/>
 	 	      	<assigned_regions>Entire Domain</assigned_regions>
 	   	</material>
	</materials>

	<initial_conditions>
		<initial_condition name="Initial Condition">
			<comments>Uniform pressure over entire domain</comments>
			<assigned_regions>Entire Domain</assigned_regions>
			<liquid_phase name = "water">
				<liquid_component name = "water">
					<uniform_pressure value = "5001325.0"/>
				</liquid_component>
			</liquid_phase>
         	</initial_condition>
	</initial_conditions>

	<boundary_conditions>
		<comments></comments>
		<boundary_condition name = "Upstream BC">
			<comments>Upstream boundary condition</comments>
			<assigned_regions>Upstream</assigned_regions>
			<liquid_phase name = "water">
				<liquid_component name = "water">
					<inward_mass_flux start="0.0" function= "constant" value="1.95e-2"/>
				</liquid_component>
			</liquid_phase>
		</boundary_condition>

		<boundary_condition name = "Downstream BC">
			<assigned_regions>Downstream</assigned_regions>
			<liquid_phase name = "water">
				<liquid_component name = "water">
					<hydrostatic start="0.0" function= "constant" value="120.0"/>
				</liquid_component>
			</liquid_phase>
		</boundary_condition>
	</boundary_conditions>

	<output>
		<vis>
			<base_filename>steady-flow</base_filename>
		  	<num_digits>5</num_digits>
		  	<time_macro>Steady State</time_macro>
		</vis>
	</output>

</amanzi_input>

