An approximate analytic expression for neoclassical toroidal plasma viscosity in tokamaks that have error fields or magnetohydrodynamic activities is presented. The expression smoothly joins transport fluxes or plasma viscosity in all the known collisionality regimes derived from the solution of the bounce averaged drift kinetic equation and should be useful in modelling results of existing and future tokamak experiments. It also incorporates some of the extensions of the known expressions to include the effects of finite ∇B drift in the non-resonant transport processes. Here, B is the magnitude of the magnetic field. The toroidal momentum balance equation is a nonlinear function of the radial electric field when the neoclassical plasma viscosity is dominant. It can have bifurcated solutions for the radial electric field and may lead to better plasma confinement as a result.