  We present our final results of the charmonium spectrum in quenched QCD on
anisotropic lattices. Simulations are made with the plaquette gauge action and
a tadpole improved clover quark action employing $\xi = a_s/a_t = 3$. We
calculate the spectrum of S- and P-states and their excitation, and study the
scaling behavior of mass splittings. Comparison is made with the experiment and
previous lattice results. The issue of hyperfine splitting for different
choices of the clover coefficients obtained by Klassen is discussed.
