<ParameterList name="Main">

  <Parameter name="Amanzi Input Format Version" type="string" value="1.1.0"/>

  <!--
     Getting started on a simple 2D steady-state saturated test
  -->
  <ParameterList name="General Description">
    <Parameter name="Model ID" type="string" value="Single-Phase Saturated"/>
    <Parameter name="Author" type="string" value="Elizabeth Keating"/>
  </ParameterList>

  <!--
     Execution Control
  -->
  <ParameterList name="Execution Control">
    <Parameter name="Flow Model" type="string" value="Single Phase"/>
    <Parameter name="Transport Model" type="string" value="Off"/>
    <Parameter name="Chemistry Model" type="string" value="Off"/>

    <ParameterList name="Time Integration Mode">
        <ParameterList name="Steady">
          <Parameter name="Start" type="double" value="0.0"/>
          <Parameter name="End" type="double" value="0.0"/>
          <Parameter name="Initial Time Step" type="double" value="1000"/>
        </ParameterList>
    </ParameterList>

    <Parameter name="Verbosity" type="string" value="High"/>

    <ParameterList name="Numerical Control Parameters">
      <ParameterList name="Unstructured Algorithm">
        <Parameter name="linear solver tolerance" type="double" value="1e-18"/>
        <Parameter name="linear solver maximum iterations" type="int" value="10"/>
      </ParameterList>
    </ParameterList>
  
  </ParameterList>
      
  <ParameterList name="Domain">
    <Parameter name="Spatial Dimension" type="int" value="3"/>
  </ParameterList>

  <!--
  
    Mesh:  three-dimensional box 100m x 50m x 50m 

  -->
  <ParameterList name="Mesh">
    <ParameterList name="Unstructured">
      <ParameterList name="Expert">
	<Parameter name="Framework" type="string" value="MSTK"/>
      </ParameterList>
      <ParameterList name="Generate Mesh">
        <ParameterList name="Uniform Structured">
          <Parameter name="Number of Cells" type="Array(int)" value="{100, 5, 50}"/>
          <Parameter name="Domain Low Coordinate" type="Array(double)" value="{0.0, 0.0,  0.0}" />
          <Parameter name="Domain High Coordinate" type="Array(double)" value="{100.0, 50.0, 50.0}" />
        </ParameterList>
      </ParameterList>
    </ParameterList>
  </ParameterList>

  <!--

    Regions: 

  -->
  <ParameterList name="Regions">
    <!--
        Domains:
    -->
    <ParameterList name="Entire Domain">
      <ParameterList name="Region: Box">
        <Parameter name="Low Coordinate" type="Array(double)" value="{0.0,0.0,0.0}"/>
        <Parameter name="High Coordinate" type="Array(double)" value="{100.0,50.0,50.0}"/>
      </ParameterList>
    </ParameterList>
    <!--
        Boundaries
    -->
    <ParameterList name="Upstream">
      <ParameterList name="Region: Box">
        <Parameter name="Low Coordinate" type="Array(double)" value="{0.0,0.0,0.0}"/>
        <Parameter name="High Coordinate" type="Array(double)" value="{0.0,50.0,50.0}"/>
      </ParameterList>
    </ParameterList>
    <ParameterList name="Downstream">
      <ParameterList name="Region: Box">
        <Parameter name="Low Coordinate" type="Array(double)" value="{100.0,0.0,0.0}"/>
        <Parameter name="High Coordinate" type="Array(double)" value="{100.0,50.0,50.0}"/>
      </ParameterList>
    </ParameterList>
    <!-- 
	 Wells:
    -->
    <ParameterList name="Well 1">
      <ParameterList name="Region: Point">
	<Parameter name="Coordinate" type="Array(double)" value="{.5, 25,49.5}"/>
      </ParameterList>
    </ParameterList>
    <ParameterList name="Well 2">
      <ParameterList name="Region: Point">
	<Parameter name="Coordinate" type="Array(double)" value="{50.5, 25, 25.5}"/>
      </ParameterList>
    </ParameterList>
    <ParameterList name="Well 2t">
      <ParameterList name="Region: Point">
	<Parameter name="Coordinate" type="Array(double)" value="{50.5, 25, 49.5}"/>
      </ParameterList>
    </ParameterList>
    <ParameterList name="Well 2b">
      <ParameterList name="Region: Point">
	<Parameter name="Coordinate" type="Array(double)" value="{50.5, 25, 0.5}"/>
      </ParameterList>
    </ParameterList>
    <ParameterList name="Well 3">
      <ParameterList name="Region: Point">
	<Parameter name="Coordinate" type="Array(double)" value="{99.5, 25, 49.5}"/>
      </ParameterList>
    </ParameterList>

  </ParameterList>


  <!--

    Materials:

  -->
  <ParameterList name="Material Properties">
    <ParameterList name="Soil">
      <Parameter name="Assigned Regions" type="Array(string)" value="{Entire Domain}"/>
      <ParameterList name="Porosity: Uniform">
        <Parameter name="Value" type="double" value="0.25"/>
      </ParameterList>
      <ParameterList name="Intrinsic Permeability: Uniform">
        <Parameter name="Value" type="double" value="1.E-12"/>
      </ParameterList>
    </ParameterList>
  </ParameterList>


  <!--

    Phases:

  -->
  <ParameterList name="Phase Definitions">
    <ParameterList name="Aqueous">
      <ParameterList name="Phase Properties">
        <ParameterList name="Viscosity: Uniform">
          <Parameter name="Viscosity" type="double" value="1.002e-3 "/>
        </ParameterList>
        <ParameterList name="Density: Uniform">
          <Parameter name="Density" type="double" value="998.2 "/>
        </ParameterList>
      </ParameterList>

    </ParameterList>
  </ParameterList>

  <!--

    Initial Conditions:

  -->
  <ParameterList name="Initial Conditions">

    <ParameterList name="Initial Condition">

      <Parameter name="Assigned Regions" type="Array(string)" value="{Entire Domain}"/>
      
      <ParameterList name="IC: Uniform Pressure">
        <Parameter name="Phase" type="string" value="Aqueous"/>
        <Parameter name="Value" type="double" value="5001325.0"/>
      </ParameterList>

    </ParameterList>

  </ParameterList>

  <!--

    Boundary Conditions

  -->
  <ParameterList name="Boundary Conditions">

    <ParameterList name="Upstream BC">
      <Parameter name="Assigned Regions" type="Array(string)" value="{Upstream}"/>
      <ParameterList name="BC: Flux">
        <Parameter name="Times" type="Array(double)" value="{0.0, 1.0}"/>
        <Parameter name="Time Functions" type="Array(string)" value="{Constant}"/>
        <Parameter name="Inward Mass Flux" type="Array(double)" value="{1.95e-2, 1.95e-2}"/>
      </ParameterList>
    </ParameterList>

    <ParameterList name="Downstream BC">
      <Parameter name="Assigned Regions" type="Array(string)" value="{Downstream}"/>
      <ParameterList name="BC: Hydrostatic">
        <Parameter name="Times" type="Array(double)" value="{0.0, 1.0}"/>
        <Parameter name="Time Functions" type="Array(string)" value="{Constant}"/>
        <Parameter name="Water Table Height" type="Array(double)" value="{120, 120}"/>
      </ParameterList>
    </ParameterList>

  </ParameterList>

  <ParameterList name="Output">

    <ParameterList name="Time Macros">
      <ParameterList name="Steady State">
	<Parameter name="Values" type="Array(double)" value="{0.}"/>
      </ParameterList>
    </ParameterList>

    <ParameterList name="Observation Data">

      <Parameter name="Observation Output Filename" type="string" value="observation.out"/>
      <ParameterList name="Pressure 1">
        <Parameter name="Region" type="string" value="Well 1"/>
        <Parameter name="Functional" type="string" value="Observation Data: Point"/>
        <Parameter name="Variable" type="string" value="Aqueous pressure"/>
        <Parameter name="Time Macro" type="string" value="Steady State"/>
      </ParameterList>
      <ParameterList name="Pressure 2">
        <Parameter name="Region" type="string" value="Well 2"/>
        <Parameter name="Functional" type="string" value="Observation Data: Point"/>
        <Parameter name="Variable" type="string" value="Aqueous pressure"/>
        <Parameter name="Time Macro" type="string" value="Steady State"/>
      </ParameterList>
      <ParameterList name="Pressure 3">
        <Parameter name="Region" type="string" value="Well 2t"/>
        <Parameter name="Functional" type="string" value="Observation Data: Point"/>
        <Parameter name="Variable" type="string" value="Aqueous pressure"/>
        <Parameter name="Time Macro" type="string" value="Steady State"/>
      </ParameterList>
      <ParameterList name="Head 3">
        <Parameter name="Region" type="string" value="Well 2t"/>
        <Parameter name="Functional" type="string" value="Observation Data: Point"/>
        <Parameter name="Variable" type="string" value="Hydraulic Head"/>
        <Parameter name="Time Macro" type="string" value="Steady State"/>
      </ParameterList>
      <ParameterList name="Pressure 4">
        <Parameter name="Region" type="string" value="Well 2b"/>
        <Parameter name="Functional" type="string" value="Observation Data: Point"/>
        <Parameter name="Variable" type="string" value="Aqueous pressure"/>
        <Parameter name="Time Macro" type="string" value="Steady State"/>
      </ParameterList>
      <ParameterList name="Pressure 5">
        <Parameter name="Region" type="string" value="Well 3"/>
        <Parameter name="Functional" type="string" value="Observation Data: Point"/>
        <Parameter name="Variable" type="string" value="Aqueous pressure"/>
        <Parameter name="Time Macro" type="string" value="Steady State"/>
      </ParameterList>

    </ParameterList>

    <ParameterList name="Visualization Data">
      <Parameter name="File Name Base" type="string" value="steady-flow"/>
      <Parameter name="File Name Digits" type="string" value="5"/>
      <Parameter name="Time Macro" type="Array(string)" value="{Steady State}"/>
    </ParameterList>

  </ParameterList>

</ParameterList>

