Added backup saturation_T routine

It actually has exactly the same form as saturation_p. Interesting.  Possibility for optimization?

Signed-off-by: Ian Bell <ian.h.bell@gmail.com>