SUBROUTINE ABERRAT_XYZ(JDE,XYZ) C C FOR A JULIAN EPHEMERIS DATE, RETURN THE RECTANGULAR COMPONENTS C OF THE ANNUAL ABERRATION CORRECTION. C C B. KNAPP, 1993-06-18, 2000-09-20, 2001-04-24 (i*4) C C USES THE METHOD OF C. RON AND J. VONDRAK ["EXPANSION OF ANNUAL C ABERRATION INTO TRIGONMETRIC SERIES", BULL. ASTRON. INST. C CZECH., VOL. 37 (1986), PP. 96-103.], AS PRESENTED IN J. MEEUS, C ASTRONOMICAL ALGORITHMS, WILLMAN-BELL, 1991. C C C RCS DATA C C $Header$ C C $Log$ C C IMPLICIT NONE C C INPUT REAL*8 JDE C C OUTPUT REAL*8 XYZ(3) C C MEAN LONGITUDES OF PERTURBING SOURCES AS A FUNCTION OF TIME REAL*8 SOURCE(2,11),MLON(11) C C ARGUMENT COMBINATION COEFFICIENTS REAL*8 COEF(11,36) C C SINE AND COSINE COEFFICIENTS FOR RECTANGULAR VELOCITY C COMPONENTS DX, DY, DZ (1E-8 AU/DAY) REAL*8 SDX(2,36) REAL*8 CDX(2,36) REAL*8 SDY(2,36) REAL*8 CDY(2,36) REAL*8 SDZ(2,36) REAL*8 CDZ(2,36) C C VARIABLES REAL*8 T, ARG, SA, CA, DX, DY, DZ INTEGER*4 I, J C DATA SOURCE & / 3.1761467D0, 1021.3285546D0, !L2 & 1.7534703D0, 628.3075849D0, !L3 & 6.2034809D0, 334.0612431D0, !L4 & 0.5995465D0, 52.9690965D0, !L5 & 0.8740168D0, 21.3299095D0, !L6 & 5.4812939D0, 7.4781599D0, !L7 & 5.3118863D0, 3.8133036D0, !L8 & 3.8103444D0, 8399.6847337D0, !L' (W) & 5.1984667D0, 7771.3771486D0, !D & 2.3555559D0, 8328.6914289D0, !M' (L) & 1.6279052D0, 8433.4661601D0 !F & / C DATA COEF ! L2 L3 L4 L5 L6 L7 L8 L' D M' F No. & / 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, ! 1 & 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, ! 2 & 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, ! 3 & 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, ! 4 & 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, ! 5 & 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, ! 6 & 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, ! 7 & 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, ! 8 & 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, ! 9 & 0, 2, 0, -1, 0, 0, 0, 0, 0, 0, 0, !10 & 0, 3, -8, 3, 0, 0, 0, 0, 0, 0, 0, !11 & 0, 5, -8, 3, 0, 0, 0, 0, 0, 0, 0, !12 & 2, -1, 0, 0, 0, 0, 0, 0, 0, 0, 0, !13 & 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, !14 & 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, !15 & 0, 1, 0, -2, 0, 0, 0, 0, 0, 0, 0, !16 & 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, !17 & 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, !18 & 2, -2, 0, 0, 0, 0, 0, 0, 0, 0, 0, !19 & 0, 1, 0, -1, 0, 0, 0, 0, 0, 0, 0, !20 & 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, !21 & 0, 3, 0, -2, 0, 0, 0, 0, 0, 0, 0, !22 & 1, -2, 0, 0, 0, 0, 0, 0, 0, 0, 0, !23 & 2, -3, 0, 0, 0, 0, 0, 0, 0, 0, 0, !24 & 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, !25 & 2, -4, 0, 0, 0, 0, 0, 0, 0, 0, 0, !26 & 0, 3, -2, 0, 0, 0, 0, 0, 0, 0, 0, !27 & 0, 0, 0, 0, 0, 0, 0, 1, 2, -1, 0, !28 & 8,-12, 0, 0, 0, 0, 0, 0, 0, 0, 0, !29 & 8,-14, 0, 0, 0, 0, 0, 0, 0, 0, 0, !30 & 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, !31 & 3, -4, 0, 0, 0, 0, 0, 0, 0, 0, 0, !32 & 0, 2, 0, -2, 0, 0, 0, 0, 0, 0, 0, !33 & 3, -3, 0, 0, 0, 0, 0, 0, 0, 0, 0, !34 & 0, 2, -2, 0, 0, 0, 0, 0, 0, 0, 0, !35 & 0, 0, 0, 0, 0, 0, 0, 1, -2, 0, 0 !36 & / C DATA SDX & / -1719914.D0, -2.D0, ! 1 & 6434.D0, 141.D0, ! 2 & 715.D0, 0.D0, ! 3 & 715.D0, 0.D0, ! 4 & 486.D0, -5.D0, ! 5 & 159.D0, 0.D0, ! 6 & 0.D0, 0.D0, ! 7 & 39.D0, 0.D0, ! 8 & 33.D0, 0.D0, ! 9 & 31.D0, 0.D0, !10 & 8.D0, 0.D0, !11 & 8.D0, 0.D0, !12 & 21.D0, 0.D0, !13 & -19.D0, 0.D0, !14 & 17.D0, 0.D0, !15 & 16.D0, 0.D0, !16 & 16.D0, 0.D0, !17 & 11.D0, 0.D0, !18 & 0.D0, 0.D0, !19 & -11.D0, 0.D0, !20 & -7.D0, 0.D0, !21 & -10.D0, 0.D0, !22 & -9.D0, 0.D0, !23 & -9.D0, 0.D0, !24 & 0.D0, 0.D0, !25 & 0.D0, 0.D0, !26 & 8.D0, 0.D0, !27 & 8.D0, 0.D0, !28 & -4.D0, 0.D0, !29 & -4.D0, 0.D0, !30 & -6.D0, 0.D0, !31 & -1.D0, 0.D0, !32 & 4.D0, 0.D0, !33 & 0.D0, 0.D0, !34 & 5.D0, 0.D0, !35 & 5.D0, 0.D0 !36 & / C DATA CDX & / -25.D0, 0.D0, ! 1 & 28007.D0, -107.D0, ! 2 & 0.D0, 0.D0, ! 3 & 0.D0, 0.D0, ! 4 & -236.D0, -4.D0, ! 5 & 0.D0, 0.D0, ! 6 & 0.D0, 0.D0, ! 7 & 0.D0, 0.D0, ! 8 & -10.D0, 0.D0, ! 9 & 1.D0, 0.D0, !10 & -28.D0, 0.D0, !11 & -28.D0, 0.D0, !12 & 0.D0, 0.D0, !13 & 0.D0, 0.D0, !14 & 0.D0, 0.D0, !15 & 0.D0, 0.D0, !16 & 0.D0, 0.D0, !17 & -1.D0, 0.D0, !18 & -11.D0, 0.D0, !19 & -2.D0, 0.D0, !20 & -8.D0, 0.D0, !21 & 0.D0, 0.D0, !22 & 0.D0, 0.D0, !23 & 0.D0, 0.D0, !24 & -9.D0, 0.D0, !25 & -9.D0, 0.D0, !26 & 0.D0, 0.D0, !27 & 0.D0, 0.D0, !28 & -7.D0, 0.D0, !29 & -7.D0, 0.D0, !30 & -5.D0, 0.D0, !31 & -1.D0, 0.D0, !32 & -6.D0, 0.D0, !33 & -7.D0, 0.D0, !34 & -5.D0, 0.D0, !35 & 0.D0, 0.D0 !36 & / C DATA SDY & / 25.D0, -13.D0, ! 1 & 25697.D0, -95.D0, ! 2 & 6.D0, 0.D0, ! 3 & 0.D0, 0.D0, ! 4 & -216.D0, -4.D0, ! 5 & 2.D0, 0.D0, ! 6 & 0.D0, 0.D0, ! 7 & 0.D0, 0.D0, ! 8 & -9.D0, 0.D0, ! 9 & 1.D0, 0.D0, !10 & 25.D0, 0.D0, !11 & -25.D0, 0.D0, !12 & 0.D0, 0.D0, !13 & 0.D0, 0.D0, !14 & 0.D0, 0.D0, !15 & 0.D0, 0.D0, !16 & 1.D0, 0.D0, !17 & -1.D0, 0.D0, !18 & -10.D0, 0.D0, !19 & -2.D0, 0.D0, !20 & -8.D0, 0.D0, !21 & 0.D0, 0.D0, !22 & 0.D0, 0.D0, !23 & 0.D0, 0.D0, !24 & -8.D0, 0.D0, !25 & 8.D0, 0.D0, !26 & 0.D0, 0.D0, !27 & 0.D0, 0.D0, !28 & -6.D0, 0.D0, !29 & 6.D0, 0.D0, !30 & -4.D0, 0.D0, !31 & -2.D0, 0.D0, !32 & -5.D0, 0.D0, !33 & -6.D0, 0.D0, !34 & -4.D0, 0.D0, !35 & 0.D0, 0.D0 !36 & / C DATA CDY & / 1578089.D0, 156.D0, ! 1 & -5904.D0, -130.D0, ! 2 & -657.D0, 0.D0, ! 3 & -656.D0, 0.D0, ! 4 & -446.D0, 5.D0, ! 5 & -147.D0, 0.D0, ! 6 & 26.D0, 0.D0, ! 7 & -36.D0, 0.D0, ! 8 & -30.D0, 0.D0, ! 9 & -28.D0, 0.D0, !10 & 8.D0, 0.D0, !11 & -8.D0, 0.D0, !12 & -19.D0, 0.D0, !13 & 17.D0, 0.D0, !14 & -16.D0, 0.D0, !15 & 15.D0, 0.D0, !16 & -15.D0, 0.D0, !17 & -10.D0, 0.D0, !18 & 0.D0, 0.D0, !19 & 9.D0, 0.D0, !20 & 6.D0, 0.D0, !21 & 9.D0, 0.D0, !22 & -9.D0, 0.D0, !23 & -8.D0, 0.D0, !24 & 0.D0, 0.D0, !25 & 0.D0, 0.D0, !26 & -8.D0, 0.D0, !27 & -7.D0, 0.D0, !28 & 4.D0, 0.D0, !29 & -4.D0, 0.D0, !30 & 5.D0, 0.D0, !31 & -7.D0, 0.D0, !32 & -4.D0, 0.D0, !33 & 0.D0, 0.D0, !34 & -5.D0, 0.D0, !35 & -5.D0, 0.D0 !36 & / C DATA SDZ & / 10.D0, 32.D0, ! 1 & 11141.D0, -48.D0, ! 2 & -15.D0, 0.D0, ! 3 & 0.D0, 0.D0, ! 4 & -94.D0, 0.D0, ! 5 & -6.D0, 0.D0, ! 6 & 0.D0, 0.D0, ! 7 & 0.D0, 0.D0, ! 8 & -5.D0, 0.D0, ! 9 & 0.D0, 0.D0, !10 & 11.D0, 0.D0, !11 & -11.D0, 0.D0, !12 & 0.D0, 0.D0, !13 & 0.D0, 0.D0, !14 & 0.D0, 0.D0, !15 & 1.D0, 0.D0, !16 & -3.D0, 0.D0, !17 & -1.D0, 0.D0, !18 & -4.D0, 0.D0, !19 & -1.D0, 0.D0, !20 & -3.D0, 0.D0, !21 & 0.D0, 0.D0, !22 & 0.D0, 0.D0, !23 & 0.D0, 0.D0, !24 & -3.D0, 0.D0, !25 & 3.D0, 0.D0, !26 & 0.D0, 0.D0, !27 & 0.D0, 0.D0, !28 & -3.D0, 0.D0, !29 & 3.D0, 0.D0, !30 & -2.D0, 0.D0, !31 & 1.D0, 0.D0, !32 & -2.D0, 0.D0, !33 & -3.D0, 0.D0, !34 & -2.D0, 0.D0, !35 & 0.D0, 0.D0 !36 & / C DATA CDZ & / 684185.D0, -358.D0, ! 1 & -2559.D0, -55.D0, ! 2 & -282.D0, 0.D0, ! 3 & -285.D0, 0.D0, ! 4 & -193.D0, 0.D0, ! 5 & -61.D0, 0.D0, ! 6 & -59.D0, 0.D0, ! 7 & -16.D0, 0.D0, ! 8 & -13.D0, 0.D0, ! 9 & -12.D0, 0.D0, !10 & 3.D0, 0.D0, !11 & -3.D0, 0.D0, !12 & -8.D0, 0.D0, !13 & 8.D0, 0.D0, !14 & -7.D0, 0.D0, !15 & 7.D0, 0.D0, !16 & -6.D0, 0.D0, !17 & -5.D0, 0.D0, !18 & 0.D0, 0.D0, !19 & 4.D0, 0.D0, !20 & 3.D0, 0.D0, !21 & 4.D0, 0.D0, !22 & -4.D0, 0.D0, !23 & -4.D0, 0.D0, !24 & 0.D0, 0.D0, !25 & 0.D0, 0.D0, !26 & -3.D0, 0.D0, !27 & -3.D0, 0.D0, !28 & 2.D0, 0.D0, !29 & -2.D0, 0.D0, !30 & 2.D0, 0.D0, !31 & -4.D0, 0.D0, !32 & -2.D0, 0.D0, !33 & 0.D0, 0.D0, !34 & -2.D0, 0.D0, !35 & -2.D0, 0.D0 !36 & / C C C *** EXECUTION BEGINS HERE *** C C DETERMINE TIME ARGUMENT T (JULIAN CENTURIES SINCE J2000.0) T = (JDE-2451545.D0)/36525.D0 C C DETERMINE THE MEAN LONGITUDES (RADIANS) OF THE PERTURBING C SOURCES L2,...,F DO 10 I=1,11 MLON(I) = SOURCE(1,I)+SOURCE(2,I)*T 10 CONTINUE C C DEVELOP AND SUM THE TRIGONOMETRIC SERIES DX = 0. DY = 0. DZ = 0. DO 30 J=1,36 ARG = 0. DO 20 I=1,11 ARG = ARG+COEF(I,J)*MLON(I) 20 CONTINUE SA = SIN(ARG) CA = COS(ARG) DX = DX + (SDX(1,J)+SDX(2,J)*T)*SA + (CDX(1,J)+CDX(2,J)*T)*CA DY = DY + (SDY(1,J)+SDY(2,J)*T)*SA + (CDY(1,J)+CDY(2,J)*T)*CA DZ = DZ + (SDZ(1,J)+SDZ(2,J)*T)*SA + (CDZ(1,J)+CDZ(2,J)*T)*CA 30 CONTINUE XYZ(1) = DX XYZ(2) = DY XYZ(3) = DZ RETURN C END