Astro-PAL
view release on metacpan or search on metacpan
palsrc/palAoppa.c view on Meta::CPAN
{ 1 , 0 , 2 , 1 , 2 , 0 , 0 , 0 , 0 },
{ 1 , 0 , 2 , -1 , 2 , 0 , 0 , 0 , 0 },
{ -1 , 0 , 0 , 2 , -2 , 0 , 0 , 0 , 0 },
{ 0 , 0 , 2 , 1 , 1 , 0 , 0 , 0 , 0 },
{ -1 , 0 , 2 , 0 , -1 , 0 , 0 , 0 , 0 },
{ -1 , 0 , 2 , 4 , 1 , 0 , 0 , 0 , 0 },
{ 0 , 0 , 2 , 2 , 0 , 0 , 0 , 0 , 0 },
{ 1 , 1 , 2 , -2 , 1 , 0 , 0 , 0 , 0 },
{ 0 , 0 , 1 , 0 , 1 , 0 , 0 , 0 , 0 },
{ -1 , 0 , 2 , -1 , 1 , 0 , 0 , 0 , 0 },
{ -2 , 0 , 2 , 2 , 1 , 0 , 0 , 0 , 0 },
{ 2 , -1 , 0 , 0 , 0 , 0 , 0 , 0 , 0 },
{ 4 , 0 , 2 , 0 , 2 , 0 , 0 , 0 , 0 },
{ 2 , 1 , 2 , -2 , 2 , 0 , 0 , 0 , 0 },
{ 0 , 1 , 2 , 1 , 2 , 0 , 0 , 0 , 0 },
{ 1 , 0 , 4 , -2 , 2 , 0 , 0 , 0 , 0 },
{ 1 , 1 , 0 , 0 , -1 , 0 , 0 , 0 , 0 },
{ -2 , 0 , 2 , 4 , 1 , 0 , 0 , 0 , 0 },
{ 2 , 0 , 2 , 0 , 0 , 0 , 0 , 0 , 0 },
{ -1 , 0 , 1 , 0 , 0 , 0 , 0 , 0 , 0 },
{ 1 , 0 , 0 , 1 , 0 , 0 , 0 , 0 , 0 },
{ 0 , 1 , 0 , 2 , 1 , 0 , 0 , 0 , 0 },
{ -1 , 0 , 4 , 0 , 1 , 0 , 0 , 0 , 0 },
{ -1 , 0 , 0 , 4 , 1 , 0 , 0 , 0 , 0 },
{ 2 , 0 , 2 , 2 , 1 , 0 , 0 , 0 , 0 },
{ 2 , 1 , 0 , 0 , 0 , 0 , 0 , 0 , 0 },
{ 0 , 0 , 5 , -5 , 5 , -3 , 0 , 0 , 0 },
{ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 2 , 0 },
{ 0 , 0 , 1 , -1 , 1 , 0 , 0 , -1 , 0 },
{ 0 , 0 , -1 , 1 , -1 , 1 , 0 , 0 , 0 },
{ 0 , 0 , -1 , 1 , 0 , 0 , 2 , 0 , 0 },
{ 0 , 0 , 3 , -3 , 3 , 0 , 0 , -1 , 0 },
{ 0 , 0 , -8 , 8 , -7 , 5 , 0 , 0 , 0 },
{ 0 , 0 , -1 , 1 , -1 , 0 , 2 , 0 , 0 },
{ 0 , 0 , -2 , 2 , -2 , 2 , 0 , 0 , 0 },
{ 0 , 0 , -6 , 6 , -6 , 4 , 0 , 0 , 0 },
{ 0 , 0 , -2 , 2 , -2 , 0 , 8 , -3 , 0 },
{ 0 , 0 , 6 , -6 , 6 , 0 , -8 , 3 , 0 },
{ 0 , 0 , 4 , -4 , 4 , -2 , 0 , 0 , 0 },
{ 0 , 0 , -3 , 3 , -3 , 2 , 0 , 0 , 0 },
{ 0 , 0 , 4 , -4 , 3 , 0 , -8 , 3 , 0 },
{ 0 , 0 , -4 , 4 , -5 , 0 , 8 , -3 , 0 },
{ 0 , 0 , 0 , 0 , 0 , 2 , 0 , 0 , 0 },
{ 0 , 0 , -4 , 4 , -4 , 3 , 0 , 0 , 0 },
{ 0 , 1 , -1 , 1 , -1 , 0 , 0 , 1 , 0 },
{ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 1 , 0 },
{ 0 , 0 , 1 , -1 , 1 , 1 , 0 , 0 , 0 },
{ 0 , 0 , 2 , -2 , 2 , 0 , -2 , 0 , 0 },
{ 0 , -1 , -7 , 7 , -7 , 5 , 0 , 0 , 0 },
{ -2 , 0 , 2 , 0 , 2 , 0 , 0 , -2 , 0 },
{ -2 , 0 , 2 , 0 , 1 , 0 , 0 , -3 , 0 },
{ 0 , 0 , 2 , -2 , 2 , 0 , 0 , -2 , 0 },
{ 0 , 0 , 1 , -1 , 1 , 0 , 0 , 1 , 0 },
{ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 2 },
{ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 1 },
{ 2 , 0 , -2 , 0 , -2 , 0 , 0 , 3 , 0 },
{ 0 , 0 , 1 , -1 , 1 , 0 , 0 , -2 , 0 },
{ 0 , 0 , -7 , 7 , -7 , 5 , 0 , 0 , 0 }
};
double psi[ 194 ][4] = {
{ 3341.5000000000000 , 17206241.800000001 , 3.1000000000000001 , 17409.500000000000 },
{ -1716.8000000000000 , -1317185.3000000000 , 1.3999999999999999 , -156.80000000000001 },
{ 285.69999999999999 , -227667.00000000000 , 0.29999999999999999 , -23.500000000000000 },
{ -68.599999999999994 , -207448.00000000000 , 0.0000000000000000 , -21.399999999999999 },
{ 950.29999999999995 , 147607.89999999999 , -2.2999999999999998 , -355.00000000000000 },
{ -66.700000000000003 , -51689.099999999999 , 0.20000000000000001 , 122.59999999999999 },
{ -108.59999999999999 , 71117.600000000006 , 0.0000000000000000 , 7.0000000000000000 },
{ 35.600000000000001 , -38740.199999999997 , 0.10000000000000001 , -36.200000000000003 },
{ 85.400000000000006 , -30127.599999999999 , 0.0000000000000000 , -3.1000000000000001 },
{ 9.0000000000000000 , 21583.000000000000 , 0.10000000000000001 , -50.299999999999997 },
{ 22.100000000000001 , 12822.799999999999 , 0.0000000000000000 , 13.300000000000001 },
{ 3.3999999999999999 , 12350.799999999999 , 0.0000000000000000 , 1.3000000000000000 },
{ -21.100000000000001 , 15699.400000000000 , 0.0000000000000000 , 1.6000000000000001 },
{ 4.2000000000000002 , 6313.8000000000002 , 0.0000000000000000 , 6.2000000000000002 },
{ -22.800000000000001 , 5796.8999999999996 , 0.0000000000000000 , 6.0999999999999996 },
{ 15.699999999999999 , -5961.1000000000004 , 0.0000000000000000 , -0.59999999999999998 },
{ 13.100000000000000 , -5159.1000000000004 , 0.0000000000000000 , -4.5999999999999996 },
{ 1.8000000000000000 , 4592.6999999999998 , 0.0000000000000000 , 4.5000000000000000 },
{ -17.500000000000000 , 6336.0000000000000 , 0.0000000000000000 , 0.69999999999999996 },
{ 16.300000000000001 , -3851.0999999999999 , 0.0000000000000000 , -0.40000000000000002 },
{ -2.7999999999999998 , 4771.6999999999998 , 0.0000000000000000 , 0.50000000000000000 },
{ 13.800000000000001 , -3099.3000000000002 , 0.0000000000000000 , -0.29999999999999999 },
{ 0.20000000000000001 , 2860.3000000000002 , 0.0000000000000000 , 0.29999999999999999 },
{ 1.3999999999999999 , 2045.3000000000000 , 0.0000000000000000 , 2.0000000000000000 },
{ -8.5999999999999996 , 2922.5999999999999 , 0.0000000000000000 , 0.29999999999999999 },
{ -7.7000000000000002 , 2587.9000000000001 , 0.0000000000000000 , 0.20000000000000001 },
{ 8.8000000000000007 , -1408.0999999999999 , 0.0000000000000000 , 3.7000000000000002 },
{ 1.3999999999999999 , 1517.5000000000000 , 0.0000000000000000 , 1.5000000000000000 },
{ -1.8999999999999999 , -1579.7000000000000 , 0.0000000000000000 , 7.7000000000000002 },
{ 1.3000000000000000 , -2178.5999999999999 , 0.0000000000000000 , -0.20000000000000001 },
{ -4.7999999999999998 , 1286.8000000000000 , 0.0000000000000000 , 1.3000000000000000 },
{ 6.2999999999999998 , 1267.2000000000000 , 0.0000000000000000 , -4.0000000000000000 },
{ -1.0000000000000000 , 1669.3000000000000 , 0.0000000000000000 , -8.3000000000000007 },
{ 2.3999999999999999 , -1020.0000000000000 , 0.0000000000000000 , -0.90000000000000002 },
{ 4.5000000000000000 , -766.89999999999998 , 0.0000000000000000 , 0.0000000000000000 },
{ -1.1000000000000001 , 756.50000000000000 , 0.0000000000000000 , -1.7000000000000000 },
{ -1.3999999999999999 , -1097.3000000000000 , 0.0000000000000000 , -0.50000000000000000 },
{ 2.6000000000000001 , -663.00000000000000 , 0.0000000000000000 , -0.59999999999999998 },
{ 0.80000000000000004 , -714.10000000000002 , 0.0000000000000000 , 1.6000000000000001 },
{ 0.40000000000000002 , -629.89999999999998 , 0.0000000000000000 , -0.59999999999999998 },
{ 0.29999999999999999 , 580.39999999999998 , 0.0000000000000000 , 0.59999999999999998 },
{ -1.6000000000000001 , 577.29999999999995 , 0.0000000000000000 , 0.50000000000000000 },
{ -0.90000000000000002 , 644.39999999999998 , 0.0000000000000000 , 0.0000000000000000 },
{ 2.2000000000000002 , -534.00000000000000 , 0.0000000000000000 , -0.50000000000000000 },
{ -2.5000000000000000 , 493.30000000000001 , 0.0000000000000000 , 0.50000000000000000 },
{ -0.10000000000000001 , -477.30000000000001 , 0.0000000000000000 , -2.3999999999999999 },
{ -0.90000000000000002 , 735.00000000000000 , 0.0000000000000000 , -1.7000000000000000 },
{ 0.69999999999999996 , 406.19999999999999 , 0.0000000000000000 , 0.40000000000000002 },
{ -2.7999999999999998 , 656.89999999999998 , 0.0000000000000000 , 0.0000000000000000 },
{ 0.59999999999999998 , 358.00000000000000 , 0.0000000000000000 , 2.0000000000000000 },
{ -0.69999999999999996 , 472.50000000000000 , 0.0000000000000000 , -1.1000000000000001 },
{ -0.10000000000000001 , -300.50000000000000 , 0.0000000000000000 , 0.0000000000000000 },
{ -1.2000000000000000 , 435.10000000000002 , 0.0000000000000000 , -1.0000000000000000 },
{ 1.8000000000000000 , -289.39999999999998 , 0.0000000000000000 , 0.0000000000000000 },
{ 0.59999999999999998 , -422.60000000000002 , 0.0000000000000000 , 0.0000000000000000 },
{ 0.80000000000000004 , -287.60000000000002 , 0.0000000000000000 , 0.59999999999999998 },
{ -38.600000000000001 , -392.30000000000001 , 0.0000000000000000 , 0.0000000000000000 },
{ 0.69999999999999996 , -281.80000000000001 , 0.0000000000000000 , 0.59999999999999998 },
{ 0.59999999999999998 , -405.69999999999999 , 0.0000000000000000 , 0.0000000000000000 },
{ -1.2000000000000000 , 229.00000000000000 , 0.0000000000000000 , 0.20000000000000001 },
{ 1.1000000000000001 , -264.30000000000001 , 0.0000000000000000 , 0.50000000000000000 },
palsrc/palAoppa.c view on Meta::CPAN
{ 0.0000000000000000 , -14.400000000000000 , 0.0000000000000000 , 0.0000000000000000 },
{ 12.400000000000000 , 0.0000000000000000 , 0.0000000000000000 , 0.0000000000000000 },
{ -12.000000000000000 , 0.0000000000000000 , 0.0000000000000000 , 0.0000000000000000 },
{ -0.80000000000000004 , 0.0000000000000000 , 0.0000000000000000 , 0.0000000000000000 },
{ 10.900000000000000 , 0.10000000000000001 , 0.0000000000000000 , 0.0000000000000000 },
{ -10.800000000000001 , 0.0000000000000000 , 0.0000000000000000 , 0.0000000000000000 },
{ 10.500000000000000 , 0.0000000000000000 , 0.0000000000000000 , 0.0000000000000000 },
{ -10.400000000000000 , 0.0000000000000000 , 0.0000000000000000 , 0.0000000000000000 },
{ -11.199999999999999 , 0.0000000000000000 , 0.0000000000000000 , 0.0000000000000000 },
{ 10.500000000000000 , 0.10000000000000001 , 0.0000000000000000 , 0.0000000000000000 },
{ -1.3999999999999999 , 0.0000000000000000 , 0.0000000000000000 , 0.0000000000000000 },
{ 0.0000000000000000 , 0.10000000000000001 , 0.0000000000000000 , 0.0000000000000000 },
{ 0.69999999999999996 , 0.0000000000000000 , 0.0000000000000000 , 0.0000000000000000 },
{ -10.300000000000001 , 0.0000000000000000 , 0.0000000000000000 , 0.0000000000000000 },
{ -10.000000000000000 , 0.0000000000000000 , 0.0000000000000000 , 0.0000000000000000 },
{ 9.5999999999999996 , 0.0000000000000000 , 0.0000000000000000 , 0.0000000000000000 },
{ 9.4000000000000004 , 0.10000000000000001 , 0.0000000000000000 , 0.0000000000000000 },
{ 0.59999999999999998 , 0.0000000000000000 , 0.0000000000000000 , 0.0000000000000000 },
{ -87.700000000000003 , 4.4000000000000004 , -0.40000000000000002 , -6.2999999999999998 },
{ 46.299999999999997 , 22.399999999999999 , 0.50000000000000000 , -2.3999999999999999 },
{ 15.600000000000000 , -3.3999999999999999 , 0.10000000000000001 , 0.40000000000000002 },
{ 5.2000000000000002 , 5.7999999999999998 , 0.20000000000000001 , -0.10000000000000001 },
{ -30.100000000000001 , 26.899999999999999 , 0.69999999999999996 , 0.0000000000000000 },
{ 23.199999999999999 , -0.50000000000000000 , 0.0000000000000000 , 0.59999999999999998 },
{ 1.0000000000000000 , 23.199999999999999 , 3.3999999999999999 , 0.0000000000000000 },
{ -12.199999999999999 , -4.2999999999999998 , 0.0000000000000000 , 0.0000000000000000 },
{ -2.1000000000000001 , -3.7000000000000002 , -0.20000000000000001 , 0.10000000000000001 },
{ -18.600000000000001 , -3.7999999999999998 , -0.40000000000000002 , 1.8000000000000000 },
{ 5.5000000000000000 , -18.699999999999999 , -1.8000000000000000 , -0.50000000000000000 },
{ -5.5000000000000000 , -18.699999999999999 , 1.8000000000000000 , -0.50000000000000000 },
{ 18.399999999999999 , -3.6000000000000001 , 0.29999999999999999 , 0.90000000000000002 },
{ -0.59999999999999998 , 1.3000000000000000 , 0.0000000000000000 , 0.0000000000000000 },
{ -5.5999999999999996 , -19.500000000000000 , 1.8999999999999999 , 0.0000000000000000 },
{ 5.5000000000000000 , -19.100000000000001 , -1.8999999999999999 , 0.0000000000000000 },
{ -17.300000000000001 , -0.80000000000000004 , 0.0000000000000000 , 0.90000000000000002 },
{ -3.2000000000000002 , -8.3000000000000007 , -0.80000000000000004 , 0.29999999999999999 },
{ -0.10000000000000001 , 0.0000000000000000 , 0.0000000000000000 , 0.0000000000000000 },
{ -5.4000000000000004 , 7.7999999999999998 , -0.29999999999999999 , 0.0000000000000000 },
{ -14.800000000000001 , 1.3999999999999999 , 0.0000000000000000 , 0.29999999999999999 },
{ -3.7999999999999998 , 0.40000000000000002 , 0.0000000000000000 , -0.20000000000000001 },
{ 12.600000000000000 , 3.2000000000000002 , 0.50000000000000000 , -1.5000000000000000 },
{ 0.10000000000000001 , 0.0000000000000000 , 0.0000000000000000 , 0.0000000000000000 },
{ -13.600000000000000 , 2.3999999999999999 , -0.10000000000000001 , 0.0000000000000000 },
{ 0.90000000000000002 , 1.2000000000000000 , 0.0000000000000000 , 0.0000000000000000 },
{ -11.900000000000000 , -0.50000000000000000 , 0.0000000000000000 , 0.29999999999999999 },
{ 0.40000000000000002 , 12.000000000000000 , 0.29999999999999999 , -0.20000000000000001 },
{ 8.3000000000000007 , 6.0999999999999996 , -0.10000000000000001 , 0.10000000000000001 },
{ 0.0000000000000000 , 0.0000000000000000 , 0.0000000000000000 , 0.0000000000000000 },
{ 0.40000000000000002 , -10.800000000000001 , 0.29999999999999999 , 0.0000000000000000 },
{ 9.5999999999999996 , 2.2000000000000002 , 0.29999999999999999 , -1.2000000000000000 }
};
/* interval between fundamental epoch j2000.0 and given epoch (jc). */
t = (date-DJM0)/DJC;
/* mean anomaly of the moon. */
el = 134.96340251*PAL__DD2R+
fmod(t*(1717915923.2178+
t*( 31.8792+
t*( 0.051635+
t*( - 0.00024470)))),TURNAS)*PAL__DAS2R;
/* mean anomaly of the sun. */
elp = 357.52910918*PAL__DD2R+
fmod(t*( 129596581.0481+
t*( - 0.5532+
t*( 0.000136+
t*( - 0.00001149)))),TURNAS)*PAL__DAS2R;
/* mean argument of the latitude of the moon. */
f = 93.27209062*PAL__DD2R+
fmod(t*(1739527262.8478+
t*( - 12.7512+
t*( - 0.001037+
t*( 0.00000417)))),TURNAS)*PAL__DAS2R;
/* mean elongation of the moon from the sun. */
d = 297.85019547*PAL__DD2R+
fmod(t*(1602961601.2090+
t*( - 6.3706+
t*( 0.006539+
t*( - 0.00003169)))),TURNAS)*PAL__DAS2R;
/* mean longitude of the ascending node of the moon. */
om = 125.04455501*PAL__DD2R+
fmod(t*( - 6962890.5431+
t*( 7.4722+
t*( 0.007702+
t*( - 0.00005939)))),TURNAS)*PAL__DAS2R;
/* mean longitude of venus. */
ve = 181.97980085*PAL__DD2R+fmod(210664136.433548*t,TURNAS)*PAL__DAS2R;
/* mean longitude of mars.*/
ma = 355.43299958*PAL__DD2R+fmod( 68905077.493988*t,TURNAS)*PAL__DAS2R;
/* mean longitude of jupiter. */
ju = 34.35151874*PAL__DD2R+fmod( 10925660.377991*t,TURNAS)*PAL__DAS2R;
/* mean longitude of saturn. */
sa = 50.07744430*PAL__DD2R+fmod( 4399609.855732*t,TURNAS)*PAL__DAS2R;
/* geodesic nutation (fukushima 1991) in microarcsec. */
dp = -153.1*sin(elp)-1.9*sin(2*elp);
de = 0.0;
/* shirai & fukushima (2001) nutation series. */
for (j=NTERMS-1; j >= 0; j--) {
theta = ((double)na[j][0])*el+
((double)na[j][1])*elp+
((double)na[j][2])*f+
((double)na[j][3])*d+
((double)na[j][4])*om+
((double)na[j][5])*ve+
((double)na[j][6])*ma+
((double)na[j][7])*ju+
((double)na[j][8])*sa;
c = cos(theta);
s = sin(theta);
dp += (psi[j][0] + psi[j][2]*t)*c + (psi[j][1] + psi[j][3]*t)*s;
de += (eps[j][0] + eps[j][2]*t)*c + (eps[j][1] + eps[j][3]*t)*s;
( run in 1.259 second using v1.01-cache-2.11-cpan-ceb78f64989 )