  In my talk I gave a status update on the extension of the CoLoRFulNNLO
subtraction method for computing QCD jet cross sections with hadrons in the
initial state. The scheme has been fully worked out previously for
electron-positron collisions and recently important steps have been made
towards generalizing it to be able to deliver corrections of the same order for
LHC processes as well. In particular, the important bottleneck of regularizing
multiple real emissions has been addressed. We demonstrate the numerical
stability of the CoLoRFulNNLO method by computing the doubly real contribution
for Higgs-boson production in gluon-gluon fusion and for W production.
