;; ;; Makes a mask to zero out everything out ;; of the good range of mk4 pixels. ;; Alice 1999 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) r0=1.0018 & dr=0.005035 & mk4rad=512./( 2.* (r0 + dr*(384.-1.))) r0=0.854772 & dr=0.00537989 ;;; True for new calibration feb 17, 2000 NumPix =384. GoodPix=NumPix-10. mk4rad =512./( 2.* (r0 + dr*(GoodPix-1.))) xcen = 255.5 & ycen = 255.5 & radius=1.00 rnder = 0.5005 dpx = fix ( cos(dp) * mk4rad * radius + xcen + rnder ) ; For disk size dpy = fix ( sin(dp) * mk4rad * radius + ycen + rnder ) ; For disk size mask = bytarr(512,512) mask(*,*) = 0 ;; ensure data outside the good range is zeroed for radius = 1.10 , 2.85 , 0.001 do begin dpx = fix ( cos(dp) * mk4rad * radius + xcen + rnder ) ; For disk size dpy = fix ( sin(dp) * mk4rad * radius + ycen + rnder ) ; For disk size mask (dpx,dpy) = 1 endfor tv,mask openw,1,'rawmk4_mask' writeu,1,mask close,1 end