We present the implementation of several colorsinglet final-state processes at Next-to-Next-to Leading Order (NNLO) accuracy in QCD to the publicly available parton-level Monte Carlo program MCFM. Specifically we discuss the processesDecays of the unstable bosons are fully included, resulting in a flexible fully differential Monte Carlo code. The NNLO corrections have been calculated using the non-local N -jettiness subtraction approach. Special attention is given to the numerical aspects of running MCFM for these processes at this order. We pay particular attention to the systematic uncertainties due to the power corrections induced by the N -jettiness regularization scheme and the evaluation time needed to run the hybrid openMP/MPI version of MCFM at NNLO on multiprocessor systems.