Polytrope definition was incorrect; fixed
Added piecewise polytropes; added option not to go to zero-enthalpy surface