First thing generally is to sort the transitions by wavelength/frequency however you want to do things.
You can do this in excel, but it is more efficient to write a program to do this.
If you do this in excel, you want to create a set of x values that are uniformly spaced, with a point density that is good enough to give you a smoothness that you can live with. I.e. 3 points over your FWHM is probably too coarse, but 3,000 is way to many. You decide. Populate a column vector with these values.
You can create the spectrum of Gaussian lines in a variety of ways. The simplest (but computationally inefficient way) is to add up all of the contributions for all of the lines at each X grid point.
e.g. if x = 8.0000, y = SUM(Intensity_i*(exp[-(8.0000-x_i)^2/sigma] ) Where "Intensity_i" are your y-values in your original spreadsheet, the "x_i" are the x values in your original spreadsheet, and sigma is obtained from the Gaussian that you are assuming. If the different transitions have different effective Gaussian widths, you could have a third column in your original spreadsheet. In excel, you may get a floating point error if you are too far away from the line (i.e. if x is too far away from x_i), in which case you may have to play with the limits of the sum. In a program, you could fix these with something that checks to see if x is more than 4 sigma away from the xi, and not bother to calculate the intensity if the line falls outside of this range.