flow123d_version: 4.0.0a01
problem: !Coupling_Sequential
  description: Test15 - 2D unsteady flow LM
  mesh:
    mesh_file: ../00_mesh/rectangle_16x1_uns.msh
  flow_equation: !Flow_Darcy_LMH
    output_specific:
      raw_flow_output: ./raw_output.txt
    input_fields:
      - region: ALL
        storativity: 1.0
      - region: domain
        anisotropy: 1
        # default
        cross_section: 1
        # default
        water_source_density: !FieldFormula
          value: -X[0]
        init_piezo_head: 0
      - region: .left
        bc_type: dirichlet
        bc_pressure: 0
    nonlinear_solver:
      linear_solver: !Petsc
        a_tol: 1.0e-12
        r_tol: 1.0e-12
    output:
      fields:
        - pressure_p0
        - velocity_p0
    balance:
      cumulative: true
    time:
      end_time: 4.0
      max_dt: 0.4
      min_dt: 0.4
    output_stream:
      file: ./flow.pvd
      format: !vtk
        variant: ascii
      times:
        - step: 0.4
      precision: 15
  solute_equation: !Coupling_OperatorSplitting
    transport: !Solute_Advection_FV
      input_fields:
        - region: domain
          init_conc: 0
        - region: .left
          bc_conc: 10
          time: 0.0
        - region: .left
          bc_conc: 0
          time: 0.1
    output_stream:
      file: ./transport.pvd
      format: !vtk
        variant: ascii
      times:
        - step: 0.4
      precision: 15
    substances:
      - name: conc
        molar_mass: 1.0
    reaction_term: !Sorption
      solvent_density: 1.0
      substances:
        - conc
      solubility:
        - 1.0
      input_fields:
        - region: ALL
          rock_density: 1.0
          sorption_type: linear
          distribution_coefficient: 0.6
          isotherm_other: 0.4
      output:
        fields: []
    time:
      end_time: 4.0
    balance:
      cumulative: true
