  The spectral energy distribution (SED) of a non-spherical star could differ
significantly from the SED of a spherical star with the same average
temperature and luminosity. Calculation of the SED of a deformed star is often
approximated as a composite of several spectra, each produced by a plane
parallel model of given effective temperature and gravity. The weighting of
these spectra over the stellar surface, and hence the inferred effective
temperature and luminosity, will be dependent on the inclination of the
rotation axis of the star with respect to the observer, as well as the
temperature and gravity distribution on the stellar surface. Here we calculate
the surface conditions of rapidly rotating stars with a 2D stellar structure
and evolution code and compare the effective temperature distribution to that
predicted by von Zeipel's law. We calculate the composite spectrum for a
deformed star by interpolating within a grid of intensity spectra of plane
parallel model atmospheres and integrating over the surface of the star. Using
this method, we find that the deduced variation of effective temperature with
inclination can be as much as 3000 K for an early B star, depending on the
details of the underlying model.
