search

rodluger / starry

Tools for mapping stars and planets.
https://starry.readthedocs.io
MIT License
138 stars 32 forks source link

Issue with computing flux when total=False in starry oblate #295

Open shashankdholakia opened 2 years ago

shashankdholakia commented 2 years ago 
map33 = starry.Map(udeg=2, gdeg=2,oblate=True)
map33.omega=0.2
map33.beta=0.23
map33.wav=510
map33.tpole=7500
map33.f = (1-2/(0.2**2 + 2))
map33[1] = 0.2
map33[2] = 0.2

star = starry.Primary(map33, m=1, r=1)

planet = starry.kepler.Secondary(map=starry.Map(ydeg=1),
    m=0,  # mass in solar masses
    r=0.1,  # radius in solar radii
    porb=1,  # orbital period in days
)

system = starry.System(star, planet)
t = np.linspace(0,1,100)
system.flux(t,total=False)

This works when the oblate flag is set to False.