On the deviation from Maxwellian of the ion velocity distribution functions in the turbulent magnetosheath