The radial profiles of the Gauss–Hermite representations (up to h ₄) with 1σ errors as derived from the spectroscopic data (gray) and the corresponding fit of the best dynamical orbit model we found for each galaxy (colored). The 1σ errors are quantified by fitting Gauss–Hermites to 100 Monte Carlo realizations on the velocity-binned LOSVDs. The models fit the observed kinematics quite well considering deviations from axisymmetry in the data cannot be emulated by axisymmetric models.