  ESTER two-dimensional models solve the steady state structure of fast
rotating early-type stars including the large scale flows associated with the
baroclinicity of the radiative zones. Models are compared successfully to the
fundamental parameters of the two main components of the triple system $\delta$
Velorum that have been derived from interferometric and orbit measurements.
Testing the models on the Be star Achernar ($\alpha$ Eri), we cannot reproduce
the data and conclude that this star has left the main sequence and is likely
crossing the Herzsprung gap. Computing main sequence evolution of fast rotating
stars at constant angular momentum shows that their criticality increases with
time suggesting that the Be phenomenon and the ensuing mass ejections is the
result of evolution.
