A widely-used procedure for the solution of high frequency vibroacoustic problems is statistical energy analysis (SEA). Statistical energy analysis is a lumped parameter approach that treats a buildup system as a collection of subsystems with average properties. While the SEA has found wide acceptance, the approach does not provide any information about the spatial variation of the dynamic response within a subsystem. Integral equations are well suited as a post processing tool to evaluate sound field within a domain using energy density and intensity at the surface of the domain. In the present work, energy flow integral equations are extended to post-process SEA results to predict interior sound field at desired locations.