vSpec = spectrum() in vacuum vSpec = spectrum(fRefractiveIndex) 1D-T (XY Gauss assumed) 2D-TX (Y Gauss assumed) = int int dx dy |E(omega, x, y)|^2 [J s]