;+
; NAME:
;	CvSky_GSEQ
; PURPOSE:
;	Euler angles for Ecliptic to GSEQ conversion
; CATEGORY:
;	smei/gen/idl/ephem
; CALLING SEQUENCE:
	FUNCTION CvSky_GSEQ, UT, degrees=degrees
; INPUTS:
;	UT	array[n]; type: time structure
; OPTIONAL INPUT PARAMETERS:
;	/degrees    if set all angles are in degrees (default: radians)
; OUTPUTS:
;	Result		array[3,n]; type: float
;					Euler angles
; INCLUDE:
	@compile_opt.pro				; On error, return to caller
; CALLS:
;	big_eph, SubArray, SuperArray
; PROCEDURE:
;	Geocentric Solar Equatorial (GSEQ) coordinates:
;	x-axis points from Earth to Sun
;	y-axis is parallel to Sun's equatorial plane pointing
;		east of the Sun
;	z-axis is pointing north and completes righthanded
;		coordinate frame
;
;	The Sun's rotation axis lies in the x-z plane.
; MODIFICATION HISTORY:
;	OCT-2006, Paul Hick (UCSD/CASS; pphick@ucsd.edu)
;-
message, 'not finished'

A = big_eph(UT			, $					
	body	= jpl_body(/sun  ,/string)	, $
	center	= jpl_body(/earth,/string)	, $
	/to_ecliptic		, $
	/to_sphere			, $
	/precess			, $
	degrees	= degrees	, $
	/onebody			, $
	/silent)

A = SubArray(A,element=0)	; Geocentric ecliptic longitude of Sun

; If we zero out rows 2 & 3, we have the Euler Angles.
; Column 1 is longitude.

A = SuperArray(A,3,/lead)
A = SubArray(A,element=1,replace=0.0d0)
A = SubArray(A,element=2,replace=0.0d0)

RETURN, A  &  END