I implemented some changes to the isochrone.py file which now instead of calculating a lower, nominal, and upper value when calculating each parameters, creates a normal distribution from the given x values in the evaluate function. Then, it takes 1,000 samples from which then I pass through the interpolate function in a for loop. Then, I take all these values and adopt the average and standard deviation as the value and uncertainty on the parameter which is then provided to the results. I made these changes because in making more SEDs, there have been high uncertainties calculated with the previous code which implied negative values of certain parameters (like mass, for example).
Also, I added a .txt file containing the ATMO2020 evolutionary track models using Strong non-equilibrium chemistry for T-Y transition dwarfs and made it an option to choose from in the model dictionary in the isochrone.py file since the hybrid-solar age stops working super well for these cooler objects .
The modifications work with my installation of sedkit but I haven't been able to get pytest working in my dev version. Curious to see if the tests pass in the repo! Let me know what you think.
isochrone.py
file which now instead of calculating a lower, nominal, and upper value when calculating each parameters, creates a normal distribution from the given x values in the evaluate function. Then, it takes 1,000 samples from which then I pass through the interpolate function in a for loop. Then, I take all these values and adopt the average and standard deviation as the value and uncertainty on the parameter which is then provided to the results. I made these changes because in making more SEDs, there have been high uncertainties calculated with the previous code which implied negative values of certain parameters (like mass, for example).isochrone.py
file since the hybrid-solar age stops working super well for these cooler objects .The modifications work with my installation of sedkit but I haven't been able to get pytest working in my dev version. Curious to see if the tests pass in the repo! Let me know what you think.