pro printa,ac,ilog
common mpsf,nx,ny,np,psfm,wp,xp,yp,gp
common filen,fname,logfil
common spher,csph
common camer,camera,lam
o=np*2-2
keep =ac
if(csph ne 0.)then ac(7+o)=csph
ac=ac*lam/0.547
;
; Takes Zernike RMs coefficients and prints out a summary of peak to peak and
; RMs coefficients.
;
for is=0,np-1 do begin
   print,'x   =r*cos(t)    b2  ',$
      ac(0+2*is),1.8992573*ac(0+2*is)-0.673796*8.3345629*ac(3+o),$
      ' star number:',is
   print,'y   =r*sin(t)    b3  ',$
      ac(1+2*is),1.8992573*ac(1+2*is)-0.673796*8.3345629*ac(4+o)
endfor
print,' '
if(n_elements(ac+o) ge 8+o)then begin
   print,'fo  =r^2         b4 ',$
      ac(2+o),3.8874443*ac(2+o)-1.108900*16.895979*ac(7+o)
endif else begin
   print,'fo  =r^2         b4  ',ac(2+o ),3.8874443*ac(2+o)
endelse
print,'cx  =r^3*cos(t)  b7  ',ac(3+o ),ac(3+o)*8.3345629
print,'cy  =r^3*sin(t)  b8  ',ac(4+o ),ac(4+o)*8.3345629
if(n_elements(ac) ge 8+o) then begin 
   if(n_elements(ac) ge 12+o)then begin
      print,'ax  =r^2*cos(2t) b5  ',$
         ac(5+o ),2.3137662*ac(5+o)-0.750864*12.033645*ac(10+o)
      print,'ay  =r^2*sin(2t) b6  ',$
         ac(6+o ),2.3137662*ac(6+o)-0.750864*12.033645*ac(11+o)
   endif else begin
      print,'ax  =r^2*cos(2t) b5  ',ac(5+o ),2.3137662*ac(5+o)
      print,'ay  =r^2*sin(2t) b6  ',ac(6+o ),2.3137662*ac(6+o)
   endelse
   print,'sa  =r^4         b11 ',ac(7+o ),ac(7+o)*16.895979
endif
if(n_elements(ac) ge 10+o) then begin 
   print,'clx =r^3*cos(3t) b9  ',ac(8+o ),ac(8+o)*2.6701691 
   print,'cly =r^3*sin(3t) b10 ',ac(9+o ),ac(9+o)*2.6701691 
endif
if(n_elements(ac) ge 14+o) then begin 
   print,'sax =r^4*cos(2t) b12 ',ac(10+o),ac(10+o)*12.033645 
   print,'say =r^4*sin(2t) b13 ',ac(11+o),ac(11+o)*12.033645 
   print,'ashx=r^4*cos(4t) b14 ',ac(12+o),ac(12+o)*29.851527 
   print,'ashy=r^4*sin(4t) b15 ',ac(13+o),ac(13+o)*29.851527 
endif       
;
; Do logfile
;
if(ilog eq 1) then begin
for is=0,np-1 do begin
   printf,logfil,'x   =r*cos(t)    b2  ',$
      ac(0+2*is),1.8992573*ac(0+2*is)-0.673796*8.3345629*ac(3+o),$
      ' star number:',is
   printf,logfil,'y   =r*sin(t)    b3  ',$
      ac(1+2*is),1.8992573*ac(1+2*is)-0.673796*8.3345629*ac(4+o)
endfor
printf,logfil,' '
if(n_elements(ac+o) ge 8+o)then begin
   printf,logfil,'fo  =r^2         b4 ',$
      ac(2+o),3.8874443*ac(2+o)-1.108900*16.895979*ac(7+o)
endif else begin
   printf,logfil,'fo  =r^2         b4  ',ac(2+o ),3.8874443*ac(2+o)
endelse
printf,logfil,'cx  =r^3*cos(t)  b7  ',ac(3+o ),ac(3+o)*8.3345629
printf,logfil,'cy  =r^3*sin(t)  b8  ',ac(4+o ),ac(4+o)*8.3345629
if(n_elements(ac) ge 8+o) then begin 
   if(n_elements(ac) ge 12+o)then begin
      printf,logfil,'ax  =r^2*cos(2t) b5  ',$
         ac(5+o ),2.3137662*ac(5+o)-0.750864*12.033645*ac(10+o)
      printf,logfil,'ay  =r^2*sin(2t) b6  ',$
         ac(6+o ),2.3137662*ac(6+o)-0.750864*12.033645*ac(11+o)
   endif else begin
      printf,logfil,'ax  =r^2*cos(2t) b5  ',ac(5+o ),2.3137662*ac(5+o)
      printf,logfil,'ay  =r^2*sin(2t) b6  ',ac(6+o ),2.3137662*ac(6+o)
   endelse
   printf,logfil,'sa  =r^4         b11 ',ac(7+o ),ac(7+o)*16.895979
endif
if(n_elements(ac) ge 10+o) then begin 
   printf,logfil,'clx =r^3*cos(3t) b9  ',ac(8+o ),ac(8+o)*2.6701691 
   printf,logfil,'cly =r^3*sin(3t) b10 ',ac(9+o ),ac(9+o)*2.6701691 
endif
if(n_elements(ac) ge 14+o) then begin 
   printf,logfil,'sax =r^4*cos(2t) b12 ',ac(10+o),ac(10+o)*12.033645 
   printf,logfil,'say =r^4*sin(2t) b13 ',ac(11+o),ac(11+o)*12.033645 
   printf,logfil,'ashx=r^4*cos(4t) b14 ',ac(12+o),ac(12+o)*29.851527 
   printf,logfil,'ashy=r^4*sin(4t) b15 ',ac(13+o),ac(13+o)*29.851527 
endif       
endif
ac=keep
return
end