The driving force dependence of the surface velocity and the average height of faceted merged steps, the terrace-surface slope, and the elementary step velocity are studied using the Monte Carlo method in the nonequilibrium steady state. The Monte Carlo study is based on a lattice model, the restricted solid-on-solid model with point-contact-type step-step attraction (p-RSOS model). The main focus of this paper is a change of the “kink density” on the vicinal surface. The temperature is selected to be in the step-faceting zone [N. Akutsu, AIP Adv. 6, 035301 (2016)] where the vicinal surface is surrounded by the (001) terrace and the (111) faceted step at equilibrium. Long time simulations are performed at this temperature to obtain steady states for the different driving forces that influence the growth/recession of the surface. A Wulff figure of the p-RSOS model is produced through the anomalous surface tension calculated using the density-matrix renormalization group method. The characteristics of the faceted macrostep profile at equilibrium are classified with respect to the connectivity of the surface tension. This surface tension connectivity also leads to a faceting diagram, where the separated areas are, respectively, classified as a Gruber-Mullins-Pokrovsky-Talapov zone, step droplet zone, and step-faceting zone. Although the p-RSOS model is a simplified model, the model shows a wide variety of dynamics in the step-faceting zone. There are four characteristic driving forces: , and . For the absolute value of the driving force, is smaller than , the step attachment-detachments are inhibited, and the vicinal surface consists of (001) terraces and the (111) side surfaces of the faceted macrosteps. For , the surface grows/recedes intermittently through the two-dimensional (2D) heterogeneous nucleation at the facet edge of the macrostep. For , the surface grows/recedes with the successive attachment-detachment of steps to/from a macrostep. When exceeds , the macrostep vanishes and the surface roughens kinetically. Classical 2D heterogeneous multinucleation was determined to be valid with slight modifications based on the Monte Carlo results of the step velocity and the change in the surface slope of the “terrace.” The finite-size effects were also determined to be distinctive near equilibrium.