a b s t r a c tStimulated by recent THz measurements of the methanol spectrum in one of our laboratories, undertaken in support of NASA programs related to the Herschel Space Observatory (HSO) and the Atacama Large Millimeter Array (ALMA), we have carried out a global analysis of available microwave and high-resolution infrared data for the first three torsional states (m t = 0, 1, 2), and for J values up to 30. This global fit of approximately 5600 frequency measurements and 19 000 Fourier transform far infrared (FTFIR) wavenumber measurements to 119 parameters reaches the estimated experimental measurement accuracy for the FTFIR transitions, and about twice the estimated experimental measurement accuracy for the microwave, submillimeter-wave, and terahertz transitions. The present fit is essentially a continuation of our earlier work, but we have greatly expanded our previous data set and have added a large number of new torsion-rotation interaction terms to the Hamiltonian in our previously used computer program. The results, together with a number of calculated (but unmeasured) transitions, including their line strength, estimated uncertainty, and lower state energy, are made available in the supplementary material as a database formatted to be useful for astronomical searches. Some discussion of several open spectroscopic problems, e.g., (i) an improved notation for the numerous parameters in the torsion-rotation Hamiltonian, (ii) possible causes of the failure to fit frequency measurements to the estimated measurement uncertainty, and (iii) pitfalls to be avoided when intercomparing apparently identical parameters from the internal axis method and the rho axis method are also given.