amaurea / enlib

5 stars 12 forks source link

error in aberration if modulation == False #43

Closed ajvanengelen closed 6 years ago

ajvanengelen commented 6 years ago

Hi, I am trying to run the aberration routine with modulation turned off, and am finding the following:

/global/u2/e/engelen/enlib/aberration.py in aberrate(imap, dir, beta, mode, order, recenter, modulation) 50 # import pdb 51 # pdb.set_trace() ---> 52 omap *= pos[2+pol,None] 53 return omap 54

IndexError: index 3 is out of bounds for axis 0 with size 3

amaurea commented 6 years ago

I've added an argument to make modulation optional now. Does this help?

ajvanengelen commented 6 years ago

Looks like it worked! Not sure why it didn't work when I tried to do the same thing, but am glad this is working now. Thanks!