The neutral diffusion approximation in a consistent tokamak edge plasma-neutral computation