The pressure distributions on the entire computation domain are obtained through solving a unified Reynolds equation system without identifying hydrodynamic or asperity contact regions.