  We present the implementation of several color-singlet 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
processes $pp\rightarrow H$, $pp\rightarrow Z$, $pp\rightarrow W$,
$pp\rightarrow HZ$, $pp\rightarrow HW$ and $pp\rightarrow\gamma\gamma$. Decays
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
multi-processor systems.
