;+
; NAME:
;	CvSky_GSE
; PURPOSE:
;	Euler angles for Ecliptic to GSE conversion
; CATEGORY:
;	smei/gen/idl/ephem
; CALLING SEQUENCE:
	FUNCTION CvSky_GSE, 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 Ecliptic (GSE) coordinates:
;	x-axis points from Earth to Sun
;	y-axis is in the ecliptic plane opposing orbital motion
;	z-axis  points to ecliptic north
; MODIFICATION HISTORY:
;	AUG-2002, Paul Hick (UCSD/CASS)
;	OCT-2006, Pual Hick (UCSD/CASS; pphick@ucsd.edu)
;		Replaced NewcombSun by big_eph
;-

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

; Column 1 now contains the geocentric ecliptic longitude of the Sun
; If we zero out columns 2 & 3, we have the Euler Angles.

A = SubArray(A,element=1,replace=0.0d0)
A = SubArray(A,element=2,replace=0.0d0)

RETURN, A  &  END