“…where μ = 1.265 is the mean atomic mass (assuming solar metal abundances; Anders & Grevesse 1989), m H is the mass of the hydrogen atom, n H is the hydrogen number density, p t is the total pressure, that is, the sum of the thermal pressure and the magnetic pressure (the factor 1/ √ 4π is absorbed in the definition of B), E is the total energy density, that is, the sum of the thermal energy density (ρ ), the kinetic energy density and the magnetic energy, u is the plasma velocity, g is the solar gravity, Λ(T ) is the radiative loss function for optically thin plasma, F c is the conductive flux, H is a heating function whose only role is to keep the unperturbed atmosphere in energy equilibrium, c s is the sound speed for an isothermal plasma, Φ is a free parameter (<1, Giuliani 1984) that determines the degree of saturation of the thermal conduction; we set Φ = 0.9, which corresponds to quite an efficient conduction. The radiative losses were computed according to version 7 of the CHIANTI code (Landi et al 2012), assuming a density of 10 9 cm −3 and ionization equilibrium according to Dere (2009).…”