Python toolkit for calculating stellar limb darkening profiles and model-specific coefficients using the stellar atmosphere spectrum library by Husser et al. (2013). Described in Parviainen & Aigrain, MNRAS 453, 3821–3826 (2015).
Hi! After a recent update of the cached files, I get inconsistent results retrieving the non-linear limb-darkening law (I removed the cache and downloaded again). I'm comparing the intensity difference of the profiles for a given temperature (4250K, 5800K, 7500K). A few months ago the cumulative difference was at most 2.5% and now is more than 10x larger. Funny enough for the third temperature, nothing changed. I'm using python 3.8 to run the code:
Hi! After a recent update of the cached files, I get inconsistent results retrieving the non-linear limb-darkening law (I removed the cache and downloaded again). I'm comparing the intensity difference of the profiles for a given temperature (4250K, 5800K, 7500K). A few months ago the cumulative difference was at most 2.5% and now is more than 10x larger. Funny enough for the third temperature, nothing changed. I'm using python 3.8 to run the code:
filters=[BoxcarFilter("w",380,800)] sc = LDPSetCreator(teff=(parvar[k], sig_steff),logg=(slogg, sig_slogg),z=(sz, sig_sz),filters=filters, cache=os.path.abspath(os.path.join("yourpath" ,os.pardir))+"/cache/") ps = sc.create_profiles(nsamples=1000) ldcn,qe = ps.coeffs_qd(do_mc=True) ldcl,qel= ps.coeffs_ln(do_mc=True) ldcnl,qenl= ps.coeffs_nl(do_mc=True)
Limb-darkening coefficients:
4250K [[0.9454834]] [[0.64804808 0.3428646 ]] [[-0.70811569 5.80458496 -7.32818848 3.22719344]]
5800K [[0.97699707]] [[0.31915583 0.68739933]] [[-1.02702287 7.47210631 -9.45970715 3.98864554]]
7500K [[0.58144531]] [[0.50178719 0.12041476]] [[-0.11783424 1.00176812 -0.1979341 -0.0744706 ]]
Best regards, Eduardo Cristo.