Fixed error in setting the mass of the particles and in reading in the radii. Changed the bulk density initialization to be a function of the number of particles that I want.