Cut of results of our two simulations RunN (left panels) and RunS (right panels) along the subsolar line (x-axis). Panels a–b: proton and electron densities. Panels c–d: magnetic field components and amplitude. Panels e–f: current density along Y_MSO of protons, electrons, and total. Panels g–h: magnetic P_mag ≡ B²/2μ₀, ram , and total pressures. These cuts were performed at the same time as Fig. 1 (t = 15.9 R/V_x). We highlight the bow shock and magnetopause stand-off distances obtained from our model (solid red vertical lines) and the ones obtained using the average model of Winslow et al. (2013) constrained by MESSENGER observations (dashed red vertical lines).