Changed from numpy matrices to arrays.

Bug in PREOS.calcLnPhi() failed to return all the results (premature return)

Bug in calcZ() returned all the real roots of the EOS, including the non-physical middle root (when 3 real solutions exist). This caused problems with the saturation line detection.