  The use of variational nuclear motion programs to compute line lists of
transition frequencies and intensities is now a standard procedure. The ExoMol
project has used this technique to generate line lists for studies of hot
bodies such as the atmospheres of exoplanets and cool stars. The resulting line
list can be huge: many contain 10 billion or more transitions. This software
update considers changes made to our programs during the course of the project
to allow for such calculations. This update considers three programs: {\sc Duo}
which computed vibronic spectra for diatomics, {\sc DVR3D} which computes
rotation-vibration spectra for triatomics, and {\sc TROVE} which computes
rotation-vibration spectra for general polyatomic systems. Important updates in
functionality include the calculation of quasibound (resonance) states and
Land\'e $g$-factors by {\sc Duo} and the calculation of resonance states by
{\sc DVR3D}. Significant algorithmic improvements are reported for both {\sc
DVR3D} and {\sc TROVE}. All three programs are publically available from
ccpforge.cse.rl.ac.uk. Future developments are also considered.
