A full-F isothermal gyro-fluid model and code (which is based on the full distribution function F compared to only small fluctuations) is extended to handle self-consistent coupling of multiple quasi-neutral ion species via the polarisation equation in the long wavelength approximation. The numerical model is used to determine two-dimensional interchange driven ‘blob’ transport in a plasma with ...