Removed an extra kernel_igamma from the rho_dh calculation

There was an extra kernel_igamma in the calculation of rho_sh which was causing an offset to equivalent Gadget-2 runs.