;+
; vignet.pro
; * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
; Determines SMM C/P vignetting correction for radius r (arcsec units).
; * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
; SYNTAX:
;	v = cpvignet (as)
;	as:	radial distance from sun center in arcsecond units.
;
; SUBROUTINES USED:
;	lint.pro
;
; HISTORY:
;	Andrew L. Stanger   HAO/NCAR   SMM C/P   1 March 1984
;	17 June 1998: Conversion from C to IDL.
;	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 cpvignet, as		; as = distance in arcsec from sun center.

;--- 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 vignetting 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]

;--- Vignetting ratios as a fraction of transmission from 0.0 Rsun
;    out to 8.4 Rsun in increments of 0.1 solar radius.
;
;    vignetting estimate :  Dale Gary   7 July 1982   (smvig2).

vigtab =								    $
     [1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, $
      1.000, 1.000, 1.000, 1.000, 1.000, 35.34, 19.31, 12.53, 9.083, 7.003, $
      5.656, 4.706, 4.034, 3.519, 3.113, 2.789, 2.536, 2.327, 2.160, 2.017, $
      1.895, 1.787, 1.685, 1.590, 1.498, 1.419, 1.347, 1.281, 1.222, 1.174, $
      1.135, 1.099, 1.069, 1.044, 1.023, 1.009, 1.001, 1.000, 1.000, 1.000, $
      1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, $
      1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, $
      1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000, $
      1.000, 1.000, 1.000, 1.000, 1.000, 1.000, 1.000]

;-------------------------------------------------------------------------------
   size_vigtab = SIZE (vigtab)
   ntab = size_vigtab (1)

   v = lint (as, astab, vigtab, ntab)

   RETURN, v
END