pro seta,ac
common mpsf,nx,ny,np,psfm,wp,xp,yp,gp
common camer,camera,lam
;case 2
a=[-2.92788,3.20713,6.06396,-2.99452,3.77422,0.113594,0.602555,-4.38223]
; case 3
a=[-2.80042,2.29146,6.30618,-1.81163,3.76222,0.876647,1.48907,-5.04596,$
 -0.309251,-4.998904e-02,-1.02750,-1.65498,-2.624347e-02,-5.816614e-02]
; case 4
a=[-3.37238,1.86156,-6.99383,-1.49684,5.65420,-1.15969,-1.40797,6.78794,$
-0.380177,-4.747110e-02,0.923253,1.50126,2.724005e-02, 8.259022e-02]
; case 1
a=[0.538994,0.704896,2.79230,-0.230412,0.580961]
;
; Convert peak to peak to RMS zernike coefficients.
;
convert,a,ac
;ac=[0.489898,0.289400,0.718287,0.250000,-2.764536e-02]
; Perkin Elmer's fit
;ac=[-2.6453,0.08973,0.52430,-0.10045,0.0789783,-0.029463,-0.1011,0.1038]
;CASE 4
;wavefront min,max     -1.23589      1.87896
;f=sqrt(total((psfc-psfm)^2)) 1.092274e-02
ac=[-0.203860,0.253063,0.122993,0.281560,-7.085917e-02,-2.482805e-02,$
-7.304214e-02,0.560844,-2.480943e-02,-1.245739e-02,-2.388870e-02,$
2.130217e-02,-8.261522e-04, 2.061806e-03]
;Case 3
;wavefront min,max     -3.40792      1.46799
;f=sqrt(total((psfc-psfm)^2)) 1.075737e-02
;               
;CASE2
ac=[-0.20411096,0.27579070,7.53120926e-02,0.28348956, -5.80129640e-02,$
2.85533387e-02,-9.29477786e-02,0.55248076]
ac=[-0.214913,0.393141,-4.605331e-02,0.455995,-0.154643,-7.156378e-02,$
   -0.161557,-0.503294,7.402786e-02,8.667597e-02,-5.219442e-02,$
    -0.243497, 4.446048e-03, 4.697140e-04]
;ac=[-0.18247018,0.4084,
ac=[0.,0.,-9.1725e-02,0.247,-0.2200,-3.8962e-02,-0.1475,-0.4905]
;clx,cly      0.20518609  4.11406771e-02
;sax,say,ashx,ashy -8.80002056e-02     -0.18024235 -1.15860552e-02
;  1.95001895e-03
ac=[0.,0.,-0.111,-0.0484,-0.295,0.0442,0.108,-0.517] 
ac= [0.0,0.0,-0.162158, 0.073461 ,-0.195927 , 0.0147374, 0.06178 ,-0.511532]
AC(2)=2.3
; fit to outer ring
ac=[-0.0,0.0,2.7253,0.19016,-0.0123893,0.0190495,-0.0221012,-0.50772]
ac(2)=0.
;AC(2)=AC(2)-0.3*6.0643/0.547
;ac=[0.,0.,5.093902e-02,-3.690926e-02,-2.349226e-02,0.119262,-0.127774,-0.656758]
;FOC first light mosaic
ac=[0.01732279,-0.01171299,0.01202544,-0.01431259,0.01679647,-0.02462909,$
    0.00578140,-0.00548071,0.01443863,-0.01414735,0.01532835,-0.03182312,$
    0.05480345,-0.01544212,-0.0374712,0.09616800 ,-0.100272 ,-0.695656]
; FOC TOT
ac=[0.,0.,0.05480345,-0.01544212,-0.0374712,0.09616800 ,-0.100272 ,-0.5]
ac=ac-ac
;ac=[0.01732279,-0.01171299,$
;    0.05480345,-0.01544212,-0.0374712,0.09616800 ,-0.100272 ,-0.695656]
;Hughes
;ac=[0.,0.,0.3551,-0.0138,-0.069,0.038,-0.037,0.143]
;ac=[-0.228207, 1.936953e-02,-0.269009, 7.902537e-02,5.275906e-02,-0.252481,$
;-7.494647e-02,0.138270, 5.636291e-02,-7.605404e-02,-9.896935e-02,-0.5]
;ac=[0.,0.,0.,0.,0.,0., 2.72538,0.190166,-1.238933e-02,$
;    1.904951e-02, -2.210126e-02, -0.507721]
;ac(6)=ac(6)-0.3*6.06/0.547     
;FOc Second light - first fit (128x128 no spiders)
ac=[7.96141250e-02,0.12260372,-0.16995673,-1.59251541e-03,-8.86293248e-03,$
3.72349790e-02,5.76951428e-02,-0.36726184]
ac=[6.092761e-02,0.140581,-0.546648,-1.739560e-02,1.229592e-02,1.444851e-02,$
0.117893,-0.636687]
ac=(550./170.)*[0.,0.,-0.9,0.,0.,0.,0.,-0.5]
ac=[-5.127561e-02,-0.181187,-2.81419, 3.768681e-02,-9.022825e-02,9.378795e-02,$
     0.284081,-1.95970]
;F288 at 170nm
ac=[7.10385240e-02,-0.10325683,-2.5467824,0.96858429,-0.30889215,$
    1.32183229e-02,0.11858188 ,-2.1387963]*170./210.
;f/96 at 210
ac=[0.285155,0.126159,-2.44335,0.112262,-0.275273,1.346612e-02,$
    2.913026e-02,-1.82885]
ac=ac-ac &ac(7)=-0.5*547./358. &ac(2)=-0.9*547./358.
ac=[-5.441917e-02,-2.817061e-02,-1.61515, 5.439053e-02,$
    -2.001082e-02,-9.495627e-02,-8.306878e-02,    -0.688733]
ac=[0.,0.,-1.,0.,0.,0.,0.,-0.5]*0.547/lam
temp=fltarr(n_elements(ac)+np*2-2)
for ip=0,np-1 do begin
    temp(0+2*ip)=ac(0)
    temp(1+2*ip)=ac(1) 
endfor
temp(np*2)=ac(2:*)
ac=temp
return
end