First row: left panel: rotational velocity (dashed blue curve) and circular velocity (solid black curve) at the final simulation time. Middle panel: m = 2 Fourier amplitudes, A₂/A₀, as a function of radius estimated from the stellar density. Curves of different colors present the time evolution of A₂/A₀. To better see the evolution of the bar strength with time, in the right panel we show the amplitude averaged over the bar maximum. Second row: face-on density maps of the stellar component for different times, as indicated. Third row: corresponding edge-on view. Contour spacing is logarithmic. Fourth row: changes in angular momentum, ΔL, as a function of radius, estimated in a time window of 0.52 Gyr, centered on the times of the snapshots shown above. Both axes are divided by the circular velocity, thus units are kpc (galactic radius). Strong variations are seen with cosmic time due to satellite perturbations and increase in bar strength.