nray=720*5 nraym1=nray-1 acirc = !PI * 2.0 / float(nray) dp = findgen(nray) * acirc + acirc * (0.25 * (float(nray))) dpx = intarr(nray) dpy = intarr(nray) xcen = 255.5 & ycen = 255.5 & mk3rad = 104.4435 & radius=1.00 rnder = 0.5005 dpx = fix ( cos(dp) * mk3rad * radius + xcen + rnder ) ; For disk size dpy = fix ( sin(dp) * mk3rad * radius + ycen + rnder ) ; For disk size mask = bytarr(512,512) mask(*,*) = 0 for radius = 1.1224, 2.45 , 0.001 do begin dpx = fix ( cos(dp) * mk3rad * radius + xcen + rnder ) ; For disk size dpy = fix ( sin(dp) * mk3rad * radius + ycen + rnder ) ; For disk size mask (dpx,dpy) = 1 endfor tv,mask openw,1,'rawmk3_mask' writeu,1,mask close,1 end