codes for the polytrope paper