rodluger / starry

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

Oblate maps unexpected amplitude #305

Open shashankdholakia opened 1 year ago

shashankdholakia commented 1 year ago

I'm running into a strange issue with amplitudes of gravity darkened maps that is leading to unexpected behavior when performing fits. Here is an attached MWE of the issue:

map = starry.Map(gdeg=4, oblate=True) #ydeg = 2*order_approx udeg=2

map.omega=0.225
map.beta = 0.08
map.wav=786.5 #TESS center bandpass
map.tpole=17300
map.f = 1-2/(0.225**2 + 2)

map.amp = 1.0
print(map.amp.eval())

The last line returns 6.112358459501673, I think it should return 1.0