The absorption spectra of Tm3+ in LiTmF4 have been measured at 2, 10, 30, and 50 K in the spectral interval 4000-25 000 cm-1. The energy levels of the ground-state configuration were calculated by diagonalizing the Hamiltonian of the electron-electron interaction, the spin-orbit coupling, and the crystal field in a basis of the whole configuration. The electrostatic parameter F2, the spin-orbit parameter ζ, and the crystal-field parameters Bkq were varied to obtain the best agreement with the experimentally observed levels. As the model does not account for configuration mixing and minor magnetic effects, it was necessary after optimizing F2 and ζ to match the centers of gravity for the multiplets before the final adjustment of the B parameters. When this was done, the standard deviation was lowered from 170 to 12 cm-1. The B parameters obtained for Tm3+ have been compared to those of Tb3+, Ho3+, and Er3+ in LiLnF4, and they follow a common trend. The intensities of the transitions from the ground state were calculated in the Judd-Ofelt scheme, fitting six complex intensity parameters A(kqλ) for best agreement with the experimentally observed intensities. The model was only able to give a rough estimate of the transition probabilities. The obtained relative standard deviation was 1.1. Contrary to what was found in the case of the energy calculations, it was important for the intensity calculations that the B parameters were allowed to take complex values. The imaginary part of the A parameters was not important to the intensities.