Calculation of ΔSsurr at constant P