Added Psat as output, optimised polynomials through pass-by-reference (thanks Ian), wrote brine documentation