Todo List

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.