This dataset includes the parallel application and the main results supporting the research article "Modeling snow saltation: the effect of grain size and interparticle cohesion" published at the Journal of Geophysical Research: Atmospheres. The code is a flow solver based on the Large Eddy Simulation (LES) technique coupled with a Lagrangian Stochastic Model (LSM). The interaction of snow particles with the bed is modeled with statistical and physically-based models for aerodynamic entrainment, rebound and splash, following the works of Groot Zwaaftink et al. (2014), Comola and Lehning (2017) and Sharma et al. (2018). This algorithm was also used by Sigmund et al. (2021) to model snow sublimation.