Dimensionless stresses as a function of time for simulations sg-mhd-tau20 (τ_cΩ₀ = 20), sg-mhd-tau10 (τ_cΩ₀ = 10), sg-mhd-2 (τ_cΩ₀ = 10, Q₀ ∼ 0.75), and MRI-compare. Shown in blue is the Reynolds contribution, α_r, in red the gravitational contribution, α_g, and in black the Maxwell part, α_m. All stresses were calculated using Eqs. (8) and (9), and the factor γ = 1.64 was also used for the isothermal MRI, allowing for an easier comparison. The curves for sg-mhd-tau20 and 10 start at tΩ₀ = 120, at which point the simulations were restarted from the pure-GI state of GI072 (see Löhnert & Peeters 2022). Then, the cooling time was set to τ_cΩ₀ = 20, and background heating was turned off (only cooling). At tΩ₀ = 200, either one of two things happens, depending on the simulation. In sg-mhd-tau20, the ZNF field was introduced directly at tΩ₀ = 200, highlighted by the solid vertical line. In sg-mhd-tau10, the pure-GI phase was prolonged from tΩ₀ = 200 (vertical dashed line) to tΩ₀ = 400 (vertical solid line), but the cooling time was reduced to τ_cΩ₀ = 10. At tΩ₀ = 400, the ZNF field was introduced. For both sg-mhd-tau20 and 10, the values of α_g and α_r were smoothed by convolving the respective time series with a Gaussian function (standard deviation σ = 3).