We present the results of axisymmetric, two-dimentional magnetohydrodynamic (MHD) simulations of weakly ionized gas torus threaded by large scale vertical magnetic fields. The gas torus corresponds to the 100pc scale circumnuclear torus observed by HST in nearby AGNs (e.g. NGC4261) or 1010M⊙ circumnuclear gas found by CO observations in luminous IR galaxies and quasars (e.g. Scoville et al. 1991). The initial state of simulation is a constant angular momentum polytropic torus threaded by uniform vertical magnetic fields. The torus is assumed to be rotating in a static, spherical hot halo. The model parameters are Eth = vs02/(γvk02) = 5 ×10−3 and Eth = vA02/vK02 = 6.6×10−6 where γ is the adiabatic index and vs0 and va0 are the sound speed and the Alfvén speed at r = r0 respectively.