  We make an extensive numerical study of a two dimensional nonconservative
model proposed by Olami-Feder-Christensen to describe earthquake behavior. By
analyzing the distribution of earthquake sizes using a multiscaling method, we
find evidence that the model is critical, with no characteristic length scale
other than the system size, in agreement with previous results. However, in
contrast to previous claims, we find convergence to universal behaviour as the
system size increases, over a range of values of the dissipation parameter,
$\alpha$. We also find that both ``free'' and ``open'' boundary conditions tend
to the same result.
  Our analysis indicates that, as $L$ increases, the behaviour slowly converges
toward a power law distribution of earthquake sizes
  $P(s) \sim s^{-\tau}$ with exponent $\tau \simeq 1.8$. The universal value of
$\tau$ we find numerically agrees quantitatively with the empirical value
($\tau=B+1$) associated with the Gutenberg-Richter law.
