Calculation of fugacity coefficient according to (Duan 2003)
function fugacity_H2O input SI.Pressure p; input SI.Temp_K T; output Real phi; end fugacity_H2O;