- `c0 = numpy.zeros((numpoints), dtype=float)` - set_type should set c0 for each particle of that type