function jouleThomsonCoefficient input ThermodynamicState state "Thermodynamic state"; output Real my(unit = "K/Pa") "Isothermal throttling coefficient"; end jouleThomsonCoefficient;