Abstract. In hybrid models, which combine hydrodynamical and transport approaches to describe different stages of heavy-ion collisions, conversion of fluid to individual particles, particlization, is a non-trivial technical problem. We describe in detail how to find the particlization hypersurface in a 3+1 dimensional model, and how to sample the particle distributions evaluated using the Cooper-Frye procedure to create an ensemble of particles as an initial state for the transport stage. We also discuss the role and magnitude of the negative contributions in the Cooper-Frye procedure.