  It is argued that the use of the initial Gaussian energy density profile for
hydrodynamics leads to much better uniform description of the RHIC heavy-ion
data than the use of the standard initial condition obtained from the Glauber
model. With the modified Gaussian initial conditions we successfully reproduce
the transverse-momentum spectra, v2, and the pionic HBT radii (including their
azimuthal dependence). The emerging consistent picture of hadron production
hints that a solution of the long standing RHIC HBT puzzle has been found.
