dynamic update of particle numbers working
need to parallelize