  We present molecular dynamics simulation measurements of the viscosities of
the Gay-Berne phenomenological model of liquid crystals in the nematic and
isotropic phases. The temperature dependence of the rotational and shear
viscosities, including the nonmonotonic behavior of one shear viscosity are in
good agreement with experimental data. The bulk viscosities are significantly
larger than the shear viscosities, again in agreement with experiment.
