We present our first numerical results of axisymmetric magnetohydrodynamic simula- tions for neutrino-cooled accretion tori around rotating black holes in general relativity. We consider tori of mass ∼ 0.1-0.4Maround a black hole of mass M =4 Mand spin a =0 - 0.9M ; such systems are candidates for the central engines of gamma-ray bursts (GRBs) formed after the collapse of massive rotating stellar cores and the merger of a black hole and a neutron star. In this paper, we consider the short-term evolution of a torus for a duration of ≈ 60 ms, focusing on short-hard GRBs. Simulations were performed with a plausible microphysical equation of state that takes into account neutronization, the nuclear statistical equilibrium of a gas of free nucleons and α-particles, black body radiation, and a relativistic Fermi gas (neutrinos, electrons, and positrons). Neutrino-emission processes, such as e ± capture onto free nucleons, e ± pair annihilation, plasmon decay, and nucleon- nucleon bremsstrahlung are taken into account as cooling processes. Magnetic braking and the magnetorotational instability in the accretion tori play a role in angular momentum redistribution, which causes turbulent motion, resultant shock heating, and mass accretion onto the black hole. The mass accretion rate is found to be ˙ M∗ ∼ 1-10M� /s, and the shock heating increases the temperature to ∼ 10 11 K. This results in a maximum neutrino emission rate of Lν = several ×10 53 ergs/s and a conversion efficiency Lν / ˙ M∗c 2 on the order of a few percent for tori with mass Mt ≈ 0.1-0.4Mand for moderately high black hole spins. These results are similar to previous results in which the phenomenological α-viscosity prescription with the α-parameter of αv =0 .01-0.1 is used. It is also found that the neutrino luminosity can be enhanced by the black hole spin, in particular for large spins, i.e., a & 0.75M ;i f the accretion flow is optically thin with respect to neutrinos, the conversion efficiency may be & 10% for a & 0.9M. Angular momentum transport, and the resulting shock heating caused by magnetic stress induce time-varying neutrino luminosity, which is a favorable property for explaining the variability of the luminosity curve of GRBs.