Overview of large-scale structure and magnetosphere boundaries in our simulations. Panels a–f show the results of RunN. Panels a–c show dipolar Y_MSO = 0 cuts of the ion density n_i (a), magnetic field amplitude |B| (b), and plasma current density along the y-axis J_y (c). Panels d–f show equatorial Z_MSO = 0 cuts of the same quantities for RunN. Panels g–l show the results of RunS in the same format as those of RunN. All plots represent a time of t = 15.9 R/V_x. The average bow shock (dashed black line) and magnetopause (solid black line) profiles found by Winslow et al. (2013) using MESSENGER observations are added. The white region around the planet corresponds to grid cells with zero macro-particles.