Monte Carlo algorithms of calculating the mathematical expectation and the variance of solar reflected intensity in a given direction, developed for a Poissonian model of broken cloudiness, are extended to the case of a three-layer cloudy-aerosol atmosphere located over a Lambertian reflecting underlying surface.