“…The formulation can be gyroaveraged, for instance, for limiting the computational cost for electrons For instance, since the variation of the gyroradius among the particle species varies as 1∶40∶160 (e−:H+:O+), it is necessary to keep the gyrotrajectory when computing particle loss of heavy ions through the magnetopause (e.g., Sorathia et al, ). Global coupled MHD/test particle codes are well adapted, for instance, for azimuthal transport, that is, solving for particle gradient‐curvature drift motion, for rapid particle energization occurring during interplanetary shocks on the front end of coronal mass ejections (e.g., Hudson et al, ; Kress et al, , ), for drift‐orbit bifurcation trajectory (Ukhorskiy et al, ), for acceleration at dipolarization fronts (Ukhorskiy et al, in this collection), for solar wind ion entering the magnetosphere (Sorathia et al, ), for energetic particle injections in the inner magnetosphere during substorms (e.g., Gkioulidou et al, ), or O+ ion outflow directly injected within the radiation belts (Gkioulidou et al, ), or for the sudden depletion (e.g., Ukhorskiy et al, ) and rapid recovery of the outer belt (e.g., Sorathia et al, , in this collection). These codes can also be used to generate diffusion coefficients (e.g., Ukhorskiy & Sitnov, ).…”