calculates saturation pressure of water, with the equation given in Duan2003
function p_sat_H2O_Duan2003 input SI.Temp_K T; output SI.Pressure p_sat; end p_sat_H2O_Duan2003;