Panel (a): Circularity plot showing the distribution of stellar orbits in the best-fit dynamical model of NGC 3521. Horizontal lines denote the separations between cold, warm, hot, and counter-rotating orbits, with default cut values of λ_z < 0.8; 0.8 > λ_z > 0.25. Panel (b): Enclosed mass profiles for the stellar (mass-follows-light) component, dark matter, and their sum. Panel (c): Intrinsic axis ratios q and p as functions of distance from the centre; T quantifies triaxiality. Vertical dashed lines indicate the bulge radius.