  We present results of Monte Carlo simulations of the gauge glass model in
three dimensions using exchange Monte Carlo. We show for the first time clear
evidence of the vortex glass ordered phase at finite temperature. Using finite
size scaling we obtain estimates for the correlation length exponent, nu = 1.39
+/- 0.20, the correlation function exponent, eta = -0.47 +/- 0.07, and the
dynamic exponent z = 4.2 +/- 0.6. Using our values for z and nu we calculate
the resistivity exponent to be s = 4.5 +/- 1.1. Finally, we provide a plausible
lower bound on the the zero-temperature stiffness exponent, theta >= 0.18.
