FIG. 2: Radial tokamak edge pedestal profiles of the plasma pressure before and after an ELM blowout event across the separatrix into the scrape-off layer (SOL). The area shaded in grey shows the present radial numerical simulation domain. spatial range from the global scale to smaller than the ion gyroradius scale, with a ratio δ = ρ$_{s}$/a = 0 . 001875 between ion gyroradius ρ$_{s}$ and minor torus radius a . In the SOL region ( r > a ) the parallel boundary condition is replaced by a Debye sheath model whose treatment is given in Ref. [44]. The initial conditions are thus prescribed and are based on experimentally diagnosed radial temperature and density pedestal profiles T ( r ) and n ( r ) for each species (electrons and ions). A consistent electrostatic potential φ ( r ) is derived by numerically solving the neoclassical balances (parallel dynamics, toroidal drifts, and collisional dissipation, but not transport) in a pre-processing step with a modified (zonally frozen) GEMR setup, resulting in a time-steady 2D dissipative solution. The s − α MHD equilibrium in turn is set internally in GEMR by computing the Shafranov shift, and in the present setup would not be consistently described by prescribed in terms of external or coupled shaped equilibrium solvers [18]. The parallel and perpendicular electron and ion temperatures, T$_{e}$$_{‖}$ , T$_{e}$$_{⊥}$ , T$_{i}$$_{‖}$ and T$_{i}$$_{⊥}$ , are directly adopted and fixed from experimentally derived values by filtering the zonal component out of the total time derivative ∂T/∂t = S$_{T}$ . Their values are zonally frozen by setting where 〈 S$_{T}$ 〉 is the zonal (parallel and perpendicular) average of all right-hand-side terms S$_{T}$ 17 ∂T ∂t = S$_{T}$ − 〈 S$_{T}$ 〉 (52) V. PROFILE PRE-EQUILIBRATION