A two‐dimensional, time‐dependent model has been constructed to study the zonally averaged structure of the middle atmosphere (16–116 km) allowing interaction among dynamics, radiation and photochemistry. The zonally averaged dynamics are governed by a stream function equation for the residual Eulerian meridional circulation wherein the effects of wave transience and dissipation have been neglected. The resulting circulation is thus driven by diabatic heating and cooling, with Rayleigh friction introduced to balance the momentum budget. The temperature structure is computed from the zonal mean thermodynamic equation, while the appropriate continuity equations are solved to determine the distribution of the various chemical species and families. In the chemical continuity equations, two types of eddy transport processes are present in addition to the transport by the residual Eulerian circulation. Small‐scale disturbances are assumed to produce turbulent mixing which is modeled in terms of horizontal and vertical diffusion coefficients, while steady state planetary waves give rise to fluxes which are expressed in terms of observed wave structure and the photochemical lifetime of the species transported. The circulation, temperature structure, and distribution of chemical constituents obtained with the model are for the most part in satisfactory agreement with observations. In particular, the direct, equator‐to‐pole circulation computed for the stratosphere seems to be capable of explaining the major features of the distribution of trace constituents there, although the effect of the Fickian diffusion included in the model is important in determining the magnitude of gradients.