Refactored p(T,rho) to use delta instead of rho

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