;
;
xrange = [-60,60]
yrange = xrange
xc = cos(2*!pi*findgen(500)/499.)*32
yc = sin(2*!pi*findgen(500)/499.)*32
old_device=!d.(0)
set_plot,'null
plot, xc, yc, xrange=xrange, yrange=yrange
set_plot,old_device
mk_square2, /qdebug
;
for im=0,2 do begin
    rarr = [24, 12]
    case im of
	0: begin & tit='No adjustment' & end
	1: begin & tit='M1 adjustment (69.9/100)' & rarr=rarr*69.9/100 & end
	2: begin & tit='M2 adjustment (121/100)' & rarr=rarr*121./100 & end
    endcase
    plot, xc, yc, xrange=xrange, yrange=xrange, tit=tit
    for irad=0,n_elements(rarr)-1 do begin
	for i=1,8 do begin
	    x0 = cos(!pi/180*360*i/8.)*rarr(irad)
	    y0 = sin(!pi/180*360*i/8.)*rarr(irad)
	    oplot, xc+x0, yc+y0, linestyle=irad
	end
    end
    pause
end
;
end