#=========================== regression =======================================
REGRESSION
  CELLS_PER_PROCESS 2
  CELL_IDS
    1059
    2781
    5349
    9241
    15201
  /
/

#=========================== numerical methods ================================
NUMERICAL_METHODS FLOW
  TIMESTEPPER
    MAXIMUM_NUMBER_OF_TIMESTEPS 0
  /
END

#=========================== discretization ===================================
GRID
  TYPE STRUCTURED
  NXYZ 20 20 40
  BOUNDS
     0.  0.  0.
    20. 20. 10.
  /
END

#=========================== datasets =========================================
DATASET layer1_dataset
  HDF5_DATASET_NAME top_surface
  FILENAME strata_dataset.h5
END

DATASET layer2_dataset
  HDF5_DATASET_NAME middle_surface
  FILENAME strata_dataset.h5
END

DATASET layer3_dataset
  HDF5_DATASET_NAME bottom_surface
  FILENAME strata_dataset.h5
END

DATASET layer4_dataset
  HDF5_DATASET_NAME undulating_top_surface
  FILENAME strata_dataset.h5
END

DATASET layer5_dataset
  HDF5_DATASET_NAME undulating_middle_surface
  FILENAME strata_dataset.h5
END

#=========================== material properties ==============================
MATERIAL_PROPERTY soil1
  ID 1
  POROSITY 0.25
  TORTUOSITY 1.
  CHARACTERISTIC_CURVES sf1
  PERMEABILITY
    PERM_ISO 1.d-20
  /
END

MATERIAL_PROPERTY soil2
  ID 2
  POROSITY 0.25
  TORTUOSITY 1.
  CHARACTERISTIC_CURVES sf1
  PERMEABILITY
    PERM_ISO 1.d-20
  /
END

MATERIAL_PROPERTY soil3
  ID 3
  POROSITY 0.25
  TORTUOSITY 1.
  CHARACTERISTIC_CURVES sf1
  PERMEABILITY
    PERM_ISO 1.d-20
  /
END

#=========================== characteristic curves ============================
CHARACTERISTIC_CURVES sf1
  SATURATION_FUNCTION VAN_GENUCHTEN
    M 0.5
    ALPHA  1.d-4
    LIQUID_RESIDUAL_SATURATION 0.1
  /
  PERMEABILITY_FUNCTION MUALEM_VG_LIQ
    M 0.5
    LIQUID_RESIDUAL_SATURATION 0.1
  /
END

#=========================== output options ===================================
OUTPUT
  FORMAT HDF5
  VARIABLES
    MATERIAL_ID_KLUDGE_FOR_VISIT
  /
END

#=========================== times ============================================
TIME
  FINAL_TIME 1.d0 d
  INITIAL_TIMESTEP_SIZE 1.d-1 d
  MAXIMUM_TIMESTEP_SIZE 1.d-1 d
END

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

REGION lower_half_y
  COORDINATES
    -1.d20 -1.d20 -1.d20
     1.d20  10.    1.d20
  /
END

REGION upper_half_y_lower_4m_z
  COORDINATES
    -1.d20 -1.d20 -1.d20
     1.d20  1.d20  4.
  /
END

#=========================== flow conditions ==================================
FLOW_CONDITION initial
  TYPE
    LIQUID_PRESSURE HYDROSTATIC
  /
  DATUM 0. 0. 10.
  LIQUID_PRESSURE 101325.
END

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

