This paper describes the ecology of drift wave turbulence and mean flows in the coupled drift-ion acoustic wave plasma of a CSDX linear device. A 1D reduced model that studies the spatiotemporal evolution of plasma mean density n, and mean flows v y and v z , in addition to fluctuation intensity e, is presented. Here, e ¼ hñ 2 þ ðr ?/ Þ 2 þṽ 2 z i is the conserved energy field. The model uses a mixing length l mix inversely proportional to both axial and azimuthal flow shear. This form of l mix closes the loop on total energy. The model self-consistently describes variations in plasma profiles, including mean flows and turbulent stresses. It investigates the energy exchange between the fluctuation intensity and mean profiles via particle flux hñṽ x i and Reynolds stresses hṽ xṽy i and hṽ xṽz i. Acoustic coupling breaks parallel symmetry and generates a parallel residual stress P res xz. The model uses a set of equations to explain the acceleration of v y and v z via P res xy / r n and P res xy / r n. Flow dynamics in the parallel direction are related to those in the perpendicular direction through an empirical coupling constant r VT. This constant measures the degree of symmetry breaking in the hk m k z i correlator and determines the efficiency of r n in driving v z. The model also establishes a relation between r v y and r v z , via the ratio of the stresses P res xy and P res xz. When parallel to perpendicular flow coupling is weak, axial Reynolds power P Re xz ¼ Àhṽ xṽz ir v z is less than the azimuthal Reynolds power P Re xy ¼ Àhṽ xṽy ir v y. The model is then reduced to a 2-field predator/prey model where v z is parasitic to the system and fluctuations evolve self-consistently. Finally, turbulent diffusion in CSDX follows the scaling: D CSDX ¼ D B q 0:6 ? , where D B is the Bohm diffusion coefficient and q ? is the ion gyroradius normalized to the density gradient jr n= nj À1 .