A computational study of thin liquid films over a solid surface is reported. The lubrication equation is numerically solved using an in-house code, which implements the finite volume method. Small slope approximation is abandoned, and a more accurate model for capillary pressure estimation is presented, allowing us to correctly investigate higher contact angles, when compared to the maximum value allowed by small slope approximation. Disjoining pressure is used for modeling substrate wettability. The in-house solver is first validated: a 1D flowing film driven by gravity is simulated and the disjoining pressure model is verified for contact angles up to 60°; replicating literature experimental investigations, a uniform film covering an inclined plate is perturbed, inducing the generation of a large dry patch; rivulet buildup is simulated; and the numerical results are compared with fully 3D computations found in the literature and verified with analytical evidences. Then, a film flowing over an inclined plate bounded by lateral walls, which is a complex configuration commonly used for studying liquid behavior in structured packing, is investigated and relevant parameters are reported. [ABSTRACT FROM AUTHOR]