A numerical model for calculation of the incoherent component of the field scattered from random rough surfaces is described. This model is based on the point scattering approach, where the mean scatterer amplitudes are calculated from deterministic models. These amplitudes are then scaled by a complex circular Gaussian random variable to simulate scattering from a surface with minimal coherence length. The resulting simulated fields are shown to agree with theory for the mean field, mean square field, statistical distribution, and the spatial coherence length.