Time series of the snapshots of neutral and charged densities for the simulations L1-WN (top left group of panels), L1-WN-NN (top right group of panels), L1-WN-NN-B (bottom left group of panels), P-WN-NN (bottom right group of panels). The in plane (x and z components) velocities of neutrals and charged particles are represented by orange arrows in the plots which show the densities of neutrals and charged particles, respectively. For the simulations with sheared magnetic field (L1-WN, L1-WN-NN, L1-WN-NN-B), the iso-contours of the absolute value of the y-component of the magnetic field potential, |A_y|, are shown at the bottom panels. The iso-contours are equally spaced between 0.6 and , where is the maximum value of |A_y|.