;+
; NAME:
;	cv
; CALLING SEQUENCE:
	PRO cv
; INCLUDE:
	@compile_opt.pro
;-

print, txt_read('smei_sgp4_orbits.txt',d)

whatis, d
n = n_elements(d)

FOR i=0L,n-1 DO BEGIN
	d[i] = strmid(d[i],0,66)
	tt = strmid(d[i],8,18)
	tt = TimeSet(tt, format='YYYY_DOY_hhmmssmss')

	; Calculate the angle from the orbit normal to the Sun

	rsun = big_eph(tt,		$		; Direction to sun
		body	= 'sun'		, $
		center	= 'earth'	, $
		/to_sphere			, $
		/onebody			, $
		/to_equatorial		, $
		/degrees			, $
		/location			)

		; Direction orbit normal

	rnorm = sgp4_orbit_axis(tt,/degrees,/to_sphere,/precess)

	; Elongation orbit normal

	elo_orbnorm = sphere_distance(rsun,rnorm,/degrees)

	ra_sun		= ra_fictitious_sun(tt,/degrees)
	dra_orbnorm	= AngleRange(rnorm[0]-ra_sun,/degrees,/pi)
	dec_orbnorm	= rnorm[1]

	d[i] += string( elo_orbnorm, format='(F10.4)') +	$
			string( ra_sun     , format='(F10.4)') +	$
			string( dra_orbnorm, format='(F10.4)') +	$
			string( dec_orbnorm, format='(F10.4)')

	;print, d[i]
ENDFOR

openw, /get_lun, iu, 'smei_sgp4_normal.txt'
FOR i=0L,n-1 DO printf, iu, d[i]
free_lun, iu

RETURN  &  END