Need to be able to set per particle density

[BryanRumsey]

• rho0 = numpy.zeros((numpoints), dtype=float)
• model.set_type should set rho0 for each particle of that type