Numerical experiments are carried out for a plane wave propagating in the atmospheric turbulence for a weak to strong fluctuation condition, i.e., the Rytov index being in a large range of 2\times {10}^{-3} to 20. Mainly two categories of propagation events are explored for the same range of Rytov index. In one category the propagation distance and also the Fresnel length are kept fixed with the turbulence strength changing. In the other the turbulence strength is kept fixed with the distance changing. The statistical characteristics of the scintillation index, the maximum and minimum of the intensity, the fractal dimension of the intensity image, and the number density of the phase singularity are analyzed. The behaviors of the fractal dimension and the density of the phase singularity present obvious differences for the two categories of propagation. The fractal dimension depends both on the Rytov index and the Fresnel length. In both weak and strong fluctuation conditions the dimension generally increases with the Rytov index, but is at minimum at the onset region. The phase singularity density is coincident with the theoretical results under a weak fluctuation condition, and has a slowly increasing manner with the Rytov index in the strong fluctuation condition. The dependence on the Fresnel size is confident and there is no saturation for the phase singularity.