;+
; NAME		sunray.pro
;
; PURPOSE	Draw radial rays.
;
; SYNTAX	sunray, rmin, rmax, anginc, xcen, ycen, scroll, pixrs, cindex, $
;                       xmin, xmax, ymin, ymax
;		rmin, rmax:	start/stop radii [Rsun]
;		rinc:		increment for radius.
;		anginc:		angular increment [degrees].
;		dotsize:	dot size [1, 3, 5, 7, ...] [pixels]
;		xcen, ycen:	sun center.
;		pixrs:		#pixels/Rsun.
;		scroll:		spacecraft roll.
;		cindex:		color index.
;		xmin, ymin:	Lower left  corner of window.
;		xmax, ymax:	Upper right corner of window.
;
; AUTHOR	Andrew L. Stanger   HAO/NCAR   10 Jan 2006
;		Adapted from "glcp_suncir.c".
;-

PRO sunray, rmin, rmax, rinc, anginc, dotsize, $
	    xcen, ycen, pixrs, scroll, cindex, $
	    xmin, xmax, ymin, ymax

   nres = 1			; Low resolution only.
   cproll = 43.0		; C/P axis is 43 deg CCW from +Y axis.
   p = FIX (dotsize) / 2

   ;*** Draw radial scans (r=r1min -> r1max) every ang1inc degrees. ***

   FOR th = 0.0, 360.0, anginc DO	$
   BEGIN ;{
;      IF (th EQ 180.0) THEN rmin = 0.2 ELSE rmin = 0.4

      FOR r = rmin, rmax, rinc DO	$
      BEGIN ;{
	 ierr = rcoord (r, th, xx, yy, 1, scroll+cproll, xcen, ycen, pixrs)
	 xg = ((xx / nres) + 0.5) + xmin
	 yg = ((yy / nres) + 0.5) + ymin

	 FOR yp = yg - p, yg + p DO		$
	 BEGIN ;{
	 IF (xg-p GE xmin AND xg+p LE xmax  AND		$
	     yp   GE ymin AND yp   LE ymax) THEN	$
	     BEGIN ;{
		wvec, xg-p, yp, xg+p, yp, cindex
	     END   ;}
         END   ;}

      END   ;}

   END   ;}

RETURN
END