Wave propagations are best described by three major approaches: the spherical wave approach represented by Rayleigh integral, the planar wave approach represented by spatial Fourier analysis, and the Gaussian beam approach. Spherical and planar wave approaches, both in form of double integrals, are analytically inconvenient and computationally expensive. In contrast, the Gaussian beam approach, with huge advantages in analytical operability and computational efficiency, has been quietly but firmly winning popularity. Dr. Breazeale’s article: A diffraction beam field expressed as the superposition of Gaussian beams, co-authored with Wen, is not only a most representative work for Gaussian beam approach but also a most cited academic work in sound field analysis and computation. As a memory of Dr. Breazeale’s pioneer role in the raise of Gaussian beam method, this article presents a new enhancement to the existing Gaussian beam method. A modified form of Gaussian beam function introduced as the base function set for wave field analysis greatly reduces the error rooted from the parabolic approximation, the number of Gaussian beams needed, and the run-time computation burden. The process for calculating the beam coefficients and beam parameters is also greatly simplified.