problem: !SequentialCoupling
  description: Test12 - narrow radioaoctive decay chain, with transport switched enabled
    but with zero pressure gradient
  mesh:
    mesh_file: ./input/sit_decay.msh
  primary_equation: !Steady_MH
    input_fields:
      - rid: 1
        anisotropy: 0.01
      - r_set: BOUNDARY
        bc_type: dirichlet
        bc_piezo_head: 0
    output:
      output_stream: 
        file: ./decay.pvd
        format: !vtk
          variant: ascii
        name: flow_output_stream
      output_fields:
        - pressure_p0
        - pressure_p1
        - velocity_p0
    balance: true
    solver: !Petsc
      a_tol: 1.0e-07
  secondary_equation: !TransportOperatorSplitting
    input_fields:
      - rid: 1
        init_conc:
          - 0.01
          - 0.02
          - 0.03
          - 0.04
          - 0.05
          - 0.06
      - r_set: BOUNDARY
        bc_conc:
          - 0.01
          - 0.02
          - 0.03
          - 0.04
          - 0.05
          - 0.06
    output_stream: 
      file: ./transport_decay.pvd
      format: !vtk
        variant: ascii
      time_step: 0.5
      name: transport_output_stream
    substances:
      - A
      - B
      - C
      - D
      - E
      - F
    time:
      end_time: 10.0
    reaction_term: !RadioactiveDecay
      decays:
        - radionuclide: D
          half_life: 2.5
          products:
            - F
        - radionuclide: F
          half_life: 2.5
          products:
            - B
    balance:
      cumulative: true