Upper panels, left: Horizontal and temporal average of the unsigned vertical component of the magnetic Poynting flux in different runs as a function of height. For better visualization, the Poynting flux has been normalized over the density. Middle: Emerging Poynting flux, vertical component, computed as $\langle S_{\mathrm{EM}}^{\mathrm{emer}}\rangle =\langle |v_z(B_x^2+B_y^2)/{\mu }_0|\rangle$. Right: Shear Poynting flux, vertical component, computed as $\langle S_{\mathrm{EM}}^{\mathrm{shear}}\rangle =\langle |-B_z(v_x{B}_x+v_y{B}_y)/{\mu }_0|\rangle$. Black lines are for the Dyn model, red for the Bz50 model, and blue for the Bz200 model. The solid thick line is for 5 km resolution, the dashed medium-thick line is for 10 km resolution, and the dotted thin line is for 20 km resolution. Bottom panels represent the same but without dividing by density and taking into account the sign while averaging. Negative values mean downward propagating flux. The apparent discontinuity in the curves is due to representation, which is caused by the change of sign.