A Monte-Carlo based technique and computer algorithm are presented for calculating the pulse response functions of a multilayer spherical cloudy atmosphere (optical thicknesses from 1 to 100) to isotropic pulses from a subcloud point source.
Depending on the optical thickness of the cloud layer the relative systematic rms error of the algorithm varies from 2.5% to 25%.