Abstract.  Toroidal rotation plays an important role in transport and stability of burning plasmas, and clarification of the mechanism of rotations spontaneously generated by perpendicular NBI and ICRF can help evaluate the rotations in the plasmas quantitatively, contributing to accurate control of them. In order to study the properties of plasma rotations, one-dimensional multi-fluid transport code TASK/TX has been developed, which is the first code to self-consistently solve a tokamak plasma evolution including a plasma rotation and a radial electric field E_r. TASK/TX automatically generates a return current j_bulk due to the charge separation originated from the difference in electron and ion banana orbit widths, resulting in a × torque on a bulk plasma. Through the simulations, we confirm the generated torque drives a toroidal rotation accompanied by the modification of a radial electric field without having to add extra torque input term in equations of motion, while a poloidal rotation almost remains unchanged. It is found that the toroidal rotation decreases with a weaker dependence than inversely proportional as the turbulent viscosity increases.