We can have multiple populations, each with its charge.
But the density used in Ohm and in Electrons is the total ion particledensity, which is wrong if populations have different charges.
We should allow for calculating the total charge density $\rhoc = \sum{pop} Z{pop}n{pop}$
We should changerho_ to be the charge density, since the total particle density will not be used anyway.
We can have multiple populations, each with its charge. But the density used in Ohm and in Electrons is the total ion particledensity, which is wrong if populations have different charges. We should allow for calculating the total charge density $\rhoc = \sum{pop} Z{pop}n{pop}$
We should change
rho_
to be the charge density, since the total particle density will not be used anyway.