Subprogram electrode_species_bc (box)

Set Neumann boundary conditions in the actual flux computation

Subprogram m_chemistry::get_rates (fields, rates, n_cells, energy_eV)

These reactions do not take into account a variable gas_temperature

Module m_dielectric

Make sure species densities are initially zero inside the dielectric

Use special prolongation for multigrid when there are surface charges

Subprogram m_gas::gas_initialize (tree, cfg)

compute gas molecular weight from gas composition

Subprogram m_photoi_helmh::photoi_helmh_bc (box, nb, iv, coords, bc_val, bc_type)

Think about good (and efficient) boundary conditions for Helmholtz equations, see also Bourdon et al. PSST 2017. Setting a Dirichlet zero b.c. is inaccurate if the streamer gets close to that boundary, otherwise it should be quite reasonable.