Work on slice_constant_charge; still with wrong particle numbers.