;+
; NAME:
;	CvSky_Geographic
; PURPOSE:
;	Euler angles for Ecliptic to Geographic conversion
; CATEGORY:
;	smei/gen/idl/ephem
; CALLING SEQUENCE:
	FUNCTION CvSky_Geographic, UT, degrees=Degrees
; INPUTS:
;	UT	array[n]; type: time structure
; OPTIONAL INPUT PARAMETERS:
;	/degrees    if set all angles are in degrees (default: radians)
; OUTPUTS:
;	abc	array[3,n]; type: float
;		Euler angles
; INCLUDE:
	@compile_opt.pro				; On error, return to caller
; CALLS:
;	ToDegrees, SubArray, SuperArray, TimeGet
; PROCEDURE:
; MODIFICATION HISTORY:
;	AUG-2002, Paul Hick (UCSD/CASS; pphick@ucsd.edu)
;-

; 0h GMST on Jan 1, 1996, 17h17m24.7013s UT

;A = TimeGet(UT, /mjd, /scalar)-50084.220424783d0
A = TimeGet(UT, /njd, /scalar)+1460.779575217d0

B = 0.25d0+(A-!earth.siderealp*long(A/!earth.siderealp))/!earth.siderealp

A = SuperArray(0*A,3,/lead)

A = SubArray(A, element=0, add=-90d0)
A = SubArray(A, element=1, add=-23.43929d0)
A = SubArray(A, element=2, add=360d0*(B-long(B)))

A = A/ToDegrees(degrees=Degrees)

RETURN, A  &  END