Hi Jan,

Here are the calibration matrices as promised. They are contained in the idl
save file attached here called 20110311_cal.sav. If you restore this file in idl, you will
have three variables, o_matrix[4,4,9], d_matrix[4,4,9] and regions[9]. Regions contains
the 9 wavelengths. The o and d matrix variables contain the 4 by 4 matrices for the 9
wavelength regions. The d_matrix is the one to use to demodulate the polarization states.
For example, if you print out regions you will get:

IDL> print,regions
 530.300      587.600      637.500      656.300      789.200 854.200      1074.70      1079.80      1083.00

and you can see that h-alpha is the fourth one (index 3).
The you can print d_matrix[*,*,3] to get the h-alpha demodulation matrix and see:

IDL> print,d_matrix[*,*,3]
     0.321168    0.0297936     0.402418     0.246620
    -0.488591     0.495348     0.407239    -0.413996
     0.567988   -0.0557511     0.205606    -0.717844
    0.0487035    -0.837206     0.509547     0.278956

You you will obtain I,Q,U,V by performing a matrix multiply of the 4 measured intensities with the
d_matrix as above.

Please let me know if you have any questions. I am very curious to see what you get when you
apply these. The matrices at 5303 and 10830 might have some problems. We will need to repeat
these calibrations at the telescope when things settle down. I have also attached the idl routines
used to compute these. The comps_cal is the main program. Eventually these routines need to
be modified to read in the fits files from cals obtained at the telescope.

Best Wishes,       Steve

Steven Tomczyk
Assistant Director
High Altitude Observatory (HAO)
National Center for Atmospheric Research (NCAR)
3090 Center Green Dr.               303-497-1579
Boulder, CO  80301                    tomczyk@ucar.edu  
