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

	<model_description name="unconfined_head_head_no_recharge_2d">
		<comments>Steady-State Richards: unconfined_head_head_no_recharge_2d</comments>
		<model_id>1SSUncFlow1/PORFLOW 4.2.1</model_id>
		<author>Greg Flach</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>3.0e+9</time>
			</time_macro>
		</macros>	
	</definitions>

	<process_kernels>
		<comments>Variably saturated flow</comments>
      		<flow state = "on" model = "richards"/>
      		<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="3.0e+9" init_dt="0.01" mode = "steady" method="bdf1" />
	</execution_controls>

	<numerical_controls>
	        <linear_solver>
			<max_iterations>10</max_iterations>
			<tolerance>1.0e-18</tolerance>
		</linear_solver>
	</numerical_controls>

        <mesh framework="mstk">
		<comments>Two-dimensional box 30.48m x 18.288m (100ft x 60ft)</comments>
		<dimension>2</dimension>
		<generate>
			<number_of_cells nx = "100" ny = "60"/>
			<box low_coordinates = "0.0, 0.0" high_coordinates = "30.48, 18.288"/>
		</generate>
	</mesh>

	<regions>
  		<comments></comments>
		<region name ="Aquifer">
			<comments>One region comprising the entire domain</comments>
			<box low_coordinates= "0.0, 0.0"   high_coordinates="30.48, 18.288"/>
		</region>
		<region name ="Left">
			<comments>Left boundary</comments>
			<box low_coordinates= "0.0, 0.0"   high_coordinates="-1.0, 0.0"/>
		</region>
		<region name ="Right">
			<comments>Right boundary</comments>
			<box low_coordinates= "30.48, 0.0" high_coordinates="1.0, 0.0"/>
		</region>
		<point name="Point0"   coordinate=" 0.1524, 0.1524"/>
		<point name="Point10"  coordinate=" 3.2004, 0.1524"/>
		<point name="Point20"  coordinate=" 6.2484, 0.1524"/>
		<point name="Point30"  coordinate=" 9.2964, 0.1524"/>
		<point name="Point40"  coordinate="12.3444, 0.1524"/>
		<point name="Point50"  coordinate="15.3924, 0.1524"/>
		<point name="Point60"  coordinate="18.4404, 0.1524"/>
		<point name="Point70"  coordinate="21.4884, 0.1524"/>
		<point name="Point80"  coordinate="24.5364, 0.1524"/>
		<point name="Point90"  coordinate="27.5844, 0.1524"/>
		<point name="Point100" coordinate="30.3276, 0.1524"/>
	</regions>
	
	<materials>
		<material name="Aquifer">
		       	<comments>Aquifer</comments>
               		<mechanical_properties>
                         	<porosity value="0.43"/>
 	 	      	</mechanical_properties>
 	 	        <permeability x="3.1198e-11" y="3.1198e-10"/>
			<cap_pressure model="van_genuchten">
				<parameters alpha="1.46e-3" sr="0.052" m="0.314"/>
			</cap_pressure>
			<rel_perm model="mualem">
			</rel_perm>
 	 	      	<assigned_regions>Aquifer</assigned_regions>
 	   	</material>
	</materials>

	<initial_conditions>
		<initial_condition name="Initial Condition">
			<comments>Aquifer</comments>
			<assigned_regions>Aquifer</assigned_regions>
			<liquid_phase name = "water">
				<liquid_component name = "water">
					<pressure value = "101325.0" function="uniform"/>
				</liquid_component>
			</liquid_phase>
         	</initial_condition>
	</initial_conditions>

	<boundary_conditions>
		<comments></comments>
		<boundary_condition name = "LeftBC">
			<comments>Boundary condition at x=0</comments>
			<assigned_regions>Left</assigned_regions>
			<liquid_phase name = "water">
				<liquid_component name = "water">
					<hydrostatic start="0.0" function= "constant" value="12.192"/>
				</liquid_component>
			</liquid_phase>
		</boundary_condition>

		<boundary_condition name = "RightBC">
			<comments>Boundary condition at x=L</comments>
			<assigned_regions>Right</assigned_regions>
			<liquid_phase name = "water">
				<liquid_component name = "water">
					<hydrostatic start="0.0" function= "constant" value="6.096"/>
				</liquid_component>
			</liquid_phase>
		</boundary_condition>
	</boundary_conditions>

	<output>
	        <observations>
		        <filename>observation.out</filename>
			<liquid_phase name="water">
			        <hydraulic_head>
				        <assigned_regions>Point0</assigned_regions>
				        <functional>point</functional>
				        <time_macro>Steady State</time_macro>
				</hydraulic_head>
			        <hydraulic_head>
				        <assigned_regions>Point10</assigned_regions>
				        <functional>point</functional>
				        <time_macro>Steady State</time_macro>
				</hydraulic_head>
			        <hydraulic_head>
				        <assigned_regions>Point20</assigned_regions>
				        <functional>point</functional>
				        <time_macro>Steady State</time_macro>
				</hydraulic_head>
			        <hydraulic_head>
				        <assigned_regions>Point30</assigned_regions>
				        <functional>point</functional>
				        <time_macro>Steady State</time_macro>
				</hydraulic_head>
			        <hydraulic_head>
				        <assigned_regions>Point40</assigned_regions>
				        <functional>point</functional>
				        <time_macro>Steady State</time_macro>
				</hydraulic_head>
			        <hydraulic_head>
				        <assigned_regions>Point50</assigned_regions>
				        <functional>point</functional>
				        <time_macro>Steady State</time_macro>
				</hydraulic_head>
			        <hydraulic_head>
				        <assigned_regions>Point60</assigned_regions>
				        <functional>point</functional>
				        <time_macro>Steady State</time_macro>
				</hydraulic_head>
			        <hydraulic_head>
				        <assigned_regions>Point70</assigned_regions>
				        <functional>point</functional>
				        <time_macro>Steady State</time_macro>
				</hydraulic_head>
			        <hydraulic_head>
				        <assigned_regions>Point80</assigned_regions>
				        <functional>point</functional>
				        <time_macro>Steady State</time_macro>
				</hydraulic_head>
			        <hydraulic_head>
				        <assigned_regions>Point90</assigned_regions>
				        <functional>point</functional>
				        <time_macro>Steady State</time_macro>
				</hydraulic_head>
			        <hydraulic_head>
				        <assigned_regions>Point100</assigned_regions>
				        <functional>point</functional>
				        <time_macro>Steady State</time_macro>
				</hydraulic_head>

			        <aqueous_pressure>
				        <assigned_regions>Point0</assigned_regions>
				        <functional>point</functional>
				        <time_macro>Steady State</time_macro>
				</aqueous_pressure>
			        <aqueous_pressure>
				        <assigned_regions>Point10</assigned_regions>
				        <functional>point</functional>
				        <time_macro>Steady State</time_macro>
				</aqueous_pressure>
			        <aqueous_pressure>
				        <assigned_regions>Point20</assigned_regions>
				        <functional>point</functional>
				        <time_macro>Steady State</time_macro>
				</aqueous_pressure>
			        <aqueous_pressure>
				        <assigned_regions>Point30</assigned_regions>
				        <functional>point</functional>
				        <time_macro>Steady State</time_macro>
				</aqueous_pressure>
			        <aqueous_pressure>
				        <assigned_regions>Point40</assigned_regions>
				        <functional>point</functional>
				        <time_macro>Steady State</time_macro>
				</aqueous_pressure>
			        <aqueous_pressure>
				        <assigned_regions>Point50</assigned_regions>
				        <functional>point</functional>
				        <time_macro>Steady State</time_macro>
				</aqueous_pressure>
			        <aqueous_pressure>
				        <assigned_regions>Point60</assigned_regions>
				        <functional>point</functional>
				        <time_macro>Steady State</time_macro>
				</aqueous_pressure>
			        <aqueous_pressure>
				        <assigned_regions>Point70</assigned_regions>
				        <functional>point</functional>
				        <time_macro>Steady State</time_macro>
				</aqueous_pressure>
			        <aqueous_pressure>
				        <assigned_regions>Point80</assigned_regions>
				        <functional>point</functional>
				        <time_macro>Steady State</time_macro>
				</aqueous_pressure>
			        <aqueous_pressure>
				        <assigned_regions>Point90</assigned_regions>
				        <functional>point</functional>
				        <time_macro>Steady State</time_macro>
				</aqueous_pressure>
			        <aqueous_pressure>
				        <assigned_regions>Point100</assigned_regions>
				        <functional>point</functional>
				        <time_macro>Steady State</time_macro>
				</aqueous_pressure>
			</liquid_phase>
		</observations>

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

</amanzi_input>

