  With (non-barotropic) equations of state valid even when the neutron, proton
and electron content of neutron star cores is not in beta equilibrium, we study
inertial and composition gravity modes of relativistic rotating neutron stars.
We solve the relativistic Euler equations in the time domain with a three
dimensional numerical code based on spectral methods, in the slow rotation,
relativistic Cowling and anelastic approximations. Principally, after a short
description of the gravity modes due to smooth composition gradients, we focus
our analysis on the question of how the inertial modes are affected by
non-barotropicity of the nuclear matter. In our study, the deviation with
respect to barotropicity results from the frozen composition of non-superfluid
matter composed of neutrons, protons and electrons, when beta equilibrium is
broken by millisecond oscillations. We show that already for moderatly fast
rotating stars the increasing coupling between polar and axial modes makes
those two cases less different than for very slowly rotating stars. In
addition, as we directly solve the Euler equations, without coupling only a few
number of spherical harmonics, we always found, for the models that we use, a
discrete spectrum for the $l = m = 2$ inertial mode. Finally, we find that, for
non-barotropic stars, the frequency of this mode, which is our main focus,
decreases in a non-negligible way, whereas the time dependence of the energy
transfer between polar and axial modes is substantially different due to the
existence of low-frequencies gravity modes.
