  Combining the hydrodynamic equations for the QGP evolution and the transport
equation for the primordially produced $J/\psi$ in the QGP, we investigate the
$J/\psi$ transverse momentum distribution as well as its suppression in the
$\sqrt{s}$=17.3A GeV Pb-Pb collisions at SPS and $\sqrt{s}$=200A GeV Au-Au
collisions at RHIC. The two sets of equations are connected by the $J/\psi$
anomalous suppression induced by its inelastic scattering with gluons in the
QGP. The calculated centrality dependence of $J/\psi$ suppression and average
transverse momentum square agree well with the SPS data. From the comparison
with the coalescence model where charm quark is fully thermalized, our
calculated elliptical flow of the primordially produced $J/\psi$ is much
smaller. This may be helpful to differentiate the $J/\psi$ production
mechanisms in relativistic heavy ion collisions.
