Add tests to stop iterating when converged to roundoff in dz_to_thickness_EOS
when in fully non-Boussinesq mode.  Also modified the code to keep track of the
layer thickness directly in non-Boussinesq mode rather than setting the layer
thickness based on a difference between interface pressures.  Boussinesq
answers are bitwise identical, but non-Boussinesq answers change.