Propagation of magnetic island caused by drift-tearing instability is numerically investigated based on a reduced set of two-fluid equations. It is found that the island propagates into the ion diamagnetic direction when the island growth is saturated, and the propagation velocity becomes small as the viscosity increases. The island propagates with zonal flow generated by drift-tearing instability, because flattening of pressure inside the island is enhanced by parallel ion velocity. The mechanism of zonal flow generation depends on the viscosity. When the viscosity is small, the flow driven by the Reynolds stress is diminished by the Maxwell stress, and the small difference between them is overcome by the ion diamagnetic stress, and thus the zonal flow directs toward the ion diamagnetic direction. When the viscosity is large, the viscous stress is counteracted by the Reynolds and the Maxwell stresses, and the small difference between them is overcome by the ion diamagnetic stress, and thus zonal flow directs toward the ion diamagnetic direction.