Non-potential magnetic field lines and magnetic isosurfaces calculated numerically by applying an affine shear transformation on the EPs with s_y = 0.6 and s_x = s_y/4. Gas pressure and gravity are included. In this model, the minimum plasma beta is around 1.4 × 10⁻³, while the maximum is around 2.5, which is well above one. In this example, we used a bipole with d = 0.25 and z₀ = 0.5 in the domain x_min = −2, x_max = 2, y_min = −2, y_max = 2, z_min = 0, and z_max = 4.