First, the simulation has a laminar or a Sweet-Parker reconnection phase, where only one elongated current sheet is formed and develops inflow and outflow velocity fields corresponding to the \(y\)  and \(z\)  direction, respectively [1]}, [2]}, [3]} due to the initial structure of the magnetic field. This can be clearly identified for both MFMS and TF in the first column of Fig. REF , as well as in the animation corresponding to this figure. After this first phase, we have a plasmoid instability phase where many plasmoids are formed and are advected from the middle of the reconnection region towards the outflow region along the \(z\)  axis (the last column of Fig. REF ). During this phase, the presence of plasmoids breaks the current sheet into multiple thinner current sheets, which reconnect faster than in the Sweet-Parker regime [4]}. This is also in agreement with previous 2D [5]} and 3D [6]} simulations.
