Neoclassical polarization drift is known to play critical role in the dynamical behavior of a sheared radial electric field E-r in a toroidal confinement device. However, basic studies on the effect of radial electric shear on neoclassical polarization drift have not yet appeared in the literature. In the present report, the neoclassical polarization drift speed V-NP of collisionless single ions is studied using a guiding-center code in a time-varying, spatially sheared E-r in a realistic tokamak geometry. It is found numerically that the V-NP for single ions is not only a function of the time derivative partial derivative E-r/partial derivative t, but also a strong function of the radial shear Delta r partial derivative E-r/partial derivative r if the shear length is on the same order as the ion banana width Delta r. Comparison with an analytic investigation reveals that this effect is simply due to the finite banana modification to the orbital average E-r. An approximate analytic formula has been presented for collisionless single banana ions in a conventional tokamak magnetic geometry. The trapped-passing boundary layer physics is not treated. (c) 2006 American Institute of Physics.