This algorithm derives estimates of broadband fluxes and heating rates consistent with liquid and ice water content estimates from the CloudSat Profiling Radar (CPR), CALIPSO, and MODIS. For each radar profile, a broadband radiative transfer model is used to calculate upwelling and downwelling longwave and shortwave fluxes at each CPR range gate from the surface to the lower stratosphere. Profiles of cloud ice and liquid water content and cloud particle effective radii are defined based on the CloudSat 2B-LWC and 2B-IWC products while precipitation properties are defined using the CloudSat 2C-PRECIP-COLUMN dataset. The properties of thin or low clouds that are not detected by CloudSat are assigned based on CALIOP backscatter, CALIPSO’s CAL_LID_L2_05kmCLay, and the CloudSat 2B-TAU product while aerosol location and optical properties are specified using CALIPSO’s CAL_LID_L2_05kmALay product. Ancillary atmospheric state variables are interpolated from ECMWF analyses and surface albedos are assigned based on seasonally-varying maps of surface reflectance properties in combination with daily snow and sea ice cover maps from passive microwave instruments. Two additional sets of flux calculations are performed with all clouds and all aerosols removed, respectively. Corresponding profiles of atmospheric heating are inferred from the vertical derivative of these fluxes.