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

  <model_description name="SR Test Problem">
    <comments>Steady-State Saturated: linear pressure profile</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>

  <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>
      <density>998.2</density>
      <viscosity>1.002E-03</viscosity>
    </liquid_phase>
  </phases>

  <execution_controls>
    <verbosity level="high"/>
    <execution_control_defaults mode="steady" method="bdf1" />
    <execution_control  start="0.0" end="0.0" init_dt="1000" mode = "steady" method="bdf1" />
  </execution_controls>
  
  <numerical_controls>
    <linear_solver>
      <max_iterations>30</max_iterations>
      <tolerance>1.0e-18</tolerance>
    </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, 25.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, 25.5)"/>
  </regions>

  <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">
	  <pressure value = "101325.0" function="uniform"/>
	</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>
  
  <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>

  <definitions>
    <macros>
      <time_macro name="Steady State">
	<time>0.0</time>
      </time_macro>
    </macros>	
  </definitions>
  
  <output>
    <observations>
      <filename>observation.out</filename>
      <liquid_phase name="water">
	<hydraulic_head>
	  <assigned_regions>Well 1</assigned_regions>
	  <functional>point</functional>
	  <time_macro>Steady State</time_macro>
	</hydraulic_head>
	<hydraulic_head>
	  <assigned_regions>Well 2</assigned_regions>
	  <functional>point</functional>
	  <time_macro>Steady State</time_macro>
	</hydraulic_head>
	<hydraulic_head>
	  <assigned_regions>Well 2t</assigned_regions>
	  <functional>point</functional>
	  <time_macro>Steady State</time_macro>
	</hydraulic_head>
	<hydraulic_head>
	  <assigned_regions>Well 2b</assigned_regions>
	  <functional>point</functional>
	  <time_macro>Steady State</time_macro>
	</hydraulic_head>
	<hydraulic_head>
	  <assigned_regions>Well 3</assigned_regions>
	  <functional>point</functional>
	  <time_macro>Steady State</time_macro>
	</hydraulic_head>
	
	<aqueous_pressure>
	  <assigned_regions>Well 1</assigned_regions>
	  <functional>point</functional>
	  <time_macro>Steady State</time_macro>
	</aqueous_pressure>
	<aqueous_pressure>
	  <assigned_regions>Well 2</assigned_regions>
	  <functional>point</functional>
	  <time_macro>Steady State</time_macro>
	</aqueous_pressure>
	<aqueous_pressure>
	  <assigned_regions>Well 2t</assigned_regions>
	  <functional>point</functional>
	  <time_macro>Steady State</time_macro>
	</aqueous_pressure>
	<aqueous_pressure>
	  <assigned_regions>Well 2b</assigned_regions>
	  <functional>point</functional>
	  <time_macro>Steady State</time_macro>
	</aqueous_pressure>
	<aqueous_pressure>
	  <assigned_regions>Well 3</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>
