The reason I'm having trouble writing a working interpolation script for the Li & Han grids is that for color-magnitude isochrones, you can organize your data by mass, and plot from point to point along the curve of the isochrone. Then all you have to do is interpolate to smooth. But isochrone in Teff log(g) space double back on themselves such that there is no way to organize them along the length of the curve so that a point-to-point plot can be created. The curve it creates is not a real function. The only other way to do it would be to plot each point to its nearest neighbor, but that script would run quite slowly in python, I'm afraid.
No comments:
Post a Comment