Statistics of scalar concentration in a double scalar mixing layer are calculated using a Lagrangian stochastic model coupled to the interaction by exchange with the conditional mean micromixing model. Excellent agreement is obtained with recent direct numerical simulation results. The model accurately reproduces a wide variety of behavior across both physical and concentration space and over both unconditional and conditional scalar statistics. (c) 2006 American Institute of Physics.