We study numerically magnetosheath ion motion in the magnetopause current sheet in the presence of magnetic fluctuations. Test particles are injected at the magnetosheath side only. The magnetic field turbulence is modelled as a power law spectrum in phase space, which reaches maximum intensity in the center of the magnetopause current sheet, and decreases towards the magnetosheath and magnetosphere boundaries, creating magnetic islands. The number of particles entering the magnetosphere, reflected from the magnetopause and flowing away from the flanks is computed, as a function of the fluctuation level of the turbulence and of the magnetic field and flow geometry parameters. All those quantities appear to be strongly dependent on the fluctuation level, with the number of particles entering the magnetopause increasing with the fluctuation level. We calculate particle density, bulk velocity and temperature and other distribution function moments. Temperature and bulk velocity are increasing in the direction normal to the magnetopause, while density is decreasing, in a qualitative accordance with the observations.