compute the saturated liquid specific h(p)
function hl_p = ThermoSysPro.Properties.WaterSteam.BaseIF97.Regions.hl_p;