function density_phX input Modelica.SIunits.Pressure p; input SpecificEnthalpy h; input MassFraction X[:] = fill(0, 0) "mass fraction m_NaCl/m_Sol"; input FixedPhase phase = 0 "2 for two-phase, 1 for one-phase, 0 if not known"; output Modelica.SIunits.Density d = Modelica.Media.Water.WaterIF97_pT.density_ph(p, h); end density_phX;