This function calculates the entropy of saturated vapor by taking the limit as the pressure approaches the saturation pressure from below.
However, sSatVapor[0.85]
returns Limit[obj$3516393[m], m -> 1, Direction -> 1].
Thanks.
sSatVapor[t_] := Module[{psat, obj},
psat = ThermodynamicData["CarbonDioxide",
"LiquidVaporPhaseBoundary", {"Temperature" ->
Quantity[t, "DegreesCelsius"]}];
obj[\[Mu]_?NumericQ] :=
QuantityMagnitude[
ThermodynamicData["CarbonDioxide",
"Entropy", {"Temperature" -> Quantity[t, "DegreesCelsius"],
"Pressure" -> \[Mu]*psat}]];
Limit[obj[m], m -> 1, Direction -> 1]]