We present the computational methodology that allows rigorous and efficient nine-dimensional (9D) quantum calculations of the intermolecular vibrational states of noncovalently bound trimers of diatomic molecules, with the monomers treated as rigid. The full 9D vibrational Hamiltonian of the trimer is partitioned into a 3D "frame" (or stretching) Hamiltonian and a 6D "bend" Hamiltonian. These two Hamiltonians are diagonalized separately, and a certain number of their lowest-energy eigenstates is included in the final 9D product contracted basis in which the full 9D intermolecular vibrational Hamiltonian is diagonalized. This methodology is applied to the 9D calculations of the intermolecular vibrational levels of (HF)3, a prototypical hydrogen-bonded trimer, on the rigid-monomer version of an ab initio calculated potential energy surface (PES). They are the first to include fully the stretch-bend coupling present in the trimer. The frequencies of all bending fundamentals considered from the present 9D calculations are about 10% lower than those from the earlier quantum 6D calculations that considered only the bending modes of the HF trimer. This means that the stretch-bend coupling is strong, and it is imperative to include it in any accurate treatment of the (HF)3 vibrations aiming to assess the accuracy of the PES employed. Moreover, the 9D results are in better agreement with the limited available spectroscopic data that those from the 6D calculations. In addition, the 9D results show sensitivity to the value of the HF bond length, equilibrium or vibrationally averaged, used in the calculations. The implication is that full-dimensional 12D quantum calculations will be required to obtain definitive vibrational excitation energies for a given PES. Our study also demonstrates that the nonadditive three-body interactions are very significant in (HF)3 and have to be included in order to obtain accurate intermolecular vibrational energy levels of the trimer.