(a–c) Distributions of the relative angle θ_ρ, θ_v, and θ_B between the local magnetic field and the rotated (by 90^∘) density gradient ∇ρ, velocity gradient ∇v, and magnetic field strength gradient ∇B, respectively. The relative angles are evaluated separately for the WNM (T  >  5000 K), UNM (200 K <T <  5000 K), and CNM (T  <  200 K). Here a relative angle of 0^∘ corresponds to parallel alignment, while 90^∘ indicates perpendicularity. (d) Distribution of the sonic Mach number M_s = v/c_s in each phase. (e) Distribution of the thermal plasma beta, ﹩{\beta }_{{\rm{th}}}=(\rho {c}_{s}^{2}/\gamma )/({B}^{2}/8\pi )﹩. (f) Distribution of the turbulent plasma beta, β_tur = ρv²/(B²/8π). All quantities are phase separated based on temperature and computed in three dimensions.