A solar radiation routine has been developed for use in climate studies. It includes the absorption and scattering due to ozone, water vapor, oxygen, carbon dioxide, clouds, and aerosols. Rayleigh scattering is also included. The UV and visible region (λ < 0.69 μm) is grouped into four bands. An effective coefficient for ozone absorption and an effective cross section for Rayleigh scattering are computed for each band. In the near-infrared region (λ > 0.69 μm), the broadband parameterization is used to compute the absorption by water vapor in a clear atmosphere, and the k-distribution method is applied to compute fluxes in a scattering atmosphere. The reflectivity and transmissivity of a scattering layer are computed analytically using the delta-four-stream discrete-ordinate approximation. The two-stream adding method is then applied to compute fluxes for a composite of clear and scattering layers. Compared to the results of high spectral resolution and detailed multiple-scattering calculations, ... Abstract A solar radiation routine has been developed for use in climate studies. It includes the absorption and scattering due to ozone, water vapor, oxygen, carbon dioxide, clouds, and aerosols. Rayleigh scattering is also included. The UV and visible region (λ < 0.69 μm) is grouped into four bands. An effective coefficient for ozone absorption and an effective cross section for Rayleigh scattering are computed for each band. In the near-infrared region (λ > 0.69 μm), the broadband parameterization is used to compute the absorption by water vapor in a clear atmosphere, and the k-distribution method is applied to compute fluxes in a scattering atmosphere. The reflectivity and transmissivity of a scattering layer are computed analytically using the delta-four-stream discrete-ordinate approximation. The two-stream adding method is then applied to compute fluxes for a composite of clear and scattering layers. Compared to the results of high spectral resolution and detailed multiple-scattering calculations, ...