  We study the oscillations of relativistic stars, incorporating key physics
associated with internal composition, thermal gradients and crust elasticity.
Our aim is to develop a formalism which is able to account for the
state-of-the-art understanding of the complex physics associated with these
systems. As a first step, we build models using a modern equation of state
including composition gradients and density discontinuities associated with
internal phase-transitions (like the crust-core transition and the point where
muons first appear in the core). In order to understand the nature of the
oscillation spectrum, we carry out cooling simulations to provide realistic
snapshots of the temperature distribution in the interior as the star evolves
through adolescence. The associated thermal pressure is incorporated in the
perturbation analysis, and we discuss the presence of $g$-modes arising as a
result of thermal effects. We also consider interface modes due to
phase-transitions and the gradual formation of the star's crust and the
emergence of a set of shear modes.
