The simulation of the nonlinear effects of laser atmospheric propagation is reported in this paper. A set of atmospheric hydrodynamics equations and nonlinear electromagnetic wave equation which interacts the atmosphere are given. The wave equation describing the stimulated thermal Rayleigh Scattering(STRS) induced small perturbation instability is derived. And turbulence is theoretically treated. Basing on thais equations, we per formed a set of 4D programs(TURBLM) to simulate thermal blooming and phase compensation for it, small perturbation instability, and the interation of turbulence with thermal blooming. The simulated results accord with theories and experiments.