function elems2007, p ;+ ; NAME: ; elems ; PURPOSE: (one line) ; Return orbital elements to pass to ephem ; DESCRIPTION: ; Return orbital elements to pass to ephem ; CATEGORY: ; Astronomy ; CALLING SEQUENCE: ; e = elems(p) ; INPUTS: ; p - planet number. 1 = mercury, etc. ; OPTIONAL INPUT PARAMETERS: ; none ; KEYWORD INPUT PARAMETERS: ; none ; KEYWORD OUTPUT PARAMETERS: ; none ; OUTPUTS: ; elem - vector of 7 elements ; e Eccentricity (degrees) ; i Inclination w.r.t xy-plane (degrees) ; LAN Longitude of Ascending Node (degrees) ; APF Argument of Perifocus (degrees) ; tau Julian Date of last periapsis ; n Mean motion (au/day) ; A semimajor axis. ; COMMON BLOCKS: ; None ; SIDE EFFECTS: ; RESTRICTIONS: ; None ; PROCEDURE: ; MODIFICATION HISTORY: ; Written 2000, by Leslie Young, SwRI ;- JDTTBD=2454101.500000000 ; = A.D. 2007-Jan-01 00:00:00.0000 (CT) case p of 1: begin EC= 2.056392449938306E-01 & QR= 3.074957814359186E-01 & IN= 2.855252549704066E+01 OM= 1.098578532395651E+01 & W = 6.757671484727892E+01 & Tp= 2.454141367679792E+06 N = 4.092343092329704E+00 & MA= 1.968477759955049E+02 & TA=-1.686185651770218E+02 A = 3.870984052246268E-01 & AD= 4.667010290133350E-01 & PR= 8.796916384521852E+01 end 2: begin EC= 6.776557216350436E-03 & QR= 7.184264488491972E-01 & IN= 2.443396698899500E+01 OM= 8.006577478148969E+00 & W = 1.239093809145851E+02 & Tp= 2.454209731815406E+06 N = 1.602144115493335E+00 & MA= 1.865970338382263E+02 & TA=-1.734914343833213E+02 A = 7.233281232626823E-01 & AD= 7.282297976761672E-01 & PR= 2.246988872715412E+02 end 3: begin EC= 1.670361885202152E-02 & QR= 9.832911829533065E-01 & IN= 2.343839018865727E+01 OM= 1.488741110702458E-04 & W = 1.029571077694497E+02 & Tp= 2.454104323523186E+06 N = 9.856169826206077E-01 & MA= 3.572170875972404E+02 & TA=-2.877822874913971E+00 A = 9.999947135016749E-01 & AD= 1.016698244050043E+00 & PR= 3.652534466713571E+02 end 4: begin EC= 9.337028465207686E-02 & QR= 1.381478526956870E+00 & IN= 2.467715211406103E+01 OM= 3.371230448954941E+00 & W = 3.329992583271292E+02 & Tp= 2.454255975314428E+06 N = 5.240018973952790E-01 & MA= 2.790546421392854E+02 & TA=-9.164634916379779E+01 A = 1.523751652488824E+00 & AD= 1.666024778020778E+00 & PR= 6.870204130738771E+02 end 5: begin EC= 4.848896937605290E-02 & QR= 4.950379421754955E+00 & IN= 2.323562753034352E+01 OM= 3.253376899252890E+00 & W = 1.271929661993569E+01 & Tp= 2.451320059060060E+06 N = 8.309488982554245E-02 & MA= 3.139189760025633E+02 & TA=-5.025505318278805E+01 A = 5.202650586729170E+00 & AD= 5.454921751703385E+00 & PR= 4.332396381484100E+03 ; EC= 4.893098616120922E-02 & QR= 4.947549130289248E+00 & IN= 2.323487974285201E+01 ; OM= 3.250916099996487E+00 & W = 1.169400819148276E+01 & Tp= 2.455640368127158E+06 ; N = 8.310825827427457E-02 & MA= 2.321073502380617E+02 & TA=-1.321532111065472E+02 ; A = 5.202092653948951E+00 & AD= 5.456636177608655E+00 & PR= 4.331699490223042E+03 end 6: begin EC= 5.427704645703578E-02 & QR= 9.029829559784201E+00 & IN= 2.255116475108028E+01 OM= 5.950726048622594E+00 & W = 8.797537943954090E+01 & Tp= 2.452837272592614E+06 N = 3.341124890185893E-02 & MA= 4.223941657672530E+01 & TA= 4.663675273934493E+01 A = 9.548070633112719E+00 & AD= 1.006631170644124E+01 & PR= 1.077481422671304E+04 end 7: begin EC= 4.764424259840767E-02 & QR= 1.827078788124309E+01 & IN= 2.366332843571104E+01 OM= 1.848854023276101E+00 & W = 1.709612677648794E+02 & Tp= 2.439574830355044E+06 N = 1.172941787724339E-02 & MA= 1.703893786303578E+02 & TA= 1.712501829189402E+02 A = 1.918483480489803E+01 & AD= 2.009888172855297E+01 & PR= 3.069206023416109E+04 end 8: begin EC= 7.100668586882949E-03 & QR= 2.990389549535517E+01 & IN= 2.229635100780807E+01 OM= 3.481801345063034E+00 & W = 2.799006594819591E+01 & Tp= 2.466049476861602E+06 N = 5.963228315789381E-03 & MA= 2.887514860624985E+02 & TA=-7.202118313380058E+01 A = 3.011775166853548E+01 & AD= 3.033160784171579E+01 & PR= 6.036998433328392E+04 end 9: begin EC= 2.543748319688926E-01 & QR= 2.967363028012835E+01 & IN= 2.343924587525028E+01 OM= 4.395974009773578E+01 & W = 1.839936675515462E+02 & Tp= 2.447899596660452E+06 N = 3.925806519630605E-03 & MA= 2.434747256451740E+01 & TA= 4.091371916521526E+01 A = 3.979698051030697E+01 & AD= 4.992033074048560E+01 & PR= 9.170090227316496E+04 end endcase elem = [EC, IN, OM, W, TP, N, A] return, elem end