home *** CD-ROM | disk | FTP | other *** search
/ Usenet 1994 January / usenetsourcesnewsgroupsinfomagicjanuary1994.iso / sources / sun / volume1 / calentool / patch4b / patches04b
Encoding:
Text File  |  1989-10-26  |  51.5 KB  |  2,035 lines
  1. *** /tmp/,RCSt1a05555    Wed Jul 19 20:46:02 1989
  2. --- moon.c    Wed Jul 19 20:29:00 1989
  3. ***************
  4. *** 1,5
  5.   /*
  6. !  * $Header: moon.c,v 2.1 89/05/09 14:19:22 billr Exp $
  7.    */
  8.   /*
  9.    * moon.c
  10.  
  11. --- 1,5 -----
  12.   /*
  13. !  * $Header: moon.c,v 2.2 89/07/19 20:28:12 billr Exp $
  14.    */
  15.   /*
  16.    * moon.c
  17. ***************
  18. *** 9,15
  19.    * Bill Randle, Tektronix, Inc. to interface to the calentool
  20.    * program.
  21.    */
  22. ! #include "ct.h"        /* for the NO_SUN_MOON #define */
  23.   #ifndef NO_SUN_MOON
  24.   
  25.   /*
  26.  
  27. --- 9,15 -----
  28.    * Bill Randle, Tektronix, Inc. to interface to the calentool
  29.    * program.
  30.    */
  31. ! #include "ct.h"     /* for the NO_SUN_MOON #define */
  32.   #ifndef NO_SUN_MOON
  33.   
  34.   /*
  35. ***************
  36. *** 16,22
  37.   
  38.       A Moon for the Sun
  39.   
  40. !     Release 2.0
  41.   
  42.       Designed and implemented by John Walker in December 1987,
  43.       revised and updated in February of 1988.
  44.  
  45. --- 16,22 -----
  46.   
  47.       A Moon for the Sun
  48.   
  49. !     Release 2.3
  50.   
  51.       Designed and implemented by John Walker in December 1987,
  52.       revised and updated in February of 1988.
  53. ***************
  54. *** 20,25
  55.   
  56.       Designed and implemented by John Walker in December 1987,
  57.       revised and updated in February of 1988.
  58.   
  59.       cc -O moontool.c -o moontool -lm -lsuntool -lsunwindow -lpixrect
  60.   
  61.  
  62. --- 20,26 -----
  63.   
  64.       Designed and implemented by John Walker in December 1987,
  65.       revised and updated in February of 1988.
  66. +     revised and updated again in June of 1988 by Ron Hitchens
  67.   
  68.       Make with:
  69.   
  70. ***************
  71. *** 21,26
  72.       Designed and implemented by John Walker in December 1987,
  73.       revised and updated in February of 1988.
  74.   
  75.       cc -O moontool.c -o moontool -lm -lsuntool -lsunwindow -lpixrect
  76.   
  77.       Adding  appropriate  floating  point  options  to your hardware.  This
  78.  
  79. --- 22,29 -----
  80.       revised and updated in February of 1988.
  81.       revised and updated again in June of 1988 by Ron Hitchens
  82.   
  83. +     Make with:
  84. + 
  85.       cc -O moontool.c -o moontool -lm -lsuntool -lsunwindow -lpixrect
  86.   
  87.       Adding  appropriate  floating  point  options  to your hardware.  This
  88. ***************
  89. *** 75,80
  90.       program.  Please leave the original attribution information intact    so
  91.       that credit and blame may be properly apportioned.
  92.   
  93.   */
  94.   
  95.   /*  Astronomical constants  */
  96.  
  97. --- 78,147 -----
  98.       program.  Please leave the original attribution information intact    so
  99.       that credit and blame may be properly apportioned.
  100.   
  101. +    Revision history:
  102. + 
  103. +         1.0  11/5/87   First version.  Only displayed icon, no
  104. +                        open window information display.
  105. + 
  106. +         2.0  3/27/88   First posting to comp.sources.unix.
  107. + 
  108. +         2.1  6/16/88   Bug fix.  Table of phases didn't update
  109. +                        at the moment of the new moon.  Call on
  110. +                        phasehunt didn't convert civil Julian date
  111. +                        to astronomical Julian date.  Reported by
  112. +                        Dag Bruck (dag@control.lth.se).
  113. + 
  114. +         2.2  2/27/89   Michael McClary  (michael@xanadu.COM)
  115. +                        Added moon map, derived from the sun "fullmoon"
  116. +                        image file.  (It has not been checked for
  117. +                        rotation from the correct orientation.  Also, a
  118. +                        fixed icon doesn't model the librations of the
  119. +                        moon as viewed from the earth.) Also: tweaked
  120. +                        corners of icon.
  121. + 
  122. +         2.3  6/7/89    Bug fix.  Table of phases skipped the phases
  123. +                        for July 1989.  This occurred due to sloppy
  124. +                        maintenance of the synodic month index in the
  125. +                        interchange of information between phasehunt()
  126. +                        and meanphase().  I simplified and corrected
  127. +                        the handling of the month index as phasehunt()
  128. +                        steps along and removed unneeded code from
  129. +                        meanphase().  Reported by Bill Randle of
  130. +                        Tektronix, Inc. (billr@saab.CNA.TEK.COM).
  131. + 
  132. + Additional History:
  133. + 
  134. +     June 1988    Modified by Ron Hitchens to produce version 2.1
  135. +             modified icon generation to show surface texture
  136. +              on visible moon face.  Eliminated "illegal" direct
  137. +              modification of icon image memory.
  138. +             added a menu to allow switching in and out of
  139. +              test mode, for entertainment value mostly.
  140. +             reworked timer behaviour so that process doesn't
  141. +              wake up unnecessarily.
  142. +             trap sigwinch signals to notice more easily when the
  143. +              tool opens and closes.
  144. +             modified layout of information in open window display
  145. +              to reduce the amount of pixels modified in each
  146. +              update.  Batched pixwin updates so that only one
  147. +              screen rasterop per cycle is done.
  148. +             changed open window to display white-on-black for a
  149. +              more aesthetic look, and to suggest the effect of
  150. +              looking at the moon in the nighttime sky.
  151. +             setup default tool and canvas colors to be the same
  152. +              as B&W monochrome, for those us lucky enough to have
  153. +              color monitors and who have the default monochrome
  154. +              colors set to something other than B&W (I like white
  155. +              on dark blue myself)
  156. +             various code reformatting and pretty-printing to suit
  157. +              my own coding taste (I got a bit carried away).
  158. +             code tweaking to make lint happy.
  159. +             returned my hacked version to John.
  160. + 
  161. +             Ron Hitchens
  162. +                 ronbo@vixen.uucp
  163. +                 ...!uunet!cs.utexas.edu!vixen!ronbo
  164. +                 hitchens@cs.utexas.edu
  165.   */
  166.   
  167.   
  168. ***************
  169. *** 77,82
  170.   
  171.   */
  172.   
  173.   /*  Astronomical constants  */
  174.   
  175.   #define epoch        2444238.5       /* 1980 January 0.0 */
  176.  
  177. --- 144,167 -----
  178.                   hitchens@cs.utexas.edu
  179.   */
  180.   
  181. + 
  182. + #include <stdio.h>
  183. + #include <math.h>
  184. + #include <sys/time.h>
  185. + 
  186. + #include <suntool/sunview.h>
  187. + #include <suntool/canvas.h>
  188. + 
  189. + #define TINYFONT    "/usr/lib/fonts/fixedwidthfonts/screen.r.7"
  190. + 
  191. + /*
  192. +  * define standard B/W monochrome colors as defaults in case we're running
  193. +  * on a color system with the monochrome colors set differently
  194. +  */
  195. + #define FG_DEFAULT    { 0, 0, 0 }        /* black */
  196. + #define BG_DEFAULT    { 255, 255, 255 }    /* white */
  197. + 
  198. + 
  199.   /*  Astronomical constants  */
  200.   
  201.   #define epoch        2444238.5       /* 1980 January 0.0 */
  202. ***************
  203. *** 80,86
  204.   /*  Astronomical constants  */
  205.   
  206.   #define epoch        2444238.5       /* 1980 January 0.0 */
  207. ! #define    J1970        2440587.5       /* VAX clock Epoch 1970 Jan 1 (0h UT) */
  208.   
  209.   /*  Constants defining the Sun's apparent orbit  */
  210.   
  211.  
  212. --- 165,171 -----
  213.   /*  Astronomical constants  */
  214.   
  215.   #define epoch        2444238.5       /* 1980 January 0.0 */
  216. ! #define J1970       2440587.5      /* VAX clock Epoch 1970 Jan 1 (0h UT) */
  217.   
  218.   /*  Constants defining the Sun's apparent orbit  */
  219.   
  220. ***************
  221. *** 114,122
  222.   
  223.   #define earthrad    6378.16       /* Radius of Earth in kilometres */
  224.   
  225. - #include <stdio.h>
  226. - #include <math.h>
  227. - #include <sys/time.h>
  228.   
  229.   #include <suntool/sunview.h>
  230.   #include <suntool/canvas.h>
  231.  
  232. --- 199,204 -----
  233.   
  234.   #define earthrad    6378.16       /* Radius of Earth in kilometres */
  235.   
  236.   
  237.   #define PI 3.14159265358979323846  /* Assume not near black hole nor in
  238.                         Tennessee */
  239. ***************
  240. *** 118,174
  241.   #include <math.h>
  242.   #include <sys/time.h>
  243.   
  244. - #include <suntool/sunview.h>
  245. - #include <suntool/canvas.h>
  246. - 
  247. - /* Icon definition.  This is just a black field with rounded corners
  248. -    that blend into the root desktop pattern. The image of the moon and
  249. -    the text are added by the program later. */
  250. - 
  251. - static short moon_img[64][4] = {
  252. - /* Format_version=1, Width=64, Height=64, Depth=1, Valid_bits_per_item=16 */
  253. -     0x8FFF,0xFFFF,0xFFFF,0xFFE8,0x9FFF,0xFFFF,0xFFFF,0xFFF8,
  254. -     0x3FFF,0xFFFF,0xFFFF,0xFFFE,0x7FFF,0xFFFF,0xFFFF,0xFFFE,
  255. -     0xFFFF,0xFFFF,0xFFFF,0xFFFE,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  256. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  257. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  258. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  259. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  260. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  261. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  262. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  263. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  264. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  265. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  266. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  267. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  268. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  269. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  270. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  271. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  272. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  273. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  274. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  275. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  276. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  277. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  278. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  279. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  280. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  281. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,0xFFFF,
  282. -     0xFFFF,0xFFFF,0xFFFF,0xFFFF,0x7FFF,0xFFFF,0xFFFF,0xFFFE,
  283. -     0xFFFF,0xFFFF,0xFFFF,0xFFFE,0xBFFF,0xFFFF,0xFFFF,0xFFFC,
  284. -     0x3FFF,0xFFFF,0xFFFF,0xFFFA,0x27FF,0xFFFF,0xFFFF,0xFFE2
  285. - };
  286. - DEFINE_ICON_FROM_IMAGE(moon_icon, moon_img);
  287. - mpr_static(icon_mpr, 64, 64, 1, moon_img);
  288. - 
  289. - extern Canvas mcanvas;
  290. - extern Pixfont *font;
  291. - extern char *monthnames[];
  292. - static Pixwin *cpw;
  293. - static int charhgt, charwid;
  294. - 
  295.   #define PI 3.14159265358979323846  /* Assume not near black hole nor in
  296.                         Tennessee */
  297.   
  298.  
  299. --- 200,205 -----
  300.   #define earthrad    6378.16       /* Radius of Earth in kilometres */
  301.   
  302.   
  303.   #define PI 3.14159265358979323846  /* Assume not near black hole nor in
  304.                         Tennessee */
  305.   
  306. ***************
  307. *** 182,188
  308.   #define dsin(x) (sin(torad((x))))              /* Sin from deg */
  309.   #define dcos(x) (cos(torad((x))))              /* Cos from deg */
  310.   
  311. - /*  Forward functions  */
  312.   
  313.   double jtime(), phase();
  314.   void phasehunt();
  315.  
  316. --- 213,218 -----
  317.   #define dsin(x) (sin(torad((x))))              /* Sin from deg */
  318.   #define dcos(x) (cos(torad((x))))              /* Cos from deg */
  319.   
  320.   
  321.   /*
  322.    * Moon image.  This is a standard icon-sized picture of the moon's face.
  323. ***************
  324. *** 184,192
  325.   
  326.   /*  Forward functions  */
  327.   
  328. ! double jtime(), phase();
  329. ! void phasehunt();
  330. ! void drawmoon(), jyear(), jhms();
  331.   
  332.   
  333.   /*  DRAWMOON  --  Construct icon for moon, given phase of moon.  */
  334.  
  335. --- 214,224 -----
  336.   #define dcos(x) (cos(torad((x))))              /* Cos from deg */
  337.   
  338.   
  339. ! /*
  340. !  * Moon image.  This is a standard icon-sized picture of the moon's face.
  341. !  * The visible part, as calculated by the current time, is extracted from
  342. !  * this image to create the displayed image.
  343. !  */
  344.   
  345.   static short moon_img [] = {
  346.   #include "moon.icon"
  347. ***************
  348. *** 188,193
  349.   void phasehunt();
  350.   void drawmoon(), jyear(), jhms();
  351.   
  352.   
  353.   /*  DRAWMOON  --  Construct icon for moon, given phase of moon.  */
  354.   
  355.  
  356. --- 220,229 -----
  357.    * this image to create the displayed image.
  358.    */
  359.   
  360. + static short moon_img [] = {
  361. + #include "moon.icon"
  362. + };
  363. + mpr_static(moon_mpr, 64, 64, 1, moon_img);
  364.   
  365.   extern Pixfont    *font;            /* pointer to regular icon font */
  366.   extern Pixfont    *sfont;            /* pointer to tiny icon font */
  367. ***************
  368. *** 189,195
  369.   void drawmoon(), jyear(), jhms();
  370.   
  371.   
  372. ! /*  DRAWMOON  --  Construct icon for moon, given phase of moon.  */
  373.   
  374.   static void drawmoon(ph)
  375.   double ph;
  376.  
  377. --- 225,239 -----
  378.   };
  379.   mpr_static(moon_mpr, 64, 64, 1, moon_img);
  380.   
  381. ! extern Pixfont    *font;            /* pointer to regular icon font */
  382. ! extern Pixfont    *sfont;            /* pointer to tiny icon font */
  383. ! extern Canvas    mcanvas;            /* handle for the canvas */
  384. ! extern char *monthnames[];
  385. ! static Pixrect    *icon_mpr;        /* actual displayed pixrect */
  386. ! static Pixwin    *cpw;            /* pointer to the canvas' pixwin */
  387. ! static int    charhgt, charwid;    /* default std font height/width */
  388. ! static int info_col;
  389. ! static struct timezone tzp;
  390.   
  391.   static char *labels [] = {
  392.       "Moon phase:",
  393. ***************
  394. *** 191,198
  395.   
  396.   /*  DRAWMOON  --  Construct icon for moon, given phase of moon.  */
  397.   
  398. ! static void drawmoon(ph)
  399. ! double ph;
  400.   {
  401.       int i, j, lx, rx;
  402.       int lb[4];
  403.  
  404. --- 235,265 -----
  405.   static int info_col;
  406.   static struct timezone tzp;
  407.   
  408. ! static char *labels [] = {
  409. !     "Moon phase:",
  410. !     "Age of moon:",
  411. !     "Moon's distance:",
  412. !     "Moon subtends:",
  413. !     "Last new moon:",
  414. !     "First quarter:",
  415. !     "Full moon:",
  416. !     "Last quarter:",
  417. !     "Next new moon:"
  418. ! };
  419. ! 
  420. ! /*  Forward functions  */
  421. ! static double        jtime(), phase();
  422. ! static void        phasehunt(), set_mode();
  423. ! static void        drawmoon(), jyear(), jhms();
  424. ! 
  425. ! 
  426. ! /*
  427. !  * DRAWMOON  --  Construct icon for moon, given phase of moon.
  428. !  */
  429. ! 
  430. ! static void drawmoon (ph, src_pr, dst_pr)
  431. ! double        ph;
  432. ! Pixrect        *src_pr, *dst_pr;
  433.   {
  434.       register int    i, lx, rx;
  435.       register double    cp, xscale;
  436. ***************
  437. *** 194,202
  438.   static void drawmoon(ph)
  439.   double ph;
  440.   {
  441. !     int i, j, lx, rx;
  442. !     int lb[4];
  443. !     double cp, xscale;
  444.   
  445.       xscale = cos(2 * PI * ph);
  446.       for (i = 0; i < 24; i++) {
  447.  
  448. --- 261,272 -----
  449.   double        ph;
  450.   Pixrect        *src_pr, *dst_pr;
  451.   {
  452. !     register int    i, lx, rx;
  453. !     register double    cp, xscale;
  454. ! #define RADIUS        27.0
  455. ! #define IRADIUS        27
  456. ! #define OFFSET        28
  457. ! #define CENTER        32
  458.   
  459.       /* Clear the destination pixrect to all one-bits (black) */
  460.       pr_rop (dst_pr, 0, 0, 64, 64, PIX_SET, (Pixrect *)0, 0, 0);
  461. ***************
  462. *** 198,219
  463.       int lb[4];
  464.       double cp, xscale;
  465.   
  466. !     xscale = cos(2 * PI * ph);
  467. !     for (i = 0; i < 24; i++) {
  468. !        lb[0] = lb[1] = lb[2] = lb[3] = 0xFFFF;
  469. !        cp = 24.0 * cos(asin(i / 24.0));
  470. !        if (ph < 0.5) {
  471. !           rx = 32 + cp;
  472. !           lx = 32 + xscale * cp;
  473. !        } else {
  474. !           lx = 33 - cp;
  475. !           rx = 33 - xscale * cp;
  476. !        }
  477. !        for (j = lx; j <= rx; j++) {
  478. !           lb[j >> 4] &= (0x8000 >> (j & 0xF)) ^ 0xFFFF;
  479. !        }
  480. !        for (j = 0; j < 4; j++)
  481. !           moon_img[28 + i][j] = moon_img[28 - i][j] = lb[j];
  482.       }
  483.   }
  484.   
  485.  
  486. --- 268,298 -----
  487.   #define OFFSET        28
  488.   #define CENTER        32
  489.   
  490. !     /* Clear the destination pixrect to all one-bits (black) */
  491. !     pr_rop (dst_pr, 0, 0, 64, 64, PIX_SET, (Pixrect *)0, 0, 0);
  492. ! 
  493. !     xscale = cos (2 * PI * ph);
  494. !     for (i = 0; i < IRADIUS; i++) {
  495. !         cp = RADIUS * cos (asin (i / RADIUS));
  496. !         if (ph < 0.5) {
  497. !             rx = CENTER + cp;
  498. !             lx = CENTER + xscale * cp;
  499. !         } else {
  500. !             lx = CENTER - cp;
  501. !             rx = CENTER - xscale * cp;
  502. !         }
  503. ! 
  504. !         /*
  505. !          * We now know the left and right endpoints of the scan line
  506. !          * for this y coordinate.  We raster-op the corresponding
  507. !          * scanlines from the source pixrect to the destination
  508. !          * pixrect, offsetting to properly place it in the pixrect and
  509. !          * reflecting vertically.
  510. !          */
  511. !         pr_rop (dst_pr, lx, OFFSET + i, (rx - lx) + 1, 1, PIX_SRC,
  512. !             src_pr, lx, OFFSET + i);
  513. !         pr_rop (dst_pr, lx, OFFSET - i, (rx - lx) + 1, 1, PIX_SRC,
  514. !             src_pr, lx, OFFSET - i);
  515.       }
  516.   }
  517.   
  518. ***************
  519. *** 217,222
  520.       }
  521.   }
  522.   
  523.   /*  MOON_DATA  -- print useful info about the moon */
  524.   
  525.   moon_data(seconds)
  526.  
  527. --- 296,307 -----
  528.       }
  529.   }
  530.   
  531. + #define prt(y) pw_text(cpw, info_col, charhgt*(y), PIX_SRC, font, tbuf)
  532. + #define prtxy(x,y) pw_text(cpw, charwid*(y+1), charhgt*(x), PIX_SRC,\
  533. +     font,tbuf)
  534. + #define EPL(x) (x), (x) == 1 ? "" : "s"
  535. + #define APOS(x) (x + 11)
  536. + 
  537.   /*  MOON_DATA  -- print useful info about the moon */
  538.   
  539.   moon_data(seconds)
  540. ***************
  541. *** 222,238
  542.   moon_data(seconds)
  543.   long seconds;
  544.   {
  545. !     int lunation, wclosed;
  546. !     long t;
  547. !     double jd, p, aom, cphase, cdist, cangdia, csund, csuang, lptime;
  548. !     double phasar[5];
  549. !     struct pr_prpos tloc;
  550. !     char amsg[12], tbuf[80];
  551. !     static double faketime = 0.0;
  552. !     static short moonilast[64][4] = {0};
  553. !     int yy, mm, dd, hh, mmm, ss;
  554. !     Pixfont *pfont;
  555. !     struct tm *gm, *ltm, *localtime();
  556.       struct timeval tvp;
  557.       struct timezone tzp;
  558.       long clock;
  559.  
  560. --- 307,320 -----
  561.   moon_data(seconds)
  562.   long seconds;
  563.   {
  564. !     int        lunation, wclosed;
  565. !     int        i, yy, mm, dd, hh, mmm, ss;
  566. !     int        aom_d, aom_h, aom_m;
  567. !     long        t;
  568. !     double        jd, p, aom, cphase, cdist, cangdia, csund, csuang;
  569. !     double        phasar [5];
  570. !     char        tbuf[80];
  571. !     struct tm    *gm, *ltm, *localtime();
  572.       struct timeval tvp;
  573.       struct pr_prpos    tloc;
  574.       struct pr_size    txt_size;
  575. ***************
  576. *** 234,240
  577.       Pixfont *pfont;
  578.       struct tm *gm, *ltm, *localtime();
  579.       struct timeval tvp;
  580. !     struct timezone tzp;
  581.       long clock;
  582.       char *atp, *asctime();
  583.   
  584.  
  585. --- 316,323 -----
  586.       char        tbuf[80];
  587.       struct tm    *gm, *ltm, *localtime();
  588.       struct timeval tvp;
  589. !     struct pr_prpos    tloc;
  590. !     struct pr_size    txt_size;
  591.       long clock;
  592.       char *atp, *asctime();
  593.   
  594. ***************
  595. *** 238,245
  596.       long clock;
  597.       char *atp, *asctime();
  598.   
  599. ! 
  600. !     pfont = pf_open("/usr/lib/fonts/fixedwidthfonts/screen.r.7");
  601.       cpw = canvas_pixwin(mcanvas);
  602.       charwid = font->pf_defaultsize.x;
  603.       charhgt = font->pf_defaultsize.y;
  604.  
  605. --- 321,328 -----
  606.       long clock;
  607.       char *atp, *asctime();
  608.   
  609. !     icon_mpr = mem_create (64, 64, 1);
  610. !     pr_rop (icon_mpr, 0, 0, 64, 64, PIX_SRC, &moon_mpr, 0, 0);
  611.       cpw = canvas_pixwin(mcanvas);
  612.       charwid = font->pf_defaultsize.x;
  613.       charhgt = font->pf_defaultsize.y;
  614. ***************
  615. *** 243,249
  616.       cpw = canvas_pixwin(mcanvas);
  617.       charwid = font->pf_defaultsize.x;
  618.       charhgt = font->pf_defaultsize.y;
  619. !     jd = jtime((gm = gmtime(&seconds)));
  620.       gettimeofday(&tvp, &tzp); /* for timezone info */
  621.       p = phase(jd, &cphase, &aom, &cdist, &cangdia, &csund, &csuang);
  622.       drawmoon(p);
  623.  
  624. --- 326,333 -----
  625.       cpw = canvas_pixwin(mcanvas);
  626.       charwid = font->pf_defaultsize.x;
  627.       charhgt = font->pf_defaultsize.y;
  628. !     info_col = charwid * 20;
  629. !     jd = jtime ((gm = gmtime (&seconds)));
  630.       gettimeofday(&tvp, &tzp); /* for timezone info */
  631.       p = phase (jd, &cphase, &aom, &cdist, &cangdia, &csund, &csuang);
  632.       aom_d = (int) aom;
  633. ***************
  634. *** 245,258
  635.       charhgt = font->pf_defaultsize.y;
  636.       jd = jtime((gm = gmtime(&seconds)));
  637.       gettimeofday(&tvp, &tzp); /* for timezone info */
  638. !     p = phase(jd, &cphase, &aom, &cdist, &cangdia, &csund, &csuang);
  639. !     drawmoon(p);
  640. !         sprintf(amsg, " %dd %dh  ",
  641. !        (int) aom, ((int) (24 * (aom - floor(aom)))));
  642. !     tloc.pr = (Pixrect *) icon_get(&moon_icon, ICON_IMAGE);
  643. !     tloc.pos.x = 2;
  644. !     tloc.pos.y = 62;
  645. !     pf_text(tloc, PIX_NOT(PIX_SRC), pfont, amsg);
  646.   
  647.       /* Only update icon if it changed (this eliminates gratuitous
  648.          flashing of the icon on-screen). */
  649.  
  650. --- 329,338 -----
  651.       info_col = charwid * 20;
  652.       jd = jtime ((gm = gmtime (&seconds)));
  653.       gettimeofday(&tvp, &tzp); /* for timezone info */
  654. !     p = phase (jd, &cphase, &aom, &cdist, &cangdia, &csund, &csuang);
  655. !     aom_d = (int) aom;
  656. !     aom_h = (int) (24 * (aom - floor(aom)));
  657. !     aom_m = (int) (1440 * (aom - floor(aom))) % 60;
  658.   
  659.       drawmoon(p, &moon_mpr, icon_mpr);
  660.   
  661. ***************
  662. *** 254,261
  663.       tloc.pos.y = 62;
  664.       pf_text(tloc, PIX_NOT(PIX_SRC), pfont, amsg);
  665.   
  666. !     /* Only update icon if it changed (this eliminates gratuitous
  667. !        flashing of the icon on-screen). */
  668.   
  669.       if (bcmp(moonilast, moon_img, sizeof moon_img) != 0)
  670.          bcopy(moon_img, moonilast, sizeof moon_img);
  671.  
  672. --- 334,340 -----
  673.       aom_h = (int) (24 * (aom - floor(aom)));
  674.       aom_m = (int) (1440 * (aom - floor(aom))) % 60;
  675.   
  676. !     drawmoon(p, &moon_mpr, icon_mpr);
  677.   
  678.       if (aom_d == 0) {
  679.           (void)sprintf(tbuf, "%dh %dm", aom_h, aom_m);
  680. ***************
  681. *** 257,264
  682.       /* Only update icon if it changed (this eliminates gratuitous
  683.          flashing of the icon on-screen). */
  684.   
  685. !     if (bcmp(moonilast, moon_img, sizeof moon_img) != 0)
  686. !        bcopy(moon_img, moonilast, sizeof moon_img);
  687.   
  688.       /* Update textual information for open window */
  689.   
  690.  
  691. --- 336,351 -----
  692.   
  693.       drawmoon(p, &moon_mpr, icon_mpr);
  694.   
  695. !     if (aom_d == 0) {
  696. !         (void)sprintf(tbuf, "%dh %dm", aom_h, aom_m);
  697. !     } else {
  698. !         (void)sprintf(tbuf, "%dd %dh", aom_d, aom_h);
  699. !     }
  700. !     txt_size = pf_textwidth (strlen (tbuf), sfont, tbuf);
  701. !     tloc.pos.x = (64 - txt_size.x) / 2;
  702. !     tloc.pos.y = 63;
  703. !     tloc.pr = icon_mpr;
  704. !     pf_text (tloc, PIX_NOT(PIX_SRC), sfont, tbuf);
  705.   
  706.       /* Update textual information for open window */
  707.   
  708. ***************
  709. *** 262,269
  710.   
  711.       /* Update textual information for open window */
  712.   
  713. ! #define prt(x) pw_text(cpw, charwid, charhgt * (x), PIX_SRC, font, tbuf)
  714. ! #define prtxy(x,y) pw_text(cpw,charwid*(y+1),charhgt*(x),PIX_SRC,font,tbuf)
  715.   
  716.           sprintf(tbuf, "Moon phase:      %d%%   [0%% = New, 100%% = Full]  ",
  717.          (int) (cphase * 100));
  718.  
  719. --- 349,356 -----
  720.   
  721.       /* Update textual information for open window */
  722.   
  723. !     /* start batching updates to the pixwin */
  724. !     pw_batch_on (cpw);
  725.   
  726.       paint_labels();
  727.       /* moon phase */
  728. ***************
  729. *** 265,272
  730.   #define prt(x) pw_text(cpw, charwid, charhgt * (x), PIX_SRC, font, tbuf)
  731.   #define prtxy(x,y) pw_text(cpw,charwid*(y+1),charhgt*(x),PIX_SRC,font,tbuf)
  732.   
  733. !         sprintf(tbuf, "Moon phase:      %d%%   [0%% = New, 100%% = Full]  ",
  734. !        (int) (cphase * 100));
  735.       prt(1);
  736.   
  737.       /* Information about the Moon */
  738.  
  739. --- 352,361 -----
  740.       /* start batching updates to the pixwin */
  741.       pw_batch_on (cpw);
  742.   
  743. !     paint_labels();
  744. !     /* moon phase */
  745. !     (void)sprintf (tbuf, "%d%%   [0%% = New, 100%% = Full]  ",
  746. !         (int) (cphase * 100));
  747.       prt(1);
  748.   
  749.       /* Information about the Moon */
  750. ***************
  751. *** 271,281
  752.   
  753.       /* Information about the Moon */
  754.   
  755. ! #define EPL(x) (x), (x) == 1 ? "" : "s"
  756. !     sprintf(tbuf,
  757. !            "Age of moon:     %d day%s, %d hour%s, %d minute%s.       ",
  758. !        EPL((int) aom), EPL(((int) (24 * (aom - floor(aom))))),
  759. !        EPL(((int) (1440 * (aom - floor(aom)))) % 60));
  760.       prt(2);
  761.       sprintf(tbuf,
  762.              "Moon's distance: %ld kilometres, %.1f Earth radii.  ",
  763.  
  764. --- 360,368 -----
  765.   
  766.       /* Information about the Moon */
  767.   
  768. !     /* age of moon */
  769. !     (void)sprintf (tbuf, "%d day%s, %d hour%s, %d minute%s.       ",
  770. !         EPL(aom_d), EPL(aom_h), EPL(aom_m));
  771.       prt(2);
  772.   
  773.       /* moon distance */
  774. ***************
  775. *** 277,285
  776.          EPL((int) aom), EPL(((int) (24 * (aom - floor(aom))))),
  777.          EPL(((int) (1440 * (aom - floor(aom)))) % 60));
  778.       prt(2);
  779. !     sprintf(tbuf,
  780. !            "Moon's distance: %ld kilometres, %.1f Earth radii.  ",
  781. !        (long) cdist, cdist / earthrad);
  782.       prt(3);
  783.       sprintf(tbuf,
  784.              "Moon subtends:   %.4f degrees.       ", cangdia);
  785.  
  786. --- 364,373 -----
  787.       (void)sprintf (tbuf, "%d day%s, %d hour%s, %d minute%s.       ",
  788.           EPL(aom_d), EPL(aom_h), EPL(aom_m));
  789.       prt(2);
  790. ! 
  791. !     /* moon distance */
  792. !     (void)sprintf (tbuf, "%ld kilometres, %.1f Earth radii.  ",
  793. !         (long) cdist, cdist / earthrad);
  794.       prt(3);
  795.   
  796.       /* moon subtends */
  797. ***************
  798. *** 281,288
  799.              "Moon's distance: %ld kilometres, %.1f Earth radii.  ",
  800.          (long) cdist, cdist / earthrad);
  801.       prt(3);
  802. !     sprintf(tbuf,
  803. !            "Moon subtends:   %.4f degrees.       ", cangdia);
  804.       prt(4);
  805.       /* Draw the moon icon in the text window */
  806.       pw_rop(cpw, 60 * charwid, charhgt, 64, 64, PIX_SRC,
  807.  
  808. --- 369,377 -----
  809.       (void)sprintf (tbuf, "%ld kilometres, %.1f Earth radii.  ",
  810.           (long) cdist, cdist / earthrad);
  811.       prt(3);
  812. ! 
  813. !     /* moon subtends */
  814. !     (void)sprintf (tbuf, "%.4f degrees.       ", cangdia);
  815.       prt(4);
  816.   
  817.       /* paint the moon's image in the upper right of the canvas */
  818. ***************
  819. *** 284,292
  820.       sprintf(tbuf,
  821.              "Moon subtends:   %.4f degrees.       ", cangdia);
  822.       prt(4);
  823. -     /* Draw the moon icon in the text window */
  824. -     pw_rop(cpw, 60 * charwid, charhgt, 64, 64, PIX_SRC,
  825. -        &icon_mpr, 0, 0);
  826.   
  827.   
  828.       /* Calculate times of phases of this lunation. */
  829.  
  830. --- 373,378 -----
  831.       /* moon subtends */
  832.       (void)sprintf (tbuf, "%.4f degrees.       ", cangdia);
  833.       prt(4);
  834.   
  835.       /* paint the moon's image in the upper right of the canvas */
  836.       pw_rop (cpw, 60 * charwid, charhgt, 64, 64, PIX_SRC, icon_mpr, 0, 0);
  837. ***************
  838. *** 288,293
  839.       pw_rop(cpw, 60 * charwid, charhgt, 64, 64, PIX_SRC,
  840.          &icon_mpr, 0, 0);
  841.   
  842.   
  843.       /* Calculate times of phases of this lunation. */
  844.   
  845.  
  846. --- 374,381 -----
  847.       (void)sprintf (tbuf, "%.4f degrees.       ", cangdia);
  848.       prt(4);
  849.   
  850. +     /* paint the moon's image in the upper right of the canvas */
  851. +     pw_rop (cpw, 60 * charwid, charhgt, 64, 64, PIX_SRC, icon_mpr, 0, 0);
  852.   
  853.       /*
  854.        * Calculate times of phases of this lunation.
  855. ***************
  856. *** 289,295
  857.          &icon_mpr, 0, 0);
  858.   
  859.   
  860. !     /* Calculate times of phases of this lunation. */
  861.   
  862.   #define APOS(x) (x + 11)
  863.       phasehunt(jd, phasar);
  864.  
  865. --- 377,386 -----
  866.       /* paint the moon's image in the upper right of the canvas */
  867.       pw_rop (cpw, 60 * charwid, charhgt, 64, 64, PIX_SRC, icon_mpr, 0, 0);
  868.   
  869. !     /*
  870. !      * Calculate times of phases of this lunation.
  871. !      */
  872. !     phasehunt(jd + 0.5, phasar);
  873.   
  874.       for (i = 0; i < 5; i++) {
  875.           fmt_phase_time (phasar[i], tbuf);
  876. ***************
  877. *** 291,310
  878.   
  879.       /* Calculate times of phases of this lunation. */
  880.   
  881. ! #define APOS(x) (x + 11)
  882. !     phasehunt(jd, phasar);
  883. !     lptime = phasar[0];
  884. !     lunation = floor(((lptime + 7) - lunatbase) / synmonth) + 1;
  885. !     /* convert to seconds local time, so we can use localtime()
  886. !        to handle dst calculations */
  887. !     clock = (long)((lptime - J1970) * 24. * 3600.);
  888. !     ltm = localtime(&clock);
  889. !     sprintf(tbuf,
  890. !            "Last new moon:   %02d:%02d %s %2d %s %d             ",
  891. !        ltm->tm_hour, ltm->tm_min, timezone(tzp.tz_minuteswest, ltm->tm_isdst), ltm->tm_mday, monthnames[ltm->tm_mon], ltm->tm_year+1900);
  892. !     prt(APOS(0));
  893. !         sprintf(tbuf, "Lunation %d    ", lunation);
  894. !     prtxy(APOS(0), 52);
  895.   
  896.       lptime = phasar[1];
  897.       /* convert to seconds local time, so we can use localtime()
  898.  
  899. --- 382,391 -----
  900.        */
  901.       phasehunt(jd + 0.5, phasar);
  902.   
  903. !     for (i = 0; i < 5; i++) {
  904. !         fmt_phase_time (phasar[i], tbuf);
  905. !         prt(APOS(i));
  906. !     }
  907.   
  908.       lunation = floor(((phasar[0] + 7) - lunatbase) / synmonth) + 1;
  909.       (void)sprintf(tbuf, "Lunation %d ", lunation);
  910. ***************
  911. *** 306,320
  912.           sprintf(tbuf, "Lunation %d    ", lunation);
  913.       prtxy(APOS(0), 52);
  914.   
  915. !     lptime = phasar[1];
  916. !     /* convert to seconds local time, so we can use localtime()
  917. !        to handle dst calculations */
  918. !     clock = (long)((lptime - J1970) * 24. * 3600.);
  919. !     ltm = localtime(&clock);
  920. !     sprintf(tbuf,
  921. !            "First quarter:   %02d:%02d %s %2d %s %d              ",
  922. !        ltm->tm_hour, ltm->tm_min, timezone(tzp.tz_minuteswest, ltm->tm_isdst), ltm->tm_mday, monthnames[ltm->tm_mon], ltm->tm_year+1900);
  923. !     prt(APOS(1));
  924.   
  925.       lptime = phasar[2];
  926.       /* convert to seconds local time, so we can use localtime()
  927.  
  928. --- 387,397 -----
  929.           prt(APOS(i));
  930.       }
  931.   
  932. !     lunation = floor(((phasar[0] + 7) - lunatbase) / synmonth) + 1;
  933. !     (void)sprintf(tbuf, "Lunation %d ", lunation);
  934. !     prtxy(APOS(0), 49);                  
  935. !     (void)sprintf(tbuf, "Lunation %d ", lunation + 1);
  936. !     prtxy(APOS(4), 49);
  937.   
  938.       /* flush the pixwin updates to the screen */
  939.       pw_batch_off (cpw);
  940. ***************
  941. *** 316,330
  942.          ltm->tm_hour, ltm->tm_min, timezone(tzp.tz_minuteswest, ltm->tm_isdst), ltm->tm_mday, monthnames[ltm->tm_mon], ltm->tm_year+1900);
  943.       prt(APOS(1));
  944.   
  945. !     lptime = phasar[2];
  946. !     /* convert to seconds local time, so we can use localtime()
  947. !        to handle dst calculations */
  948. !     clock = (long)((lptime - J1970) * 24. * 3600.);
  949. !     ltm = localtime(&clock);
  950. !     sprintf(tbuf,
  951. !            "Full moon:       %02d:%02d %s %2d %s %d              ",
  952. !        ltm->tm_hour, ltm->tm_min, timezone(tzp.tz_minuteswest, ltm->tm_isdst), ltm->tm_mday, monthnames[ltm->tm_mon], ltm->tm_year+1900);
  953. !     prt(APOS(2));
  954.   
  955.       lptime = phasar[3];
  956.       /* convert to seconds local time, so we can use localtime()
  957.  
  958. --- 393,400 -----
  959.       (void)sprintf(tbuf, "Lunation %d ", lunation + 1);
  960.       prtxy(APOS(4), 49);
  961.   
  962. !     /* flush the pixwin updates to the screen */
  963. !     pw_batch_off (cpw);
  964.   
  965.       return;
  966.   }
  967. ***************
  968. *** 326,340
  969.          ltm->tm_hour, ltm->tm_min, timezone(tzp.tz_minuteswest, ltm->tm_isdst), ltm->tm_mday, monthnames[ltm->tm_mon], ltm->tm_year+1900);
  970.       prt(APOS(2));
  971.   
  972. !     lptime = phasar[3];
  973. !     /* convert to seconds local time, so we can use localtime()
  974. !        to handle dst calculations */
  975. !     clock = (long)((lptime - J1970) * 24. * 3600.);
  976. !     ltm = localtime(&clock);
  977. !     sprintf(tbuf,
  978. !            "Last quarter:    %02d:%02d %s %2d %s %d              ",
  979. !        ltm->tm_hour, ltm->tm_min, timezone(tzp.tz_minuteswest, ltm->tm_isdst), ltm->tm_mday, monthnames[ltm->tm_mon], ltm->tm_year+1900);
  980. !     prt(APOS(3));
  981.   
  982.       lptime = phasar[4];
  983.       /* convert to seconds local time, so we can use localtime()
  984.  
  985. --- 396,404 -----
  986.       /* flush the pixwin updates to the screen */
  987.       pw_batch_off (cpw);
  988.   
  989. !     return;
  990. ! }
  991. ! #undef APOS
  992.   
  993.   
  994.   /*
  995. ***************
  996. *** 336,342
  997.          ltm->tm_hour, ltm->tm_min, timezone(tzp.tz_minuteswest, ltm->tm_isdst), ltm->tm_mday, monthnames[ltm->tm_mon], ltm->tm_year+1900);
  998.       prt(APOS(3));
  999.   
  1000. !     lptime = phasar[4];
  1001.       /* convert to seconds local time, so we can use localtime()
  1002.          to handle dst calculations */
  1003.       clock = (long)rint((lptime - J1970) * 24. * 3600.);
  1004.  
  1005. --- 400,418 -----
  1006.   }
  1007.   #undef APOS
  1008.   
  1009. ! 
  1010. ! /*
  1011. !  * FMT_PHASE_TIME -- Format the provided julian date into the provided buffer
  1012. !  *        in local time format for screen display
  1013. !  */
  1014. ! 
  1015. ! fmt_phase_time (utime, buf)
  1016. !     double    utime;
  1017. !     char    *buf;
  1018. ! {
  1019. !     long    clock;
  1020. !     struct tm    *ltm, *localtime();
  1021. ! 
  1022.       /* convert to seconds local time, so we can use localtime()
  1023.          to handle dst calculations */
  1024.       clock = (long)((utime - J1970) * 24. * 3600.);
  1025. ***************
  1026. *** 339,345
  1027.       lptime = phasar[4];
  1028.       /* convert to seconds local time, so we can use localtime()
  1029.          to handle dst calculations */
  1030. !     clock = (long)rint((lptime - J1970) * 24. * 3600.);
  1031.       ltm = localtime(&clock);
  1032.       sprintf(tbuf,
  1033.              "Next new moon:   %02d:%02d %s %2d %s %d              ",
  1034.  
  1035. --- 415,421 -----
  1036.   
  1037.       /* convert to seconds local time, so we can use localtime()
  1038.          to handle dst calculations */
  1039. !     clock = (long)((utime - J1970) * 24. * 3600.);
  1040.       ltm = localtime(&clock);
  1041.       (void)sprintf (buf, " %02d:%02d %s %2d %s %d             ",
  1042.       ltm->tm_hour, ltm->tm_min, timezone(tzp.tz_minuteswest, ltm->tm_isdst), ltm->tm_mday, monthnames[ltm->tm_mon], ltm->tm_year+1900);
  1043. ***************
  1044. *** 341,353
  1045.          to handle dst calculations */
  1046.       clock = (long)rint((lptime - J1970) * 24. * 3600.);
  1047.       ltm = localtime(&clock);
  1048. !     sprintf(tbuf,
  1049. !            "Next new moon:   %02d:%02d %s %2d %s %d              ",
  1050. !        ltm->tm_hour, ltm->tm_min, timezone(tzp.tz_minuteswest, ltm->tm_isdst), ltm->tm_mday, monthnames[ltm->tm_mon], ltm->tm_year+1900);
  1051. !     prt(APOS(4));
  1052. !         sprintf(tbuf, "Lunation %d    ", lunation + 1);
  1053. !     prtxy(APOS(4), 52);
  1054. ! #undef APOS
  1055.   }
  1056.   
  1057.   /*  JDATE  --  Convert internal GMT date and time to Julian day
  1058.  
  1059. --- 417,424 -----
  1060.          to handle dst calculations */
  1061.       clock = (long)((utime - J1970) * 24. * 3600.);
  1062.       ltm = localtime(&clock);
  1063. !     (void)sprintf (buf, " %02d:%02d %s %2d %s %d             ",
  1064. !     ltm->tm_hour, ltm->tm_min, timezone(tzp.tz_minuteswest, ltm->tm_isdst), ltm->tm_mday, monthnames[ltm->tm_mon], ltm->tm_year+1900);
  1065.   }
  1066.   
  1067.   
  1068. ***************
  1069. *** 350,357
  1070.   #undef APOS
  1071.   }
  1072.   
  1073. - /*  JDATE  --  Convert internal GMT date and time to Julian day
  1074. -            and fraction.  */
  1075.   
  1076.   static long jdate(t)
  1077.   struct tm *t;
  1078.  
  1079. --- 421,426 -----
  1080.       ltm->tm_hour, ltm->tm_min, timezone(tzp.tz_minuteswest, ltm->tm_isdst), ltm->tm_mday, monthnames[ltm->tm_mon], ltm->tm_year+1900);
  1081.   }
  1082.   
  1083.   
  1084.   /*
  1085.    * PAINT_LABELS -- Draw the labels into the canvas (open) window and right
  1086. ***************
  1087. *** 353,360
  1088.   /*  JDATE  --  Convert internal GMT date and time to Julian day
  1089.              and fraction.  */
  1090.   
  1091. ! static long jdate(t)
  1092. ! struct tm *t;
  1093.   {
  1094.       long c, m, y;
  1095.   
  1096.  
  1097. --- 422,435 -----
  1098.   }
  1099.   
  1100.   
  1101. ! /*
  1102. !  * PAINT_LABELS -- Draw the labels into the canvas (open) window and right
  1103. !  *        justify them.  Done once at startup.  We paint the labels
  1104. !  *        separately to minimize the amount of screen real-estate
  1105. !  *        being modified on each update.
  1106. !  */
  1107. ! 
  1108. ! paint_labels ()
  1109.   {
  1110.       int    i;
  1111.   
  1112. ***************
  1113. *** 356,362
  1114.   static long jdate(t)
  1115.   struct tm *t;
  1116.   {
  1117. !     long c, m, y;
  1118.   
  1119.       y = t->tm_year + 1900;
  1120.       m = t->tm_mon + 1;
  1121.  
  1122. --- 431,437 -----
  1123.   
  1124.   paint_labels ()
  1125.   {
  1126. !     int    i;
  1127.   
  1128.       for (i = 0; i < 4; i++) {
  1129.           pw_text (cpw, charwid * (17 - strlen (labels [i])),
  1130. ***************
  1131. *** 358,363
  1132.   {
  1133.       long c, m, y;
  1134.   
  1135.       y = t->tm_year + 1900;
  1136.       m = t->tm_mon + 1;
  1137.       if (m > 2)
  1138.  
  1139. --- 433,463 -----
  1140.   {
  1141.       int    i;
  1142.   
  1143. +     for (i = 0; i < 4; i++) {
  1144. +         pw_text (cpw, charwid * (17 - strlen (labels [i])),
  1145. +             charhgt * (i + 1), PIX_SRC,
  1146. +             font, labels[i]);
  1147. +     }
  1148. +     for (i = 4; i < 9; i++) {
  1149. +         pw_text (cpw, charwid * (17 - strlen (labels [i])),
  1150. +             charhgt * (i + 7), PIX_SRC,
  1151. +             font, labels[i]);
  1152. +     }
  1153. + }
  1154. + 
  1155. + 
  1156. + /*
  1157. +  * JDATE  --  Convert internal GMT date and time to Julian day
  1158. +  *           and fraction.
  1159. +  */
  1160. + 
  1161. + static
  1162. + long
  1163. + jdate (t)
  1164. +     struct tm    *t;
  1165. + {
  1166. +     long        c, m, y;
  1167. + 
  1168.       y = t->tm_year + 1900;
  1169.       m = t->tm_mon + 1;
  1170.       if (m > 2) {
  1171. ***************
  1172. *** 360,370
  1173.   
  1174.       y = t->tm_year + 1900;
  1175.       m = t->tm_mon + 1;
  1176. !     if (m > 2)
  1177. !        m = m - 3;
  1178. !     else {
  1179. !        m = m + 9;
  1180. !        y--;
  1181.       }
  1182.       c = y / 100L;           /* Compute century */
  1183.       y -= 100L * c;
  1184.  
  1185. --- 460,470 -----
  1186.   
  1187.       y = t->tm_year + 1900;
  1188.       m = t->tm_mon + 1;
  1189. !     if (m > 2) {
  1190. !         m = m - 3;
  1191. !     } else {
  1192. !         m = m + 9;
  1193. !         y--;
  1194.       }
  1195.       c = y / 100L;           /* Compute century */
  1196.       y -= 100L * c;
  1197. ***************
  1198. *** 368,375
  1199.       }
  1200.       c = y / 100L;           /* Compute century */
  1201.       y -= 100L * c;
  1202. !     return t->tm_mday + (c * 146097L) / 4 + (y * 1461L) / 4 +
  1203. !         (m * 153L + 2) / 5 + 1721119L;
  1204.   }
  1205.   
  1206.   /* JTIME --    Convert internal GMT date and time to astronomical Julian
  1207.  
  1208. --- 468,475 -----
  1209.       }
  1210.       c = y / 100L;           /* Compute century */
  1211.       y -= 100L * c;
  1212. !     return (t->tm_mday + (c * 146097L) / 4 + (y * 1461L) / 4 +
  1213. !         (m * 153L + 2) / 5 + 1721119L);
  1214.   }
  1215.   
  1216.   
  1217. ***************
  1218. *** 372,380
  1219.           (m * 153L + 2) / 5 + 1721119L;
  1220.   }
  1221.   
  1222. - /* JTIME --    Convert internal GMT date and time to astronomical Julian
  1223. -            time (i.e. Julian date plus day fraction, expressed as
  1224. -            a double).  */
  1225.   
  1226.   static double jtime(t)
  1227.   struct tm *t;
  1228.  
  1229. --- 472,477 -----
  1230.           (m * 153L + 2) / 5 + 1721119L);
  1231.   }
  1232.   
  1233.   
  1234.   /*
  1235.    * JTIME --    Convert internal GMT date and time to astronomical Julian
  1236. ***************
  1237. *** 376,383
  1238.              time (i.e. Julian date plus day fraction, expressed as
  1239.              a double).  */
  1240.   
  1241. ! static double jtime(t)
  1242. ! struct tm *t;
  1243.   {
  1244.       return (jdate(t) - 0.5) + 
  1245.          (t->tm_sec + 60 * (t->tm_min + 60 * t->tm_hour)) / 86400.0;
  1246.  
  1247. --- 473,488 -----
  1248.   }
  1249.   
  1250.   
  1251. ! /*
  1252. !  * JTIME --    Convert internal GMT date and time to astronomical Julian
  1253. !  *           time (i.e. Julian date plus day fraction, expressed as
  1254. !  *           a double).
  1255. !  */
  1256. ! 
  1257. ! static
  1258. ! double
  1259. ! jtime (t)
  1260. !     struct tm *t;
  1261.   {
  1262.       return (jdate (t) - 0.5) + 
  1263.          (t->tm_sec + 60 * (t->tm_min + 60 * t->tm_hour)) / 86400.0;
  1264. ***************
  1265. *** 379,385
  1266.   static double jtime(t)
  1267.   struct tm *t;
  1268.   {
  1269. !     return (jdate(t) - 0.5) + 
  1270.          (t->tm_sec + 60 * (t->tm_min + 60 * t->tm_hour)) / 86400.0;
  1271.   }
  1272.   
  1273.  
  1274. --- 484,490 -----
  1275.   jtime (t)
  1276.       struct tm *t;
  1277.   {
  1278. !     return (jdate (t) - 0.5) + 
  1279.          (t->tm_sec + 60 * (t->tm_min + 60 * t->tm_hour)) / 86400.0;
  1280.   }
  1281.   
  1282. ***************
  1283. *** 383,390
  1284.          (t->tm_sec + 60 * (t->tm_min + 60 * t->tm_hour)) / 86400.0;
  1285.   }
  1286.   
  1287. - /*  JYEAR  --  Convert Julian date to year, month, day, which are
  1288. -            returned via integer pointers to integers.  */
  1289.   
  1290.   static void jyear(td, yy, mm, dd)
  1291.   double td;
  1292.  
  1293. --- 488,493 -----
  1294.          (t->tm_sec + 60 * (t->tm_min + 60 * t->tm_hour)) / 86400.0;
  1295.   }
  1296.   
  1297.   
  1298.   /*
  1299.    * JYEAR  --  Convert Julian date to year, month, day, which are
  1300. ***************
  1301. *** 386,394
  1302.   /*  JYEAR  --  Convert Julian date to year, month, day, which are
  1303.              returned via integer pointers to integers.  */
  1304.   
  1305. ! static void jyear(td, yy, mm, dd)
  1306. ! double td;
  1307. ! int *yy, *mm, *dd;
  1308.   {
  1309.       double j, d, y, m;
  1310.   
  1311.  
  1312. --- 489,504 -----
  1313.   }
  1314.   
  1315.   
  1316. ! /*
  1317. !  * JYEAR  --  Convert Julian date to year, month, day, which are
  1318. !  *           returned via integer pointers to integers.  
  1319. !  */
  1320. ! 
  1321. ! static
  1322. ! void
  1323. ! jyear (td, yy, mm, dd)
  1324. !     double    td;
  1325. !     int    *yy, *mm, *dd;
  1326.   {
  1327.       double j, d, y, m;
  1328.   
  1329. ***************
  1330. *** 392,398
  1331.   {
  1332.       double j, d, y, m;
  1333.   
  1334. !     td += 0.5;           /* Astronomical to civil */
  1335.       j = floor(td);
  1336.       j = j - 1721119.0;
  1337.       y = floor(((4 * j) - 1) / 146097.0);
  1338.  
  1339. --- 502,508 -----
  1340.   {
  1341.       double j, d, y, m;
  1342.   
  1343. !     td += 0.5;                /* Astronomical to civil */
  1344.       j = floor(td);
  1345.       j = j - 1721119.0;
  1346.       y = floor(((4 * j) - 1) / 146097.0);
  1347. ***************
  1348. *** 406,412
  1349.       d = floor((d + 5.0) / 5.0);
  1350.       y = (100.0 * y) + j;
  1351.       if (m < 10.0)
  1352. !        m = m + 3;
  1353.       else {
  1354.          m = m - 9;
  1355.          y = y + 1;
  1356.  
  1357. --- 516,522 -----
  1358.       d = floor((d + 5.0) / 5.0);
  1359.       y = (100.0 * y) + j;
  1360.       if (m < 10.0)
  1361. !         m = m + 3;
  1362.       else {
  1363.           m = m - 9;
  1364.           y = y + 1;
  1365. ***************
  1366. *** 408,415
  1367.       if (m < 10.0)
  1368.          m = m + 3;
  1369.       else {
  1370. !        m = m - 9;
  1371. !        y = y + 1;
  1372.       }
  1373.       *yy = y;
  1374.       *mm = m;
  1375.  
  1376. --- 518,525 -----
  1377.       if (m < 10.0)
  1378.           m = m + 3;
  1379.       else {
  1380. !         m = m - 9;
  1381. !         y = y + 1;
  1382.       }
  1383.       *yy = y;
  1384.       *mm = m;
  1385. ***************
  1386. *** 416,422
  1387.       *dd = d;
  1388.   }
  1389.   
  1390. - /*  JHMS  --  Convert Julian time to hour, minutes, and seconds.  */
  1391.   
  1392.   static void jhms(j, h, m, s)
  1393.   double j;
  1394.  
  1395. --- 526,531 -----
  1396.       *dd = d;
  1397.   }
  1398.   
  1399.   
  1400.   /*
  1401.    * JHMS  --  Convert Julian time to hour, minutes, and seconds.
  1402. ***************
  1403. *** 418,426
  1404.   
  1405.   /*  JHMS  --  Convert Julian time to hour, minutes, and seconds.  */
  1406.   
  1407. ! static void jhms(j, h, m, s)
  1408. ! double j;
  1409. ! int *h, *m, *s;
  1410.   {
  1411.       long ij;
  1412.   
  1413.  
  1414. --- 527,541 -----
  1415.   }
  1416.   
  1417.   
  1418. ! /*
  1419. !  * JHMS  --  Convert Julian time to hour, minutes, and seconds.
  1420. !  */
  1421. ! 
  1422. ! static
  1423. ! void
  1424. ! jhms(j, h, m, s)
  1425. !     double j;
  1426. !     int *h, *m, *s;
  1427.   {
  1428.       long ij;
  1429.   
  1430. ***************
  1431. *** 424,430
  1432.   {
  1433.       long ij;
  1434.   
  1435. !     j += 0.5;           /* Astronomical to civil */
  1436.       ij = (j - floor(j)) * 86400.0;
  1437.       *h = ij / 3600L;
  1438.       *m = (ij / 60L) % 60L;
  1439.  
  1440. --- 539,545 -----
  1441.   {
  1442.       long ij;
  1443.   
  1444. !     j += 0.5;                /* Astronomical to civil */
  1445.       ij = (j - floor(j)) * 86400.0;
  1446.       *h = ij / 3600L;
  1447.       *m = (ij / 60L) % 60L;
  1448. ***************
  1449. *** 431,447
  1450.       *s = ij % 60L;
  1451.   }
  1452.   
  1453. - /*  MEANPHASE  --  Calculates mean phase of the Moon for a given
  1454. -            base date and desired phase:
  1455. -                0.0   New Moon
  1456. -                0.25  First quarter
  1457. -                0.5   Full moon
  1458. -                0.75  Last quarter
  1459. -             Beware!!!  This routine returns meaningless
  1460. -                     results for any other phase arguments.  Don't
  1461. -             attempt to generalise it without understanding
  1462. -             that the motion of the moon is far more complicated
  1463. -             that this calculation reveals. */
  1464.   
  1465.   static double meanphase(sdate, phase, usek)
  1466.   double sdate, phase;
  1467.  
  1468. --- 546,551 -----
  1469.       *s = ij % 60L;
  1470.   }
  1471.   
  1472.   
  1473.   /*
  1474.    * MEANPHASE  --  Calculates mean phase of the Moon for a given
  1475. ***************
  1476. *** 443,451
  1477.               that the motion of the moon is far more complicated
  1478.               that this calculation reveals. */
  1479.   
  1480. ! static double meanphase(sdate, phase, usek)
  1481. ! double sdate, phase;
  1482. ! double *usek;
  1483.   {
  1484.       int yy, mm, dd;
  1485.       double k, t, t2, t3, nt1;
  1486.  
  1487. --- 547,567 -----
  1488.   }
  1489.   
  1490.   
  1491. ! /*
  1492. !  * MEANPHASE  --  Calculates mean phase of the Moon for a given
  1493. !  *        base date.  This argument K to this function is
  1494. !  *        the precomputed synodic month index, given by:
  1495. !  *
  1496. !  *            K = (year - 1900) * 12.3685
  1497. !  *
  1498. !  *        where year is expressed as a year and fractional
  1499. !  *        year.
  1500. !  */
  1501. ! 
  1502. ! static
  1503. ! double
  1504. ! meanphase (sdate, k)
  1505. !     double    sdate, k;
  1506.   {
  1507.       double    t, t2, t3, nt1;
  1508.   
  1509. ***************
  1510. *** 447,454
  1511.   double sdate, phase;
  1512.   double *usek;
  1513.   {
  1514. !     int yy, mm, dd;
  1515. !     double k, t, t2, t3, nt1;
  1516.   
  1517.       jyear(sdate, &yy, &mm, &dd);
  1518.   
  1519.  
  1520. --- 563,569 -----
  1521.   meanphase (sdate, k)
  1522.       double    sdate, k;
  1523.   {
  1524. !     double    t, t2, t3, nt1;
  1525.   
  1526.       /* Time in Julian centuries from 1900 January 0.5 */
  1527.       t = (sdate - 2415020.0) / 36525;
  1528. ***************
  1529. *** 450,459
  1530.       int yy, mm, dd;
  1531.       double k, t, t2, t3, nt1;
  1532.   
  1533. -     jyear(sdate, &yy, &mm, &dd);
  1534. - 
  1535. -     k = (yy + ((mm - 1) * (1.0 / 12.0)) - 1900) * 12.3685;
  1536. - 
  1537.       /* Time in Julian centuries from 1900 January 0.5 */
  1538.       t = (sdate - 2415020.0) / 36525;
  1539.       t2 = t * t;           /* Square for frequent use */
  1540.  
  1541. --- 565,570 -----
  1542.   {
  1543.       double    t, t2, t3, nt1;
  1544.   
  1545.       /* Time in Julian centuries from 1900 January 0.5 */
  1546.       t = (sdate - 2415020.0) / 36525;
  1547.       t2 = t * t;           /* Square for frequent use */
  1548. ***************
  1549. *** 459,465
  1550.       t2 = t * t;           /* Square for frequent use */
  1551.       t3 = t2 * t;           /* Cube for frequent use */
  1552.   
  1553. -     *usek = k = floor(k) + phase;
  1554.       nt1 = 2415020.75933 + synmonth * k
  1555.             + 0.0001178 * t2
  1556.             - 0.000000155 * t3
  1557.  
  1558. --- 570,575 -----
  1559.       t2 = t * t;           /* Square for frequent use */
  1560.       t3 = t2 * t;           /* Cube for frequent use */
  1561.   
  1562.       nt1 = 2415020.75933 + synmonth * k
  1563.           + 0.0001178 * t2
  1564.           - 0.000000155 * t3
  1565. ***************
  1566. *** 461,469
  1567.   
  1568.       *usek = k = floor(k) + phase;
  1569.       nt1 = 2415020.75933 + synmonth * k
  1570. !           + 0.0001178 * t2
  1571. !           - 0.000000155 * t3
  1572. !           + 0.00033 * dsin(166.56 + 132.87 * t - 0.009173 * t2);
  1573.   
  1574.       return nt1;
  1575.   }
  1576.  
  1577. --- 571,579 -----
  1578.       t3 = t2 * t;           /* Cube for frequent use */
  1579.   
  1580.       nt1 = 2415020.75933 + synmonth * k
  1581. !         + 0.0001178 * t2
  1582. !         - 0.000000155 * t3
  1583. !         + 0.00033 * dsin(166.56 + 132.87 * t - 0.009173 * t2);
  1584.   
  1585.       return nt1;
  1586.   }
  1587. ***************
  1588. *** 468,477
  1589.       return nt1;
  1590.   }
  1591.   
  1592. - /*  TRUEPHASE  --  Given a K value used to determine the
  1593. -            mean phase of the new moon, and a phase
  1594. -            selector (0.0, 0.25, 0.5, 0.75), obtain
  1595. -            the true, corrected phase time.  */
  1596.   
  1597.   static double truephase(k, phase)
  1598.   double k, phase;
  1599.  
  1600. --- 578,583 -----
  1601.       return nt1;
  1602.   }
  1603.   
  1604.   
  1605.   /*
  1606.    * TRUEPHASE  --  Given a K value used to determine the
  1607. ***************
  1608. *** 473,480
  1609.              selector (0.0, 0.25, 0.5, 0.75), obtain
  1610.              the true, corrected phase time.  */
  1611.   
  1612. ! static double truephase(k, phase)
  1613. ! double k, phase;
  1614.   {
  1615.       double t, t2, t3, pt, m, mprime, f;
  1616.       int apcor = FALSE;
  1617.  
  1618. --- 579,595 -----
  1619.   }
  1620.   
  1621.   
  1622. ! /*
  1623. !  * TRUEPHASE  --  Given a K value used to determine the
  1624. !  *        mean phase of the new moon, and a phase
  1625. !  *        selector (0.0, 0.25, 0.5, 0.75), obtain
  1626. !  *        the true, corrected phase time.
  1627. !  */
  1628. ! 
  1629. ! static
  1630. ! double
  1631. ! truephase(k, phase)
  1632. !     double k, phase;
  1633.   {
  1634.       double t, t2, t3, pt, m, mprime, f;
  1635.       int apcor = FALSE;
  1636. ***************
  1637. *** 545,551
  1638.          apcor = TRUE;
  1639.       }
  1640.       if (!apcor) {
  1641. !            fprintf(stderr, "TRUEPHASE called with invalid phase selector.\n");
  1642.          abort();
  1643.       }
  1644.       return pt;
  1645.  
  1646. --- 660,667 -----
  1647.          apcor = TRUE;
  1648.       }
  1649.       if (!apcor) {
  1650. !            (void)fprintf (stderr,
  1651. !         "TRUEPHASE called with invalid phase selector.\n");
  1652.          abort();
  1653.       }
  1654.       return pt;
  1655. ***************
  1656. *** 551,560
  1657.       return pt;
  1658.   }
  1659.   
  1660. - /*  PHASEHUNT  --  Find time of phases of the moon which surround
  1661. -            the current date.  Five phases are found, starting
  1662. -            and ending with the new moons which bound the
  1663. -            current lunation.  */
  1664.   
  1665.   static void phasehunt(sdate, phases)
  1666.   double sdate;
  1667.  
  1668. --- 667,672 -----
  1669.       return pt;
  1670.   }
  1671.   
  1672.   
  1673.   /*
  1674.    * PHASEHUNT  --  Find time of phases of the moon which surround
  1675. ***************
  1676. *** 556,564
  1677.              and ending with the new moons which bound the
  1678.              current lunation.  */
  1679.   
  1680. ! static void phasehunt(sdate, phases)
  1681. ! double sdate;
  1682. ! double phases[5];
  1683.   {
  1684.       double adate, k1, k2, nt1, nt2;
  1685.   
  1686.  
  1687. --- 668,685 -----
  1688.   }
  1689.   
  1690.   
  1691. ! /*
  1692. !  * PHASEHUNT  --  Find time of phases of the moon which surround
  1693. !  *        the current date.  Five phases are found, starting
  1694. !  *        and ending with the new moons which bound the
  1695. !  *        current lunation.
  1696. !  */
  1697. ! 
  1698. ! static
  1699. ! void
  1700. ! phasehunt (sdate, phases)
  1701. !     double    sdate;
  1702. !     double    phases [5];
  1703.   {
  1704.       int    yy, mm, dd;
  1705.       double    adate, k1, k2, nt1, nt2;
  1706. ***************
  1707. *** 560,566
  1708.   double sdate;
  1709.   double phases[5];
  1710.   {
  1711. !     double adate, k1, k2, nt1, nt2;
  1712.   
  1713.       adate = sdate - 45;
  1714.       nt1 = meanphase(adate, 0.0, &k1);
  1715.  
  1716. --- 681,688 -----
  1717.       double    sdate;
  1718.       double    phases [5];
  1719.   {
  1720. !     int    yy, mm, dd;
  1721. !     double    adate, k1, k2, nt1, nt2;
  1722.   
  1723.       adate = sdate - 45;
  1724.       jyear(adate, &yy, &mm, &dd);
  1725. ***************
  1726. *** 563,569
  1727.       double adate, k1, k2, nt1, nt2;
  1728.   
  1729.       adate = sdate - 45;
  1730. !     nt1 = meanphase(adate, 0.0, &k1);
  1731.       while (TRUE) {
  1732.          adate += synmonth;
  1733.          nt2 = meanphase(adate, 0.0, &k2);
  1734.  
  1735. --- 685,694 -----
  1736.       double    adate, k1, k2, nt1, nt2;
  1737.   
  1738.       adate = sdate - 45;
  1739. !     jyear(adate, &yy, &mm, &dd);
  1740. !     k1 = floor((yy + ((mm - 1) * (1.0 / 12.0)) - 1900) * 12.3685);
  1741. ! 
  1742. !     adate = nt1 = meanphase(adate, k1);
  1743.       while (TRUE) {
  1744.           adate += synmonth;
  1745.           k2 = k1 + 1;
  1746. ***************
  1747. *** 565,576
  1748.       adate = sdate - 45;
  1749.       nt1 = meanphase(adate, 0.0, &k1);
  1750.       while (TRUE) {
  1751. !        adate += synmonth;
  1752. !        nt2 = meanphase(adate, 0.0, &k2);
  1753. !        if (nt1 <= sdate && nt2 > sdate)
  1754. !           break;
  1755. !        nt1 = nt2;
  1756. !        k1 = k2;
  1757.       }
  1758.       phases[0] = truephase(k1, 0.0);
  1759.       phases[1] = truephase(k1, 0.25);
  1760.  
  1761. --- 690,702 -----
  1762.   
  1763.       adate = nt1 = meanphase(adate, k1);
  1764.       while (TRUE) {
  1765. !         adate += synmonth;
  1766. !         k2 = k1 + 1;
  1767. !         nt2 = meanphase(adate, k2);
  1768. !         if (nt1 <= sdate && nt2 > sdate)
  1769. !             break;
  1770. !         nt1 = nt2;
  1771. !         k1 = k2;
  1772.       }
  1773.       phases[0] = truephase(k1, 0.0);
  1774.       phases[1] = truephase(k1, 0.25);
  1775. ***************
  1776. *** 579,585
  1777.       phases[4] = truephase(k2, 0.0);
  1778.   }
  1779.   
  1780. - /*  KEPLER  --    Solve the equation of Kepler.  */
  1781.   
  1782.   static double kepler(m, ecc)
  1783.   double m, ecc;
  1784.  
  1785. --- 705,710 -----
  1786.       phases[4] = truephase(k2, 0.0);
  1787.   }
  1788.   
  1789.   
  1790.   /*
  1791.    * KEPLER  --    Solve the equation of Kepler.
  1792. ***************
  1793. *** 581,588
  1794.   
  1795.   /*  KEPLER  --    Solve the equation of Kepler.  */
  1796.   
  1797. ! static double kepler(m, ecc)
  1798. ! double m, ecc;
  1799.   {
  1800.       double e, delta;
  1801.   #define EPSILON 1E-6
  1802.  
  1803. --- 706,719 -----
  1804.   }
  1805.   
  1806.   
  1807. ! /*
  1808. !  * KEPLER  --    Solve the equation of Kepler.
  1809. !  */
  1810. ! 
  1811. ! static
  1812. ! double
  1813. ! kepler(m, ecc)
  1814. !     double m, ecc;
  1815.   {
  1816.       double e, delta;
  1817.   #define EPSILON 1E-6
  1818. ***************
  1819. *** 589,597
  1820.   
  1821.       e = m = torad(m);
  1822.       do {
  1823. !        delta = e - ecc * sin(e) - m;
  1824. !        e -= delta / (1 - ecc * cos(e));
  1825. !     } while (abs(delta) > EPSILON);
  1826.       return e;
  1827.   }
  1828.   
  1829.  
  1830. --- 720,728 -----
  1831.   
  1832.       e = m = torad(m);
  1833.       do {
  1834. !         delta = e - ecc * sin(e) - m;
  1835. !         e -= delta / (1 - ecc * cos(e));
  1836. !     } while (abs (delta) > EPSILON);
  1837.       return e;
  1838.   }
  1839.   
  1840. ***************
  1841. *** 595,601
  1842.       return e;
  1843.   }
  1844.   
  1845. - /*  PHASE  --  Calculate phase of moon as a fraction:
  1846.   
  1847.       The argument is the time for which the phase is requested,
  1848.       expressed as a Julian date and fraction.  Returns the terminator
  1849.  
  1850. --- 726,731 -----
  1851.       return e;
  1852.   }
  1853.   
  1854.   
  1855.   /*
  1856.    * PHASE  --  Calculate phase of moon as a fraction:
  1857. ***************
  1858. *** 597,610
  1859.   
  1860.   /*  PHASE  --  Calculate phase of moon as a fraction:
  1861.   
  1862. !     The argument is the time for which the phase is requested,
  1863. !     expressed as a Julian date and fraction.  Returns the terminator
  1864. !     phase angle as a percentage of a full circle (i.e., 0 to 1),
  1865. !     and stores into pointer arguments the illuminated fraction of
  1866. !         the Moon's disc, the Moon's age in days and fraction, the
  1867. !     distance of the Moon from the centre of the Earth, and the
  1868. !     angular diameter subtended by the Moon as seen by an observer
  1869. !     at the centre of the Earth.
  1870.   
  1871.   */
  1872.   
  1873.  
  1874. --- 727,744 -----
  1875.   }
  1876.   
  1877.   
  1878. ! /*
  1879. !  * PHASE  --  Calculate phase of moon as a fraction:
  1880. !  *
  1881. !  *    The argument is the time for which the phase is requested,
  1882. !  *    expressed as a Julian date and fraction.  Returns the terminator
  1883. !  *    phase angle as a percentage of a full circle (i.e., 0 to 1),
  1884. !  *    and stores into pointer arguments the illuminated fraction of
  1885. !  *    the Moon's disc, the Moon's age in days and fraction, the
  1886. !  *    distance of the Moon from the centre of the Earth, and the
  1887. !  *    angular diameter subtended by the Moon as seen by an observer
  1888. !  *    at the centre of the Earth.
  1889. !  */
  1890.   
  1891.   static
  1892.   double
  1893. ***************
  1894. *** 606,621
  1895.       angular diameter subtended by the Moon as seen by an observer
  1896.       at the centre of the Earth.
  1897.   
  1898. ! */
  1899. ! 
  1900. ! static double phase(pdate, pphase, mage, dist, angdia, sudist, suangdia)
  1901. ! double pdate;
  1902. ! double *pphase;            /* Illuminated fraction */
  1903. ! double *mage;               /* Age of moon in days */
  1904. ! double *dist;               /* Distance in kilometres */
  1905. ! double *angdia;            /* Angular diameter in degrees */
  1906. ! double *sudist;            /* Distance to Sun */
  1907. ! double *suangdia;                  /* Sun's angular diameter */
  1908.   {
  1909.   
  1910.       double Day, N, M, Ec, Lambdasun, ml, MM, MN, Ev, Ae, A3, MmP,
  1911.  
  1912. --- 740,755 -----
  1913.    *    at the centre of the Earth.
  1914.    */
  1915.   
  1916. ! static
  1917. ! double
  1918. ! phase (pdate, pphase, mage, dist, angdia, sudist, suangdia)
  1919. !     double    pdate;
  1920. !     double    *pphase;        /* Illuminated fraction */
  1921. !     double    *mage;            /* Age of moon in days */
  1922. !     double    *dist;            /* Distance in kilometres */
  1923. !     double    *angdia;        /* Angular diameter in degrees */
  1924. !     double    *sudist;        /* Distance to Sun */
  1925. !     double    *suangdia;        /* Sun's angular diameter */
  1926.   {
  1927.   
  1928.       double    Day, N, M, Ec, Lambdasun, ml, MM, MN, Ev, Ae, A3, MmP,
  1929. ***************
  1930. *** 618,628
  1931.   double *suangdia;                  /* Sun's angular diameter */
  1932.   {
  1933.   
  1934. !     double Day, N, M, Ec, Lambdasun, ml, MM, MN, Ev, Ae, A3, MmP,
  1935. !            mEc, A4, lP, V, lPP, NP, y, x, Lambdamoon, BetaM,
  1936. !            MoonAge, MoonPhase,
  1937. !            MoonDist, MoonDFrac, MoonAng, MoonPar,
  1938. !            F, SunDist, SunAng;
  1939.   
  1940.           /* Calculation of the Sun's position */
  1941.   
  1942.  
  1943. --- 752,762 -----
  1944.       double    *suangdia;        /* Sun's angular diameter */
  1945.   {
  1946.   
  1947. !     double    Day, N, M, Ec, Lambdasun, ml, MM, MN, Ev, Ae, A3, MmP,
  1948. !         mEc, A4, lP, V, lPP, NP, y, x, Lambdamoon, BetaM,
  1949. !         MoonAge, MoonPhase,
  1950. !         MoonDist, MoonDFrac, MoonAng, MoonPar,
  1951. !         F, SunDist, SunAng;
  1952.   
  1953.           /* Calculation of the Sun's position */
  1954.   
  1955. ***************
  1956. *** 626,636
  1957.   
  1958.           /* Calculation of the Sun's position */
  1959.   
  1960. !     Day = pdate - epoch;        /* Date within epoch */
  1961. !     N = fixangle((360 / 365.2422) * Day); /* Mean anomaly of the Sun */
  1962. !     M = fixangle(N + elonge - elongp);    /* Convert from perigee
  1963. !                        co-ordinates to epoch 1980.0 */
  1964. !     Ec = kepler(M, eccent);     /* Solve equation of Kepler */
  1965.       Ec = sqrt((1 + eccent) / (1 - eccent)) * tan(Ec / 2);
  1966.       Ec = 2 * todeg(atan(Ec));   /* True anomaly */
  1967.           Lambdasun = fixangle(Ec + elongp);  /* Sun's geocentric ecliptic
  1968.  
  1969. --- 760,770 -----
  1970.   
  1971.           /* Calculation of the Sun's position */
  1972.   
  1973. !     Day = pdate - epoch;            /* Date within epoch */
  1974. !     N = fixangle((360 / 365.2422) * Day);    /* Mean anomaly of the Sun */
  1975. !     M = fixangle(N + elonge - elongp);    /* Convert from perigee
  1976. !                          co-ordinates to epoch 1980.0 */
  1977. !     Ec = kepler(M, eccent);            /* Solve equation of Kepler */
  1978.       Ec = sqrt((1 + eccent) / (1 - eccent)) * tan(Ec / 2);
  1979.       Ec = 2 * todeg(atan(Ec));        /* True anomaly */
  1980.           Lambdasun = fixangle(Ec + elongp);    /* Sun's geocentric ecliptic
  1981. ***************
  1982. *** 632,640
  1983.                          co-ordinates to epoch 1980.0 */
  1984.       Ec = kepler(M, eccent);     /* Solve equation of Kepler */
  1985.       Ec = sqrt((1 + eccent) / (1 - eccent)) * tan(Ec / 2);
  1986. !     Ec = 2 * todeg(atan(Ec));   /* True anomaly */
  1987. !         Lambdasun = fixangle(Ec + elongp);  /* Sun's geocentric ecliptic
  1988. !                            longitude */
  1989.       /* Orbital distance factor */
  1990.       F = ((1 + eccent * cos(torad(Ec))) / (1 - eccent * eccent));
  1991.       SunDist = sunsmax / F;        /* Distance to Sun in km */
  1992.  
  1993. --- 766,774 -----
  1994.                            co-ordinates to epoch 1980.0 */
  1995.       Ec = kepler(M, eccent);            /* Solve equation of Kepler */
  1996.       Ec = sqrt((1 + eccent) / (1 - eccent)) * tan(Ec / 2);
  1997. !     Ec = 2 * todeg(atan(Ec));        /* True anomaly */
  1998. !         Lambdasun = fixangle(Ec + elongp);    /* Sun's geocentric ecliptic
  1999. !                             longitude */
  2000.       /* Orbital distance factor */
  2001.       F = ((1 + eccent * cos(torad(Ec))) / (1 - eccent * eccent));
  2002.       SunDist = sunsmax / F;            /* Distance to Sun in km */
  2003. ***************
  2004. *** 637,644
  2005.                              longitude */
  2006.       /* Orbital distance factor */
  2007.       F = ((1 + eccent * cos(torad(Ec))) / (1 - eccent * eccent));
  2008. !     SunDist = sunsmax / F;        /* Distance to Sun in km */
  2009. !         SunAng = F * sunangsiz;     /* Sun's angular size in degrees */
  2010.   
  2011.   
  2012.           /* Calculation of the Moon's position */
  2013.  
  2014. --- 771,778 -----
  2015.                               longitude */
  2016.       /* Orbital distance factor */
  2017.       F = ((1 + eccent * cos(torad(Ec))) / (1 - eccent * eccent));
  2018. !     SunDist = sunsmax / F;            /* Distance to Sun in km */
  2019. !         SunAng = F * sunangsiz;        /* Sun's angular size in degrees */
  2020.   
  2021.   
  2022.           /* Calculation of the Moon's position */
  2023. ***************
  2024. *** 725,728
  2025.       *suangdia = SunAng;
  2026.       return fixangle(MoonAge) / 360.0;
  2027.   }
  2028. ! #endif    /* NO_SUN_MOON */
  2029.  
  2030. --- 859,862 -----
  2031.       *suangdia = SunAng;
  2032.       return fixangle(MoonAge) / 360.0;
  2033.   }
  2034. ! #endif  /* NO_SUN_MOON */
  2035.