#=========================== regression testing output ========================
REGRESSION
  CELL_IDS
    1
    5
    9
  /
  VARIABLES
    LIQUID_SATURATION
    SOLUTE_CONCENTRATION
    PERMEABILITY_X
    PERMEABILITY_Z
  /
END

#=========================== numerical methods ================================
NUMERICAL_METHODS FLOW

  TIMESTEPPER
  /

  NEWTON_SOLVER
    MAX_ALLOW_LIQ_PRES_CHANGE_NI 1.d0
    MAXIMUM_NUMBER_OF_ITERATIONS 12
  /

  LINEAR_SOLVER
  /

END

#=========================== discretization ===================================
GRID
  TYPE STRUCTURED
  NXYZ 3 1 3
  BOUNDS
    0.d0 0.d0 0.d0
    3.d0 1.d0 3.d0
  /
END

#=========================== material properties ==============================
MATERIAL_PROPERTY backfill
  ID 1
  POROSITY 0.2
  ROCK_DENSITY 2.60 g/cm^3
  EXTERNAL_FILE soil_compressibility.txt
  CHARACTERISTIC_CURVES cc1
  PERMEABILITY
#    PERM_HORIZONTAL 4.7615d-11
    PERM_HORIZONTAL 1.d-11
    VERTICAL_ANISOTROPY_RATIO 0.1d0
  /
END

MATERIAL_PROPERTY Hf1
  ID 2
  POROSITY 0.25d0
  ROCK_DENSITY 2.60 g/cm^3
  EXTERNAL_FILE soil_compressibility.txt
  CHARACTERISTIC_CURVES cc2
  PERMEABILITY
#    PERM_HORIZONTAL 4.7615d-11
    PERM_HORIZONTAL 1.d-12
    VERTICAL_ANISOTROPY_RATIO 0.1d0
  /
END

MATERIAL_PROPERTY Hf2
  ID 3
  POROSITY 0.35d0
  ROCK_DENSITY 2.71 g/cm^3
  EXTERNAL_FILE soil_compressibility.txt
  CHARACTERISTIC_CURVES cc3
  PERMEABILITY
#    PERM_HORIZONTAL 6.3160d-12
    PERM_HORIZONTAL 5.d-11
    VERTICAL_ANISOTROPY_RATIO 0.1d0
  /
END

#=========================== characteristic curves ============================
CHARACTERISTIC_CURVES cc1
  SATURATION_FUNCTION VAN_GENUCHTEN
    ALPHA 8.d-4
    M 0.25d0
    LIQUID_RESIDUAL_SATURATION 0.02d0
    MAX_CAPILLARY_PRESSURE 1.d8
  /
  PERMEABILITY_FUNCTION MUALEM_VG_LIQ
    SMOOTH
    M 0.25d0
    LIQUID_RESIDUAL_SATURATION 0.02d0
  /
END

CHARACTERISTIC_CURVES cc2
  SATURATION_FUNCTION VAN_GENUCHTEN
    ALPHA 6.d-4 # 0.06419 1/cm
    M 0.4d0
    LIQUID_RESIDUAL_SATURATION 0.08d0
    MAX_CAPILLARY_PRESSURE 1.d8
  /
  PERMEABILITY_FUNCTION MUALEM_VG_LIQ
    SMOOTH
    M 0.4d0
    LIQUID_RESIDUAL_SATURATION 0.08d0
  /
END

CHARACTERISTIC_CURVES cc3
  SATURATION_FUNCTION VAN_GENUCHTEN
    ALPHA 7.d-5
    M 0.5d0
    LIQUID_RESIDUAL_SATURATION 0.13d0
    MAX_CAPILLARY_PRESSURE 1.d8
  /
  PERMEABILITY_FUNCTION MUALEM_VG_LIQ
    SMOOTH
    M 0.5d0
    LIQUID_RESIDUAL_SATURATION 0.13d0
  /
END

#=========================== regions ==========================================
REGION all
  COORDINATES
    -1.d20 -1.d20 -1.d20
    1.d20 1.d20 1.d20
  /
END

REGION top
  FACE TOP
  COORDINATES
    -1.d20 -1.d20 3.d0
    1.d20 1.d20 3.d0
  /
END

REGION bottom
  FACE BOTTOM
  COORDINATES
    -1.d20 -1.d20 0.d0
    1.d20 1.d20 0.d0
  /
END

REGION backfill
  COORDINATES
    -1.d20 -1.d20 2.d0
    1.d20 1.d20 3.d0
  /
END

REGION Hf1
  COORDINATES
    -1.d20 -1.d20 1.d0
    1.d20 1.d20 2.d0
  /
END

REGION Hf2
  COORDINATES
    -1.d20 -1.d20 0.d0
    1.d20 1.d20 1.d0
  /
END

#=========================== flow conditions ==================================
FLOW_CONDITION initial
  TYPE
    LIQUID_PRESSURE HYDROSTATIC
    CONCENTRATION DIRICHLET
  /
  DATUM 0.d0 0.d0 0.1d0
  LIQUID_PRESSURE 101325.d0
  CONCENTRATION 1.d-10
END

FLOW_CONDITION recharge
  TYPE
    LIQUID_FLUX NEUMANN
    CONCENTRATION DIRICHLET
  /
  LIQUID_FLUX 4.d0 mm/y
  CONCENTRATION 1.d-10
END

#=========================== condition couplers ===============================
# initial condition
INITIAL_CONDITION
  FLOW_CONDITION initial
  REGION all
END

BOUNDARY_CONDITION bottom
  FLOW_CONDITION initial
  REGION bottom
END

BOUNDARY_CONDITION land_surface
  FLOW_CONDITION recharge
  REGION top
END

#=========================== stratigraphy couplers ============================
STRATA
  REGION backfill
  MATERIAL backfill
END

STRATA
  REGION Hf1
  MATERIAL Hf1
END

STRATA
  REGION Hf2
  MATERIAL Hf2
END

