______________________________________________________________________________
RF
REFCENT is a module which corrects a two-dimensional spectrum for the
effects of differential refraction in the atmosphere. Each column of data is
shifted vertically by means of a tapered sinc interpolation. The amount of
shift is calculated from the zenith distance of the observation, the
orientation of the slit relative to the horizon, and a model for the
refractive index (7th order fit to parameters in Allen's "Astrophysical
Quantities"). Because this module requires a lot of cpu time, it is best
to run it as a batch job.
______________________________________________________________________________
Enter an option: 21-MAR-1989 11:32:06.52
OPTION=RF
REFCENT, Version 4
Open the file of spectral data for which atmospheric dispersion effects
are to be removed:
Enter the filename: (<CR> for default)
FILE=LONG
There is one map in the file
Reading from the file
10 rows of 576 pixels in R4 format
Enter the effective wavelength at which the guiding was done:
(units of Angstroms)
WAVELENGTH=5700
Enter the scale along the slit: (arc seconds/pixel)
SCALE=.63
Enter the position angle of the slit on the sky: (degrees E of N)
POSITION_ANGLE=15
Is the horizon toward the low row side of the slit? (default = Y)
LOW=
T
Enter a choice:
H = Latitude, declination, and hour angle
from SAD header information
M = Manual input of above quantities (default)
CHOICE=M
Enter the latitude in degrees, minutes, and seconds:
LATITUDE=32 0 0
Enter the declination in degrees, minutes, and seconds:
DECLINATION=47 22 01
Enter the hour angle in hours, minutes, and seconds: (East < 0)
HOUR_ANGLE=-2 03 40
Open a file for the shifted data
NMAPS= 1 NPLNS= 1 NROWS= 10 NELS= 576
Enter the filename: (<CR> for default)
FILE=LONGNEW
Creating new file
Do you wish to copy the header from the input to output file? (default = Y)
SAME=
T
Maximum shift was: -2.74 pixels at 4500.00 Angstroms
If the file contains several maps, after doing the first
map, use option DO, e.g., OPTION=DO($=2:11).RF(MAPI=$,MAPO=$)
______________________________________________________________________________