The study deals with the development of a thermohydrodynamic (THD) computational procedure for evaluating the pressure, temperature and velocity distributions in fluid films with very rough geometry. A parametric investigation is performed to predict the bearing behaviors in the lubricating film having the absorbed layers and their interfaces determined by the rough surfaces with Gaussian distribution. The layers are expressed as functions of the standard deviations of each surface to characterize flow patterns between both the rough sur-faces. The velocity variations and the heat generation are assumed to occur in the central (shear) zone with the same bearing length and width. The coupled effect of surface roughness and shear zone dependency on hydrodynamic pressure and temperature has been found in non-contact mode. The procedure confirms the numerically determined relationship between the pressure and film gap on condition that its roughness magnitude is smaller than the fluid film thickness.