Code, data, and analysis scripts for: Relativistic runaway electron avalanches: unified density-dependent scaling and transport

This dataset contains simulation results and analysis products used to investigate relativistic runaway electron avalanches in the atmosphere and to constrain key physical relations governing avalanche development and particle transport. The work combines two complementary components: the determination of the avalanche e-folding length through optimization of the Dwyer model parameter, and the quantitative description of particle propagation beyond the electric-field region using energy spectra and bootstrap-based analysis. Simulations were performed using the CORSIKA Monte Carlo code (version 7.7500), with uniform vertical electric fields applied over a 2000 m atmospheric layer above the detector level. The archive includes representative raw simulation outputs, processed datasets (CSV), analysis scripts, energy spectra, bootstrap results, and figures and tables used in the manuscript. Due to file-size limitations, only representative subsets of raw simulation outputs are included. The complete dataset is available from the author upon request. CORSIKA is not included due to licensing restrictions.