# COMMENT_Stop_time:"1275",
problem: !SequentialCoupling
  description: Steady flow + transport with source
  mesh:
    mesh_file: ./input/test16.msh
  primary_equation: !Steady_MH
    n_schurs: 2
    input_fields:
      - r_set: BOUNDARY
        bc_type: dirichlet
        bc_pressure: !FieldFormula
          value: y
      - r_set: BULK
        cross_section: 1
        conductivity: 1
    output:
      output_stream: 
        file: ./flow_test16.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-12
      r_tol: 1.0e-12
  secondary_equation: !SoluteTransport_DG
    # dual_porosity:false,
    input_fields:
      - r_set: BOUNDARY
        bc_conc: 0
      - r_set: BULK
        init_conc: 0
        porosity: 1
        diff_m: 1.0e-09
        disp_l: 0.01
        disp_t: 0.01
        dg_penalty: 10.0
        sources_density: !FieldElementwise
          gmsh_file: input/test16_sources.msh
          field_name: density
        sources_sigma: !FieldElementwise
          gmsh_file: input/test16_sources.msh
          field_name: sigma
        sources_conc: !FieldElementwise
          gmsh_file: input/test16_sources.msh
          field_name: concentration
    dg_variant: symmetric
    output_stream: 
      file: ./transport_test16.pvd
      format: !vtk
        variant: ascii
      time_step: 0.8
      name: transport_output_stream
    sorption_enable: false
    sources_file: ./input/test16.tso
    substances:
      - conc
    time:
      end_time: 5.7
      max_dt: 1.0
      min_dt: 1.0
    solver: !Petsc
      a_tol: 1.0e-12
      r_tol: 1.0e-12
    balance:
      cumulative: true