re-use functions
what happens to derivatives along the saturation line at the critical point??
what happens to derivatives along the saturation line at the critical point??
also output dTp and dpT.

notice: the derivatives seem to be wrong at the critical point
also output dTp and dpT
what happens to derivatives along the saturation line at the critical point??