SIMULATION
  SIMULATION_TYPE SUBSURFACE
  PROCESS_MODELS
    SUBSURFACE_FLOW flow
      MODE TH
    /
  /
END

SUBSURFACE

GRID
  TYPE structured
  NXYZ 10 1 1
  DXYZ
   10.0d0
   10.0d0
   10.0d0
  END
END

REGION all
  COORDINATES
    0.0d0 0.0d0 0.0d0
    100.0d0 10.0d0 10.0d0
  /
END

REGION left
  COORDINATES
    0.d0   0.d0   0.d0
    40.0d0 10.0d0 10.0d0
  /
END

REGION right
  COORDINATES
    40.0d0 0.0d0 0.0d0
    100.0d0 10.0d0 10.0d0
  /
END

REGION left_face
  FACE WEST
  COORDINATES
    0.d0 0.d0 0.d0
    0.d0 10.0d0 10.0d0
  /
END

REGION right_face
  FACE EAST
  COORDINATES
    100.0d0 0.d0 0.d0
    100.0d0 10.0d0 10.0d0
  /
END

MATERIAL_PROPERTY beam_left
  ID 1
  POROSITY 1.d-5
  TORTUOSITY 1.d0
  ROCK_DENSITY 2.8E3
  SPECIFIC_HEAT 1.d-3
  THERMAL_CONDUCTIVITY_DRY 100 W/m-C
  THERMAL_CONDUCTIVITY_WET 100 W/m-C
  CHARACTERISTIC_CURVES cc1
  PERMEABILITY
    PERM_X 1.d-20
    PERM_Y 1.d-20
    PERM_Z 1.d-20
  /
END

MATERIAL_PROPERTY beam_right
  ID 2
  POROSITY 1.d-5
  TORTUOSITY 1.d0
  ROCK_DENSITY 2.8E3
  SPECIFIC_HEAT 1.d-3
  THERMAL_CONDUCTIVITY_DRY 300 W/m-C
  THERMAL_CONDUCTIVITY_WET 300 W/m-C
  CHARACTERISTIC_CURVES cc1
  PERMEABILITY
    PERM_X 1.d-20
    PERM_Y 1.d-20
    PERM_Z 1.d-20
  /
END

CHARACTERISTIC_CURVES cc1
  SATURATION_FUNCTION VAN_GENUCHTEN
    LIQUID_RESIDUAL_SATURATION 0.5d-1
    M 0.75
    ALPHA 1.d-3
  /
  PERMEABILITY_FUNCTION MUALEM
    LIQUID_RESIDUAL_SATURATION 0.1d0
    M 0.5d0
  /
  PERMEABILITY_FUNCTION MUALEM_VG_GAS
    LIQUID_RESIDUAL_SATURATION 0.1d0
    GAS_RESIDUAL_SATURATION 0.1d0
    M 0.5d0
  /
END

STRATA
  REGION left
  MATERIAL beam_left
END

STRATA
  REGION right
  MATERIAL beam_right
END

TIME
  FINAL_TIME 100 yr
  INITIAL_TIMESTEP_SIZE 1.d-4 day
  MAXIMUM_TIMESTEP_SIZE 2.d0 yr at 0.d0 yr
END

OUTPUT
  SNAPSHOT_FILE
    NO_PRINT_INITIAL
    #NO_PRINT_FINAL
    FORMAT HDF5
  /
END

FLOW_CONDITION initial
  TYPE
    TEMPERATURE dirichlet
    LIQUID_PRESSURE dirichlet
  /
  TEMPERATURE 1.0d0 C
  LIQUID_PRESSURE 101325 Pa
END

FLOW_CONDITION left_face
  TYPE
    TEMPERATURE dirichlet
    LIQUID_PRESSURE dirichlet
  /
  TEMPERATURE 1.d0 C
  LIQUID_PRESSURE 101325 Pa
END

FLOW_CONDITION right_face
  TYPE
    ENERGY_FLUX neumann
    LIQUID_PRESSURE dirichlet
  /
  ENERGY_FLUX +1.5d0 W/m^2
  LIQUID_PRESSURE 101325 Pa
END

INITIAL_CONDITION initial_left
  REGION left
  FLOW_CONDITION initial
END
INITIAL_CONDITION initial_right
  REGION right
  FLOW_CONDITION initial
END

BOUNDARY_CONDITION left_face
  REGION left_face
  FLOW_CONDITION left_face
END

BOUNDARY_CONDITION right_face
  REGION right_face
  FLOW_CONDITION right_face
END

END_SUBSURFACE