We develop a Monte Carlo (MC) tool incorporated with the three-subband approximation model to investigate the in-plane spln-polarized transport in GaAs/GaAlAs quantum well. Using the tool, the effects of the electron occupation of higher subbands and the intersuhband scattering on the spin dephasing have been studied. Compared with the corresponding results of the simple one-snbband approximation model, the spin dephasing length is reduced four times under 0.125 kV/cm of driving electric field at 300K by the MC tool incorporated with the three-subband approximation model, indicating that the three-subbarld approximation model predicts significantly shorter spin dephasing length with temperature increasing. Our simulation results suggest that the effects of the electron occupation of higher subbands and the intersubband scattering on the spln-dependent transport of GaAs 2-dhuensional electron gas need to be considered when the driving electric field exceeds the moderate value and the lattice temperature is above 100K. The simulation by using the MC tool incorporated with the three-subband approximation model also indicates that, under a eertain driving electric field and lattice temperature, larger channel widths cause spins to be depolarized faster. Ranges of the three components of the spins are different for three different injected spin polarizations due to the anisotropy of spin-orbit interaction.