Gibbs function for region 1, 2 or 5: g(p,T,region)
function gibbs extends Modelica.Icons.Function; input SI.Pressure p "pressure"; input SI.Temperature T "temperature (K)"; input Integer region "IF97 region, 1, 2 or 5"; output Real g "dimensionless Gibbs funcion"; end gibbs;