proper polarized analytic beam

[bhazelton]

This code is from a very old branch of Zac's that is impossible to integrate now, but we may want this code in the future:

def analytic_beam_jones(nu, za, az, sigma=0.3):
    B = np.exp(-np.tan(za/2.)**2. / 2. / sigma**2.)
    J = np.array([[np.cos(az)*np.cos(az), -np.sin(az)],
                 [np.sin(az)*np.cos(za), np.cos(az)]])
     return B * J

[jpober]

@bhazelton I think this is just hanging around for future use? Or do you have near term plans for it?

[bhazelton]

I don't have near-term plans but we didn't want to lose this snippet, so we stuck it in here. I'm ok with closing, the issue will still exist.

[jpober]

This code is used by tests in pyradiosky. In test_skymodel and is called analytic_beam_jones(). Could be useful as an analytic beam option in pyuvsim.

[bhazelton]

This has been implemented in pyuvdata as part of moving the analytic beams to pyuvdata.