NAME:
	TMO_skymotion
 PURPOSE:
	Corrects positions on CCD for skymotion during transit
 CATEGORY:
	SMEI prototype camera: TMO
 CALLING SEQUENCE:
	rr_out = TMO_skymotion(rr [, p_axis, frame, /degrees, to_zenith=to_zenith,  $
	    fov=fov, vv=vv, plotx=plotx
 INPUTS:
	rr		array[2,n]; type: float
			    positions on ccd as pairs off azimuthal angle (along long
			    dimension of fov), and radius
	p_axis		array[2]; type: float; default: [0,0]
			    position of optical axis on ccd: azimuthal angle and radius
			    Currently only p_axis[0] is used.
	frame		array[m]; type: any; default: 1
			    list of frame numbers (see PROCEDURE)
 OPTIONAL INPUT PARAMETERS:
	/south		if set, the south position (19.5 south of zenith) is assumed
	/north		if set, the north position (10.0 north of zenith) is assumed
	/degrees	if set, all angles are in degrees (default: radians)
	/plotx		plot results
 OUTPUTS:
	rr_out		array[2,n,m]; type: float
			    positions rr, corrected for motion across the sky
 OPTIONAL OUTPUT PARAMETERS:
	fov=fov 	array[*]; CCD azimuth angle (-30,..,+30) of points along
			    long dimension of fov
	vv=vv		array[*]; type: float
			    skymotion at positions 'fov' in pixels/minute
 INCLUDE:
	@compile_opt.pro	; On error, return to caller
 CALLED BY:
	TMO_tracksky
 RESTRICTIONS:
	No check is made whether the locations rr are actually inside the fov
 CALLS: ***
	BadValue, INTERPOL, InitVar, MEAN, SubArray, SuperArray, SyncDims, ToRadians, gridgen
 PROCEDURE:
 >	The SMEI fov is assumed to be along the local meridian
 > The motion of stars during transit are calculated in pixels/minute
	under the assumption that the 3 degree fov is covered by 27.16 pixels, and
	that the geographic latitude of TMO is 34.382 degree.
 > The 'frame' array is used to define the time axis in units of minutes.
	(TMO data were taken at a 1 minute cadence). Usually 'frame' will be an
	integer array of frame numbers, but the effective cadence time can be adjusted
	by multiplication of the 'frame' array with the cadence time in minutes.
 >	For each location 'rr' on the CCD inside the fov, the corrected position is
	the position of 'rr' after a time 'frame' minutes has elapsed.
 >	Only the radial component r[1,*] is corrected, i.e. stars are assumed to
	transit on the CCD in the radial direction.
 MODIFICATION HISTORY:
	JUL-2001, Paul Hick (UCSD/CASS; pphick@ucsd.edu)

 NAME:
	TMO_tracksky
 PURPOSE:

 CATEGORY:
	SMEI prototype camera: TMO
 CALLING SEQUENCE:
 INPUTS:
	box		array[2,2]; type: integer; default: defined in smei_camera
			    defines two corners of box in the form [ [x1,y1], [x2,y2] ]
			    relative to the lowerleft corner of a full image
			    The box contains nx by ny pixels (nx=x2-x1+1, ny=y2-y1+1)
 OPTIONAL INPUT PARAMETERS:
	center=center	array[2]; type: float; default: defined in smei_camera
			    'center' must be defined relative to the same origin as 'box'
	p_axis=p_axis	array[2]; type: float; default: defined in smei_camera
	frame=frame	array[m]; type: integer; default: 1
	/plotx
 OUTPUTS:
	rfov		array[2,nx,ny,m]; type: float
			    rectangular coordinates of all pixels in box corrected for skymotion
			    relative to the same origin as 'box' and 'center'
 OPTIONAL OUTPUT PARAMETERS:
	pfov		array[2,nx,ny,m]; type: float
			    polar coordinates of all pixels in box corrected for skymotion
 INCLUDE:
	@compile_opt.pro	; On error, return to caller
 CALLS: ***
	BadValue, CV_COORD, Inside_Wedge, SuperArray, TMO_skymotion, gridgen, smei_camera
	view
 CALLED BY:
	qView_TMO_tracksky
 PROCEDURE:
	For each of the pixels in 'box' the position in each frame is calculated after
	correcting it for skymotion.
 MODIFICATION HISTORY:
	JULY-2001, Paul Hick (UCSD/CASS; pphick@ucsd.edu)