Аннотация:We discuss the development of the large-eddy simulation (LES) model of the atmospheric boundary layer with embedded two-moment bulk cloud microphysics scheme well-suited for massively-parallel heterogeneous supercomputers based on GPU (Graphics Processing Units) architecture. To evaluate the LES model and its computational efficiency, we consider the numerical setup corresponding to the development of an intense Arctic cold-air outbreak case. It is shown that the dynamic closure approach for calculation of subgrid scale fluxes, applied to both heat and moisture transport, allows to correctly reproduce moist convective boundary layers with mixed-phased clouds even with coarse grid resolution. Implementation of state-of-the-art microphysics scheme for GPU systems not only led to significant speedup of the computations, but in general improved the multi-GPU scaling of the model.