The paper presents diagnostic formulas for calculation of cloud fraction in ICM RAS GCM, together with the algorithm of calculating radiative heat influxes to atmospheric layers containing multilayer broken clouds. Model results are compared with satellite observations.