  We present a Monte-Carlo simulation of energy deposition process in
relativistic heavy-ion collisions based on a new realization of the
Interacting-Gluon-Model (IGM) for high energy $N-N$ collisions. In particular
we show results for proton spectra from collisions of $E_{lab}=200 \ GeV/N$
$^{32}$S beam incident on $^{32}$S target and analyze the energy and mass
dependence of nuclear stopping power predicted by our model. Theoretical
predictions for proton rapidity distributions of both $^{208}$Pb + $^{208}$Pb
collisions at $E_{lab}=160 \ GeV/N$ CERN SPS and $^{197}$Au + $^{197}$Au at
$\sqrt{s_{NN}}=200 \ GeV$ BNL RHIC are given.
