More work on Rho-I saturation solver

Worked on (and completed most) of the rho-i RRND solver with several
changes to other various functions used.

Had to readjust the Helmholtz coefficients for zero-ed internal energy
and entropy at the triple point.

Clausius-Clapeyron functions now abide by non-dimensional naming
standards.