Simulation output and analysis code for Allen et al.’s ‘Numerical Simulations of Supercell Storms Employing the Thin-Boundary-Layer Equations’

This dataset contains the simulation output and analysis code used by Allen et al. (2026) for their article, "Numerical Simulations of Supercell Storms Employing the Thin-Boundary-Layer Equations."  Three relatively high-resolution (75-m horizontal grid spacing), large-eddy-simulation (LES) ensembles of tornadic supercell storms are used to investigate the effect of a more realistic treatment of the lower boundary condition on supercell storm simulations. These ensembles vary only in the parameterization of near-surface turbulence: one ensemble uses the semislip scheme, whereas the other two employ versions of the "thin-boundary-layer equations'' (TBLE) approach.  Each simulated storm is a long-lived supercell, with an intense mesocyclone and at least one tornadolike vortex (TLV).  The primary differences observed across the ensembles are (i) the amplitude of near-surface turbulent eddies in the environmental boundary layer and (ii) the intensity of the TLVs, with both TBLE schemes resulting in stronger turbulent eddies in the environment and stronger TLVs within the simulated storms.  Although substantial variations in precipitation distribution and cold pool strength are also found among the simulations, the differences in these larger-scale storm attributes from one ensemble to another are not statistically significant.