  A set of equations are developed that extend the macroscale magnetic
reconnection simulation model kglobal to include particle ions. The extension
from earlier versions of kglobal, which included only particle electrons,
requires the inclusion of the inertia of particle ions in the fluid momentum
equation, which was not required in the electron-only model. The new equations
will facilitate the exploration of the simultaneous non-thermal energization of
ions and electrons during magnetic reconnection in macroscale systems.
Numerical tests of the propagation of Alfv\'en waves in a plasma with
anisotropic electron and ion pressure are presented to benchmark the new model.
