Particle conservation in numerical models of the tokamak plasma edge