Recent derived procedure of fitting the Rankine-Hugoniot shock jump relations for fast MHD shocks (Lin et al., 2006) is here applied to slow-mode shocks. The present paper is a continuation of Lin et al. (2006). We modify the method to be applicable to slow shocks. We have tested the new method via a simulated slow shock. The result shows that the method can find a best fit solution that is close to the "true answer" as long as the systematic errors in the synthetic data are not too large. We also apply the method to one interplanetary slow shock and two slow shocks observed in the magnetotail. The results have been compared to the previous studies. The new method of shock fitting is demonstrated successfully for identifying slow shocks.