CO rotational diagrams of all objects. We plot the decadic logarithm of the column densities N_u divided by the statistical weight g_u, as function of the upper level energy of the transition E_u in Kelvin. SPIRE data are fitted with black, PACS data are fitted with blue and red. As the plots contain a logarithmic scale, the errors are actually not symmetric, but we use it here as an approximation.