A new gyrokinetic toroidal particle code has been developed to study the ion temperature gradient (ITG) driven turbulence in reactor relevant tokamak parameters. We use a new method based on a canonical Maxwellian distribution FCM(Pφ,ε,μ), which is defined by three constants of motion in the axisymmetric toroidal system—the canonical angular momentum Pφ, the energy ε, and the magnetic moment μ. A quasi-ballooning representation enables linear and nonlinear high-m,n global calculations to be carried out, with a good numerical convergence. Conservation properties are improved by using optimized particle loading. From comprehensive linear global analyses over a wide range of unstable toroidal mode numbers (n = 0–100) in large tokamak parameters (a/ρti = 320–460), it is found that the reversed shear configuration produces an effective stabilizing effect on the ITG mode in the q min region through global effects. In the nonlinear simulation, it is found that the new method based on FCM can simulate a zonal flow damping correctly; and spurious zonal flow oscillations, which are observed in a conventional method based on a local Maxwellian distribution FLM(ψ,ε,μ), do not appear in the nonlinear regime.