Abstract The effect of magnetic field ripple on tokamak plasma without turbulence is studied numerically and augmented with a reduced analytical model that includes neoclassical processes in the presence non-axisymmetric perturbation stochastic transport. For this study, has been implemented GYSELA gyrokinetic code. This implementation verified thanks to test toroidal angular momentum conservation. code was then successfully benchmarked against NEO code, which solves drift kinetic equation, collisionality range <mml:math xmlns:mml=“http://www.w3.org/1998/Math/MathML" overflow=“scroll”> mml:msup mml:miν</mml:mi> mml:mo⋆</mml:mo> </mml:msup> mml:mo∈</mml:mo> </mml:math> [0.05–0.5] for several amplitudes ripple. An observation, shared by model, simulations thermal drive mean poloidal velocity—measured $k_{V_P}$?> mml:msub mml:mik</mml:mi> mml:mrow mml:miV</mml:mi> mml:miP</mml:mi> </mml:msub> </mml:mrow> coefficient—decreases sharply large yet experimentally relevant amplitudes, may even change sign.