Phase-space scatter plots of the normalized local integrated luminosity L_J/L_J,max for each H₂ line. Each point corresponds to a computational cell in the simulation, with point size scaled by L_J/L_J,max to reduce overlap and emphasize more emissive regions. The reddish areas highlight the densities and temperatures ranges that most strongly contribute to the total emission. Solid black lines mark approximate contours of poloidal wind velocity (v_p = 1, 3, and 5 km s⁻¹), and the green contour indicates the isothermal sonic surface. Dashed blue lines represent the best-fit rotational temperatures (Table 4).