  The vast separation dividing the characteristic times of energy confinement
and turbulence in the core of toroidal plasmas makes first-principles
prediction on long timescales extremely challenging. Here we report the
demonstration of a multiple-timescale method that enables coupling global
gyrokinetic simulations with a transport solver to calculate the evolution of
the self-consistent temperature profile. This method, which exhibits resiliency
to the intrinsic fluctuations arising in turbulence simulations, holds
potential for integrating nonlocal gyrokinetic turbulence simulations into
predictive, whole-device models.
