High Energy Laser (HEL) propagation through turbulent atmosphere is examined via numerical simulation. The beam propagation is modeled with the paraxial equation, which in turn is written as a system of equations for a quantum fluid, via the Madelung transform. A finite volume solver is applied to the quantum fluid equations, which supports sharp gradients in beam intensity. The atmosphere is modeled via a coupled advection-diffusion equation whose initial data have Kolmogorov spectrum. In this model the combined effects of thermal blooming, beam slewing, and deep turbulence are simulated