We present results from the simultaneous application of planar laser-induced fluorescence (PLIF), particle image velocimetry (PIV) and particle tracking velocimetry (PTV), complemented by direct numerical simulations, aimed at the detailed hydrodynamic characterization of harmonically excited liquid-film flows falling under the action of gravity. The experimental campaign comprises four different aqueous-glycerol solutions corresponding to four Kapitza numbers (Ka = 14, 85, 350, 1800), spanning the Reynolds number range Re = 2.3-320, and with forcing frequencies f w = 7 and 10 Hz. PLIF was employed to generate spatiotemporally resolved film-height measurements, and PIV and PTV to generate two-dimensional velocity-vector maps of the flow field underneath the wavy film interface. The latter allows for instantaneous, highly localized velocity-profile, bulk-velocity, and flow-rate data to be retrieved, based on which the effect of local film topology on the flow field underneath the waves is studied in detail. Temporal sequences of instantaneous and local film height and bulk velocity are generated and combined into bulk flow-rate time series. The time-mean flow rates are then decomposed into steady and unsteady components, the former represented by the product of the mean film height and mean bulk velocity and the latter by the covariance of the film-height and bulk-velocity fluctuations. The steady terms are found to vary linearly with the flow Re, with the best-fit gradients approximated closely by the kinematic viscosities of the three examined liquids. The unsteady terms, typically amounting to 5%-10% of the mean and peaking at approximately 20%, are found to scale linearly with the film-height variance. And, interestingly, the instantaneous flow rate is found to vary linearly with the instantaneous film height. Both experimental and numerical flow-rate data are closely approximated by a simple analytical relationship with only minor deviations. This relationship includes terms such as the wave speed c and mean flow rate Q, which can be obtained by relatively simple and inexpensive methods, thus allowing for spatiotemporally resolved flow-rate predictions to be made without requiring explicit knowledge of the full flow-field information underneath the wavy interface.