SIMULATION
  SIMULATION_TYPE SUBSURFACE
  PROCESS_MODELS
    SUBSURFACE_FLOW flow
      MODE RICHARDS
    /
  /
END

SUBSURFACE

GRID
  TYPE structured
  NXYZ 10 1 1
  DXYZ
   2.0d0
   1.0d0
   1.0d0
  END
END

REGION all
  COORDINATES
    0.0d0 0.0d0 0.0d0
    20.0d0 1.0d0 1.0d0
  /
END

REGION left_end
  FACE WEST
  COORDINATES
    0.d0 0.d0 0.d0
    0.d0 1.0d0 1.0d0
  /
END

REGION right_end
  FACE EAST
  COORDINATES
    20.0d0 0.d0 0.d0
    20.0d0 1.0d0 1.0d0
  /
END

MATERIAL_PROPERTY beam
  ID 1
  POROSITY 0.25
  TORTUOSITY 1.d0
  ROCK_DENSITY 2.5d3 kg/m^3
  CHARACTERISTIC_CURVES cc1
  PERMEABILITY
    PERM_X 1.d-14
    PERM_Y 1.d-14
    PERM_Z 1.d-14
  /
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

EOS WATER
  #ref. density [kg/m3], ref. pressure [Pa], water compressibility [1/Pa?]
  DENSITY EXPONENTIAL 1000. 0. 1.d-8
  VISCOSITY CONSTANT 1.728d-3 Pa-s
END

FLUID_PROPERTY
  DIFFUSION_COEFFICIENT 1.d-9
END

STRATA
  REGION all
  MATERIAL beam
END

TIME
  FINAL_TIME 0.75 day
  INITIAL_TIMESTEP_SIZE 1.d-4 day
  MAXIMUM_TIMESTEP_SIZE 0.0001 day at 0.d0 day
END

OUTPUT
  SNAPSHOT_FILE
    TIMES day 0.10 0.25 0.50
    #NO_PRINT_INITIAL
    FORMAT HDF5
  /
END

FLOW_CONDITION initial
  TYPE
    LIQUID_PRESSURE dirichlet
  /
  LIQUID_PRESSURE 0.25 MPa
END

FLOW_CONDITION left_end
  TYPE
    LIQUID_PRESSURE dirichlet
  /
  LIQUID_PRESSURE LIST
    # p = p_b*t + offset; p_b=2e6 Pa/day; offset = 0.25 MPa
    TIME_UNITS day
    DATA_UNITS Pa
    INTERPOLATION LINEAR
    #time  #pressure
    0.00d0 0.25d6
    0.25d0 0.75d6
    0.50d0 1.25d6
    1.00d0 2.25d6
  /
END

FLOW_CONDITION right_end
  TYPE
    LIQUID_PRESSURE dirichlet
  /
  LIQUID_PRESSURE LIST
    # p = p_b*t + offset; p_b=2e6 Pa/day; offset = 0.25 MPa
    TIME_UNITS day
    DATA_UNITS Pa
    INTERPOLATION LINEAR
    #time  #pressure
    0.00d0 0.25d6
    0.25d0 0.75d6
    0.50d0 1.25d6
    1.00d0 2.25d6
  /
END

INITIAL_CONDITION initial
  REGION all
  FLOW_CONDITION initial
END

BOUNDARY_CONDITION left_end
  REGION left_end
  FLOW_CONDITION left_end
END

BOUNDARY_CONDITION right_end
  FLOW_CONDITION right_end
  REGION right_end
END

END_SUBSURFACE