Reference liquids can be useful in the new quantitative approach to elastohydrodynamic lubrication calculations only if the properties of these liquids are accurately modeled. The pressures in the contact sometimes exceed the pressures of the available viscometer measurements, requiring extrapolation. The previtreous piezoviscous response is always faster-than-exponential with an apparent divergence pressure. Therefore, any model used for extrapolation must possess these characteristics and up to now published work in EHL, where extrapolation of the viscosity of squalane is performed, has not. To date, the only model that has been shown to accurately extrapolate to high pressure is the Hybrid model. This expression was applied to six reference liquids for use in quantitative investigations in EHL. Viscosity measurements for trioctyl trimellitate were extended to 1.1 GPa. The results obtained with high pressure viscometers are preferable over molecular dynamics results if they are different. Molecular dynamics calculations are not currently capable of reliable results for viscosity at very high pressures.