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, ν = 1.39 ± 0.20, the correlation function exponent, η = −0.47 ± 0.07, and the dynamic exponent z = 4.2 ± 0.6. Using our values for z and ν 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, θ ≥ 0.18.