A method for calculations of rotational-vibrational states of triatomic molecules up to dissociation threshold (and scattering resonances above it) is devised, that combines hyper-spherical coordinates, sequential diagonalization-truncation procedure, optimized grid DVR, and complex absorbing potential. Efficiency and accuracy of the method and new code are tested by computing the spectrum of ozone up to dissociation threshold, using two different potential energy surfaces. In both cases good agreement with results of previous studies is obtained for the lower energy states localized in the deep (∼10 000 cm −1 ) covalent well. Upper part of the bound state spectrum, within 600 cm −1 below dissociation threshold, is also computed and is analyzed in detail. It is found that long progressions of symmetric-stretching and bending states (up to 8 and 11 quanta, respectively) survive up to dissociation threshold and even above it, whereas excitations of the asymmetric-stretching overtones couple to the local vibration modes, making assignments difficult. Within 140 cm −1 below dissociation threshold, large-amplitude vibrational states of a floppy complex O· · · O 2 are formed over the shallow van der Waals plateau. These are assigned using two local modes: the rocking-motion and the dissociative-motion progressions, up to 6 quanta in each, both with frequency ∼20 cm −1 . Many of these plateau states are mixed with states of the covalent well. Interestingly, excitation of the rocking-motion helps keeping these states localized within the plateau region, by raising the effective barrier. Published by AIP Publishing. [http://dx