In the atmosphere, convection can organize from smaller scale updrafts into more coherent structures on various scales. In bulk-plume cumulus convection parameterizations, this type of organization has to be represented in terms of how the resolved flow would “feel” convection if more coherent structures were present on the subgrid. This type of subgrid organization acts as building blocks for larger scale tropical convective organization known to modulate local and remote weather. In this work a parameterization for subgrid (and cross-grid) organization in a bulk-plume convection scheme is proposed using the stochastic, self-organizing, properties of cellular automata (CA). We investigate the effects of using a CA which can interact with three different components of the bulk-plume scheme that modulate convective activity: entrainment, triggering, and closure. The impacts of the revised schemes are studied in terms of the model's ability to organize convectively coupled equatorial waves (CCEWs). The differing impacts of adopting the stochastic CA scheme, as compared to the widely used Stochastically Perturbed Physics Tendency (SPPT) scheme, are also assessed. Results show that with the CA scheme, precipitation is more spatially and temporally organized, and there is a systematic shift in equatorial wave phase speed not seen with SPPT. Previous studies have noted a linear relationship between Gross Moist Stability (GMS) and Kelvin wave phase speed. Analysis of GMS in this study shows an increase in Kelvin wave phase speed and an increase in GMS with the CA scheme, which is tied to a shift from large-scale precipitation to convective precipitation.