Oblate maps unexpected amplitude

[shashankdholakia]

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