;+
; cpstray.pro
; * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
; Determines SMM C/P stray light for radius r (arcsec units).
; * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
;
; SYNTAX:
;	s = cpstray (as)
;	as:	radial distance from sun center in arcsecond units.

; SUBROUTINES USED:
;	lint.pro
;
; HISTORY:
;	Andrew L. Stanger   SMM C/P   HAO/NCAR   1 March 1984
;	10 December 1985: Apply Munro calibration correction.
;	19 June 1990: Apply Munro correction factor ONLY if lcc = 1.
;	17 June 1998: Conversion from C to IDL.
;                     Use cc = 2.0 (Jupiter correction).
;	1 April 1999: Original rtab values were determined
;                     using 960 arcsec/Rsun.
;                     Created new radial distance table (astab) in arcsec units.
;                     Input distance is now in arcsec units.
; * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
;-

FUNCTION cpstray, as 		; as = distance in arcsec from sun center.

Cexp = 1.0e-11			; Scale factor for strtab table.
cc   = 2.0			; Jupiter calibration correction value.

;-- Radial values in units of R/Rsun corresponding to stray light table.

rtab = 								   $
     [ 0.00, 0.10, 0.20, 0.30, 0.40, 0.50, 0.60, 0.70, 0.80, 0.90, $
       1.00, 1.10, 1.20, 1.30, 1.40, 1.50, 1.60, 1.70, 1.80, 1.90, $
       2.00, 2.10, 2.20, 2.30, 2.40, 2.50, 2.60, 2.70, 2.80, 2.90, $
       3.00, 3.10, 3.20, 3.30, 3.40, 3.50, 3.60, 3.70, 3.80, 3.90, $
       4.00, 4.10, 4.20, 4.30, 4.40, 4.50, 4.60, 4.70, 4.80, 4.90, $
       5.00, 5.10, 5.20, 5.30, 5.40, 5.50, 5.60, 5.70, 5.80, 5.90, $
       6.00, 6.10, 6.20, 6.30, 6.40, 6.50, 6.60, 6.70, 6.80, 6.90, $
       7.00, 7.10, 7.20, 7.30, 7.40, 7.50, 7.60, 7.70, 7.80, 7.90, $
       8.00, 8.10, 8.20, 8.30, 8.40, 8.50, 8.60]

;-- Radial values in units of arcseconds corresponding to stray light table.

astab =									       $
[  0.0,   96.0,  192.0,  288.0,  384.0,  480.0,  576.0,  672.0,  768.0,  864.0,$
 960.0, 1056.0, 1152.0, 1248.0, 1344.0, 1440.0, 1536.0, 1632.0, 1728.0, 1824.0,$
1920.0, 2016.0, 2112.0, 2208.0, 2304.0, 2400.0, 2496.0, 2592.0, 2688.0, 2784.0,$
2880.0, 2976.0, 3072.0, 3168.0, 3264.0, 3360.0, 3456.0, 3552.0, 3648.0, 3744.0,$
3840.0, 3936.0, 4032.0, 4128.0, 4224.0, 4320.0, 4416.0, 4512.0, 4608.0, 4704.0,$
4800.0, 2896.0, 4992.0, 5088.0, 5184.0, 5280.0, 5376.0, 5472.0, 5568.0, 5664.0,$
5760.0, 5856.0, 5952.0, 6048.0, 6144.0, 6240.0, 6336.0, 6432.0, 6528.0, 6624.0,$
6720.0, 6816.0, 6912.0, 7008.0, 7104.0, 7200.0, 7296.0, 7392.0, 7488.0, 7584.0,$
7680.0, 7776.0, 7872.0, 7968.0, 8064.0, 8160.0, 8256.0]

; Stray light values in radiance units (10**-11 B/Bsun)
; in increments of 0.1 solar radius, beginning with 0.0 Rsun and
; continuing out to 8.4 Rsun.
;
; stray light estimate : William Wagner 24 July 1982 graph (smstr2).
;

strtab =							 	    $
     [ 0.00,  0.00,  0.00,  0.00,  0.00,  0.00,  0.00,  0.00,  0.00,  0.00, $
       0.00,  0.00,  0.00,  0.00,  0.00, 90.00, 190.0, 110.0, 67.00, 45.00, $
      36.00, 29.00, 25.50, 23.00, 21.50, 20.50, 20.20, 20.10, 20.03, 20.00, $
      20.03, 20.10, 20.20, 20.50, 21.00, 23.00, 28.50, 29.50, 29.80, 29.90, $
      30.00, 30.05, 30.00, 30.00, 29.00, 27.50, 25.50, 23.50, 21.50, 20.70, $
      20.00, 20.00, 20.00, 20.00, 20.00, 20.00, 20.00, 20.00, 20.00, 20.00, $
      20.00, 20.00, 20.00, 20.00, 20.00, 20.00, 20.00, 20.00, 20.00, 20.00, $
      20.00, 20.00, 20.00, 20.00, 20.00, 20.00, 20.00, 20.00, 20.00, 20.00, $
      20.00, 20.00, 20.00, 20.00, 20.00, 20.00, 20.00]

;-------------------------------------------------------------------------------
   size_strtab = SIZE (strtab)
   ntab = size_strtab (1)

   s = lint (as, astab, strtab, ntab) * Cexp * cc

   RETURN, s

END