Automatically estimate zmax if zmax=None.  Estimated as twice the scale height at r_max