A scheme is described that provides an integrated description of turbulent transport in free convective boundary layers with shallow cumulus. The scheme uses a mass-flux formulation, as is commonly found in cumulus schemes, and a 1.5-order closure, involving turbulent kinetic energy and eddy diffusivity. Both components are active in both the subcloud and cloud layers. The scheme is called “mass flux–diffusion.” In the subcloud layer, the mass-flux component provides nonlocal transport. The scheme combines elements from schemes that are conceptually similar but differ in detail. An entraining plume model is used to find the mass flux. The mass flux is continuous through the cloud base. The lateral fractional entrainment rate is constant with height, while the detrainment-rate profile reduces the mass flux smoothly to zero at the cloud top. The eddy diffusivity comes from a turbulent kinetic energy–length scale formulation. The scheme has been implemented in a simple one-dimensional (single column) model. Results of simulations of a standard case that has been used for other model intercomparisons [Atmospheric Radiation Measurement (ARM), 21 June 1997] are shown and indicate that the scheme provides good results. The model also simulates the profile of a conserved scalar; this capability is applied to a case from the 1999 Southern Oxidants Study Nashville (Tennessee) experiment, where it produces good simulations of vertical profiles of carbon monoxide in a cloud-topped boundary layer.