The charmonium states with their different binding energies and radii dissolve at different temperatures of the medium produced in relativistic heavy-ion collisions. Relative yields of charmonium and thus their survival have potential to map the properties of Quark Gluon Plasma. In this study, we estimate the combined effect of color screening, gluon-induced dissociation and recombination on charmonium production in heavy-ion collisions (Pb+Pb ions) at centre of mass energy ($sqrt{s_{rm NN}}$) = 5.02 TeV. The rate equations of dissociation and recombination are solved separately with a 2-dimensional accelerated expansion of fireball volume. To solve the recombination rate equation, we have used an approach of Bateman solution which ensures the dissociation of the recombined charmonium in the QGP medium. The modifications of charmonium states are estimated in an expanding QGP with the conditions relevant for Pb+Pb collisions at LHC.