Change in variable name from dn to rho
Allow for saturation at surface in pressure calculation