problem: !SequentialCoupling
  description: Test2 - Transport 1D-2D (convection, dual porosity, sorption)
  mesh:
    mesh_file: ${INPUT}/compatible.msh
    sets:
      - name: subset
        region_labels:
          - 2d
          - 1d
      - name: .sides
        region_labels:
          - .2d
          - .1d
  primary_equation: !Steady_MH
    input_fields:
      - region: 1d
        anisotropy: 10
        cross_section: 1
        sigma: 0.05
      - region: 2d
        anisotropy: 1
        cross_section: 1
      - r_set: .sides
        bc_type: dirichlet
        bc_pressure: !FieldFormula
          value: x
    balance: true
    output:
      output_stream:
        file: ./test2.pvd
        format: !vtk
          variant: ascii
      output_fields:
        - pressure_p0
        - pressure_p1
        - velocity_p0
        - region_id
    solver: !Petsc
      a_tol: 1.0e-12
      r_tol: 1.0e-12
  secondary_equation: !TransportOperatorSplitting
    input_fields:
      - r_set: .sides
        bc_conc: !FieldFormula
          value: x/12.75
      - r_set: subset
        init_conc: 0
        porosity: 0.25
      - r_set: ALL
        init_conc: !FieldElementwise
          gmsh_file: ${INPUT}/test18_init.msh
          field_name: concentration
        diff_m: 0.1
        dg_penalty: 100
    output_stream:
      file: ./transport.pvd
      format: !vtk
        variant: ascii
      time_step: 0.5
    output_fields:
      - conc
      - region_id
    substances:
      - A
      - B
    time:
      end_time: 5.0
    balance:
      cumulative: true
    reaction_term: !DualPorosity
      input_fields:
        - r_set: subset
          diffusion_rate_immobile:
            - 0.01
            - 0.01
          porosity_immobile: 0.25
          init_conc_immobile:
            - 0.0
            - 0.0
      output_fields: []
      reaction_mobile: !SorptionMobile
        solvent_density: 1.0
        substances:
          - A
          - B
        solubility:
          - 1.0
          - 1.0
        input_fields:
          - r_set: subset
            rock_density: 1.0
            sorption_type:
              - linear
              - freundlich
            isotherm_mult: 0.02
            isotherm_other:
              - 0
              - 0.5
        output_fields: []
      reaction_immobile: !SorptionImmobile
        solvent_density: 1.0
        substances:
          - A
          - B
        solubility:
          - 1.0
          - 1.0
        input_fields:
          - r_set: subset
            rock_density: 1.0
            sorption_type:
              - linear
              - freundlich
            isotherm_mult: 0.02
            isotherm_other:
              - 0
              - 0.5
        output_fields: []