home *** CD-ROM | disk | FTP | other *** search
/ GEMini Atari / GEMini_Atari_CD-ROM_Walnut_Creek_December_1993.iso / files / astronmy / strchart / patch / planet.dif < prev    next >
Encoding:
Text File  |  1989-03-23  |  17.6 KB  |  593 lines
  1. *** ../starchart/planet.c    Mon Sep 19 09:11:43 1988
  2. --- planet.c    Mon Sep 26 08:00:31 1988
  3. ***************
  4. *** 45,57 ****
  5.   #include <stdio.h>
  6.   #include <math.h>
  7.   
  8.   #ifndef SYSV
  9.   #include <sys/time.h>    /* for getting current GMT (generic Unix) */
  10.   #else
  11.   #include <time.h>    /* for getting current GMT (sysV version) */
  12.   #endif
  13.   
  14. ! #define CURYEAR    1987    /* default year -- needs to be maintained */
  15.   
  16.   #ifndef PLANETFILE
  17.   #define PLANETFILE "./planet.star"
  18. --- 45,62 ----
  19.   #include <stdio.h>
  20.   #include <math.h>
  21.   
  22. + #ifdef ST_MWC /* Mark Williams C for the Atari ST */
  23. + #define SYSV /* MWC is almost SYSV compatible */
  24. + #include <osbind.h> /* Needed for Cconws() and Crawcin() */
  25. + #endif /* ST_MWC */
  26. + 
  27.   #ifndef SYSV
  28.   #include <sys/time.h>    /* for getting current GMT (generic Unix) */
  29.   #else
  30.   #include <time.h>    /* for getting current GMT (sysV version) */
  31.   #endif
  32.   
  33. ! #define CURYEAR    1989    /* default year -- needs to be maintained */
  34.   
  35.   #ifndef PLANETFILE
  36.   #define PLANETFILE "./planet.star"
  37. ***************
  38. *** 68,73 ****
  39. --- 73,84 ----
  40.   #define MAGURA 590
  41.   #define MAGNEP 800
  42.   
  43. + #ifdef ST_MWC
  44. + #define daylight tt.tm_isdst
  45. + extern char *getenv(); /* used to see if we're running from the desktop */
  46. + struct tm_t *tt;
  47. + #endif /* ST_MWC */
  48. + 
  49.   double pie, rad;
  50.   double htod(), atof(), kepler(), truean();
  51.   double longi(), lati(), poly(), aint(), range();
  52. ***************
  53. *** 119,125 ****
  54. --- 130,138 ----
  55.       gettimeofday(&tv, &tz);
  56.       return(GMT1970 + tv.tv_sec/SECSPERDAY);
  57.   #else
  58. + #ifndef ST_MWC
  59.       long time();
  60. + #endif /* ST_MWC */
  61.       return(GMT1970 + (double)time((long *)0)/SECSPERDAY);
  62.   #endif
  63.       }
  64. ***************
  65. *** 201,207 ****
  66.       aa1 = (int) (year/100);
  67.       bb1 = 2 - aa1 + (int)(aa1/4);
  68.       jd = aint(365.25*year) + aint(30.6001*(month + 1));
  69. !     jd = jd  + day + 1720994.5;
  70.       if((year + month/100) > 1582.10) jd = jd + bb1;
  71.       return(jd);
  72.       }
  73. --- 214,220 ----
  74.       aa1 = (int) (year/100);
  75.       bb1 = 2 - aa1 + (int)(aa1/4);
  76.       jd = aint(365.25*year) + aint(30.6001*(month + 1));
  77. !     jd += day + 1720994.5;
  78.       if((year + month/100) > 1582.10) jd = jd + bb1;
  79.       return(jd);
  80.       }
  81. ***************
  82. *** 218,224 ****
  83. --- 231,243 ----
  84.       double N,D,epli,thapp,omeg;
  85.       double nu2,P,Q,S,V,W,ze,l1pert,epert,w1pert,apert;
  86.       double psi,H,G,eta,th;
  87. + #ifdef ST_MWC
  88. +     time_t now;
  89.   
  90. +     now = time((long) 0);
  91. +     tt = localtime(&now);
  92. + #endif /* ST_MWC */
  93. + 
  94.   #define WRITEMODE "w"
  95.   #define OPENFAIL 0
  96.       progname = argv[0];
  97. ***************
  98. *** 321,327 ****
  99.           RA = atan2(N,D)/rad;
  100.           DEC = asin(sin(epli*rad)*sin(thapp*rad))/rad;
  101.           speak(RA,DEC,Sr, MAGSOL, "PS", "Sol");
  102. ! /* tansformation of coordinates on Mercury and output */
  103.           trans(r,b,ll,Stheta,Sr,epli, MAGMER, "PM", "Mercury");
  104.           
  105.   /* Now start on Venus */
  106. --- 340,346 ----
  107.           RA = atan2(N,D)/rad;
  108.           DEC = asin(sin(epli*rad)*sin(thapp*rad))/rad;
  109.           speak(RA,DEC,Sr, MAGSOL, "PS", "Sol");
  110. ! /* transformation of coordinates on Mercury and output */
  111.           trans(r,b,ll,Stheta,Sr,epli, MAGMER, "PM", "Mercury");
  112.           
  113.   /* Now start on Venus */
  114. ***************
  115. *** 373,380 ****
  116.       radpert = 0.000022501 * cos((2*M-2*M2-58.208)*rad)
  117.            +0.000019045 * cos((3*M-3*M2+92.577)*rad)
  118.            +0.000006887 * cos((M5-M2-118.090)*rad)
  119. !          +0.000005172 * cos((M-M2-29.110)*rad)
  120. !          +0.000003620 * cos((5*M-4*M2-104.208)*rad)
  121.            +0.000003283 * cos((4*M-4*M2+63.513)*rad)
  122.            +0.000003074 * cos((2*M5-2*M2-55.167)*rad);
  123.       r = r + radpert;
  124. --- 392,399 ----
  125.       radpert = 0.000022501 * cos((2*M-2*M2-58.208)*rad)
  126.            +0.000019045 * cos((3*M-3*M2+92.577)*rad)
  127.            +0.000006887 * cos((M5-M2-118.090)*rad)
  128. !          +0.000005172 * cos((M-M2-29.110)*rad);
  129. !     radpert+= 0.000003620 * cos((5*M-4*M2-104.208)*rad)
  130.            +0.000003283 * cos((4*M-4*M2+63.513)*rad)
  131.            +0.000003074 * cos((2*M5-2*M2-55.167)*rad);
  132.       r = r + radpert;
  133. ***************
  134. *** 425,432 ****
  135.       lonpert = 0.00705*cos((M5-M4-48.958)*rad)
  136.            +0.00607*cos((2*M5-M4-188.350)*rad)
  137.            +0.00445*cos((2*M5-2*M4-191.897)*rad)
  138. !          +0.00388*cos((M-2*M4+20.495)*rad)
  139. !          +0.00238*cos((M-M4+35.097)*rad)
  140.            +0.00204*cos((2*M-3*M4+158.638)*rad)
  141.            +0.00177*cos((3*M4-M2-57.602)*rad)
  142.            +0.00136*cos((2*M-4*M4+154.093)*rad)
  143. --- 444,451 ----
  144.       lonpert = 0.00705*cos((M5-M4-48.958)*rad)
  145.            +0.00607*cos((2*M5-M4-188.350)*rad)
  146.            +0.00445*cos((2*M5-2*M4-191.897)*rad)
  147. !          +0.00388*cos((M-2*M4+20.495)*rad);
  148. !     lonpert+= 0.00238*cos((M-M4+35.097)*rad)
  149.            +0.00204*cos((2*M-3*M4+158.638)*rad)
  150.            +0.00177*cos((3*M4-M2-57.602)*rad)
  151.            +0.00136*cos((2*M-4*M4+154.093)*rad)
  152. ***************
  153. *** 445,452 ****
  154.       radpert+= 0.000007914*cos((3*M5-2*M4-139.737)*rad)
  155.            +0.000007004*cos((2*M5-3*M4-102.888)*rad)
  156.            +0.000006620*cos((M-2*M4+113.202)*rad)
  157. !          +0.000004930*cos((3*M5-3*M4-76.243)*rad)
  158. !          +0.000004693*cos((3*M-5*M4+190.603)*rad)
  159.            +0.000004571*cos((2*M-4*M4+244.702)*rad)
  160.            +0.000004409*cos((3*M5-M4-115.828)*rad);
  161.       r = r + radpert;
  162. --- 464,471 ----
  163.       radpert+= 0.000007914*cos((3*M5-2*M4-139.737)*rad)
  164.            +0.000007004*cos((2*M5-3*M4-102.888)*rad)
  165.            +0.000006620*cos((M-2*M4+113.202)*rad)
  166. !          +0.000004930*cos((3*M5-3*M4-76.243)*rad);
  167. !     radpert+= 0.000004693*cos((3*M-5*M4+190.603)*rad)
  168.            +0.000004571*cos((2*M-4*M4+244.702)*rad)
  169.            +0.000004409*cos((3*M5-M4-115.828)*rad);
  170.       r = r + radpert;
  171. ***************
  172. *** 507,513 ****
  173.       -(0.003083 + 0.000275*nu2 - 0.000489*nu2*nu2)*sin(2*V)
  174.       +0.002472*sin(W)
  175.       +0.013619*sin(ze)
  176. !     +0.018472*sin(2*ze)
  177.       +0.006717*sin(3*ze)
  178.       +0.002775*sin(4*ze)
  179.       +(0.007275 - 0.001253*nu2)*sin(ze)*sin(Q)
  180. --- 526,533 ----
  181.       -(0.003083 + 0.000275*nu2 - 0.000489*nu2*nu2)*sin(2*V)
  182.       +0.002472*sin(W)
  183.       +0.013619*sin(ze)
  184. !     +0.018472*sin(2*ze);
  185. ! l1pert = l1pert
  186.       +0.006717*sin(3*ze)
  187.       +0.002775*sin(4*ze)
  188.       +(0.007275 - 0.001253*nu2)*sin(ze)*sin(Q)
  189. ***************
  190. *** 519,525 ****
  191.       -(0.035681 + 0.001208*nu2)*sin(ze)*cos(Q)
  192.       -0.004261*sin(2*ze)*cos(Q)
  193.       +0.002178*cos(Q)
  194. !     +(-0.006333 + 0.001161*nu2)*cos(ze)*cos(Q)
  195.       -0.006675*cos(2*ze)*cos(Q)
  196.       -0.002664*cos(3*ze)*cos(Q)
  197.       -0.002572*sin(ze)*sin(2*Q)
  198. --- 539,546 ----
  199.       -(0.035681 + 0.001208*nu2)*sin(ze)*cos(Q)
  200.       -0.004261*sin(2*ze)*cos(Q)
  201.       +0.002178*cos(Q)
  202. !     +(-0.006333 + 0.001161*nu2)*cos(ze)*cos(Q);
  203. ! l1pert = l1pert
  204.       -0.006675*cos(2*ze)*cos(Q)
  205.       -0.002664*cos(3*ze)*cos(Q)
  206.       -0.002572*sin(ze)*sin(2*Q)
  207. ***************
  208. *** 532,538 ****
  209.       -.0006764*sin(ze)*sin(Q)
  210.       -.0001110*sin(2*ze)*sin(Q)
  211.       -.0000224*sin(3*ze)*sin(Q)
  212. !     -.0000204*sin(Q)
  213.       +(.0001284 + .0000116*nu2)*cos(ze)*sin(Q)
  214.       +.0000188*cos(2*ze)*sin(Q)
  215.       +(.0001460 + .0000130*nu2)*sin(ze)*cos(Q)
  216. --- 553,560 ----
  217.       -.0006764*sin(ze)*sin(Q)
  218.       -.0001110*sin(2*ze)*sin(Q)
  219.       -.0000224*sin(3*ze)*sin(Q)
  220. !     -.0000204*sin(Q);
  221. ! epert = epert
  222.       +(.0001284 + .0000116*nu2)*cos(ze)*sin(Q)
  223.       +.0000188*cos(2*ze)*sin(Q)
  224.       +(.0001460 + .0000130*nu2)*sin(ze)*cos(Q)
  225. ***************
  226. *** 544,550 ****
  227.       +.0000508*cos(3*ze)*cos(Q)
  228.       +.0000230*cos(4*ze)*cos(Q)
  229.       +.0000108*cos(5*ze)*cos(Q)
  230. !     -(.0000956 + .0000073*nu2)*sin(ze)*sin(2*Q)
  231.       +.0000448*sin(2*ze)*sin(2*Q)
  232.       +.0000137*sin(3*ze)*sin(2*Q)
  233.       +(-.0000997 + .0000108*nu2)*cos(ze)*sin(2*Q)
  234. --- 566,573 ----
  235.       +.0000508*cos(3*ze)*cos(Q)
  236.       +.0000230*cos(4*ze)*cos(Q)
  237.       +.0000108*cos(5*ze)*cos(Q)
  238. !     -(.0000956 + .0000073*nu2)*sin(ze)*sin(2*Q);
  239. ! epert = epert
  240.       +.0000448*sin(2*ze)*sin(2*Q)
  241.       +.0000137*sin(3*ze)*sin(2*Q)
  242.       +(-.0000997 + .0000108*nu2)*cos(ze)*sin(2*Q)
  243. ***************
  244. *** 554,560 ****
  245.       +(-.0000956 +.0000099*nu2)*sin(ze)*cos(2*Q)
  246.       +.0000490*sin(2*ze)*cos(2*Q)
  247.       +.0000158*sin(3*ze)*cos(2*Q)
  248. !     +.0000179*cos(2*Q)
  249.       +(.0001024 + .0000075*nu2)*cos(ze)*cos(2*Q)
  250.       -.0000437*cos(2*ze)*cos(2*Q)
  251.       -.0000132*cos(3*ze)*cos(2*Q);
  252. --- 577,584 ----
  253.       +(-.0000956 +.0000099*nu2)*sin(ze)*cos(2*Q)
  254.       +.0000490*sin(2*ze)*cos(2*Q)
  255.       +.0000158*sin(3*ze)*cos(2*Q)
  256. !     +.0000179*cos(2*Q);
  257. ! epert = epert
  258.       +(.0001024 + .0000075*nu2)*cos(ze)*cos(2*Q)
  259.       -.0000437*cos(2*ze)*cos(2*Q)
  260.       -.0000132*cos(3*ze)*cos(2*Q);
  261. ***************
  262. *** 562,568 ****
  263.   w1pert = (0.007192 - 0.003147*nu2)*sin(V)
  264.       +(-0.020428 - 0.000675*nu2 + 0.000197*nu2*nu2)*cos(V)
  265.       +(0.007269 + 0.000672*nu2)*sin(ze)*sin(Q)
  266. !     -0.004344*sin(Q)
  267.       +0.034036*cos(ze)*sin(Q)
  268.       +0.005614*cos(2*ze)*sin(Q)
  269.       +0.002964*cos(3*ze)*sin(Q)
  270. --- 586,593 ----
  271.   w1pert = (0.007192 - 0.003147*nu2)*sin(V)
  272.       +(-0.020428 - 0.000675*nu2 + 0.000197*nu2*nu2)*cos(V)
  273.       +(0.007269 + 0.000672*nu2)*sin(ze)*sin(Q)
  274. !     -0.004344*sin(Q);
  275. ! w1pert = w1pert
  276.       +0.034036*cos(ze)*sin(Q)
  277.       +0.005614*cos(2*ze)*sin(Q)
  278.       +0.002964*cos(3*ze)*sin(Q)
  279. ***************
  280. *** 573,579 ****
  281.       -0.006603*cos(ze)*cos(Q)
  282.       -0.005356*sin(ze)*sin(2*Q)
  283.       +0.002722*sin(2*ze)*sin(2*Q)
  284. !     +0.004483*cos(ze)*sin(2*Q)
  285.       -0.002642*cos(2*ze)*sin(2*Q)
  286.       +0.004403*sin(ze)*cos(2*Q)
  287.       -0.002536*sin(2*ze)*cos(2*Q)
  288. --- 598,605 ----
  289.       -0.006603*cos(ze)*cos(Q)
  290.       -0.005356*sin(ze)*sin(2*Q)
  291.       +0.002722*sin(2*ze)*sin(2*Q)
  292. !     +0.004483*cos(ze)*sin(2*Q);
  293. ! w1pert = w1pert
  294.       -0.002642*cos(2*ze)*sin(2*①)
  295.       +0.004403*sin(ze)*cos(2*Q)
  296.       -0.002536*sin(2*ze)*cos(2*Q)
  297. ***************
  298. *** 591,597 ****
  299.       +.000693*cos(2*ze)
  300.       +.000312*cos(3*ze)
  301.       +.000147*cos(4*ze)
  302. !     +.000299*sin(ze)*sin(Q)
  303.       +.000      +.000      +.000      +.000      +.000-
  304.   #define MAGURA 590
  305.   #define MAGNEP 800
  306.   
  307. + #ifdef ST_MWC
  308. + #define daylight tt.tm_isdst
  309. + extern char *getenv(); /* used to see if we're running from the desktop */
  310. + struct tm_t *tt;
  311. + #endif /* ST_MWC */
  312. + 
  313.   double pie, rad;
  314.   double htod(), atof(), kepler(), truean();
  315.   double longi(), lati(), poly(), aint(), range();
  316. ***************
  317. *** 119,125 ****
  318. --- 130,138 ----
  319.       gettimeofday(&tv, &tz);
  320.       return(GMT1970 + tv.tv_sec/SECSPERDAY);
  321.   #else
  322. + #ifndef ST_MWC
  323.       long time();
  324. + #endif /* ST_MWC */
  325.       return(GMT1970 + (double)time((long *)0)/SECSPERDAY);
  326.   #endif
  327.       }
  328. ***************
  329. *** 201,207 ****
  330.       aa1 = (int) (year/100);
  331.       bb1 = 2 - aa1 + (int)(aa1/4);
  332.       jd = aint(365.25*year) + aint(30.6001*(month + 1));
  333. !     jd = jd  + day + 1720994.5;
  334.       if((year + month/100) > 1582.10) jd = jd + bb1;
  335.       return(jd);
  336.       }
  337. --- 214,220 ----
  338.       aa1 = (int) (year/100);
  339.       bb1 = 2 - aa1 + (int)(aa1/4);
  340.       jd = aint(365.25*year) + aint(30.6001*(month + 1));
  341. !     jd += day + 1720994.5;
  342.       if((year + month/100) > 1582.10) jd = jd + bb1;
  343.       return(jd);
  344.       }
  345. ***************
  346. *** 218,224 ****
  347. --- 231,243 ----
  348.       double N,D,epli,thapp,omeg;
  349.       double nu2,P,Q,S,V,W,ze,l1pert,epert,w1pert,apert;
  350.       double psi,H,G,eta,th;
  351. + #ifdef ST_MWC
  352. +     time_t now;
  353.   
  354. +     now = time((long) 0);
  355. +     tt = localtime(&now);
  356. + #endif /* ST_MWC */
  357. + 
  358.   #define WRITEMODE "w"
  359.   #define OPENFAIL 0
  360.       progname = argv[0];
  361. ***************
  362. *** 321,327 ****
  363.           RA = atan2(N,D)/rad;
  364.           DEC = asin(sin(epli*rad)*sin(thapp*rad))/rad;
  365.           speak(RA,DEC,Sr, MAGSOL, "PS", "Sol");
  366. ! /* tansformation of coordinates on Mercury and output */
  367.           trans(r,b,ll,Stheta,Sr,epli, MAGMER, "PM", "Mercury");
  368.           
  369.   /* Now start on Venus */
  370. --- 340,346 ----
  371.           RA = atan2(N,D)/rad;
  372.           DEC = asin(sin(epli*rad)*sin(thapp*rad))/rad;
  373.           speak(RA,DEC,Sr, MAGSOL, "PS", "Sol");
  374. ! /* transformation of coordinates on Mercury and output */
  375.           trans(r,b,ll,Stheta,Sr,epli, MAGMER, "PM", "Mercury");
  376.           
  377.   /* Now start on Venus */
  378. ***************
  379. *** 373,380 ****
  380.       radpert = 0.000022501 * cos((2*M-2*M2-58.208)*rad)
  381.            +0.000019045 * cos((3*M-3*M2+92.577)*rad)
  382.            +0.000006887 * cos((M5-M2-118.090)*rad)
  383. !          +0.000005172 * cos((M-M2-29.110)*rad)
  384. !          +0.000003620 * cos((5*M-4*M2-104.208)*rad)
  385.            +0.000003283 * cos((4*M-4*M2+63.513)*rad)
  386.            +0.000003074 * cos((2*M5-2*M2-55.167)*rad);
  387.       r = r + radpert;
  388. --- 392,399 ----
  389.       radpert = 0.000022501 * cos((2*M-2*M2-58.208)*rad)
  390.            +0.000019045 * cos((3*M-3*M2+92.577)*rad)
  391.            +0.000006887 * cos((M5-M2-118.090)*rad)
  392. !          +0.000005172 * cos((M-M2-29.110)*rad);
  393. !     radpert+= 0.000003620 * cos((5*M-4*M2-104.208)*rad)
  394.            +0.000003283 * cos((4*M-4*M2+63.513)*rad)
  395.            +0.000003074 * cos((2*M5-2*M2-55.167)*rad);
  396.       r = r + radpert;
  397. ***************
  398. *** 425,432 ****
  399.       lonpert = 0.00705*cos((M5-M4-48.958)*rad)
  400.            +0.00607*cos((2*M5-M4-188.350)*rad)
  401.            +0.00445*cos((2*M5-2*M4-191.897)*rad)
  402. !          +0.00388*cos((M-2*M4+20.495)*rad)
  403. !          +0.00238*cos((M-M4+35.097)*rad)
  404.            +0.00204*cos((2*M-3*M4+158.638)*rad)
  405.            +0.00177*cos((3*M4-M2-57.602)*rad)
  406.            +0.00136*cos((2*M-4*M4+154.093)*rad)
  407. --- 444,451 ----
  408.       lonpert = 0.00705*cos((M5-M4-48.958)*rad)
  409.            +0.00607*cos((2*M5-M4-188.350)*rad)
  410.            +0.00445*cos((2*M5-2*M4-191.897)*rad)
  411. !          +0.00388*cos((M-2*M4+20.495)*rad);
  412. !     lonpert+= 0.00238*cos((M-M4+35.097)*rad)
  413.            +0.00204*cos((2*M-3*M4+158.638)*rad)
  414.            +0.00177*cos((3*M4-M2-57.602)*rad)
  415.            +0.00136*cos((2*M-4*M4+154.093)*rad)
  416. ***************
  417. *** 445,452sin(4*ze)*sin(2*Q)
  418. ***************
  419. *** 707,713 ****
  420.       +.0000467*cos(3*psi)*sin(2*Q)
  421.       -.0000490*cos(2*Q)
  422.       -(.0002842 + .0000279*nu2)*sin(ze)*cos(2*Q)
  423. !     +(.0000128 + .0000226*nu2)*sin(2*ze)*cos(2*Q)
  424.       +.0000224*sin(3*ze)*cos(2*Q)
  425.       +(-.0001594 + .0000282*nu2)*cos(ze)*cos(2*Q)
  426.       +(.0002162 - .0000207*nu2)*cos(2*ze)*cos(2*Q)
  427. --- 740,747 ----
  428.       +.0000467*cos(3*psi)*sin(2*Q)
  429.       -.0000490*cos(2*Q)
  430.       -(.0002842 + .0000279*nu2)*sin(ze)*cos(2*Q)
  431. !     +(.0000128 + .0000226*nu2)*sin(2*ze)*cos(2*Q);
  432. ! epert = epert
  433.       +.0000224*sin(3*ze)*cos(2*Q)
  434.       +(-.0001594 + .0000282*nu2)*cos(ze)*cos(2*Q)
  435.       +(.0002162 - .0000207*nu2)*cos(2*ze)*cos(2*Q)
  436. ***************
  437. *** 716,722 ****
  438.       +.0000343*cos(4*ze)*cos(2*Q)
  439.       +.0000469*sin(3*psi)*cos(2*Q)
  440.       -.0000242*cos(3*psi)*cos(2*Q)
  441. !     -.0000205*sin(ze)*sin(3*Q)
  442.       +.0000262*sin(3*ze)*sin(3*Q)
  443.       +.0000208*cos(ze)*cos(3*Q)
  444.       -.0000271*cos(3*ze)*cos(3*Q)
  445. --- 750,757 ----
  446.       +.0000343*cos(4*ze)*cos(2*Q)
  447.       +.0000469*sin(3*psi)*cos(2*Q)
  448.       -.0000242*cos(3*psi)*cos(2*Q)
  449. !     -.0000205*sin(ze)*sin(3*Q);
  450. ! epert = epert
  451.       +.0000262*sin(3*ze)*sin(3*Q)
  452.       +.0000208*cos(ze)*cos(3*Q)
  453.       -.0000271*cos(3*ze)*cos(3*Q)
  454. ***************
  455. *** 728,734 ****
  456.       -0.007075*sin(ze)
  457.       -0.075825*sin(ze)*sin(Q)
  458.       -0.024839*sin(2*ze)*sin(Q)
  459. !     -0.008631*sin(3*ze)*sin(Q)
  460.       -0.072586*cos(Q)
  461.       -0.150383*cos(ze)*cos(Q)
  462.       +0.026897*cos(2*ze)*cos(Q)
  463. --- 763,770 ----
  464.       -0.007075*sin(ze)
  465.       -0.075825*sin(ze)*sin(Q)
  466.       -0.024839*sin(2*ze)*sin(Q)
  467. !     -0.008631*sin(3*ze)*sin(Q);
  468. ! w1pert = w1pert
  469.       -0.072586*cos(Q)
  470.       -0.150383*cos(ze)*cos(Q)
  471.       +0.026897*cos(2*ze)*cos(Q)
  472. ***************
  473. *** 737,743 ****
  474.       -(0.013597 +0.001719*nu2)*sin(ze)*sin(2*Q)
  475.       +(-0.007742 + 0.001517*nu2)*cos(ze)*sin(2*Q)
  476.       +(0.013586 - 0.001375*nu2)*cos(2*ze)*sin(2*Q)
  477. !     +(-0.013667 + 0.001239*nu2)*sin(ze)*cos(2*Q)
  478.       +0.011981*sin(2*ze)*cos(2*Q)
  479.       +(0.014861 + 0.001136*nu2)*cos(ze)*cos(2*Q)
  480.       -(0.013064 + 0.001628*nu2)*cos(2*ze)*cos(2*Q);
  481. --- 773,780 ----
  482.       -(0.013597 +0.001719*nu2)*sin(ze)*sin(2*Q)
  483.       +(-0.007742 + 0.001517*nu2)*cos(ze)*sin(2*Q)
  484.       +(0.013586 - 0.001375*nu2)*cos(2*ze)*sin(2*Q)
  485. !     +(-0.013667 + 0.001239*nu2)*sin(ze)*cos(2*Q);
  486. ! w1pert = w1pert
  487.       +0.011981*sin(2*ze)*cos(2*Q)
  488.       +(0.014861 + 0.001136*nu2)*cos(ze)*cos(2*Q)
  489.       -(0.013064 + 0.001628*nu2)*cos(2*ze)*cos(2*Q);
  490. ***************
  491. *** 751,757 ****
  492.   apert = .000572*sin(V) -.001590*sin(2*ze)*cos(Q)
  493.       +.002933*cos(V) -.000647*sin(3*ze)*cos(Q)
  494.       +.033629*cos(ze) -.000344*sin(4*ze)*cos(Q)
  495. !     -.003081*cos(2*ze) +.002885*cos(ze)*cos(Q)
  496.       -.001423*cos(3*ze) +(.002172 + .000102*nu2)*cos(2*ze)*cos(Q)
  497.       -.000671*cos(4*ze) +.000296*cos(3*ze)*cos(Q)
  498.       -.000320*cos(5*ze) -.000267*sin(2*ze)*sin(2*Q);
  499. --- 788,795 ----
  500.   apert = .000572*sin(V) -.001590*sin(2*ze)*cos(Q)
  501.       +.002933*cos(V) -.000647*sin(3*ze)*cos(Q)
  502.       +.033629*cos(ze) -.000344*sin(4*ze)*cos(Q)
  503. !     -.003081*cos(2*ze) +.002885*cos(ze)*cos(Q);
  504. ! apert = apert
  505.       -.001423*cos(3*ze) +(.002172 + .000102*nu2)*cos(2*ze)*cos(Q)
  506.       -.000671*cos(4*ze) +.000296*cos(3*ze)*cos(Q)
  507.       -.000320*cos(5*ze) -.000267*sin(2*ze)*sin(2*Q);
  508. ***************
  509. *** 763,769 ****
  510.       -.000393*sin(3*ze)*sin(Q)
  511.       -.000228*sin(4*ze)*sin(Q)
  512.       +.002138*cos(ze)*sin(Q)
  513. !     -.000999*cos(2*ze)*sin(Q)
  514.       -.000642*cos(3*ze)*sin(Q)
  515.       -.000325*cos(4*ze)*sin(Q)
  516.       -.000890*cos(Q)
  517. --- 801,808 ----
  518.       -.000393*sin(3*ze)*sin(Q)
  519.       -.000228*sin(4*ze)*sin(Q)
  520.       +.002138*cos(ze)*sin(Q)
  521. !     -.000999*cos(2*ze)*sin(Q);
  522. ! apert = apert
  523.       -.000642*cos(3*ze)*sin(Q)
  524.       -.000325*cos(4*ze)*sin(Q)
  525.       -.000890*cos(Q)
  526. ***************
  527. *** 772,778 ****
  528.       -.000856*sin(ze)*cos(2*Q)
  529.       +.000441*sin(2*ze)*cos(2*Q)
  530.       +.000296*cos(2*ze)*cos(2*Q)
  531. !     +.000211*cos(3*ze)*cos(2*Q)
  532.       -.000427*sin(ze)*sin(3*Q)
  533.       +.000398*sin(3*ze)*sin(3*Q)
  534.       +.000344*cos(ze)*cos(3*Q)
  535. --- 811,818 ----
  536.       -.000856*sin(ze)*cos(2*Q)
  537.       +.000441*sin(2*ze)*cos(2*Q)
  538.       +.000296*cos(2*ze)*cos(2*Q)
  539. !     +.000211*cos(3*ze)*cos(2*Q);
  540. ! apert = apert
  541.       -.000427*sin(ze)*sin(3*Q)
  542.       +.000398*sin(3*ze)*sin(3*Q)
  543.       +.000344*cos(ze)*cos(3*Q)
  544. ***************
  545. *** 1002,1009 ****
  546.       +.001051*cos(2*th);
  547.   trans(r,b,ll,Stheta,Sr,epli, MAGNEP, "PN", "Neptune");
  548.   
  549.   putchar('\n');
  550. ! close(logfile);
  551.   exit(0);
  552.   } /* end of program main */
  553.   
  554. --- 1042,1061 ----
  555.       +.001051*cos(2*th);
  556.   trans(r,b,ll,Stheta,Sr,epli, MAGNEP, "PN", "Neptune");
  557.   
  558. + #ifdef ST_MWC
  559. + Cconws("\r\n");
  560. + #else /* !ST_MWC */
  561.   putchar('\n');
  562. ! #endif /* ST_MWC */
  563. ! fclose(logfile);
  564. ! #ifdef ST_MWC /* We want to hold the screen if invoked from the GEM desktop; */
  565. !            /* the desktop doesn't usually set any environment variables, */
  566. !            /* and if getenv() does find anything, it's probably empty */
  567. ! if (((getenv("PATH")) == 0) || (strlen(getenv("PATH")) == 0)) {
  568. !      Cconws("press any key to continue: "); /* Hold screen if desktop */
  569. !      i = Crawcin(); /* Read raw character input */
  570. ! }
  571. ! #endif /* ST_MWC */
  572.   exit(0);
  573.   } /* end of program main */
  574.   
  575. ***************
  576. *** 1023,1031 ****
  577. --- 1075,1089 ----
  578.   double aint(z)
  579.       double z;
  580.   {
  581. + #ifdef ST_MWC /* In Mark Williams C ints are 16 bits long, so use longs */
  582. +     long trunk;
  583. +     
  584. +     trunk = (long)z;
  585. + #else /* !ST_MWC */
  586.       int trunk;
  587.       
  588.       trunk = (int)z;
  589. + #endif /* ST_MWC */
  590.       z = (double) trunk;
  591.       return(z);
  592.   }
  593.