Two-dimensional solitary waves at the surface of a film flow down a vertical plane are considered. When the system is subjected to inlet white noise, solitary waves are formed after an inception region and interact with each other. Using open-domain simulations of reduced equation models, we investigate numerically their late time process dynamics. Close to the instability threshold, the waves synchronize themselves into bound states. For higher values of the Reynolds number, the separation distance between the waves increases and the synchronization process at work is weaker. Performing statistics, we show that the mean characteristics of the waves correspond to the minimal value of the mean film thickness along the traveling-wave branch of solutions. In this regime, synchronization occurs through the waves tails which is associated with a change of scaling of the waves features. A similar behavior is observed performing simulations in periodic domains: the selected waves maximize the mean flow rate. (C) 2016 Elsevier Ltd. All rights reserved.