A Monte Carlo Code to Compute Energy Fluxes in Cometary Nuclei