This paper is focused on the numerical implementation of thermo-viscoplastic constitutive equation in finite strain contexts using FEM. A length scale parameter is introduced implicitly through viscosity in order to address strain localization in the boundary value problems. The constitutive law is formulated in the corotational material coordinate system. Thus, the evolution equations of all internal variables can be expressed in terms of their material time derivatives. The numerical algorithm for implementing the material model in a finite element program is also formulated, and several numerical examples are carried out. As results, they show good performances of the present framework in describing the strain localization problem.