;
; nim polarimeter retarder calibration test
;

dpath='/usr2/mpenn/02sep24/polarimetry/'

dark=swap_endian(rfits(dpath+'polar_dark.d.i.001'))
pol000=swap_endian(rfits(dpath+'polar_000.o.i.001'))-dark
pol002=swap_endian(rfits(dpath+'polar_002.o.i.001'))-dark
pol004=swap_endian(rfits(dpath+'polar_004.o.i.001'))-dark
pol006=swap_endian(rfits(dpath+'polar_006.o.i.001'))-dark
pol008=swap_endian(rfits(dpath+'polar_008.o.i.001'))-dark
pol010=swap_endian(rfits(dpath+'polar_010.o.i.001'))-dark
pol013=swap_endian(rfits(dpath+'polar_013.o.i.001'))-dark
pol016=swap_endian(rfits(dpath+'polar_016.o.i.001'))-dark
pol020=swap_endian(rfits(dpath+'polar_020.o.i.001'))-dark
pol030=swap_endian(rfits(dpath+'polar_030.o.i.001'))-dark
pol040=swap_endian(rfits(dpath+'polar_040.o.i.001'))-dark
pol050=swap_endian(rfits(dpath+'polar_050.o.i.001'))-dark
pol060=swap_endian(rfits(dpath+'polar_060.o.i.001'))-dark
pol080=swap_endian(rfits(dpath+'polar_080.o.i.001'))-dark
pol100=swap_endian(rfits(dpath+'polar_100.o.i.001'))-dark


x=[0.,.2,.4,.6,.8,1.0,1.3,1.6,2.0,3.,4.,5.,6.,8.,10]

int=[mean(pol000),mean(pol002),mean(pol004),mean(pol006),mean(pol008),mean(pol010),$
     mean(pol013),mean(pol016),mean(pol020),mean(pol030),mean(pol040),mean(pol050),$
     mean(pol060),mean(pol080),mean(pol100)]

plot,x,int,psym=1,xran=[0,10.2],xstyl=1,xtit='Volts',$
     chars=1.6,ytit='Intensity'
oplot,x,int

pol022=swap_endian(rfits(dpath+'polar_022.o.i.001'))-dark
pol024=swap_endian(rfits(dpath+'polar_024.o.i.001'))-dark
pol026=swap_endian(rfits(dpath+'polar_026.o.i.001'))-dark
pol028=swap_endian(rfits(dpath+'polar_028.o.i.001'))-dark
pol032=swap_endian(rfits(dpath+'polar_032.o.i.001'))-dark
pol034=swap_endian(rfits(dpath+'polar_034.o.i.001'))-dark
pol036=swap_endian(rfits(dpath+'polar_036.o.i.001'))-dark
pol038=swap_endian(rfits(dpath+'polar_038.o.i.001'))-dark


x2=[2.2,2.4,2.6,2.8,3.2,3.4,3.6,3.8]
y2=[mean(pol022),mean(pol024),mean(pol026),mean(pol028),$
    mean(pol032),mean(pol034),mean(pol036),mean(pol038)]

oplot,x2,y2,psym=2
oplot,x2,y2

the_end:
end