This dataset contains data from reported Molecular Dynamics (MD) simulations of shear-driven liquid flow over an infinite gas nanofilm covering a solid surface. Four computationally expensive simulations are considered, which differ from each other from the rarefaction of the gas nanofilm (i.e. the Knudsen number). In each data folder one finds the final steady-state profiles taken normal from the surface, and the time-varying properties in the gas/liquid bulk.