3D shape of our domain (x, y and z axis with units of ×10 Mm) with over-plotted, normalised number density values at t = 261 min of the reference case. The coordinate system marked in orange describes the s and t axis, along and transverse to the magnetic field lines, respectively. The v axis is perpendicular to the magnetic field lines (and does not exist in our 2D domain). It is important to note that the x, y and z axis are not to scale.