Perturbative cross sections in QCD are riddled with kinematic divergences. While the UV divergences can be treated consistently by renormalizing the couplings and fields, there are several ways to treat the IR ones originating from partons becoming soft and/or collinear. In this talk we provide an update on the status of the CoLoRFul subtraction scheme in the context of LHC processes. The latter is based on the construction of an approximation to the cross section that mimics its IR structure and has been successfully applied in the past to processes with only final-state hadrons. The main focus of the talk will be the analytic computation of the integrated subtraction terms for the production of a colour singlet final state at NNLO accuracy. The complete scheme is implemented in a Fortran code dubbed NNLOCAL, which will also briefly be discussed.