A semiclassical surface hopping Herman-Kluk (HK) propagator is developed for multisurface, multidimensional nonadiabatic problems. Its application for two-surface model problems is presented for two-dimensional collinear systems. Both noncrossing and crossing cases are investigated. The calculation on the basis of Monte Carlo method shows that the general multisurface semiclassical HK propagator method is efficient and accurate.