IAPWS95 is a pressure-explicit equation of state, i.e. properties are expressed as a function of density and temperature. However, generating the reweighted ensemble as a function of both temperature and pressure using molecular dynamics is much easier if we have properties as a function of pressure and temperature. We therefore back-calculated a density giving the desired pressure at a given temperature and used this density and temperature to estimate properties from IAPWS95 equation of state. Specifically, at each temperature, we find the density which gives P sim -P IAP W S95 (ρ,T ) = 0, using 1000 kg/m 3 as an initial guess density. This process converges to the stable minima for all densities and pressures used in this study. The root is found using the fsolve module from scipy 1. The plots show the (first column, top-to-bottom) ρ exp (T,P), G res,exp (T,P), U res,exp (T,P) property surfaces, (second column, top to bottom) V res,exp (T,P) (low pressure), V res,exp (T,P) (high pressure) and S res,exp (T,P) property surfaces and (third column, top to bottom) C Pres,exp (T,P), κ(T,P), α(T,P) property surfaces all estimated using the IAPWS95 equation of state for liquid water. V res,exp (T,P) is shown in two separate plots as the magnitude of V res is significantly different for pressures 1-100 bar and 100-5000 bar.