A new operational ice dynamics module has been implemented within the framework of the Canadian Community Ice-Ocean Model. The formulation is based on Hibler's viscous plastic rheology, a Particle-In-Cell (PIC) approach to model ice advection, and a Zhang-Hibler scheme to solve the momentum equations. The PIC approach reduces numerical diffusion and improves the accuracy of predicting ice edge locations. Discretizing the ice cover into a large number of particles additionally simplifies modelling thickness distribution in great detail, without resorting to complex assumptions concerning a thickness distribution function. The Zhang-Hibler scheme leads to relatively high computational efficiency, and guarantees that the plastic yield conditions are satisfied. The program is tested using idealized cases, which demonstrate the model's ability to predict ice edge locations, zones of high pressure, and ridging.