Returns the current phase
function phase_ps input Units.SI.Pressure p "pressure"; input Units.SI.SpecificEntropy s "specific entropy"; output Integer phase = 0 "phase: 2 for two-phase, 1 for one phase"; end phase_ps;