The neoclassical return current and parallel viscositydue to the radial electric field (non-equilibrium poloidal rotation) in tokamaks arecalculated using the 5-D (3-D in configuration space and 2-D in velocity space) Monte Carlo code ASCOT.Results are compared with the known neoclassicalnon-ambipolar return current expressions by Shaing [Phys. Fluids B 2 (1990) 2847]and Stringer [Nucl. Fusion 33 (1993) 1249] over a wide range of poloidal Mach numbers Mp = Er/(vtiBθ) and for various collisionalities.Here, vti is the ion thermal velocity, Er is the radial electric field and Bθ is the poloidal magnetic field. Quantitative differences exist due to a simpler collision operator and thevarious other approximations used in deriving the analytic expressions.