;+ ; 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