A phase field model for step dynamics on vicinal surface is presented. Using this model, time dependent, collective motions of steps were investigated. Through numerical simulations, morphological step instabilities induced by ES barrier were analyzed, and it is shown that this model could interpret various phenomena during step flow growth such as step bunching and meandering.