Source field modeling by self-consistent Gaussian beam superposition (two-dimensional)