#=========================== output options ===================================
OUTPUT
#  FORMAT TECPLOT POINT
  VARIABLES
    LIQUID_PRESSURE
    GAS_SATURATION
    TEMPERATURE
    THERMODYNAMIC_STATE
    POROSITY
    TORTUOSITY
    PERMEABILITY
  /
END

#=========================== discretization ===================================
GRID
  TYPE structured
  NXYZ 1 1 5
  BOUNDS
    0.d0 0.d0 0.d0
    1.d0 1.d0 5.d0
  /
END

#=========================== regression =======================================
REGRESSION
  CELL_IDS
    1
    3
    5
  /
END

#=========================== chemistry ========================================
CHEMISTRY
  PRIMARY_SPECIES
    Tracer
  /
  OUTPUT
    FREE_ION
    ALL
  /
END

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

  TIMESTEPPER
    # end simulation without taking a time step
  #  MAX_STEPS 0
  /

  NEWTON_SOLVER
    USE_EUCLIDEAN_NORM_CONVERGENCE
  /

  LINEAR_SOLVER
    SOLVER ITERATIVE
  /
  
END

NUMERICAL_METHODS TRANSPORT

  LINEAR_SOLVER
    SOLVER ITERATIVE
  /

END

#=========================== fluid properties =================================
FLUID_PROPERTY
  DIFFUSION_COEFFICIENT 1.d-9
END

EOS WATER
  DENSITY CONSTANT 1.d3
  ENTHALPY CONSTANT 1.8890d0
  VISCOSITY CONSTANT 8.904156d-4
END

#=========================== saturation functions =============================
CHARACTERISTIC_CURVES cc1
#  TEST
  SATURATION_FUNCTION VAN_GENUCHTEN
    LIQUID_RESIDUAL_SATURATION 0.1d0
    ALPHA 1.d-4
    M 0.5d0
    MAX_CAPILLARY_PRESSURE 1.d6

  /
  PERMEABILITY_FUNCTION MUALEM
    PHASE LIQUID
    LIQUID_RESIDUAL_SATURATION 0.1d0
    M 0.5d0
  /
  PERMEABILITY_FUNCTION MUALEM_VG_GAS
    PHASE GAS
    LIQUID_RESIDUAL_SATURATION 0.1d0
    GAS_RESIDUAL_SATURATION 1.d-40
    M 0.5d0
  /
/

#=========================== times ============================================
TIME
  FINAL_TIME 0.1d0 y
  INITIAL_TIMESTEP_SIZE 1.d-6 y
  MAXIMUM_TIMESTEP_SIZE 1.d0 y
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 5.d0
    1.d20 1.d20 5.d0
  /
END

#=========================== flow conditions ==================================
FLOW_CONDITION top
  TYPE
    GAS_PRESSURE dirichlet
    GAS_SATURATION dirichlet
    TEMPERATURE dirichlet
  /
  GAS_PRESSURE 200001.d0
  GAS_SATURATION 0.99
  TEMPERATURE 25.d0
/

FLOW_CONDITION initial
  TYPE
    LIQUID_PRESSURE dirichlet
    MOLE_FRACTION dirichlet
    TEMPERATURE dirichlet
  /
  LIQUID_PRESSURE 200001.d0
  MOLE_FRACTION 1.d-7
  TEMPERATURE 25.d0
/

#=========================== transport conditions =============================
TRANSPORT_CONDITION top
  TYPE dirichlet_zero_gradient
  CONSTRAINT top
    CONCENTRATIONS
      Tracer 1.d-3  F
    /
  /
END

TRANSPORT_CONDITION initial
  TYPE zero_gradient
  CONSTRAINT middle
    CONCENTRATIONS
      Tracer 1.d-10  F
    /
  /
END

#=========================== constraints ======================================

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

BOUNDARY_CONDITION
  FLOW_CONDITION top
  TRANSPORT_CONDITION top
  REGION top
END

#=========================== stratigraphy couplers ============================
STRATA
  REGION all
  MATERIAL soil
END

