The interaction between radiation and clouds represents a source of uncertainty in numerical weather prediction (NWP), due to both intrinsic problems of one-dimensional radiation schemes and poor representation of clouds. The underlying question addressed in this study is how large the bias is of radiative heating rates in NWP models for shallow cumulus clouds and how it scales with various parameters, such as solar zenith angle, surface albedo, cloud cover and liquid water path.