Plotting function bug in plot_number_densities_1D

In `plot.py`, `plot_number_densities_1D` still assumes that we have two dummy species (and will therefore not plot the first and last species). This has not been noticed since we removed the dummy species in #212.