Dissociation of charmonia in hot QCD is studied on a quenched anisotropic lattice using standard plaquette gauge action and O(a) improved Wilson fermion action. Simulations are carried out at temperatures in the range 0.88 TC to 2.3 TC. We calculate effective masses and Bethe–Salpeter wavefunctions for the ground state (1S, 1P) and some excited states (2S, 2P) with variational analysis. To distinguish between bound states and scattering states, we apply two methods: first, we compare effective masses for charmonium correlation functions under different spatial boundary conditions. Because the scattering states are sensitive to the boundary conditions, finite volume effects are expected when a charmonium dissociates. Second, we study if the wavefunctions show a sign of scattering states which will be characterized by a broad wavefunction. We find no clear evidence of dissociation for the charmonium states up to 2.3 TC so far.