function density_ph input AbsolutePressure p; input SpecificEnthalpy h; input FixedPhase phase = 0; output Density d; end density_ph;