“…Here, symbols n ( i ) denote the mass density of the neutrals (ions), V n (V i ) the velocity of the neutrals (ions), p n (p ie ) the neutral (ion+electron) gas pressure with p ie = p i + p e = 2p i , B the magnetic field, η the magnetic resistivity coefficient taken in its simplified form as η = ( i ν ei + n ν en )/(e 2 n 2 e ) where ν en and ν ei are respectively the electron-neutral and electron-ion collisions frequencies [47], μ the magnetic permeability of the medium, γ = 5/3 the adiabatic index, T i (T n ) is the proton (hydrogen atom) temperature, μ i = μ n = 1 are the mean masses of plasma species and neutrals, respectively, m H is the hydrogen mass, k B is Boltzmann's constant, g = [0, −g, 0] the vector of gravitational acceleration with its magnitude g = 274.78 m s −2 , acting along the y-axis in Cartesian coordinates. The reaction rates of the electron impact ionization and radiative recombination, Γ ion,rec i,n , momentum collisional, S i,n , and energy source, Q i,n , terms are described in great detail in [42]…”