home *** CD-ROM | disk | FTP | other *** search
/ Simtel MSDOS - Coast to Coast / simteldosarchivecoasttocoast2.iso / astrnomy / de118i.zip / INI118D.H < prev    next >
Text File  |  1993-01-31  |  37KB  |  1,231 lines
  1. #if UNK
  2.  
  3. DOUBLE DTR = (DOUBLE )1.74532925199432957692E-2;
  4. DOUBLE RTD = (DOUBLE )5.72957795130823208768E1;
  5. DOUBLE RTS = (DOUBLE )2.06264806247096355156E5; /* arc seconds per radian */
  6. DOUBLE STR = (DOUBLE )4.84813681109535993590E-6; /* radians per arc second */
  7. DOUBLE PI  = (DOUBLE )3.14159265358979323846;
  8. /* Standard epochs.  Note Julian epochs (J) are measured in
  9.  * years of 365.25 days.
  10.  */
  11. DOUBLE J2000 = (DOUBLE )2451545.0;    /* 2000 January 1.5 */
  12. DOUBLE B1950 = (DOUBLE )2433282.423;    /* 1950 January 0.923 Besselian epoch */
  13. DOUBLE J1900 = (DOUBLE )2415020.0;    /* 1900 January 0, 12h UT */
  14. DOUBLE B1900 = (DOUBLE )2415020.313;
  15.  
  16. DOUBLE Jcentury = (DOUBLE )3.6525e4; /* 1 Julian century */
  17. DOUBLE Jmillenium = (DOUBLE )3.6525e5;
  18. DOUBLE Jtmillenium = (DOUBLE )3.6525e6;
  19.  
  20. #if !(LDOUBLE)
  21. double Zero = 0.0;
  22. double Fourth = 0.25;
  23. double Half = 0.5;
  24. double One = 1.0;
  25. double OneandaHalf = 1.5;
  26. double Two = 2.0;
  27. double Three = 3.0;
  28. double ThreeandaHalf = 3.5;
  29. double Four = 4.0;
  30. double Five = 5.0;
  31. double Six = 6.0;
  32. double Ten = 10.0;
  33. #endif
  34.  
  35. /* Solar system barycentric velocity and position
  36.  * state of Mercury, Venus, EMB, ..., Pluto, MOON, Sun
  37.  * in the order dx/dt, x, dy/dt, y, dz/dt, z for each object.
  38.  *
  39.  * EMB is the arithmetic average of Earth and Moon, weighted by
  40.  * their masses. The coordinates of the Moon variable are the
  41.  * the solar system barycentric coordinates of the Moon
  42.  * minus those of the Earth.
  43.  */
  44.  
  45. /* Note, these particular numbers are derived from the heliocentric
  46.  * planet states and barycentric Sun state given in the
  47.  * documentation for the DE118.
  48.  */
  49.  
  50. /* June 28.0, 1969 */
  51. DOUBLE JD0 = ((DOUBLE )2440400.50);
  52.  
  53. DOUBLE yn0[6*NTOTAL] = {
  54. /* Lunar librations */
  55.  1.00015500951082782e-4,
  56.  6.03343019729665792e-3,
  57.  1.49039876740924137e-5,
  58.  3.82393007236771923e-1,
  59. -1.19434270765326686e-4,
  60.  1.28120436555489181e0,
  61. /* Mercury */
  62.  3.7085069798210382186e-003,
  63.  3.6030663368975339466e-001,
  64.  2.4854958767430945324e-002,
  65. -9.4812876741771684223e-002,
  66.  1.2929109014677844626e-002,
  67. -8.7466840117233140436e-002,
  68. /* Venus */
  69.  1.1156645711264016669e-002,
  70.  6.0786466491731583464e-001,
  71.  1.5494075513638794325e-002,
  72. -3.5518362463675619232e-001,
  73.  6.2773904546696609267e-003,
  74. -1.9824142909855515515e-001,
  75.  /* EMB */
  76.  1.6833251020051496668e-002,
  77.  1.0820747754938311664e-001,
  78.  1.5602036176919105255e-003,
  79. -9.2711110739430602933e-001,
  80.  6.7646174015847273137e-004,
  81. -4.0209347855944090112e-001,
  82. /* Mars */
  83.  1.4481919298277924969e-002,
  84. -1.2796408611369531836e-001,
  85.  8.0528538390447499843e-005,
  86. -1.3262618005333617013e+000,
  87. -3.5188931029397090065e-004,
  88. -6.0530808652523961512e-001,
  89. /* Jupiter */
  90.  1.0053452569924098185e-003,
  91. -5.3896824544609061333e+000,
  92. -6.5298425191689416643e-003,
  93. -7.7026549518616593034e-001,
  94. -2.8258787532429609536e-003,
  95. -1.9866431165907522014e-001,
  96. /* Saturn */
  97. -3.1594662504012930114e-003,
  98.  7.9527768530257360864e+000,
  99.  4.3714634278372622354e-003,
  100.  4.5078822184006553686e+000,
  101.  1.9441395169137103763e-003,
  102.  1.5201955253183338898e+000,
  103. /* Uranus */
  104.  1.7108310564806817248e-004,
  105. -1.8278236586353147533e+001,
  106. -3.7646704682815900043e-003,
  107. -9.5764572881482056433e-001,
  108. -1.6519678610257000136e-003,
  109. -1.6132190397271035415e-001,
  110. /* Neptune */
  111.  2.6225242764289213785e-003,
  112. -1.6367191358770888335e+001,
  113. -1.5277473123858904045e-003,
  114. -2.3760896725373076342e+001,
  115. -6.9183197562182804864e-004,
  116. -9.3213866179497290101e+000,
  117. /* Pluto */
  118.  2.8177758090360373050e-004,
  119. -3.0447680255169362534e+001,
  120. -3.1469590804946202045e-003,
  121. -5.3177934960261367037e-001,
  122. -1.0794238049289112837e-003,
  123.  9.0596584886274922101e+000,
  124. #if MOON
  125.  5.98752118417335956e-4,
  126. -8.35703164195952601e-4,
  127. -1.74153713527242722e-4,
  128. -1.98543915768166071e-3,
  129. -8.84771962437116281e-5,
  130. -1.08326877048661754e-3,
  131. #endif
  132. /* Sun */
  133. -2.8369446340813151639e-007,
  134.  4.5144118714356666407e-003,
  135.  5.1811944086463255444e-006,
  136.  7.2282841152065867346e-004,
  137.  2.2306588340621263489e-006,
  138.  2.4659100492567986271e-004,
  139. /* Asteroid states follow */
  140. };
  141.  
  142.  
  143.  
  144. /* Test state, 100 days later: */
  145. #if 0
  146. DOUBLE JD1 = ((DOUBLE )2440500.50);
  147.  
  148. DOUBLE yn1[6*NTOTAL] = {
  149. -3.85893930164594646e-5,
  150.  1.29694398195755468e-2,
  151.  1.14905257263582262e-5,
  152.  3.82508699238156682e-1,
  153.  4.75391383166468672e-5,
  154.  1.27428371119265801e0,
  155.  
  156. -2.43015281791706444e-002,
  157.  2.40773527079857565e-001,
  158.  1.87910251556777273e-002,
  159.  1.96954621939863219e-001,
  160.  1.25591312851252066e-002,
  161.  8.10048178490518384e-002,
  162.  
  163. -1.63440955802718280e-002,
  164. -4.21860322028048782e-001,
  165. -1.14246422850682140e-002,
  166.  5.17869026250609690e-001,
  167. -4.11622909296383714e-003,
  168.  2.60203829195766570e-001,
  169.  
  170. -3.94484632235194096e-003,
  171.  9.81177706008309405e-001,
  172.  1.53727031456690289e-002,
  173.  1.96484576859672724e-001,
  174.  6.66612671337020014e-003,
  175.  8.51297144135776147e-002,
  176.  
  177.  8.50596535330489333e-003,
  178.  1.14162585220489816e+000,
  179.  1.16329928693268469e-002,
  180. -7.02290330317807430e-001,
  181.  5.11130674190038601e-003,
  182. -3.52984793305744779e-001,
  183.  
  184.  1.97739891612302323e-003,
  185. -5.24031415190823051e+000,
  186. -6.33006771988879063e-003,
  187. -1.41424053721658884e+000,
  188. -2.76393462026492719e-003,
  189. -4.78580847895956282e-001,
  190.  
  191. -3.45070902605290058e-003,
  192.  7.62218572268596801e+000,
  193.  4.19480008161587867e-003,
  194.  4.93634138197818904e+000,
  195.  1.88368464231530789e-003,
  196.  1.71165065659011950e+000,
  197.  
  198.  2.59339002502375430e-004,
  199. -1.82567140746315272e+001,
  200. -3.75914213576224408e-003,
  201. -1.33385143258585295e+000,
  202. -1.65079191573496991e-003,
  203. -3.26466517509008289e-001,
  204.  
  205.  2.63977839373647663e-003,
  206. -1.61040739057229584e+001,
  207. -1.50241016013489575e-003,
  208. -2.39124059584972814e+001,
  209. -6.81887221476645979e-004,
  210. -9.39007319300685950e+000,
  211.  
  212.  3.09967640221318611e-004,
  213. -3.04180939220586550e+001,
  214. -3.14632190270354497e-003,
  215. -8.46445850516399332e-001,
  216. -1.08776892324236050e-003,
  217.  8.95129835863062695e+000,
  218.  
  219. #if MOON
  220. -4.22681402426694080e-4,
  221. -1.72474926563029233e-3,
  222. -3.21609581451700993e-4,
  223.  1.82647356024942241e-3,
  224. -1.82676776570974036e-4,
  225.  9.67995005426182924e-4,
  226. #endif
  227. -9.96208174079740373e-007,  /* rotated from DE200 */
  228.  4.45122038203489470e-003,
  229.  5.06057698826071977e-006,
  230.  1.23363163514498492e-003,
  231.  2.19376032512757903e-006,
  232.  4.67156762898216060e-004,
  233. };
  234. #endif /* 0 */
  235.  
  236. #if 1
  237. DOUBLE JD1 = ((DOUBLE )2440800.50);
  238.  
  239. DOUBLE yn1[6*NTOTAL] = {
  240.  1.19798509852348462e-4,
  241.  3.26749396278042651e-2,
  242. -6.85863913398886621e-5,
  243.  3.85255316277713118e-1,
  244. -9.87921606265917336e-5,
  245.  1.25594900035445747e0,
  246.  
  247.  1.1640059959070938016e-002,
  248. -3.4350297408563828055e-001,
  249. -1.7964875644231269124e-002,
  250. -2.5118910315168408963e-001,
  251. -1.0817311005210173814e-002,
  252. -9.9170891227838064646e-002,
  253.  
  254.  1.8597282203783061216e-002,
  255. -2.6914630075961837156e-001,
  256. -6.5864799228651037249e-003,
  257. -6.1655661757510340864e-001,
  258. -4.1428517490915759141e-003,
  259. -2.6070731202257053265e-001,
  260.  
  261.  1.3083581197550196016e-002,
  262.  6.4260215465460413745e-001,
  263.  9.8846107560980890857e-003,
  264. -7.2083597294516633464e-001,
  265.  4.2862797207098424858e-003,
  266. -3.1263839868006851566e-001,
  267.  
  268. -1.0289202815795734684e-002,
  269. -1.0394941428468679346e+000,
  270. -7.0723334989418045152e-003,
  271.  1.1522403836312279294e+000,
  272. -2.9703138309275424941e-003,
  273.  5.5677401157133424433e-001,
  274.  
  275.  4.6228653666202124358e-003,
  276. -4.2373615722355993645e+000,
  277. -5.0555918692615562352e-003,
  278. -3.1464733091277224306e+000,
  279. -2.2818455722493645242e-003,
  280. -1.2462241659031373286e+000,
  281.  
  282. -4.2556740044874090445e-003,
  283.  6.4633979534971069954e+000,
  284.  3.5633047172379600052e-003,
  285.  6.1037564142437140792e+000,
  286.  1.6573768541155527860e-003,
  287.  2.2444618449141959815e+000,
  288.  
  289.  5.2274776790945661154e-004,
  290. -1.8139340904007687965e+001,
  291. -3.7317457890160985748e-003,
  292. -2.4578888288247007605e+000,
  293. -1.6425078479167742140e-003,
  294. -8.2063943185823257564e-001,
  295.  
  296.  2.6898588576179062356e-003,
  297. -1.5304564931392374665e+001,
  298. -1.4254385368896515149e-003,
  299. -2.4351618762495032811e+001,
  300. -6.5161737589378315099e-004,
  301. -9.5901150170962550989e+000,
  302.  
  303.  3.9518201228036738355e-004,
  304. -3.0312344985766555864e+001,
  305. -3.1426251869952901448e-003,
  306. -1.7898553889579568107e+000,
  307. -1.1124350707669052540e-003,
  308.  8.6212535443336042115e+000,
  309.  
  310. #if MOON
  311. -4.58408893102436927e-4,
  312. -1.60130208907259920e-3,
  313. -2.78697065392994499e-4,
  314.  1.95720839700677492e-3,
  315. -1.70697259101564644e-4,
  316.  9.68483037769266028e-4,
  317. #endif
  318.  /* rotated from DE200 */
  319. /*
  320. -3.44322101864905357e-006,
  321.  3.77983010651375442e-003,
  322.  4.03605242892639870e-006,
  323.  2.63055823065286220e-003,
  324.  1.81002374725178528e-006,
  325.  1.08188006521861207e-003,
  326. */
  327. -3.4432210186490844452e-006,
  328.  3.7798301065137354516e-003,
  329.  4.0360524289263751236e-006,
  330.  2.6305582306528893588e-003,
  331.  1.8100237472517860139e-006,
  332.  1.0818800652186113487e-003,
  333. };
  334. #endif /* 1 */
  335.  
  336.  
  337. /* Speed of light.
  338.  */
  339. /* DOUBLE CLIGHT = (DOUBLE)2.99792458e5; */
  340. DOUBLE C = ((DOUBLE ) 173.144632720536344565);
  341.  
  342. /* Earth's mass divided by Moon's mass
  343.  */
  344. #define EMRATd ((DOUBLE ) 8.13005869999999999376E1)
  345. DOUBLE EMRAT = EMRATd;
  346.  
  347. /* Gaussian gravitational constant */
  348. DOUBLE KG = 0.01720209895;
  349.  
  350. /* GM's of the solar system bodies
  351.  * These are scaled such that GMsun = k^2
  352.  * (k = Gaussian gravitational constant).
  353.  */
  354. DOUBLE GMs[] = {
  355. #if LIBRAT
  356.  0.0, /*place holder for libration */
  357. #endif
  358.  4.91254745145081175785E-11, /* Mercury */
  359.  7.24345620963276523095E-10, /* Venus */
  360.  
  361. #if MOON
  362. /* (GMB*EMRAT)/(1.0+EMRAT), Earth only */
  363.  8.88769273403302327042E-10,
  364. #else
  365.  8.99701165855730774179E-10, /* EMB */
  366. #endif
  367.  9.54952894222405763492E-11, /* Mars */
  368.  2.82534210344592625472E-7, /* Jupiter */
  369.  8.45946850483065929285E-8, /* Saturn */
  370.  1.28881623813803488851E-8, /* Uranus */
  371.  1.53211248128427618918E-8, /* Neptune */
  372.  2.27624775186369921644E-12, /* Pluto */
  373. #if MOON
  374. /* (GMB)/(1.0+EMRAT),  Moon only */
  375.  1.09318924524284471369E-11,
  376. #endif
  377.  2.95912208285591102582E-4, /* Sun */
  378. #if AROIDS
  379.  1.746e-13, /* Ceres */
  380.  3.2e-14,   /* Pallas */
  381.  4.08e-14,  /* Vesta */
  382.  1.6e-15,   /* Iris */
  383.  2.6e-15    /* Bamberga */
  384. #endif
  385. };
  386.  
  387.  
  388. /* Astronomical unit */
  389. DOUBLE AU = 1.4959787066e11;
  390. /* kilometers: */
  391. DOUBLE RADS = 6.96e5;
  392. DOUBLE RADM = 1738.0;
  393. DOUBLE RADE = 6378.14;
  394.  
  395.  
  396.  
  397. DOUBLE Je[] = {
  398. 1.08263e-3, /* J2E */
  399. -2.54e-6,   /* J3E */
  400. -1.61e-6    /* J4E */
  401. };
  402.  
  403. DOUBLE Jm[] = {
  404. 2.02150907893e-4,  /* J2M */
  405. 1.21260448837e-5,  /* J3M */
  406. -1.45383007072e-7, /* J4M */
  407. };
  408.  
  409. DOUBLE Cnm[8] = {
  410.  2.230351309e-5,   /* C22 */
  411.  3.07082741328e-5, /* C31 */
  412.  4.88840471683e-6, /* C32 */
  413.  1.43603108489e-6, /* C33 */
  414. -7.17780149806e-6, /* C41 */
  415. -1.43951838385e-6, /* C42 */
  416. -8.54788154819e-8, /* C43 */
  417. -1.5490389313e-7   /* C44 */
  418. };
  419.  
  420. DOUBLE Snm[8] = {
  421.  0.0,             /* S22 */
  422.  5.61066891941e-6, /* S31 */
  423.  1.68743052295e-6, /* S32 */
  424. -3.343544677e-7,   /* S33 */
  425.  2.94743374914e-6, /* S41 */
  426. -2.8843721272e-6,  /* S42 */
  427. -7.88967312839e-7, /* S43 */
  428.  5.6404155572e-8   /* S44 */
  429. };
  430.  
  431.  
  432. DOUBLE K2M = 2.22160557154e-2;
  433. DOUBLE LOVENO  = 0.30;
  434. DOUBLE PHASE = 4.0700012e-2;
  435.  
  436. /* DE118 initial psidot: */
  437. DOUBLE PSLP1 = 2.29971502189818919e-1;
  438. /* PSLPI = 2 pi / PSLP1 = period in days */
  439. DOUBLE PSLPI =  27.32158222801637949485;
  440. /* PSLPA + PSLPB = PSLPI */
  441. DOUBLE PSLPIA =  27.3215820789337158203125; /* only top 26 bits are nonzero */
  442. DOUBLE PSLPIB =  1.490826636745416323896386e-7;
  443. DOUBLE JDEPOCH = 2440400.5;
  444.  
  445. /* (C-A)/B */
  446. #define LBETd 6.31686773468e-4
  447. DOUBLE LBET = LBETd;
  448. /* (B-A)/C */
  449. #define LGAMd 2.28022183594e-4
  450. DOUBLE LGAM = LGAMd;
  451. /* C/MR^2 */
  452. /*
  453. DOUBLE CMRSQ 0.3906895261319410091
  454. DOUBLE LGAMBET = (LGAMd - LBETd)/(1.0-LBETd*LGAMd);
  455. DOUBLE CMR2 = CMRSQ*AMd*AMd;
  456. DOUBLE BMR2 = CMRSQ*AMd*AMd*(1.0 + LGAMd)/(1.0 + LBETd);
  457. DOUBLE AMR2 = CMRSQ*AMd*AMd*(1.0 - LBETd * LGAMd)/(1.0 + LBETd);
  458. */
  459. DOUBLE LGAMBET = -4.03664648017289733947E-4;
  460. /* equatorial radius of Earth, in au */
  461. DOUBLE AE =  4.26352325064817808471E-5;
  462. /* equatorial radius of Moon, in au */
  463. #define AMd  1.16178124180819699985E-5
  464. DOUBLE AM = AMd;
  465.  
  466. DOUBLE CMRSQ = 0.3906895261319410091;
  467. DOUBLE CMR2 = 5.2732758299330413853212309E-11;
  468. DOUBLE BMR2 = 5.2711485389894113965222358E-11;
  469. DOUBLE AMR2 = 5.2699461151199766018387453E-11;
  470. #endif /* UNK */
  471. #if IBMPC
  472. short DTR[] = {0x9d39,0xa252,0xdf46,0x3f91};
  473. short RTD[] = {0xc1f8,0x1a63,0xa5dc,0x404c};
  474. short RTS[] = {0xad7b,0x7331,0x2dc6,0x4109};
  475. short STR[] = {0x8f9d,0xb2ff,0x55a5,0x3ed4};
  476. short PI[] = {0x2d18,0x5444,0x21fb,0x4009};
  477. short J2000[] = {0x0000,0x8000,0xb42c,0x4142};
  478. short B1950[] = {0xdd2f,0x3624,0x9081,0x4142};
  479. short J1900[] = {0x0000,0x0000,0x6cd6,0x4142};
  480. short B1900[] = {0x624e,0x2810,0x6cd6,0x4142};
  481. short Jcentury[] = {0x0000,0x0000,0xd5a0,0x40e1};
  482. short Jmillenium[] = {0x0000,0x0000,0x4b08,0x4116};
  483. short Jtmillenium[] = {0x0000,0x0000,0xddca,0x414b};
  484. short Zero[] = {0x0000,0x0000,0x0000,0x0000};
  485. short Fourth[] = {0x0000,0x0000,0x0000,0x3fd0};
  486. short Half[] = {0x0000,0x0000,0x0000,0x3fe0};
  487. short One[] = {0x0000,0x0000,0x0000,0x3ff0};
  488. short OneandaHalf[] = {0x0000,0x0000,0x0000,0x3ff8};
  489. short Two[] = {0x0000,0x0000,0x0000,0x4000};
  490. short Three[] = {0x0000,0x0000,0x0000,0x4008};
  491. short ThreeandaHalf[] = {0x0000,0x0000,0x0000,0x400c};
  492. short Four[] = {0x0000,0x0000,0x0000,0x4010};
  493. short Five[] = {0x0000,0x0000,0x0000,0x4014};
  494. short Six[] = {0x0000,0x0000,0x0000,0x4018};
  495. short Ten[] = {0x0000,0x0000,0x0000,0x4024};
  496.  
  497. short JD0[] = {0x0000,0x4000,0x9e68,0x4142};
  498. short yn0[4*6*NTOTAL] = {
  499. 0xa479,0x3903,0x37ed,0x3f1a,
  500. 0xd720,0x9614,0xb682,0x3f78,
  501. 0xc10e,0x77ee,0x4184,0x3eef,
  502. 0x4406,0x8513,0x7920,0x3fd8,
  503. 0x7f80,0x25d8,0x4f19,0xbf1f,
  504. 0xcf9f,0x2618,0x7fd0,0x3ff4,
  505. 0x0dd3,0x8640,0x614d,0x3f6e,
  506. 0xad8d,0x8e0e,0x0f43,0x3fd7,
  507. 0xbbc3,0x0c32,0x7394,0x3f99,
  508. 0x787c,0x1cd8,0x45a8,0xbfb8,
  509. 0x5fd5,0xa313,0x7a93,0x3f8a,
  510. 0xb8df,0x11c9,0x643a,0xbfb6,
  511. 0x18a2,0xa3b4,0xd94b,0x3f86,
  512. 0xd806,0x9906,0x73a0,0x3fe3,
  513. 0x1500,0x9ce7,0xbb5b,0x3f8f,
  514. 0xef55,0x18f8,0xbb54,0xbfd6,
  515. 0xfd65,0x2b4f,0xb652,0x3f79,
  516. 0x657d,0xa358,0x5ff9,0xbfc9,
  517. 0x467c,0x5a77,0x3cbc,0x3f91,
  518. 0xdb72,0x3941,0xb37c,0x3fbb,
  519. 0x6d8d,0xe0d6,0x8ff7,0x3f59,
  520. 0x8a59,0xe9c0,0xaae4,0xbfed,
  521. 0x60c3,0x8686,0x2a92,0x3f46,
  522. 0x4052,0x4916,0xbbe6,0xbfd9,
  523. 0xc307,0x4e27,0xa8b2,0x3f8d,
  524. 0xde62,0x8e75,0x6120,0xbfc0,
  525. 0xfa3e,0xc14f,0x1c2d,0x3f15,
  526. 0x990b,0x4b33,0x385e,0xbff5,
  527. 0xb969,0x145b,0x0fb9,0xbf37,
  528. 0x2af2,0x1074,0x5eaf,0xbfe3,
  529. 0x86d9,0x3fff,0x78b9,0x3f50,
  530. 0xf03b,0xead6,0x8f08,0xc015,
  531. 0xf544,0x4115,0xbf09,0xbf7a,
  532. 0xfa7f,0xd2e1,0xa603,0xbfe8,
  533. 0x0667,0x1a7c,0x264c,0xbf67,
  534. 0xa25a,0x08ba,0x6dd5,0xbfc9,
  535. 0x39d8,0x86fc,0xe1e1,0xbf69,
  536. 0x867d,0xbc40,0xcfa4,0x401f,
  537. 0xdfa5,0xc754,0xe7cf,0x3f71,
  538. 0x00c3,0x46b9,0x0812,0x4012,
  539. 0x635e,0xe937,0xda4f,0x3f5f,
  540. 0x48d8,0x8b0c,0x52b8,0x3ff8,
  541. 0xc0dd,0xaff3,0x6c98,0x3f26,
  542. 0xeffd,0x834e,0x473a,0xc032,
  543. 0xdb5c,0x1153,0xd716,0xbf6e,
  544. 0x3d5a,0xa7cd,0xa508,0xbfee,
  545. 0x3d83,0xfbfa,0x10da,0xbf5b,
  546. 0x7aaa,0x36d8,0xa632,0xbfc4,
  547. 0xdebd,0x0001,0x7bd5,0x3f65,
  548. 0x4b76,0x40bd,0x5e00,0xc030,
  549. 0x657e,0x2f91,0x07d6,0xbf59,
  550. 0x1c61,0x20b7,0xc2ca,0xc037,
  551. 0xebac,0xdad4,0xab81,0xbf46,
  552. 0xeea9,0xc96a,0xa48c,0xc022,
  553. 0x6d46,0x7ea6,0x7771,0x3f32,
  554. 0x0471,0x2c57,0x729b,0xc03e,
  555. 0x738a,0xcaa4,0xc7a6,0xbf69,
  556. 0x676f,0x2067,0x0456,0xbfe1,
  557. 0x52c0,0x7c34,0xaf6e,0xbf51,
  558. 0x2af1,0x8eb3,0x1e8b,0x4022,
  559. 0x7ce0,0x622d,0x9eb2,0x3f43,
  560. 0x89e4,0xe133,0x6262,0xbf4b,
  561. 0x4f4f,0x0216,0xd3a1,0xbf26,
  562. 0x9bd6,0x8803,0x43c4,0xbf60,
  563. 0x7dd2,0xa840,0x319a,0xbf17,
  564. 0x0681,0xfc4a,0xbf8e,0xbf51,
  565. 0x48a5,0x71ec,0x09d5,0xbe93,
  566. 0x6a9d,0x3593,0x7db4,0x3f72,
  567. 0x746b,0xdfe0,0xbb43,0x3ed5,
  568. 0x49f8,0x31ad,0xaf86,0x3f47,
  569. 0xfbb2,0xa9ae,0xb64d,0x3ec2,
  570. 0x8354,0x4d38,0x291c,0x3f30,
  571. };
  572. short JD1[] = {0x0000,0x4000,0x9f30,0x4142};
  573. short yn1[4*6*NTOTAL] = {
  574. 0xa4f1,0xba4a,0x678a,0x3f1f,
  575. 0x4c98,0x0a33,0xbac5,0x3fa0,
  576. 0x409d,0x3b24,0xfac1,0xbf11,
  577. 0x776e,0xea01,0xa805,0x3fd8,
  578. 0x42dd,0x655e,0xe5d4,0xbf19,
  579. 0x75ea,0xfa9f,0x185d,0x3ff4,
  580. 0x110e,0x66c8,0xd6be,0x3f87,
  581. 0xb307,0xe5f1,0xfbf3,0xbfd5,
  582. 0x8750,0x6579,0x6562,0xbf92,
  583. 0x79b6,0x75c9,0x137b,0xbfd0,
  584. 0xfba7,0xe7fb,0x2762,0xbf86,
  585. 0xeab8,0x7689,0x6343,0xbfb9,
  586. 0x5d56,0x7b70,0x0b2a,0x3f93,
  587. 0x8633,0x67e6,0x39b1,0xbfd1,
  588. 0xda89,0xbdd4,0xfa6c,0xbf7a,
  589. 0xc222,0xf193,0xbad4,0xbfe3,
  590. 0x7407,0x4c64,0xf818,0xbf70,
  591. 0xc2b1,0xb8bd,0xaf6d,0xbfd0,
  592. 0x4358,0x8add,0xcb90,0x3f8a,
  593. 0x94c2,0x64d2,0x9032,0x3fe4,
  594. 0xf4e0,0xff72,0x3e61,0x3f84,
  595. 0x8df8,0x9a32,0x1116,0xbfe7,
  596. 0x228f,0x738d,0x8e7d,0x3f71,
  597. 0x8045,0x7c73,0x0244,0xbfd4,
  598. 0x1a91,0x6cc4,0x1281,0xbf85,
  599. 0x92d8,0x9c3e,0xa1c4,0xbff0,
  600. 0x8a30,0x1159,0xf7e1,0xbf7c,
  601. 0x39d3,0x9ccd,0x6f93,0x3ff2,
  602. 0xf56c,0x1867,0x5533,0xbf68,
  603. 0xc57a,0xbb5e,0xd117,0x3fe1,
  604. 0xdec3,0xf903,0xef6c,0x3f72,
  605. 0x50f0,0xe978,0xf30e,0xc010,
  606. 0xafa7,0x1bd8,0xb52c,0xbf74,
  607. 0x8613,0x32c3,0x2bfa,0xc009,
  608. 0x4682,0x836f,0xb160,0xbf62,
  609. 0x9f24,0xc040,0xf088,0xbff3,
  610. 0x53cd,0xcabf,0x6e65,0xbf71,
  611. 0x36af,0xfe3d,0xda84,0x4019,
  612. 0x30d1,0xa73d,0x30ca,0x3f6d,
  613. 0xb53a,0x1f17,0x6a3f,0x4018,
  614. 0x3d27,0xd8ad,0x278a,0x3f5b,
  615. 0x35db,0x6968,0xf4a8,0x4001,
  616. 0xee4c,0x489a,0x2120,0x3f41,
  617. 0xb546,0xd871,0x23ab,0xc032,
  618. 0xaca3,0xc3dd,0x9209,0xbf6e,
  619. 0x0bf2,0x9e48,0xa9c1,0xc003,
  620. 0xdc66,0x5f5f,0xe92d,0xbf5a,
  621. 0x73c0,0xa034,0x42ad,0xbfea,
  622. 0xde76,0xfa60,0x090a,0x3f66,
  623. 0xad6e,0xef47,0x9bef,0xc02e,
  624. 0x5f1e,0xf97f,0x5ab8,0xbf57,
  625. 0x937f,0xafed,0x5a03,0xc038,
  626. 0xfb0f,0xa8d3,0x5a29,0xbf45,
  627. 0x9d22,0x8e36,0x2e23,0xc023,
  628. 0x162a,0xd197,0xe60d,0x3f39,
  629. 0xefdb,0xd74a,0x4ff5,0xc03e,
  630. 0xfa3f,0x0cda,0xbe90,0xbf69,
  631. 0x4cf0,0x6782,0xa33f,0xbffc,
  632. 0x0ed3,0x0fdd,0x39e4,0xbf52,
  633. 0xdfb3,0xf1ce,0x3e14,0x4021,
  634. 0xf74e,0x343d,0x0ad3,0xbf3e,
  635. 0x1bd1,0x06a1,0x3c59,0xbf5a,
  636. 0xd241,0xc806,0x43c2,0xbf32,
  637. 0xac3c,0x41cf,0x0890,0x3f60,
  638. 0x5f17,0x4a6c,0x5fa6,0xbf26,
  639. 0x6280,0x7caa,0xbc39,0x3f4f,
  640. 0x5695,0xc604,0xe242,0xbecc,
  641. 0x2615,0xd624,0xf6e0,0x3f6e,
  642. 0xd7aa,0xa4db,0xedad,0x3ed0,
  643. 0xa17d,0x3244,0x8cae,0x3f65,
  644. 0x7a70,0x281a,0x5dfe,0x3ebe,
  645. 0x486d,0xdfe4,0xb9bb,0x3f51,
  646. };
  647.  
  648. short CLIGHT[] = {0xf3b6,0xd4fd,0x4c41,0x4112};
  649. short C[] = {0x9453,0xd4cc,0xa4a0,0x4065};
  650. short EMRAT[] = {0xa693,0xd141,0x533c,0x4054};
  651. short KG[] = {0xb69a,0x51a6,0x9d6d,0x3f91};
  652. short GMs[] = {
  653. 0x0000,0x0000,0x0000,0x0000,
  654. 0x30a7,0xbfef,0x01cb,0x3dcb,
  655. 0x9305,0x5596,0xe369,0x3e08,
  656. 0x614e,0x7bf8,0x89b2,0x3e0e,
  657. 0xe19b,0x330d,0x3fe2,0x3dda,
  658. 0x594a,0x991b,0xf5e6,0x3e92,
  659. 0xa2cf,0x700f,0xb54d,0x3e76,
  660. 0xddff,0x9a62,0xad57,0x3e4b,
  661. 0xf7c8,0x4ffd,0x7370,0x3e50,
  662. 0xde3d,0x7dd0,0x05a7,0x3d84,
  663. 0x5e14,0xbc9f,0x0a1b,0x3da8,
  664. 0x95e6,0x41cf,0x6495,0x3f33,
  665. 0x5b6d,0xc1f4,0x92a0,0x3d48,
  666. 0x5616,0x9ee7,0x03af,0x3d22,
  667. 0x8dc2,0x5100,0xf7e6,0x3d26,
  668. 0x89bc,0x97d8,0xd2b2,0x3cdc,
  669. 0xefe9,0x1b5f,0x6b31,0x3ce7,
  670. };
  671. short AU[] = {0x0000,0x2d22,0x6a5d,0x4241};
  672. short RADS[] = {0x0000,0x0000,0x3d80,0x4125};
  673. short RADM[] = {0x0000,0x0000,0x2800,0x409b};
  674. short RADE[] = {0x3d71,0xd70a,0xea23,0x40b8};
  675. short AE[] = {0x85dd,0x6f26,0x5a67,0x3f06};
  676. short Je[] = {
  677. 0x107b,0x1c65,0xbce1,0x3f51,
  678. 0x2b21,0xc46c,0x4e9b,0xbec5,
  679. 0x1061,0xb881,0x02e5,0xbebb,
  680. };
  681.  
  682. short AM[] = {0xde5a,0xfaf9,0x5d43,0x3ee8};
  683.  
  684. short Jm[] = {
  685. 0x70b6,0x1396,0x7f0f,0x3f2a,
  686. 0x6b20,0xeb27,0x6e1e,0x3ee9,
  687. 0xdbeb,0x743c,0x8352,0xbe83,
  688. };
  689. short Cnm[] = {
  690. 0x2d4f,0xbfb9,0x630d,0x3ef7,
  691. 0xd5a6,0x43d2,0x1998,0x3f00,
  692. 0x9b68,0x74f4,0x80e2,0x3ed4,
  693. 0xd94e,0xe02c,0x17b4,0x3eb8,
  694. 0x0276,0x0d60,0x1b1b,0xbede,
  695. 0x98f0,0x3391,0x26af,0xbeb8,
  696. 0x6738,0x398d,0xf20f,0xbe76,
  697. 0x72ec,0x0db0,0xca75,0xbe84,
  698. };
  699. short Snm[] = {
  700. 0x0000,0x0000,0x0000,0x0000,
  701. 0x0433,0xeddb,0x8868,0x3ed7,
  702. 0xf086,0x7b25,0x4f75,0x3ebc,
  703. 0x2d26,0x8009,0x702a,0xbe96,
  704. 0x9be9,0xd6f2,0xb990,0x3ec8,
  705. 0x0bce,0x58aa,0x3224,0xbec8,
  706. 0x7ebf,0x7803,0x792d,0xbeaa,
  707. 0x4c3e,0xcc08,0x4820,0x3e6e,
  708. };
  709. short K2M[] = {0x9b68,0x42f9,0xbfce,0x3f96};
  710. short LOVENO[] = {0x3333,0x3333,0x3333,0x3fd3};
  711. short PHASE[] = {0x3ed7,0xc8f9,0xd6a1,0x3fa4};
  712. short PSLP1[] = {0x68ee,0xc875,0x6fb4,0x3fcd};
  713. short PSLPI[] = {0x4e1f,0x3680,0x5253,0x403b};
  714. short PSLPIA[] = {0x0000,0x3400,0x5253,0x403b};
  715. short PSLPIB[] = {0x397b,0xfa49,0x0270,0x3e84};
  716. short JDEPOCH[] = {0x0000,0x4000,0x9e68,0x4142};
  717. short LBET[] = {0x2ad9,0x0445,0xb2f9,0x3f44};
  718. short LGAM[] = {0x44a6,0xa481,0xe327,0x3f2d};
  719. short LGAMBET[] = {0x94f0,0x7637,0x745e,0xbf3a};
  720. short CMRSQ[] = {0x1755,0xa468,0x010e,0x3fd9};
  721. short CMR2[] = {0x0314,0xd853,0xfd79,0x3dcc};
  722. short BMR2[] = {0xc125,0x685d,0xfa7b,0x3dcc};
  723. short AMR2[] = {0xdbbe,0x3056,0xf8ca,0x3dcc};
  724. #endif
  725. #if MIEEE
  726. short DTR[] = {0x3f91,0xdf46,0xa252,0x9d39};
  727. short RTD[] = {0x404c,0xa5dc,0x1a63,0xc1f8};
  728. short RTS[] = {0x4109,0x2dc6,0x7331,0xad7b};
  729. short STR[] = {0x3ed4,0x55a5,0xb2ff,0x8f9d};
  730. short PI[] = {0x4009,0x21fb,0x5444,0x2d18};
  731. short J2000[] = {0x4142,0xb42c,0x8000,0x0000};
  732. short B1950[] = {0x4142,0x9081,0x3624,0xdd2f};
  733. short J1900[] = {0x4142,0x6cd6,0x0000,0x0000};
  734. short B1900[] = {0x4142,0x6cd6,0x2810,0x624e};
  735. short Jcentury[] = {0x40e1,0xd5a0,0x0000,0x0000};
  736. short Jmillenium[] = {0x4116,0x4b08,0x0000,0x0000};
  737. short Jtmillenium[] = {0x414b,0xddca,0x0000,0x0000};
  738. short Zero[] = {0x0000,0x0000,0x0000,0x0000};
  739. short Fourth[] = {0x3fd0,0x0000,0x0000,0x0000};
  740. short Half[] = {0x3fe0,0x0000,0x0000,0x0000};
  741. short One[] = {0x3ff0,0x0000,0x0000,0x0000};
  742. short OneandaHalf[] = {0x3ff8,0x0000,0x0000,0x0000};
  743. short Two[] = {0x4000,0x0000,0x0000,0x0000};
  744. short Three[] = {0x4008,0x0000,0x0000,0x0000};
  745. short ThreeandaHalf[] = {0x400c,0x0000,0x0000,0x0000};
  746. short Four[] = {0x4010,0x0000,0x0000,0x0000};
  747. short Five[] = {0x4014,0x0000,0x0000,0x0000};
  748. short Six[] = {0x4018,0x0000,0x0000,0x0000};
  749. short Ten[] = {0x4024,0x0000,0x0000,0x0000};
  750.  
  751. short JD0[] = {0x4142,0x9e68,0x4000,0x0000};
  752. short yn0[4*6*NTOTAL] = {
  753. 0x3f1a,0x37ed,0x3903,0xa479,
  754. 0x3f78,0xb682,0x9614,0xd720,
  755. 0x3eef,0x4184,0x77ee,0xc10e,
  756. 0x3fd8,0x7920,0x8513,0x4406,
  757. 0xbf1f,0x4f19,0x25d8,0x7f80,
  758. 0x3ff4,0x7fd0,0x2618,0xcf9f,
  759. 0x3f6e,0x614d,0x8640,0x0dd3,
  760. 0x3fd7,0x0f43,0x8e0e,0xad8d,
  761. 0x3f99,0x7394,0x0c32,0xbbc3,
  762. 0xbfb8,0x45a8,0x1cd8,0x787c,
  763. 0x3f8a,0x7a93,0xa313,0x5fd5,
  764. 0xbfb6,0x643a,0x11c9,0xb8df,
  765. 0x3f86,0xd94b,0xa3b4,0x18a2,
  766. 0x3fe3,0x73a0,0x9906,0xd806,
  767. 0x3f8f,0xbb5b,0x9ce7,0x1500,
  768. 0xbfd6,0xbb54,0x18f8,0xef55,
  769. 0x3f79,0xb652,0x2b4f,0xfd65,
  770. 0xbfc9,0x5ff9,0xa358,0x657d,
  771. 0x3f91,0x3cbc,0x5a77,0x467c,
  772. 0x3fbb,0xb37c,0x3941,0xdb72,
  773. 0x3f59,0x8ff7,0xe0d6,0x6d8d,
  774. 0xbfed,0xaae4,0xe9c0,0x8a59,
  775. 0x3f46,0x2a92,0x8686,0x60c3,
  776. 0xbfd9,0xbbe6,0x4916,0x4052,
  777. 0x3f8d,0xa8b2,0x4e27,0xc307,
  778. 0xbfc0,0x6120,0x8e75,0xde62,
  779. 0x3f15,0x1c2d,0xc14f,0xfa3e,
  780. 0xbff5,0x385e,0x4b33,0x990b,
  781. 0xbf37,0x0fb9,0x145b,0xb969,
  782. 0xbfe3,0x5eaf,0x1074,0x2af2,
  783. 0x3f50,0x78b9,0x3fff,0x86d9,
  784. 0xc015,0x8f08,0xead6,0xf03b,
  785. 0xbf7a,0xbf09,0x4115,0xf544,
  786. 0xbfe8,0xa603,0xd2e1,0xfa7f,
  787. 0xbf67,0x264c,0x1a7c,0x0667,
  788. 0xbfc9,0x6dd5,0x08ba,0xa25a,
  789. 0xbf69,0xe1e1,0x86fc,0x39d8,
  790. 0x401f,0xcfa4,0xbc40,0x867d,
  791. 0x3f71,0xe7cf,0xc754,0xdfa5,
  792. 0x4012,0x0812,0x46b9,0x00c3,
  793. 0x3f5f,0xda4f,0xe937,0x635e,
  794. 0x3ff8,0x52b8,0x8b0c,0x48d8,
  795. 0x3f26,0x6c98,0xaff3,0xc0dd,
  796. 0xc032,0x473a,0x834e,0xeffd,
  797. 0xbf6e,0xd716,0x1153,0xdb5c,
  798. 0xbfee,0xa508,0xa7cd,0x3d5a,
  799. 0xbf5b,0x10da,0xfbfa,0x3d83,
  800. 0xbfc4,0xa632,0x36d8,0x7aaa,
  801. 0x3f65,0x7bd5,0x0001,0xdebd,
  802. 0xc030,0x5e00,0x40bd,0x4b76,
  803. 0xbf59,0x07d6,0x2f91,0x657e,
  804. 0xc037,0xc2ca,0x20b7,0x1c61,
  805. 0xbf46,0xab81,0xdad4,0xebac,
  806. 0xc022,0xa48c,0xc96a,0xeea9,
  807. 0x3f32,0x7771,0x7ea6,0x6d46,
  808. 0xc03e,0x729b,0x2c57,0x0471,
  809. 0xbf69,0xc7a6,0xcaa4,0x738a,
  810. 0xbfe1,0x0456,0x2067,0x676f,
  811. 0xbf51,0xaf6e,0x7c34,0x52c0,
  812. 0x4022,0x1e8b,0x8eb3,0x2af1,
  813. 0x3f43,0x9eb2,0x622d,0x7ce0,
  814. 0xbf4b,0x6262,0xe133,0x89e4,
  815. 0xbf26,0xd3a1,0x0216,0x4f4f,
  816. 0xbf60,0x43c4,0x8803,0x9bd6,
  817. 0xbf17,0x319a,0xa840,0x7dd2,
  818. 0xbf51,0xbf8e,0xfc4a,0x0681,
  819. 0xbe93,0x09d5,0x71ec,0x48a5,
  820. 0x3f72,0x7db4,0x3593,0x6a9d,
  821. 0x3ed5,0xbb43,0xdfe0,0x746b,
  822. 0x3f47,0xaf86,0x31ad,0x49f8,
  823. 0x3ec2,0xb64d,0xa9ae,0xfbb2,
  824. 0x3f30,0x291c,0x4d38,0x8354,
  825. };
  826. short JD1[] = {0x4142,0x9f30,0x4000,0x0000};
  827. short yn1[4*6*NTOTAL] = {
  828. 0x3f1f,0x678a,0xba4a,0xa4f1,
  829. 0x3fa0,0xbac5,0x0a33,0x4c98,
  830. 0xbf11,0xfac1,0x3b24,0x409d,
  831. 0x3fd8,0xa805,0xea01,0x776e,
  832. 0xbf19,0xe5d4,0x655e,0x42dd,
  833. 0x3ff4,0x185d,0xfa9f,0x75ea,
  834. 0x3f87,0xd6be,0x66c8,0x110e,
  835. 0xbfd5,0xfbf3,0xe5f1,0xb307,
  836. 0xbf92,0x6562,0x6579,0x8750,
  837. 0xbfd0,0x137b,0x75c9,0x79b6,
  838. 0xbf86,0x2762,0xe7fb,0xfba7,
  839. 0xbfb9,0x6343,0x7689,0xeab8,
  840. 0x3f93,0x0b2a,0x7b70,0x5d56,
  841. 0xbfd1,0x39b1,0x67e6,0x8633,
  842. 0xbf7a,0xfa6c,0xbdd4,0xda89,
  843. 0xbfe3,0xbad4,0xf193,0xc222,
  844. 0xbf70,0xf818,0x4c64,0x7407,
  845. 0xbfd0,0xaf6d,0xb8bd,0xc2b1,
  846. 0x3f8a,0xcb90,0x8add,0x4358,
  847. 0x3fe4,0x9032,0x64d2,0x94c2,
  848. 0x3f84,0x3e61,0xff72,0xf4e0,
  849. 0xbfe7,0x1116,0x9a32,0x8df8,
  850. 0x3f71,0x8e7d,0x738d,0x228f,
  851. 0xbfd4,0x0244,0x7c73,0x8045,
  852. 0xbf85,0x1281,0x6cc4,0x1a91,
  853. 0xbff0,0xa1c4,0x9c3e,0x92d8,
  854. 0xbf7c,0xf7e1,0x1159,0x8a30,
  855. 0x3ff2,0x6f93,0x9ccd,0x39d3,
  856. 0xbf68,0x5533,0x1867,0xf56c,
  857. 0x3fe1,0xd117,0xbb5e,0xc57a,
  858. 0x3f72,0xef6c,0xf903,0xdec3,
  859. 0xc010,0xf30e,0xe978,0x50f0,
  860. 0xbf74,0xb52c,0x1bd8,0xafa7,
  861. 0xc009,0x2bfa,0x32c3,0x8613,
  862. 0xbf62,0xb160,0x836f,0x4682,
  863. 0xbff3,0xf088,0xc040,0x9f24,
  864. 0xbf71,0x6e65,0xcabf,0x53cd,
  865. 0x4019,0xda84,0xfe3d,0x36af,
  866. 0x3f6d,0x30ca,0xa73d,0x30d1,
  867. 0x4018,0x6a3f,0x1f17,0xb53a,
  868. 0x3f5b,0x278a,0xd8ad,0x3d27,
  869. 0x4001,0xf4a8,0x6968,0x35db,
  870. 0x3f41,0x2120,0x489a,0xee4c,
  871. 0xc032,0x23ab,0xd871,0xb546,
  872. 0xbf6e,0x9209,0xc3dd,0xaca3,
  873. 0xc003,0xa9c1,0x9e48,0x0bf2,
  874. 0xbf5a,0xe92d,0x5f5f,0xdc66,
  875. 0xbfea,0x42ad,0xa034,0x73c0,
  876. 0x3f66,0x090a,0xfa60,0xde76,
  877. 0xc02e,0x9bef,0xef47,0xad6e,
  878. 0xbf57,0x5ab8,0xf97f,0x5f1e,
  879. 0xc038,0x5a03,0xafed,0x937f,
  880. 0xbf45,0x5a29,0xa8d3,0xfb0f,
  881. 0xc023,0x2e23,0x8e36,0x9d22,
  882. 0x3f39,0xe60d,0xd197,0x162a,
  883. 0xc03e,0x4ff5,0xd74a,0xefdb,
  884. 0xbf69,0xbe90,0x0cda,0xfa3f,
  885. 0xbffc,0xa33f,0x6782,0x4cf0,
  886. 0xbf52,0x39e4,0x0fdd,0x0ed3,
  887. 0x4021,0x3e14,0xf1ce,0xdfb3,
  888. 0xbf3e,0x0ad3,0x343d,0xf74e,
  889. 0xbf5a,0x3c59,0x06a1,0x1bd1,
  890. 0xbf32,0x43c2,0xc806,0xd241,
  891. 0x3f60,0x0890,0x41cf,0xac3c,
  892. 0xbf26,0x5fa6,0x4a6c,0x5f17,
  893. 0x3f4f,0xbc39,0x7caa,0x6280,
  894. 0xbecc,0xe242,0xc604,0x5695,
  895. 0x3f6e,0xf6e0,0xd624,0x2615,
  896. 0x3ed0,0xedad,0xa4db,0xd7aa,
  897. 0x3f65,0x8cae,0x3244,0xa17d,
  898. 0x3ebe,0x5dfe,0x281a,0x7a70,
  899. 0x3f51,0xb9bb,0xdfe4,0x486d,
  900. };
  901.  
  902. short CLIGHT[] = {0x4112,0x4c41,0xd4fd,0xf3b6};
  903. short C[] = {0x4065,0xa4a0,0xd4cc,0x9453};
  904. short EMRAT[] = {0x4054,0x533c,0xd141,0xa693};
  905. short KG[] = {0x3f91,0x9d6d,0x51a6,0xb69a};
  906. short GMs[] = {
  907. 0x0000,0x0000,0x0000,0x0000,
  908. 0x3dcb,0x01cb,0xbfef,0x30a7,
  909. 0x3e08,0xe369,0x5596,0x9305,
  910. 0x3e0e,0x89b2,0x7bf8,0x614e,
  911. 0x3dda,0x3fe2,0x330d,0xe19b,
  912. 0x3e92,0xf5e6,0x991b,0x594a,
  913. 0x3e76,0xb54d,0x700f,0xa2cf,
  914. 0x3e4b,0xad57,0x9a62,0xddff,
  915. 0x3e50,0x7370,0x4ffd,0xf7c8,
  916. 0x3d84,0x05a7,0x7dd0,0xde3d,
  917. 0x3da8,0x0a1b,0xbc9f,0x5e14,
  918. 0x3f33,0x6495,0x41cf,0x95e6,
  919. 0x3d48,0x92a0,0xc1f4,0x5b6d,
  920. 0x3d22,0x03af,0x9ee7,0x5616,
  921. 0x3d26,0xf7e6,0x5100,0x8dc2,
  922. 0x3cdc,0xd2b2,0x97d8,0x89bc,
  923. 0x3ce7,0x6b31,0x1b5f,0xefe9,
  924. };
  925. short AU[] = {0x4241,0x6a5d,0x2d22,0x0000};
  926. short RADS[] = {0x4125,0x3d80,0x0000,0x0000};
  927. short RADM[] = {0x409b,0x2800,0x0000,0x0000};
  928. short RADE[] = {0x40b8,0xea23,0xd70a,0x3d71};
  929. short AE[] = {0x3f06,0x5a67,0x6f26,0x85dd};
  930. short Je[] = {
  931. 0x3f51,0xbce1,0x1c65,0x107b,
  932. 0xbec5,0x4e9b,0xc46c,0x2b21,
  933. 0xbebb,0x02e5,0xb881,0x1061,
  934. };
  935. short AM[] = {0x3ee8,0x5d43,0xfaf9,0xde5a};
  936. short Jm[] = {
  937. 0x3f2a,0x7f0f,0x1396,0x70b6,
  938. 0x3ee9,0x6e1e,0xeb27,0x6b20,
  939. 0xbe83,0x8352,0x743c,0xdbeb,
  940. };
  941. short Cnm[] = {
  942. 0x3ef7,0x630d,0xbfb9,0x2d4f,
  943. 0x3f00,0x1998,0x43d2,0xd5a6,
  944. 0x3ed4,0x80e2,0x74f4,0x9b68,
  945. 0x3eb8,0x17b4,0xe02c,0xd94e,
  946. 0xbede,0x1b1b,0x0d60,0x0276,
  947. 0xbeb8,0x26af,0x3391,0x98f0,
  948. 0xbe76,0xf20f,0x398d,0x6738,
  949. 0xbe84,0xca75,0x0db0,0x72ec,
  950. };
  951. short Snm[] = {
  952. 0x0000,0x0000,0x0000,0x0000,
  953. 0x3ed7,0x8868,0xeddb,0x0433,
  954. 0x3ebc,0x4f75,0x7b25,0xf086,
  955. 0xbe96,0x702a,0x8009,0x2d26,
  956. 0x3ec8,0xb990,0xd6f2,0x9be9,
  957. 0xbec8,0x3224,0x58aa,0x0bce,
  958. 0xbeaa,0x792d,0x7803,0x7ebf,
  959. 0x3e6e,0x4820,0xcc08,0x4c3e,
  960. };
  961. short K2M[] = {0x3f96,0xbfce,0x42f9,0x9b68};
  962. short LOVENO[] = {0x3fd3,0x3333,0x3333,0x3333};
  963. short PHASE[] = {0x3fa4,0xd6a1,0xc8f9,0x3ed7};
  964. short PSLP1[] = {0x3fcd,0x6fb4,0xc875,0x68ee};
  965. short PSLPI[] = {0x403b,0x5253,0x3680,0x4e1f};
  966. short PSLPIA[] = {0x403b,0x5253,0x3400,0x0000};
  967. short PSLPIB[] = {0x3e84,0x0270,0xfa49,0x397b};
  968. short JDEPOCH[] = {0x4142,0x9e68,0x4000,0x0000};
  969. short LBET[] = {0x3f44,0xb2f9,0x0445,0x2ad9};
  970. short LGAM[] = {0x3f2d,0xe327,0xa481,0x44a6};
  971. short LGAMBET[] = {0xbf3a,0x745e,0x7637,0x94f0};
  972. short CMRSQ[] = {0x3fd9,0x010e,0xa468,0x1755};
  973. short CMR2[] = {0x3dcc,0xfd79,0xd853,0x0314};
  974. short BMR2[] = {0x3dcc,0xfa7b,0x685d,0xc125};
  975. short AMR2[] = {0x3dcc,0xf8ca,0x3056,0xdbbe};
  976. #endif /* MIEEE */
  977. #if DEC
  978. short DTR[] = {0036616,0175065,0011224,0164711};
  979. short RTD[] = {0041545,0027340,0151436,0007676};
  980. short RTS[] = {0044511,0067063,0114615,0065726};
  981. short STR[] = {0033642,0126455,0113774,0076351};
  982. short PI[] = {0040511,0007732,0121041,0064302};
  983. short J2000[] = {0045425,0120544,0000000,0000000};
  984. short B1950[] = {0045424,0102011,0130446,0164571};
  985. short J1900[] = {0045423,0063260,0000000,0000000};
  986. short B1900[] = {0045423,0063261,0040203,0011157};
  987. short Jcentury[] = {0044016,0126400,0000000,0000000};
  988. short Jmillenium[] = {0044662,0054100,0000000,0000000};
  989. short Jtmillenium[] = {0045536,0167120,0000000,0000000};
  990. short Zero[] = {0000000,0000000,0000000,0000000};
  991. short Fourth[] = {0037600,0000000,0000000,0000000};
  992. short Half[] = {0040000,0000000,0000000,0000000};
  993. short One[] = {0040200,0000000,0000000,0000000};
  994. short OneandaHalf[] = {0040300,0000000,0000000,0000000};
  995. short Two[] = {0040400,0000000,0000000,0000000};
  996. short Three[] = {0040500,0000000,0000000,0000000};
  997. short ThreeandaHalf[] = {0040540,0000000,0000000,0000000};
  998. short Four[] = {0040600,0000000,0000000,0000000};
  999. short Five[] = {0040640,0000000,0000000,0000000};
  1000. short Six[] = {0040700,0000000,0000000,0000000};
  1001. short Ten[] = {0041040,0000000,0000000,0000000};
  1002.  
  1003. short JD0[] = {0045424,0171502,0000000,0000000};
  1004. short yn0[4*6*NTOTAL] = {
  1005. 0034721,0137551,0144035,0021713,
  1006. 0036305,0132024,0130246,0134404,
  1007. 0034172,0006043,0137566,0004160,
  1008. 0037703,0144404,0024232,0020055,
  1009. 0134772,0074311,0027303,0176000,
  1010. 0040243,0177201,0030306,0076374,
  1011. 0036163,0005154,0031000,0067224,
  1012. 0037670,0075034,0070165,0066146,
  1013. 0036713,0116240,0060625,0157031,
  1014. 0137302,0026500,0163303,0141735,
  1015. 0036523,0152235,0014232,0177246,
  1016. 0137263,0020720,0107115,0143366,
  1017. 0036466,0145135,0016640,0142422,
  1018. 0040033,0116404,0144066,0140060,
  1019. 0036575,0155334,0163470,0124001,
  1020. 0137665,0155240,0143707,0075252,
  1021. 0036315,0131221,0055177,0165446,
  1022. 0137512,0177715,0015303,0025746,
  1023. 0036611,0162742,0151672,0031742,
  1024. 0037335,0115741,0145016,0155617,
  1025. 0035714,0077677,0003263,0066153,
  1026. 0140155,0053447,0047004,0051311,
  1027. 0035461,0052224,0032063,0003026,
  1028. 0137715,0157462,0044262,0001222,
  1029. 0036555,0042622,0070476,0014073,
  1030. 0137403,0004404,0071656,0171423,
  1031. 0034650,0160556,0005177,0150761,
  1032. 0140251,0141362,0054634,0144127,
  1033. 0135270,0076710,0121335,0145512,
  1034. 0140032,0172570,0101641,0053621,
  1035. 0035603,0142711,0177774,0033314,
  1036. 0140654,0074107,0053267,0100732,
  1037. 0136325,0174112,0004257,0125042,
  1038. 0140105,0030036,0113417,0151767,
  1039. 0136071,0031140,0151740,0031466,
  1040. 0137513,0067250,0042725,0011320,
  1041. 0136117,0007414,0033741,0147276,
  1042. 0040776,0076445,0161004,0031746,
  1043. 0036217,0037176,0035246,0176452,
  1044. 0040620,0040222,0032710,0003031,
  1045. 0035776,0151177,0044673,0015363,
  1046. 0040302,0112704,0054142,0043277,
  1047. 0035063,0062305,0077636,0003353,
  1048. 0141222,0034724,0015167,0077745,
  1049. 0136166,0134260,0105236,0155337,
  1050. 0140165,0024105,0037151,0165317,
  1051. 0135730,0103327,0157721,0166030,
  1052. 0137445,0030621,0133303,0152523,
  1053. 0036053,0157250,0000016,0172753,
  1054. 0141202,0170002,0002752,0055657,
  1055. 0135710,0037261,0076213,0025756,
  1056. 0141276,0013121,0002670,0161411,
  1057. 0135465,0056016,0153247,0056536,
  1058. 0141025,0022146,0045527,0072510,
  1059. 0035223,0135613,0172463,0065057,
  1060. 0141363,0112331,0061270,0021607,
  1061. 0136116,0036466,0052443,0116123,
  1062. 0140010,0021261,0001473,0035572,
  1063. 0135615,0075563,0160642,0113003,
  1064. 0041020,0172134,0072631,0053605,
  1065. 0035434,0172623,0010553,0163377,
  1066. 0135533,0011427,0004634,0047444,
  1067. 0135066,0116410,0010262,0075172,
  1068. 0136002,0017044,0040034,0157263,
  1069. 0134671,0106325,0041003,0167217,
  1070. 0135615,0176167,0161120,0032011,
  1071. 0132630,0047253,0107542,0042451,
  1072. 0036223,0166641,0126233,0052345,
  1073. 0033655,0155036,0177403,0121527,
  1074. 0035475,0076061,0106552,0047675,
  1075. 0033425,0131155,0046567,0156615,
  1076. 0035201,0044342,0064704,0015237,
  1077. };
  1078. short JD1[] = {0045424,0174602,0000000,0000000};
  1079. short yn1[4*6*NTOTAL] = {
  1080. 0034773,0036125,0151125,0023612,
  1081. 0037005,0153050,0050632,0062276,
  1082. 0134617,0153011,0154442,0002351,
  1083. 0037705,0040057,0050013,0135554,
  1084. 0134717,0027243,0025362,0013345,
  1085. 0040240,0141357,0152373,0127515,
  1086. 0036476,0132763,0033100,0104156,
  1087. 0137657,0157637,0027615,0114067,
  1088. 0136623,0025423,0025714,0035174,
  1089. 0137600,0115733,0127113,0146662,
  1090. 0136461,0035427,0037737,0156470,
  1091. 0137313,0015033,0132117,0052700,
  1092. 0036630,0054523,0155602,0165264,
  1093. 0137611,0146613,0037464,0030630,
  1094. 0136327,0151545,0167246,0152106,
  1095. 0140035,0153247,0106236,0010422,
  1096. 0136207,0140302,0061443,0120073,
  1097. 0137605,0075555,0142756,0012611,
  1098. 0036526,0056204,0053352,0015303,
  1099. 0040044,0100623,0023224,0123021,
  1100. 0036441,0171417,0175627,0123401,
  1101. 0140070,0104264,0150624,0067701,
  1102. 0036214,0071753,0116151,0012173,
  1103. 0137640,0011043,0161634,0001044,
  1104. 0136450,0112013,0063040,0152204,
  1105. 0140205,0007044,0160764,0113275,
  1106. 0136347,0137410,0105314,0050577,
  1107. 0040223,0076234,0163151,0147234,
  1108. 0136102,0124630,0141477,0125540,
  1109. 0040016,0104275,0155366,0025722,
  1110. 0036227,0075547,0144036,0173025,
  1111. 0140607,0114167,0045702,0103601,
  1112. 0136245,0124540,0157305,0076473,
  1113. 0140511,0057721,0113034,0030231,
  1114. 0136025,0105404,0015572,0032022,
  1115. 0140237,0102106,0001004,0174444,
  1116. 0136213,0071456,0052772,0117152,
  1117. 0040716,0152047,0170751,0132573,
  1118. 0036151,0103125,0034751,0103204,
  1119. 0040703,0050770,0174275,0124723,
  1120. 0035731,0036126,0142551,0164473,
  1121. 0040417,0122503,0045501,0127324,
  1122. 0035411,0004402,0042327,0071143,
  1123. 0141221,0016536,0141615,0125056,
  1124. 0136164,0110116,0017355,0062433,
  1125. 0140435,0047014,0171100,0057615,
  1126. 0135727,0044552,0175376,0161460,
  1127. 0140122,0012555,0000643,0116774,
  1128. 0036060,0044127,0151406,0171660,
  1129. 0141164,0157577,0075075,0065563,
  1130. 0135672,0152707,0145772,0174362,
  1131. 0141302,0150035,0077554,0115770,
  1132. 0135452,0150515,0043237,0154166,
  1133. 0141031,0070434,0070664,0164414,
  1134. 0035317,0030156,0106270,0130522,
  1135. 0141362,0077656,0135127,0077334,
  1136. 0136115,0172200,0063327,0150774,
  1137. 0140345,0014773,0036022,0063604,
  1138. 0135621,0147440,0077350,0073226,
  1139. 0041011,0170247,0107166,0176634,
  1140. 0135360,0053231,0120757,0135156,
  1141. 0135721,0161310,0032410,0157214,
  1142. 0135222,0017026,0040066,0111005,
  1143. 0036000,0042202,0007175,0060741,
  1144. 0135062,0176462,0051542,0174267,
  1145. 0035575,0160713,0162523,0012000,
  1146. 0133547,0011026,0030042,0132250,
  1147. 0036167,0133406,0130441,0030251,
  1148. 0033607,0066555,0023336,0136516,
  1149. 0036054,0062561,0111045,0005747,
  1150. 0033362,0167761,0040323,0151603,
  1151. 0035615,0146736,0177442,0041544,
  1152. };
  1153.  
  1154. short CLIGHT[] = {0044622,0061016,0123757,0116662};
  1155. short C[] = {0042055,0022406,0123144,0121226};
  1156. short EMRAT[] = {0041642,0114746,0105015,0032234};
  1157. short KG[] = {0036614,0165552,0106465,0132315};
  1158. short GMs[] = {
  1159. 0000000,0000000,0000000,0000000,
  1160. 0027530,0007135,0177571,0102467,
  1161. 0030507,0015512,0126264,0114046,
  1162. 0030564,0046623,0157703,0005155,
  1163. 0027721,0177421,0114157,0006325,
  1164. 0032627,0127464,0144332,0145115,
  1165. 0032265,0125153,0100175,0013174,
  1166. 0031535,0065274,0151426,0167765,
  1167. 0031603,0115602,0077757,0137100,
  1168. 0026440,0026473,0167206,0170745,
  1169. 0027100,0050335,0162372,0170240,
  1170. 0035233,0022252,0007174,0127461,
  1171. 0025504,0112406,0007642,0155545,
  1172. 0025020,0016574,0173472,0130255,
  1173. 0025067,0137462,0104004,0067020,
  1174. 0023746,0112624,0137304,0046741,
  1175. 0024073,0054610,0155377,0077507,
  1176. };
  1177. short AU[] = {0051413,0051351,0064420,0000000};
  1178. short RADS[] = {0045051,0166000,0000000,0000000};
  1179. short RADM[] = {0042731,0040000,0000000,0000000};
  1180. short RADE[] = {0043307,0050436,0134121,0165605};
  1181. short AE[] = {0034462,0151473,0074464,0027350};
  1182. short Je[] = {
  1183. 0035615,0163410,0161450,0101730,
  1184. 0133452,0072336,0021541,0054413,
  1185. 0133330,0013455,0142010,0101406,
  1186. };
  1187.  
  1188. short AM[] = {0034102,0165037,0153716,0171321};
  1189.  
  1190. short Jm[] = {
  1191. 0035123,0174170,0116263,0102663,
  1192. 0034113,0070367,0054473,0054401,
  1193. 0132434,0015223,0120746,0157531,
  1194. };
  1195. short Cnm[] = {
  1196. 0034273,0014155,0176711,0065164,
  1197. 0034400,0146302,0017226,0126456,
  1198. 0033644,0003423,0123644,0155476,
  1199. 0033300,0136647,0000546,0145162,
  1200. 0133760,0154330,0065400,0011661,
  1201. 0133301,0032571,0116214,0143603,
  1202. 0132267,0110171,0146153,0034675,
  1203. 0132446,0051650,0066603,0113540,
  1204. };
  1205. short Snm[] = {
  1206. 0000000,0000000,0000000,0000000,
  1207. 0033674,0041507,0067330,0020625,
  1208. 0033342,0075653,0154457,0102055,
  1209. 0132663,0100524,0000111,0064455,
  1210. 0033505,0146206,0133624,0157510,
  1211. 0133501,0110442,0142520,0057164,
  1212. 0133123,0144553,0140033,0172770,
  1213. 0032162,0040406,0060102,0060761,
  1214. };
  1215. short K2M[] = {0036665,0177162,0013714,0155474};
  1216. short LOVENO[] = {0037631,0114631,0114631,0114632};
  1217. short PHASE[] = {0037046,0132416,0043711,0173267};
  1218. short PSLP1[] = {0037553,0076646,0041653,0043554};
  1219. short PSLPI[] = {0041332,0111231,0132002,0070372};
  1220. short PSLPIA[] = {0041332,0111231,0120000,0000000};
  1221. short PSLPIB[] = {0032440,0011607,0151111,0145731};
  1222. short JDEPOCH[] = {0045424,0171502,0000000,0000000};
  1223. short LBET[] = {0035445,0113710,0021051,0053304};
  1224. short LGAM[] = {0035157,0014475,0022012,0022456};
  1225. short LGAMBET[] = {0135323,0121363,0130674,0123574};
  1226. short CMRSQ[] = {0037710,0004165,0021500,0135246};
  1227. short CMR2[] = {0027547,0165716,0141230,0014243};
  1228. short BMR2[] = {0027547,0151733,0041356,0004454};
  1229. short AMR2[] = {0027547,0143121,0101266,0156762};
  1230. #endif /* DEC */
  1231.