More comments and updates
edge cases not workign; some tabled EOS NaNat surface, need to find polytrope to compare