Fixed a bug inherited from SIR in Wittmann's equation of state option
The quantity dif is defined in a way such that if p and p1 differ by
orders of magnitude then dif-1 can be very small, signaling a premature
end of iteration