john@cs.su.oz.au (John Mackin) (10/01/90)
Submitted-by: John Mackin <john@cs.su.oz.au> Posting-number: Volume 9, Issue 62 Archive-name: sunclock/part01 # This is a shell archive. Remove anything before this line, # then unpack it by saving it in a file and typing "sh file". # # Wrapped by john on Sun Sep 30 21:41:46 EST 1990 # Contents: README sunclock.man Makefile.dist Imakefile sunclock.h bitmaps.c # bitmaps.h icon.bits.h large.bits.h astro.c sunclock.c patchlevel.h echo x - README sed 's/^@//' > "README" <<'@//E*O*F README//' Sun Clock - X11 Version 1.0 X11 version by John Mackin, <john@cs.su.oz.AU>, based on a Suntools program by John Walker, <kelvin@acad.uu.NET>. This program is a clock that shows which portion of the Earth's surface is illuminated by the Sun. It is designed to be usually iconic, but can be opened for a larger display with the time updated every second and both the local timezone and UTC displayed. The Suntools version had a menu that allowed you to speed up time, show different dates, etc., but I have never implemented any of that in the X version; there is some support for it in the code, however. The program should have been written using Xt rather than raw Xlib, and I tried that at first. I couldn't get the icon window to work, though, so I abandoned that version. Hence, the program does _not_ accept most Xt-style command line options; valid options are given in the manual page. The program has been tested on a number of different clients (MIPS, Sun, VAX) and servers (MIPS colour console, Sun MIT mono, NCD-16/19), under X11R3 and R4. If you find bugs, please report them to me, john@cs.su.oz.AU. To build the program, decide whether you want to use imake or not. If you don't, link Makefile.dist to Makefile, and edit it, otherwise edit the Imakefile and use xmkmf. The program uses two fonts, one for the icon display and one for the large window display. The names of these fonts are compiled in. The font names comes from -DBIGFONT and -DSMALLFONT in the Makefile. If you are not using imake, you will need to supply a -DSYSV option if you are System V. You should supply -DNEW_CTIME if you are using the table-driven ctime and a "struct tm" contains a tm_zone field. If neither of the ways of getting the local timezone name works on your system, please let me know. The original Suntools program, in case you want it, was posted as Volume 1, Issue 79 of comp.sources.sun. This program is public domain and may be freely copied as long as the notices at the top of sunclock.c remain intact. @//E*O*F README// chmod u=rw,g=,o= README echo x - sunclock.man sed 's/^@//' > "sunclock.man" <<'@//E*O*F sunclock.man//' @.TH SUNCLOCK 1 @.SH NAME sunclock \- show illuminated portion of Earth @.SH SYNOPSIS @.B sunclock [ @.I options ] @.SH DESCRIPTION @.I sunclock is an X11 application that displays a map of the Earth and indicates the illuminated portion of the globe by drawing sunlit areas dark on light, night areas as light on dark. It also displays the current date and time. @.SH OPTIONS The program does not use the X Toolkit, and hence only those options explicitly enumerated below may be used. Neither are any resources used. @.TP @.B \-version Show program version and patch level, and then exit. @.TP @.BI "\-display " dispname Give the name of the X server to contact. @.TP @.B \-iconic Start iconic. Usually you will want to specify this option. @.TP @.BI "\-geometry " +x+y Specify the position of the window. Only the position is used; any size information given is ignored. @.TP @.BI "\-icongeometry " +x+y Specify the position of the icon. Only the position is used; any size information given is ignored. @.SH "HOW IT WORKS" @.I sunclock calculates the position of the Sun using the algorithm in chapter 18 of: @.PP @.I "Astronomical Formulae for Calculators" by Jean Meeus, Third Edition, Richmond: Willmann-Bell, 1985. @.PP and Mercator projects the illuminated area onto map image. The Sun's position is calculated to better than one arc-second in accuracy. @.SH BUGS The maps are not rescaled if you resize the open window or icon. (I don't want to have to store a vector database for the map.) The program contravenes section 4.1.9 of the ICCCM in that its icon window is a fixed size, and any WM_ICON_SIZE property of the root window is ignored. @.PP The illuminated area shown is the area which would be sunlit if the Earth had no atmosphere. The actual illuminated area is larger because of atmospheric refraction and twilight. @.SH AUTHOR John Mackin, Basser Department of Computer Science, University of Sydney, Sydney, Australia; <john@cs.su.oz.AU>. Based on a Suntools program by John Walker, Autodesk, Inc.; <kelvin@acad.uu.NET>. @//E*O*F sunclock.man// chmod u=rw,g=,o= sunclock.man echo x - Makefile.dist sed 's/^@//' > "Makefile.dist" <<'@//E*O*F Makefile.dist//' # # Compile options: # # Set -DSYSV if on System V # Set -DNEW_CTIME if using the table-driven version of ctime (i.e., if # your struct tm contains a tm_zone field) # Set -DBIGFONT= and -DSMALLFONT= whatever fonts you like # Pick -O or -g # CFLAGS=-DBIGFONT=\"9x15\" -DSMALLFONT=\"6x10\" -g OFILES=sunclock.o bitmaps.o astro.o sunclock: $(OFILES) cc $(CFLAGS) $(OFILES) -o $@ -lX11 -lm @//E*O*F Makefile.dist// chmod u=rw,g=,o= Makefile.dist echo x - Imakefile sed 's/^@//' > "Imakefile" <<'@//E*O*F Imakefile//' # # Compile options: # # Set -DNEW_CTIME if using the table-driven version of ctime (i.e., if # your struct tm contains a tm_zone field) # Set -DBIGFONT= and -DSMALLFONT= whatever fonts you like # Pick -O or -g # CDEBUGFLAGS=-DBIGFONT=\"9x15\" -DSMALLFONT=\"6x10\" -g SRCS=sunclock.c bitmaps.c astro.c OBJS=sunclock.o bitmaps.o astro.o LOCAL_LIBRARIES=$(XLIB) -lm DEPLIBS= # I don't know what I'm doing, and I hate imake! ComplexProgramTarget(sunclock) @//E*O*F Imakefile// chmod u=rw,g=,o= Imakefile echo x - sunclock.h sed 's/^@//' > "sunclock.h" <<'@//E*O*F sunclock.h//' /* * Sun clock definitions. */ #include <X11/Xlib.h> #include <X11/Xutil.h> #include <stdio.h> #include <math.h> #include <time.h> #include "patchlevel.h" #include "bitmaps.h" #define abs(x) ((x) < 0 ? (-(x)) : x) /* Absolute value */ #define sgn(x) (((x) < 0) ? -1 : ((x) > 0 ? 1 : 0)) /* Extract sign */ #define dtr(x) ((x) * (PI / 180.0)) /* Degree->Radian */ #define rtd(x) ((x) / (PI / 180.0)) /* Radian->Degree */ #define fixangle(a) ((a) - 360.0 * (floor((a) / 360.0))) /* Fix angle */ #define PI 3.14159265358979323846 #define TERMINC 100 /* Circle segments for terminator */ #define PROJINT (60 * 10) /* Frequency of seasonal recalculation in seconds. */ @//E*O*F sunclock.h// chmod u=rw,g=,o= sunclock.h echo x - bitmaps.c sed 's/^@//' > "bitmaps.c" <<'@//E*O*F bitmaps.c//' /* * bitmaps for sun clock, X11 version */ #define DEFINE_BITS #include "bitmaps.h" @//E*O*F bitmaps.c// chmod u=rw,g=,o= bitmaps.c echo x - bitmaps.h sed 's/^@//' > "bitmaps.h" <<'@//E*O*F bitmaps.h//' /* * bitmaps for sun clock, X11 version; header file */ #ifdef DEFINE_BITS #define EXTERN #else #define EXTERN extern #endif #define icon_map_width 126 #define icon_map_height 63 EXTERN char icon_map_bits[] #ifdef DEFINE_BITS #include "icon.bits.h" #else ; #endif #define large_map_width 640 #define large_map_height 320 EXTERN char large_map_bits[] #ifdef DEFINE_BITS #include "large.bits.h" #else ; #endif @//E*O*F bitmaps.h// chmod u=rw,g=,o= bitmaps.h echo x - icon.bits.h sed 's/^@//' > "icon.bits.h" <<'@//E*O*F icon.bits.h//' = { 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,0x00,0xc0,0xff,0xff,0xff,0x07,0x00,0x00,0x10,0x00,0x00, 0x00,0x00,0x00,0x00,0x00,0x00,0xfe,0xff,0x03,0xe0,0x07,0xfe,0xf0,0x3f,0xc0, 0x0f,0x00,0x00,0x00,0x00,0x00,0xf8,0xff,0xff,0x07,0x80,0x01,0x7e,0x00,0x7c, 0xc0,0x3f,0xe0,0x01,0x00,0x00,0x00,0xf8,0xff,0x1f,0x0c,0xe0,0x01,0x00,0x00, 0xff,0x7d,0xfb,0xef,0x07,0x00,0xc1,0x1f,0xfe,0xff,0xff,0x38,0xf0,0x18,0xf0, 0x03,0xef,0x0f,0xa0,0xff,0x7f,0x16,0x73,0xf0,0xff,0xff,0xff,0x3b,0x78,0x00, 0xbc,0xff,0xff,0x03,0x06,0x00,0xe6,0x1f,0x7f,0x00,0xb8,0xfb,0xff,0x33,0xcf, 0x07,0xfe,0x7f,0xc0,0x0f,0x00,0x00,0x00,0x18,0x7c,0x06,0xf0,0xf3,0xff,0xf1, 0xc1,0x83,0xf7,0x1f,0x00,0x00,0x00,0x00,0xe0,0x1c,0xf0,0xff,0xe1,0x67,0xe8, 0xe1,0x00,0xf0,0xff,0x0f,0x00,0x00,0x00,0x80,0xff,0x07,0xc0,0xc3,0x83,0xcd, 0xbd,0x03,0x00,0x70,0xff,0x0b,0x00,0x00,0x06,0xe0,0xf0,0x00,0xf8,0x00,0x07, 0x18,0x2f,0x0e,0x00,0xf8,0x3f,0xc0,0x00,0x03,0x1e,0xe0,0x71,0x00,0x1f,0x00, 0x0f,0xb8,0x8d,0x0f,0x00,0xf8,0x00,0x60,0x00,0x02,0x0c,0x80,0x31,0x00,0x00, 0x00,0x1c,0x8e,0xdf,0x1f,0x00,0x60,0x00,0x28,0x18,0x00,0x00,0xc0,0x19,0x00, 0x00,0x00,0x08,0xc8,0xef,0x1b,0x00,0x40,0x0f,0xcf,0x8d,0x0f,0x00,0xe0,0x0f, 0x00,0x00,0x00,0x08,0x08,0xff,0x01,0x00,0xf8,0x9f,0xde,0x0d,0x03,0x00,0xb8, 0x03,0x00,0x00,0x00,0x08,0x00,0x77,0x00,0x10,0xc8,0xff,0xbd,0x03,0x00,0x80, 0xcf,0x00,0x00,0x00,0x00,0x18,0x00,0x10,0x00,0x70,0xf8,0xef,0xb7,0x01,0x00, 0x80,0xef,0x00,0x00,0x00,0x00,0x70,0x00,0x1c,0x00,0x00,0x78,0xc7,0x04,0x00, 0x00,0x80,0xfd,0x00,0x00,0x00,0x00,0xc0,0x80,0x07,0x00,0x00,0x0b,0xfc,0xc7, 0x00,0x00,0x80,0x19,0x00,0x00,0x00,0x00,0xc0,0xf1,0x06,0x00,0x00,0x0e,0x00, 0xc6,0x03,0x00,0x00,0x01,0x00,0x00,0x00,0x00,0xc0,0x13,0x0c,0x00,0x00,0x06, 0x00,0x8e,0x3f,0x00,0x80,0x01,0x00,0x00,0x80,0x00,0x00,0x17,0x0f,0x00,0x00, 0x03,0x00,0x0f,0xe7,0x70,0xf8,0x01,0x00,0x00,0x80,0x01,0x00,0xf4,0xff,0x00, 0x00,0x03,0x00,0x1c,0xc6,0xd8,0x18,0x01,0x00,0x00,0x00,0x00,0x00,0xfc,0xfb, 0x01,0x00,0x02,0x00,0x18,0x83,0x8c,0x99,0x01,0x00,0x00,0x00,0x00,0x00,0xe0, 0x03,0x00,0x00,0x03,0x0c,0xf8,0x81,0xc5,0x9b,0x03,0x00,0x00,0x00,0x00,0x00, 0x01,0xf3,0x03,0x00,0x06,0x0c,0xe0,0x00,0x47,0x9f,0x03,0x40,0x02,0x00,0x00, 0x00,0x01,0xff,0x07,0x00,0x44,0x00,0x80,0x00,0x47,0x8f,0x07,0x40,0x02,0x00, 0x00,0x00,0x00,0x0c,0x1c,0x00,0xfc,0x03,0xc4,0x00,0xc6,0xef,0x06,0x00,0x00, 0x00,0x00,0x00,0x00,0x08,0x30,0x00,0x00,0x03,0x67,0x00,0x80,0xff,0x0e,0x00, 0x00,0x00,0x00,0x00,0xc0,0x0c,0xf0,0x00,0x00,0x13,0x37,0x00,0x80,0xd7,0x5f, 0x00,0x08,0x00,0x00,0x00,0x40,0x0c,0xb0,0x03,0x00,0x12,0x1b,0x00,0x00,0xff, 0xff,0x0d,0x08,0x00,0x00,0x00,0x00,0x04,0x00,0x06,0x00,0x06,0x09,0x00,0x00, 0xfe,0x61,0x1f,0x00,0x00,0x40,0x00,0x00,0x0c,0x00,0x04,0x06,0x84,0x09,0x00, 0x00,0xf8,0xc7,0x77,0x00,0x00,0x40,0x00,0x00,0x18,0x00,0x06,0x00,0x04,0xb9, 0x00,0x00,0x80,0xbb,0x07,0x00,0x00,0x00,0x00,0x00,0x50,0x00,0x02,0x00,0x04, 0xf3,0x00,0x00,0x00,0xae,0x01,0x01,0x00,0x20,0x00,0x00,0x70,0x00,0x02,0x30, 0xcc,0xff,0x00,0x00,0x00,0xe3,0x03,0x03,0x00,0x20,0x00,0x00,0x40,0x00,0x03, 0x00,0x4c,0xe4,0x02,0x00,0xc0,0x01,0x86,0x03,0x00,0x00,0x00,0x00,0x40,0xc0, 0x01,0x00,0x0c,0x64,0x02,0x00,0x60,0x00,0x0c,0x01,0x00,0x00,0x00,0x00,0x40, 0x60,0x00,0x00,0x08,0x66,0x00,0x00,0x60,0x00,0x08,0x00,0x00,0x00,0x00,0x00, 0x60,0x60,0x00,0x00,0x18,0x03,0x00,0x00,0xe0,0xe1,0x08,0x00,0x00,0x00,0x00, 0x00,0x20,0x30,0x00,0x00,0x90,0x01,0x00,0x00,0xc0,0xfe,0x08,0x00,0x00,0x00, 0x00,0x00,0x20,0x1c,0x00,0x00,0xf0,0x00,0x00,0x00,0xc0,0x63,0x0c,0x0c,0x00, 0x00,0x00,0x00,0x20,0x0c,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0xe0,0x07,0x1c, 0x00,0x00,0x00,0x00,0x20,0x07,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x03, 0x1c,0x00,0x00,0x00,0x00,0x30,0x01,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00, 0x03,0x0f,0x00,0x00,0x00,0x00,0xb0,0x01,0x00,0x00,0x00,0x00,0x00,0x00,0x00, 0x00,0x00,0x07,0x00,0x00,0x00,0x00,0xb0,0x00,0x00,0x00,0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x01,0x00,0x00,0x00,0x00,0xf0,0x06,0x00,0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0x01,0x06,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x06,0x00, 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80,0x03, 0x00,0x00,0x00,0x00,0x1f,0xfc,0xff,0xff,0x01,0x00,0x00,0x00,0x00,0x00,0xf8, 0x03,0x00,0x00,0x08,0xff,0xe1,0x07,0x00,0x00,0x7f,0x00,0x00,0x00,0x00,0xf8, 0xfb,0x02,0x00,0xf8,0xff,0x01,0x00,0x00,0x00,0x00,0xc0,0x03,0x00,0xe0,0xff, 0x3f,0x7f,0x03,0xc0,0x0f,0x00,0x00,0x00,0x00,0x00,0x00,0xc0,0x03,0x80,0x3f, 0x06,0x00,0x80,0x1f,0x7c,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0xc0,0x01,0x60, 0x0f,0x00,0x00,0x00,0xe0,0x03,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0xe0,0x01, 0x80,0x01,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0xe0, 0x1f,0x07,0x06,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00, 0x80,0x03,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00, 0x00,0x00,0x00}; @//E*O*F icon.bits.h// chmod u=rw,g=,o= icon.bits.h echo x - large.bits.h sed 's/^@//' > "l.b.h.Z.uu" <<'@//E*O*F l.b.h.Z.uu//' begin 600 l.b.h.Z M'YV0/4#L40`"!AX8,%@81*CP8,*%#QTV9`AQ8D2*$BMJS,B1(<&-&$->'&FQ M)$B2)TU*_-@194N5(F&ZC)F2)<V7*7/BW'G38T&>,X/*'*HSIDVA19$")?KR M*-.>2J$^)>DTZ=2K5FO^E)IU:5>I5;V*Y3J6ZM:H:+&6+1F6K-NT7Q.VA;LV M[M"Y:M_F%8K7KM^Z*\_N'?S7J.#"=/5J14P8L&'&D!W+/2RYL>+`D2]K-ILY ML6>VE#=;_MRW<N>[H3^?MEI:].JFJ4?+9AW[M6O8ME6O;:W[-MC:IH/_SCU[ M*>_BOCD+[TT:>'+DRI\37[P<^M7CTRMCKWYZN_3JWIE;#V^=,?GL>L]SU^Y< M/''UWT7#=\\>/?W'Z^/S;5\>/'_[@\W77WK_Y2>>@`!2I]^`T=WG(&H)&L@" M@A+2%N&"%"YH884[W8!#0V:,`:(9(Y:X$(D9/KBAAASZ=.&#*3((X7IE1&0& MB0;=:"(>.NIDAHTH_M3C03&H%"*(8DR$8XL-ROA>@2S*!Z6*]3')U!@1S;#D M<'A@J>2.0VJIEAEE[`C#EFG%^.*$4SJY6YMK*D@EBS^RF>.61>I$@XT?WCD1 M&0V)V949??*XI(XLF7&#F3.X*66<_D'ZG)I65J0>F1,U.JB(?AID`YH3V0#: M7X0V%,.2-4ZTZ).1S5#F1I@:=,.2@E:)0Y:@2KJ?KK)1>M.M(&IJ*)BY\E@F M2V3@Z>:L#0$K:['X71;KL(,6.N2IQ(9J9W!>&C3#G$WR:I>O]*5JT!AY$KDM M&:_BP>Q&8G!ZT)+89N?L08TF*B\,JP9G+H\U-"3&DCBDR^]U0EXD;+?JOG;O MD`?+Z2B,<%9J'&6E&E2&L@OUFS$>9:1;9Z<-XQ%#<Y%]_-W([AX:),E14KP1 M#A[W&W.`%=^\(GT?WY#D0O+&T*X9HCX+8J'U\BBB3="""V*[!.KD(8CRLJPS M;NM]G$,.7W:7L],[3]S0O[4R'&;1>.!PK)`<<S3&EC:#C?4-PJY<Z+O&9HJS MN#L="7.1Y(H=>':8LC1#OS=$9`/:TYK!-<DQ_,QCH6,("W'<"$<V\$0&ZR8# MKL%NF_35DP(G]-AR1VTQA@7ZS6._25MMH[PVB!PTK?(^_/*\8L/:](QDY=!N M&`;7G?IDV9']^[A?]XXUBP]+/H/D4^=H(Z`Y2H[H3S78W.ZW1NY[_*AY$UYH MZ:1>M#Q@,Q0*_KG(KX[VO3`$S#S?-WT\Y/8Y/@ZRRV.30;8@5RR6#0EO7+H- MQ,8W.(Y\;'28V]"_DA:[0FVN8XV9H/\22`/GH6\G8@!20V30KUHYSE0`-`@9 MM)="/-0J6A%A6,36H[N]K8XN]`.@33K'E!I:KVL\09OJ9L"X)2%P?&&SEPA) MMC$090EV8=A63];7$_&-1TA0XR'^PN5`:YU/BV)9($R@%A%@V<1X:GD?CRR' MO8-$L$JDHPO+\%8GIC&E#%6S7`MG>)G328Q45)09AV*PKT"JA'Z`<XX5<[0O M-';E7S!8V$5XV,";E8:05-.;CV+80@/ND8L8P>1K'FA(2S7/81%IXD(H62`9 M@A$CM?H7Q,!8O><)ZX@>3-,I=Q*#0MT+EQLA8D-NT"X;Q`\&[;I@V@SV2I?< M8)&`,=ONSI3+7E'&ASOYU[U^64JW7611B:+7H0K)%&(*;I<L$B6/(E*#?+$- MB.:4U?>V)"P<[-&1V=R56+`)@V;*B%(3/)]+7&>&>NJ3++*;R-`BDK1:YJ^1 MP$-BI'B".,K%4%@Q\-\SG:BJ9J7K;8&BEQ5+8S5\2BN+-HQC.2\B4/+!H(,Y MTF-'(87+MA1QIC<$GA#I$[LV4G,CB6N(V:`&SJT(BP9;@NDZ-=DW]64P(LF2 MHDDOLM-JFA*H+.TAX]J(`V&]\$`%NH''ENA&4T&-C-8\$5/S![I1QK!?-"CJ M1HS9+!MA5%X;70A:0;*^?\43,V+Y5P[H><AO%L>,[U3K7QAV`S!NTV`^S=)( MT:D4AHU!>Y-422V/Q!(2J@62=/$CR,B:P(@XM(RD[1?##F>JNYF6H9*M*X@` M51IN-D9_;5/=963YR<8@T*!29*CD0KG7G&)P.PR+047=1#]WSHRTQ<+F"0W" M3WPV=&@?8HDTS22RS&Z$L=W$2!FLY5R+;=<@D<LD!O_8MR6A"R;%4N,'5;J1 MDZED!E4C[';H2JKS(8Y>[J.5$:=9QG;A30S:J]M-#9*#\TD.DOA<+:C4H\ZK MN8Z:-F642ERIDEXV)*JK')M/W^@BXXHEA)PCGR=#RM3?,N6KMUW:3Y1[1Z@> M:FBH0E6_:(9"H;XQ:5XZ(U-$NUO#FJ=UQ:4?[-HE6OYV*8;I*EB/&2E%CPW7 M?B8N;6-DF*7S2=*LA5WK,>G"9;7$"T0>:]N0>,R2,K2QL1\.L%"=*K5NVO:J M%LLK8%@6Y*UXS,!;*NX,B0;$RY09SZ1K*75=(D,TL@21'18TE3X7OFW-<D<W MV""?YMR0/0VSN'K&R)T!RU:K<J3"8QZC;`U"@X5YK&Y3G8@`O27$R+I4HIII M*:7$^,,0'UEJ;76@4A(JJ[=RVE..3AA.[P5CG:#Z*8I^GDMRF)=BC6'6>+`! MXG2Y:%/YSPRVWF)P30V#^0E16&38EUPYI#]H-AIU#:D=^91G)(4U)`=U>V.S M.2+?-Z8J4=(+[[2AW,/9>?3#23(<-!$J+Y#FR)^_SC)8]:JM<]4,IW2RUG!3 MO1#+J7MLE#-8L3Q\KB4I\Z<=F3#&OF@D"Z(*XF*Q;!D3B5[WFK9JQ22Q'(TM MSXB2&V7!F5;(QE8WF.O6)4X>G6X(*DYX>ZN-E;-TT>'[XO59;<?<WHF3[S5T M@U!ZKC$4I*@OOEZ)WYK<POK4JI?9T09"\L?%M9K*&?S:>)?19NZ5M'J]F>)- M2I6CK8JVD]OG=#S(0%@GIWF+A"APGXO;1VL#:K0O!Y.JLJI\9>TUCZ(LK,$> M;6SI0B#LHHT1E2F>PSJ180D;7U]HTO714!/6OQR>^19A^]F5=)IZ0AUVG,K0 M\@<!/@P;@P,ACG[M0[(!U,DKU'3)VW##]NZ3U?*P?1F=XVEC:+'$T*ZN>GN8 M=`Y^6ATC?`925BU?Y[MB)Z>WW`<*:FC+@?5-M2]YI7M;_RJ#QYA)=!;K!)+= ME5(\<Q'=ETK"DGYT1'K81SKEESEFUQ5!I3$+>!%*A5^FYBH"$WDMTV'[@ETM MAV;^YT+S="6993@*Z#QF$$7?I3CQ$X$;*(&R]G-G]W@7\Q=D<((@H7?WLV?: MIQ9$UF6!TB_Y%S]*MV'VAEXX&!PD1VH,55R>YBF*$S3M`DG8-E^X-G'Y86VF M%5[)5E])J$"$]16B4AHR@&W-!68Z<2_:UEH".&1)B($;\81X((>LQ43I,@.L M1(-'=X%0@V48<5KCYFQ)9B-,QGH<P6P$EA>QYBTNP2Z!(D<:"`,+1R5.!DO% MQ3">QA)T$V_NA38R<'WGI(=*@5@Q57$IMU)KMT(JH4V]A6?I=2$PYV0&1%L_ M(5^C&'Y7TV?U-8A"93G4`S>!,HEM*(J.%UAU=RZ6@U(ZX61F,UR*5WU"M4>A MIW,#HG@W(8L,]6;"^%P8EW7Q45QA<(1X4`-#HXS9-EG$F$0M0BL6]7G(6(PW MD7Z1%$#O%V;4E79*I3H?LV\A`7.$%TR%4@;YZ(VXMG#M`F*;=W")XE/X5(E7 M.(,/N%0DXS,LIHE9HH-!ARKYQA2%LTM)`VE\Y8;UMBUC0(U\<X-?$A8X4#4? MA7489UU?2)`RN(.A-'G^17^FHH(<$H[;(@;%PY$W61AF:&0=D5W`88MK]UZ^ M9B7;H3)O5$LNN#C-4H4$8S-WJ&OGMR:^LHT=P7@.>!D/PY`UUHU?.6WF"`-) M.'A<R4<')3?_Z$+1I3A+DGY19Q,;U(#BYGY?89+9]W=EYQL1EDI628)29(B_ M(D0B$Y0<00,B!Y$H1U\=`EN!V)>LMFWQ5C<8R3KI6(/3]T@SL8AC`7REHAZ( M.""QQA+J="](98H]D9E8B%!-:(2=]F:%(H]`!YE6V!@.R12V29;L11;`M&>Y MU98P`7,NJ(8V@Y<'X9J3R2)DXX2$.!'R0FEZ*2&EX8(A,9"-\99D09VAP9=% M-F55<F&@^1VJ%RA89@/HR&\QB8)SB7U#A5.]N5BW224J0U7C!X=B85\S`XK5 M(IW-21=GUG$#4E5A\9;JA)2P(GV:&5BR]T;8I)0NI&DV=&C"`9Z,-A/K=AO\ MZ!+\Z4(S-Q%^"$H>VIY=<8;JV!-\*8]<R8SB.'RO\2^?V"R\J#%OMD&)^!=R M>!M*!H^]UVU_N73%!)K@DR@,JIU=\6Y>8Q!B@#2&^:.U%X]911VR>#Y/JIL1 MP9/'M151)B\(YIO<4I_](7-B9D,]@Y'I]Z$K0YQ_44L=.B87,:*3&5!IA#3, M29,9VA!:FC8;F2-W6B46>AOSN4_'B"\!JA9R>'\ZD:9KR47W4H2OX8*\YZ/@ M5V@N,4&ZL80PNAYRB&)+B0-BF:(/N1F*_LFE+^1V*JG7'I)S/5FSW-IQ7 M!(.P&AG)A9')Y4^79**<:F,3P36/IISIQ*NF*C>3&G$840-<"4G**5_4PW*B MNIWTD5P;=Z@U9W<%)YX->A,,(V7K5Y+5A)I_>C6/.G?K009(JGZ&"H%&28.F M1ZGU5:U%\6:5%9N_225<)V##2'[A"B7C&A7RD@/I*EX3R!2[J94K!Z9]I)/9 M^6+*$2:!-I.C59%;\:_(@:&.21?<:7Y`<V\#:UK$*G8#J*UB!WR@921@!*D3 M119T2)0JM4,2NX>%@38Q4(4I:[&(-JJ'R$Z_DR'"FDK1>8\@(J<YBX(,*I,G M>K2A&)'E]F%5:%(=5)W;ZDQ7*C4&.9*H^5&H5:814II2*[(QB[3N\K1AF['% M<:R\!(S<J%V#*35R!C_WIK3(>HL3(:\YM1TH*A8;.W[?,0/E^9HZ>XTG&!84 M=(D7`5%%(87'1BU16K#H@;-,^&(ANW;]E$8DZCR_!(I?.ZL[40.3"VR%)JP: M*`-(8T$NJ[,>1BF6Y57K\U<880-8B55E2QUF@W2!>[F1H:9T,:AU)"0?MQ#* M^4#A5C*Y"8%4PEL`BF>8B$^N.U`!&X8KJQO8"!.-6I:S:S*#"QRE"K1FPI?T M%JW9L:.7JK0':EI0,[RU)+Y!.&0+UX41J48Y,!-HM;E+RA.?:WA*\:&JM*X[ MNRT(5+4<HHL[(8]\N81C0%0P<3>^E"5Y>&H+][$>)&-(R(:T2KRW:ZW"<;][ MBV<-UBP;U"_I)TNARI05<[*D>KK4.TQ(LYYH.4(+/$+Y2WK5"[X\L2]H4TNP M^T]9F1_9JQ,;S!&>*U24,Y=>E24_:[;9,:C_LH8=6[0+H51+>+`'<3X%:"I@ M>R'YN#!."J],RZ9%Z28Y;!#*NG>5BUX2;#*0"SUD)29M%FPXH%29]A3`AWHF M0X%K]XJ1NYD8(<!E/,666K:5-,-OXBGN8RUG><27=ZJ\64;[`G5&5Y[EI:8+ M)\B65#'!J38*"[C7*XE<S+*%0@:XB)]SVQ\@R1%TF27^PV,$VA;EB39J],-D M06WW2LIRRR,,Z\6X*:8@<;_`6<A_W+'T2ZA]!UAJQ##R9L:L>1'8QFOB!X6% MIV4@\;=N4F9A80.WO,F.F\&RF:I9IT8HN920!,"/,A:O5,7G$GH&`\6&?'!, MJE\_P3``B,V36<KJ&Y(`$LQN$6VDR+^W@9U`6+^D\TI?1J#O^,26YE,%!$UJ M6,O%&Z8M%A6&A,^:O,OV:,']H:G9G,N4&V7Q8P;ILJ.RU\@?)JMIF$K\-\N0 MHD8Q@%&Q>VH42HP-A-&<ZQ9VJ]&U-7`"$Q$[FIBFTD9DX)^`Q747K(Y+"#%< M&0-`/<X:[6RX3!^R'+V$F;3#W*-PNY0N2'B/EM10&JA7'!R87-`F<]2:N\,O M>Y0Z2]5YS&Z\K+52!$;?.[%@W<<W88%539F<N=0=@6W)Q44Q,+QXO:G?8<)0 MS8T$S;1C""7.IQ03U$;:U+8[T<$LT6H;+<]A40;MLEKX1#^4K,M=O+B`[=2Z M1]9;BUYMI*"'6#=8EHD_H<K\C--@7<I36QP#6GMEJ!+??!`0[,2=3<-<S5SD M*]J'%1$Y`#N[TSA;\KN=R;_@JC$G?8B<-\R?;6I*:GCNZU"YW=D-!-L3S1/2 M?,^M,]JMG;P\$7HZ>7)<-"U_[<52;`-%C),&[=UZ_-M*Q+5P)#=#^!3":!-" M-)24O<.F;0/+K)@O9;WQC<0=T=T.7=C%*L],"CJ/QJ!)V+PTJL,T*-/>LB7M MVB7(G=Z=S%-IC+O)S>&Z?1,67L/$IW<2S;?_?&I+D@-;U:\%[BTL3"7^/-.I MPY471HVET=N\0BF#BDUU.,B[+>1D,0;7_1VOFM9_O7`EA'UX],MB4<^#'9%H M&W5^(F0QKN3'V]4A@;$(FX.3_9]+:35@).4+GN5GC@>.N(+N;78I+M4A1BG: MW6&<[1OS<Q\TF^%L"</][3#0&N)'WM0,KN7E]-()SIDECM>][1WGX\O!%-IZ MS-%M!]SP;6+6K.*=-M1W[1)F3A<.%9`V0MP=WN?C[=D('MM9[H@48@/MTNDV MWA^N#N>\0[B"/-V5C!Z2#M":+K-B(7^6L6!S&-&4'AQ2G.;Q<=V5*,UO7NI` M2N@%'LP5!G/;A>5IP\L<V[]>I]7%P>HS$L*([.;7Q-`7O.R7,374SJVUW,&5 M;I[IDE$4KG5H6D:!CNF10<]<3A?J;M&)3!>@G,R:E83Q"]T@#IEDX%,H/NP7 M9KOB#<T1::&[2>YU+C?Y#@-^S5`S;L4"#\BK+=P;=.KTT:0SLEK6OK2OZQ+; M"^-#3A;VA9J&WK!,3=\-S:/7?A,3#W-=YT*;3>#!8=JHGO*"_G]94I[+:EK7 M/.(ZR^/.KH0O.FWZBM(FIKL4G/20"?5.K\8C!)X77]<T4(6[CGW%+E'RR/7- MZ8$P31F)WO632?4^CW]$VQ%<:5GB//"[._,F+\H\<9#Y-4Q7ZBM_WO,BONG$ MY]@NOQ/WXI,9#_&/J//?48E'CC2WI],9+5%C,']5C_9&[T&Q5WD&SO/C;O8C M'QD[FO7LYRU@A/2WOIB(2>I_KSH5)GS??O*CG!F?#_@S0[-T$<_-ON<T',NF M-77[/NC1W17<.;#SOONXUI[9[<J]X;YE\.WF&1$@;_RJC[H,17K:&19.:2H> M_^Z#C]>L/:.7#Q+`!\M1X>6`]5>>NO;ASTNSK]22"1+7?6%MM-*OGAEB__NY MGRZ']^L<"//"?4P5[RDLB:]7_RS&;7-_Q8'GI:NWI)^"D[$+7/=/VA0'"U7; M"H/[FGB@Z<*T//!V:5B-V%I]><HFF+]T,A-R1_RH%6]I5<FUSB?^KI=46G_] M:"0MA`@X^,Y>IF-#B`]/68:XI_/47@MK=J8-._%`7^'\(MZ&H8%YBIMAJQGX M%*)6+5(*::K]&<$(D>08WFU0>P;L_7&$)>;.]-"\JTXXP&8QA9WFF")@WB*` MC"@8T3OUAS6JX(VK5(_N!J*F:".'R$8`7'=VS`NV'^`0Z_):1>MC@JTF!8K# MAN;D'GY#>9RLK@%"D!7U'!X+6H.^H?C5MR88@T#75XB`#^^=Q;"<M,U(6"$T M<+ZA7^PI!`)-YHA&T4!OZD(@OV$G"(&72/I-$="TT0^7(%^*"X09=65/C[4O M*1)M9$G'"T$##CTX08I6X0X7!@N&H:G#Z+-:UIO"@CS2=G@0^($X\K?P>L,1 M@P\+1PKB0A`U$:+->?H*Q`\BN3NPL\@RGJZ3AJ!0-UB6YI8(B9S=BQ(/3/%- MFBG%$="@DAM81,K0(;UGN`U!X#'4="XF6/D_.3;-%I`DO!KR2,]EB=Z@5.9< ML_,LA$U=A(4]M0(UWB>$@=7/%%V8-/8%*T:I.HA3L![2PU"2\YI8B.,B%9"[ MJ02*".G8H!64&U6P>K6%DB@62I6R8X50:@GE+2M8*/H@1H!]S8PC4$1=A1.Y MG^6+3,V)E]%$7J?[AAQ2LS0\I2?2!7LXX*J9,IM^Z!`C=L*0H!3EH64`B4^A M@5T&88+,<A?DFQ%J)/3!A'0E0PS3:2J*^,]BE*MNJ//4$[GIB"4MDZF%,("X MW"'&\X`.R[,)FG0U07J6<]A^T?`3FJ!FL9:6XG9#AO%A!`X4A4A\%N(-]&'@ M33\UK+#`6!RC4<R*4D\I["DC]P>AQ34D:YXQFJ5#.@<3D&!%2S\$$#6!M^:W M%BO:#;B#DU`QPL5\`N4*8W,J@L7!!FK%U_"AVF)=Q#9`\4-MDRIT/H!/4=-Q MN]'O'47\8!SAA;S87XR/Y%G&S52D?L(-8%A@D=4<H/P%WKP/@]D@UB@C&HU< MZ!+U"5!49"8Q),P/RX%WO.)?H(J/4;&YK/#(I_@<=5".TH>9D42]]QEU8V1X MCZ[-9-BM,!B#4!,Y@12SL!^F#KR([5P(4"ME=,W$Y8@&U$!JW#DTD/609G4? M6)/"$J-=JXA>AP5VI<H((C;@'[PKCDD(I;X(X>IRH%^,#$_N;]"W#/D2D8-; MRX^6P35V!1D@JSCC.UQ+'L<<YL;VN!HC1"TA`QR2,B!(1Q&8_J$'N8F+ZMQ$ M-=S&_Y2@C62*TY&+,`Q\HQ9,'Y`,4I3K(;4%:&C*KF)>7`@!L#1X'B68'6)C M=+QWQ6'_R3H4&1HM@RT\D\*,3@*&L+?Y:AIXZV-D0/`!2G\X'"W=SPMO$>*' M[;,+PR;!4DA,=FKA!>Z$'85M^IV7O),%D%%*P#@DI`JD2/R#R6EIF<DR*>X4 MBE"!8"RA$C&.,&DEV2.+\(D)<AXU2D/#O=;CBI1^RTA-4@;]],:4)(O2CJ#Q M`4JL&@DE4F.@A$L8#_LYRPAYPO;5;;`!^2@?@1$>>1HO8I^TE9PN*H;*F:"< M=B(.T2<^4/+MBX;(?_2;7LQ+H/)#]JK@:`B[EIL<E1623_[`NT49)AYZA&YD MRT2*RN`7D9BE\$-[*NE>6DHEV0[!9&.0`5#,&7DV=O@H!V;DPXIQ\$#B2?M' M?0ZE^VI:?,,NTAJF\F]"`PW4AK'O6>K+ED@"MUV]S#)OD<TQP8BX,DWEY!LF M^P+;\,AW:`SEI<5LDD^OA\4<A.D<@!JSVHHA\3*D/X_&5""F2J`!<,SB`4U_ M^"4_8,\\=#'O/E3*=/5ZXN5.4",W3ZX)D,H&:-H<AZN:]E$X($:.4.-LI*8B MBAF3=\1-FP8<K!SQ$DT8Y;B-2W08S#SD@=N&*F-M9DV&20M160*#A.I08::S M^^)8@D;"\9J#<S,%E=90*5><9DR:YW(@.IM7,@R'G`S)8A.A<\*2QD1=:,>6 M`(Y=(1-Z2V5I"9G=SKN.Y65@YCJ):#I:'NJ<<IJM$AX$$]@B"XV\B"O`T&HB M31F$-@4EZQ0V:Z5R6DOIXQ#]&)R$"8+S8LJZ<7@0F./B"GB>XJ/\SC.&WD(/ MK2R>M.]OVH9VUW@J"7\SGE9JLH%/Y<ADXH<<&@-]#U7M3C"5+O/"V(2<BW+Y MO3LS$#`-6R=3+433V2B5?!?'J(M>\S+N"`]L#9SRP]IFSA"2GM`EQLS\\.'2 M9L/<A*1R;[8**+;54,7Y6&;4XYC('G$4"(7EZL0^S7-86=!,)<A*G.ID@TE' M.=22]*<K^1<-?25H@XXA0!#Y':BE=:(,V$EY_H7Z2;BBS904F#E1*<B7\B1H MKN/\`9\]$D(*3_JD.9?D3PPM1G1C^@N&0J1:Y6(DAN["649.)IDP8Y9K;"!' M%'NIA!<J*BO1_`*B4]&9W`T`@SCYY@Y+H8OR=MZ&$CJ9?ABN1):0)]X(T9/X MS.0%5LL9>C1_HAT)F1?X8Q3--AHK7(;*0=H*,4_-I%UUS[:IR)]9]1+I5V2D M7>&/UJHV.35S)24M5*:*378KN;05@6-UPH\KU#*(T8$G"47I3EATE$&*I='3 MQ_L(Y/IIAP*)OLT>=+,,10</1*(Q2VA&"0B*$?[H32N3T&275DPJ`<]:B/XA M2V&!!L0?<?1+7.1KX&.TTND)1]L#2H<F]8$FZ@N*ZA/!\F`HYK]+?.C%;@$B M8)D?E&=,8YDSX92R2GB4KFPI+Z4/4$DW3)#UY#<'9NP\BIRO5JI,?/H'E6@2 MU'<WDLJI1WPI:D+9XO)6$9234M2(($VWS!3U:JTR#1:%])@L,UH37:.9D[J, MM1.J4"\#+-6!SDP/B<7M6"'Q$/'<C:7)%8JQC_HSP:?!7`\=%:;F2Y@CY9C- M['28-"AZ%*:&4`/"Z;=<>SM.HE##74EFZ$)150X2ZI)BMYSA/3X,W2NI+36E M0M0=&M:6I1?SH?3L)R%35*+YZ$XL!9?`;Z8B!T?Z46VI#XV!(K5ZOM6&9Z"V M`K;)=PFUQ^'/.@H2]"G\VZN7(:B.5<4C7YCI)DV0-0!:_%7W$3W+Z=F,CI5U MQ-'5/%DQ7TF[^Z"`=#TP%M)3-YIC,<*C5-*RXIKO01DSPV:E86GT%N:_A8<E M82@[:6D@TWERR<)P3)\J9'RLK]"%H,S=&CZO423=IY]J:XG`*GI9?1A7XI$; MQ)$:KK#:,464!LT:-)6L_E930T1#`_!Q4<4IM"(YC`1)?$I^<PXFBQ(]0@4G M/2,2I/$5K[6X:D&W$`NE8II$K!(E!A2]H)8'=6O?&GE6\;X$4,>063L"A6R9 MB>L0ILAM<5>_IW-PK!AGN)(NZOH7`FM$94"0T).:13A(:C(3MMDZ)K1%"$YE MZJ%HX7"=:R_S[Q5+F,4QURK^4F&M\0]^M$/)4R4*:MU5C:%^WL?V2C5S!N43 M#A"V@<10Y+J8XL>JA)<R,S@45$6:006KI`"PH:'`W@9<&BR?YEG\<NM!PE'1 M3PF0".=-79VM53@8UHUZ7X>D)`6,A[7RN0G9*E_#J]MBKI@38XZ%V8=:J6IU M-9YGEI1&-OE%ZEXJ@>5(2Z^4>K80BYN@;.Z;CG*UR9+2F.KV3.58/+$*T\F% MU&C"9%/>&9NRA#8J(#@+]:$ZZK)S@_V0DHW9X]E3?>S9FEQ2"\(>VFCY$Q0K MW7QO@T0;JE/5^DG_PCY+A4Q+[)4&*?>A=.P/=:TE4ND)6B3+)`^K>N"C+C8D MJ%I"FE=]+<$D'0/KO!HOBUALR6Q1X(]3M<W:AF4KY;S3!..AU\U<CL%)^TC) M59RU#8XTS`[,OOIFR6%\_*J%(=P&VM13@@9KQ&FPFTQT)@8-:T5GPD%T-RB, MGD8)5%%CPQRX8ZGPP9$NVV=['PXBK+2W<&0*I0Y["G&UZX0,DA0VRXB^[C=@ M[9P'(:D\U*G*QXP+`Q:N%N6*>8'?,EQ'\6/+K<?\LM5VY8($DZL[+>Z`(%9G MM.+>)@3:%<9I2*BLLU;-V@;!>5`''=1K#30WV6)0/"`&EFVV_;DD%[BB.9!+ M^,@M0ZVR;Q)23KFI)S]K)DN(?F-UWO)+A)LS-NW)I0]]$-_:66+[H<R`Q'0V M/-=UVMQG!SF!F@/-<NF4,AS9@RML2Y%0"HM1%^9JU!BX'98L(J1;>@+>/CVK M*SK$+E[ENX,WQ3)>;R%K<2H4++CL,8TR2&]!(FEFI>VQD@*H[:V@VQ&HKH5M MKL0NGR;0.>E"^%*$)+TR-^51W=9$;37#*YFEM!;T,1J#89D,AW+RIU<WYO+8 MR2LWMRZ^8FC=4OB2#REFVE:+4O"S+E?.!C09BTUU;EZ8O5(V^([%7FDRXE*@ M\*;<<.UNU2SW7L,==,N^A)7PR10LB]<&3_JD#C>@[862CBMR7Z;:M6OVQ,L* M50L#-?9K]&6/R:6;.,")^FB'[?B5BY[.,)+=!1R1[E.P);6,D-AJ5&KF'.*K M=+2C##C4W%[?<G:EK\+["A54Z$)5Z[M_>\,&_@OQ#/9*WH>9&$@L6%7!!CA* M0"V@J4,S[:!KAHK$`V-@Y3OBL,V.DHR7T84`71)\83.I&'*T]3?P.MGU(%JN M*7M%5!PU$B5A\9M8$:_6A<`ZF.[2H+<6>UN6X'6J.!A6A.",.H5]*_U5ABXS MZLU=,EH@WF/2+<([./+&X`29506J,83!5'CN/84@2X21*AM.J@1'!-,/"WQN MWY=R6*)QN`QC8;CZ7<;PA7G#BW@$'TP9&':O:ZEEP"Q"#/A?`97Y\+`9YL": MI1"#7A[\B;O6PQV2Z#=EBLJ[ZRX4\?"MO:X8@ZVD!*A//7$DGJ>M^/"^X#Z, M>??,_7I4%.K!Z=\#?#6"6'T1Q6<XY%ZOA'J=KF9%U<6;*>]^6^5KBU^Q!*82 M%<90+L]$#)0.)_N-Q=58"2^>P\>+\:MO(3T6*QRKXW),PS9O)?[&-YAWLEAJ MS(ZU[1_F$.;Q'<-A+1Q&NL^H[<5+;@Y_X>LU@,7QJ[48Y`@F`.%U_!5PKRI. MPJU!Y79@OK$SI8A#)L.[&!.#XWK\?T506QV[I#@2CV(YG(<'<G;X0<,X]H[D MP)7S&')(ML9/KXK-8SGLDBLR/R;)I9@8YZ_+)G7I<>ZZR!UY)6=DB]%D/F\R M]LDO>9Y!9(X<6]*@2\;)(ODHRSSQNI0AT'0:*([7*,/BH0R583+JO1IZ1HH% MSMP+>3>R1N;*3M5U,(ML6I63\E,VRVZY=0PQDXR4XP-07L'7&"V["3Y#TNSQ M)8[+9UD9;R;P4F6*`['YM2>R+7=E@VR#!7,8&3`3$"#G2L5+DYDR63[(D(E2 M;N.;#)<7\UL&C?O,,UMF$2R4'7-E#LQUN3.K9KULFA7S:F[-W!@V<V;6G)II M\VD&S+A9-*/FVZR;<_-KKLVR^1[WYN'\FWES<7;-MADY!^>_C)-1S$Q6SL`Y M.AOGY,QA"1:<M<0.=C=#Y^ELJRCS=O[.L5DZZRQDFY._\G$.S]PY.M,/2*Q] MP?-L7L[9F48$7..+GMVS<#[/D&+-^>7X3)RI\WM&RJ.&.>/G^^R?"72;ZD%* M-PNS03[66_]I@;X:-<\<YV40N8\([HZPJ9H'W9D)9+R50>3-A,^8B[1\3#FR M)09@L`D)\<7F!,]M@6M+LWA6"V(%&D=GP8;>S)\#IIFOY'YN!"[\4#O@>NB: MVMEY'+-];*"#`U`D):51(#^N<@I])_*/U@W[;$BGYZOA9'@T6L3*"OKI-N87 MG8:6&?%=PZ*2^QA3C]K_XJ+=Z%[I5T#[-DZ4TZRKP)7*YLO@3F7FIA?!IVDC M-,^Y/DNL;0+"@-7&(UD*NI*6&HOC2Z?<]-#3/NRLR.@I#:50#IRQ9QB9764K M'7DNH%A\P<M!A"J["1E0->:S2H;3:Y8&'EAU28F9&?]X05SZ'LN2N.NB>RC/ M6!+$HZU50XGU,4"*,/ZE_'E`LRQG,H8_M?,X+[Y56`S%66VFN6_]6#L%N5@C MZHC$1S]S&@)V#.?>_,MI')3U,)I\R,M:3;M;2<P;\>P'>IZ^HU#GN$^3]B#M MI(:6I!E4_^>^U9X3H+4@&'6C6$P\<[7'R`K$N-1+:XXT:"?=JC6=RHK*F&_M M>!^;T/9X)"$6L!*Z(_1JR=RA9W1;3DW<-$2=1Q^SHNN72I,;:,T/:^M4/1;O M,M!4JAH#N-29)8CS_*.U#2,(FD-,SM%\LC?V7^/#3C*IK%%A%S+GZS2--'29 M4:/G5+QGI"V_:<N&M]YEY81MGA^T%=:JF-;2M)`6O:7[R,(.Q4JZ+!OM->VM MQ6\`7M02!7GA:@-7'77RU&:/=7:+OLN6G:LI,=*^.4@X6Y?M5O&TZ]K004W[ M.FU'*O(ZM)TUQ[9OJ64;3K<;':G98`T))X":$=]M4O>OIQ51$K?INF'K['5M MG]T"S^82^J,0!>JU/;@%-VY:R+O$N*5D#CVA*S>W'LX'N^>F8Y;MJSOWXB[2 MFD$B?Z7'W:X9M]J>W'L9>E5K.>V=7[?E;MR4VR:W"&UGMT^WR\;=1+:&;JM_ M++=1-_"^W=*-K!0A.U)AV_9^AMW(VW7_9ZN!-K2V1$+516X+$>V>[+LI=_?V M01W(5!9NQ7V\/[>N!GLV.VICYO,-O<WWUW;?(#IZ&V_;#;\]=[^FW_=;>L]O M_?V[\;>ZWM\`O'_S;^]=OO-W`"?@_IM\)W#37<#_MP`_X.W;@#_P"0XN^L`. "4`"[ ` end @//E*O*F l.b.h.Z.uu// chmod u=rw,g=,o= l.b.h.Z.uu uudecode <l.b.h.Z.uu uncompress <l.b.h.Z >large.bits.h rm -f l.b.h.* echo x - astro.c sed 's/^@//' > "astro.c" <<'@//E*O*F astro.c//' /* * Sun clock - astronomical routines. */ #include "sunclock.h" /* JDATE -- Convert internal GMT date and time to Julian day and fraction. */ long jdate(t) struct tm *t; { long c, m, y; y = t->tm_year + 1900; m = t->tm_mon + 1; if (m > 2) m = m - 3; else { m = m + 9; y--; } c = y / 100L; /* Compute century */ y -= 100L * c; return t->tm_mday + (c * 146097L) / 4 + (y * 1461L) / 4 + (m * 153L + 2) / 5 + 1721119L; } /* JTIME -- Convert internal GMT date and time to astronomical Julian time (i.e. Julian date plus day fraction, expressed as a double). */ double jtime(t) struct tm *t; { return (jdate(t) - 0.5) + (((long) t->tm_sec) + 60L * (t->tm_min + 60L * t->tm_hour)) / 86400.0; } /* KEPLER -- Solve the equation of Kepler. */ double kepler(m, ecc) double m, ecc; { double e, delta; #define EPSILON 1E-6 e = m = dtr(m); do { delta = e - ecc * sin(e) - m; e -= delta / (1 - ecc * cos(e)); } while (abs(delta) > EPSILON); return e; } /* SUNPOS -- Calculate position of the Sun. JD is the Julian date of the instant for which the position is desired and APPARENT should be nonzero if the apparent position (corrected for nutation and aberration) is desired. The Sun's co-ordinates are returned in RA and DEC, both specified in degrees (divide RA by 15 to obtain hours). The radius vector to the Sun in astronomical units is returned in RV and the Sun's longitude (true or apparent, as desired) is returned as degrees in SLONG. */ sunpos(jd, apparent, ra, dec, rv, slong) double jd; int apparent; double *ra, *dec, *rv, *slong; { double t, t2, t3, l, m, e, ea, v, theta, omega, eps; /* Time, in Julian centuries of 36525 ephemeris days, measured from the epoch 1900 January 0.5 ET. */ t = (jd - 2415020.0) / 36525.0; t2 = t * t; t3 = t2 * t; /* Geometric mean longitude of the Sun, referred to the mean equinox of the date. */ l = fixangle(279.69668 + 36000.76892 * t + 0.0003025 * t2); /* Sun's mean anomaly. */ m = fixangle(358.47583 + 35999.04975*t - 0.000150*t2 - 0.0000033*t3); /* Eccentricity of the Earth's orbit. */ e = 0.01675104 - 0.0000418 * t - 0.000000126 * t2; /* Eccentric anomaly. */ ea = kepler(m, e); /* True anomaly */ v = fixangle(2 * rtd(atan(sqrt((1 + e) / (1 - e)) * tan(ea / 2)))); /* Sun's true longitude. */ theta = l + v - m; /* Obliquity of the ecliptic. */ eps = 23.452294 - 0.0130125 * t - 0.00000164 * t2 + 0.000000503 * t3; /* Corrections for Sun's apparent longitude, if desired. */ if (apparent) { omega = fixangle(259.18 - 1934.142 * t); theta = theta - 0.00569 - 0.00479 * sin(dtr(omega)); eps += 0.00256 * cos(dtr(omega)); } /* Return Sun's longitude and radius vector */ *slong = theta; *rv = (1.0000002 * (1 - e * e)) / (1 + e * cos(dtr(v))); /* Determine solar co-ordinates. */ *ra = fixangle(rtd(atan2(cos(dtr(eps)) * sin(dtr(theta)), cos(dtr(theta))))); *dec = rtd(asin(sin(dtr(eps)) * sin(dtr(theta)))); } /* GMST -- Calculate Greenwich Mean Siderial Time for a given instant expressed as a Julian date and fraction. */ double gmst(jd) double jd; { double t, theta0; /* Time, in Julian centuries of 36525 ephemeris days, measured from the epoch 1900 January 0.5 ET. */ t = ((floor(jd + 0.5) - 0.5) - 2415020.0) / 36525.0; theta0 = 6.6460656 + 2400.051262 * t + 0.00002581 * t * t; t = (jd + 0.5) - (floor(jd + 0.5)); theta0 += (t * 24.0) * 1.002737908; theta0 = (theta0 - 24.0 * (floor(theta0 / 24.0))); return theta0; } @//E*O*F astro.c// chmod u=rw,g=,o= astro.c echo x - sunclock.c sed 's/^@//' > "sunclock.c" <<'@//E*O*F sunclock.c//' /* * Sun clock. X11 version by John Mackin. * * This program was derived from, and is still in part identical with, the * Suntools Sun clock program whose author's comment appears immediately * below. Please preserve both notices. * * The X11R3/4 version of this program was written by John Mackin, at the * Basser Department of Computer Science, University of Sydney, Sydney, * New South Wales, Australia; <john@cs.su.oz.AU>. This program, like * the one it was derived from, is in the public domain: `Love is the * law, love under will.' */ /* Sun clock Designed and implemented by John Walker in November of 1988. Version for the Sun Workstation. The algorithm used to calculate the position of the Sun is given in Chapter 18 of: "Astronomical Formulae for Calculators" by Jean Meeus, Third Edition, Richmond: Willmann-Bell, 1985. This book can be obtained from: Willmann-Bell P.O. Box 35025 Richmond, VA 23235 USA Phone: (804) 320-7016 This program was written by: John Walker Autodesk, Inc. 2320 Marinship Way Sausalito, CA 94965 USA Fax: (415) 389-9418 Voice: (415) 332-2344 Ext. 2829 Usenet: {sun,well,uunet}!acad!kelvin or: kelvin@acad.uu.net This program is in the public domain: "Do what thou wilt shall be the whole of the law". I'd appreciate receiving any bug fixes and/or enhancements, which I'll incorporate in future versions of the program. Please leave the original attribution information intact so that credit and blame may be properly apportioned. Revision history: 1.0 12/21/89 Initial version. 8/24/89 Finally got around to submitting. */ #define FAILFONT "fixed" #define VERSION "1.0" #include "sunclock.h" struct sunclock { int s_width; /* size of pixmap */ int s_height; Window s_window; /* associated window */ Pixmap s_pixmap; /* and pixmap */ int s_flags; /* see below */ int s_noon; /* position of noon */ short * s_wtab1; /* current width table (?) */ short * s_wtab; /* previous width table (?) */ long s_increm; /* increment for fake time */ long s_time; /* time - real or fake, see flags */ GC s_gc; /* GC for writing text into window */ char * (*s_tfunc)(); /* function to return the text */ char s_text[80]; /* and the current text that's there */ int s_textx; /* where to draw the text */ int s_texty; /* where to draw the text */ long s_projtime; /* last time we projected illumination */ int s_timeout; /* time until next image update */ struct sunclock * s_next; /* pointer to next clock context */ }; /* * bits in s_flags */ #define S_FAKE 01 /* date is fake, don't use actual time */ #define S_ANIMATE 02 /* do animation based on increment */ #define S_DIRTY 04 /* pixmap -> window copy required */ #define S_ICON 010 /* this is the icon window */ char * strrchr(); long time(); double jtime(); double gmst(); char * salloc(); char * bigtprint(); char * smalltprint(); struct sunclock * makeClockContext(); Bool evpred(); char * Wdayname[] = { "Sun", "Mon", "Tue", "Wed", "Thu", "Fri", "Sat" }; char * Monname[] = { "Jan", "Feb", "Mar", "Apr", "May", "Jun", "Jul", "Aug", "Sep", "Oct", "Nov", "Dec" }; struct geom { int mask; int x; int y; }; char * Name; char * Display_name = ""; Display * dpy; int scr; unsigned long Black; unsigned long White; GC GC_store; GC GC_invert; GC GC_bigf; GC GC_smallf; XFontStruct * SmallFont; XFontStruct * BigFont; Pixmap Mappix; Pixmap Iconpix; Window Icon; Window Clock; struct sunclock * Current; int Iconic = 0; struct geom Geom = { 0, 0, 0 }; struct geom Icongeom = { 0, 0, 0 }; main(argc, argv) int argc; register char ** argv; { char * p; Name = *argv; if (p = strrchr(Name, '/')) Name = ++p; parseArgs(argc, argv); dpy = XOpenDisplay(Display_name); if (dpy == (Display *)NULL) { fprintf(stderr, "%s: can't open display `%s'\n", Display_name); exit(1); } scr = DefaultScreen(dpy); getColors(); getFonts(); makePixmaps(); makeWindows(); makeGCs(Clock, Mappix); setAllHints(argc, argv); makeClockContexts(); XSelectInput(dpy, Clock, ExposureMask); XSelectInput(dpy, Icon, ExposureMask); XMapWindow(dpy, Clock); eventLoop(); /* * eventLoop() never returns, but one day it might, if someone adds a * menu for animation or such with a "quit" option. */ shutDown(); exit(0); } parseArgs(argc, argv) register int argc; register char ** argv; { while (--argc > 0) { ++argv; if (strcmp(*argv, "-display") == 0) { needMore(argc, argv); Display_name = *++argv; --argc; } else if (strcmp(*argv, "-iconic") == 0) Iconic++; else if (strcmp(*argv, "-geometry") == 0) { needMore(argc, argv); getGeom(*++argv, &Geom); --argc; } else if (strcmp(*argv, "-icongeometry") == 0) { needMore(argc, argv); getGeom(*++argv, &Icongeom); --argc; } else if (strcmp(*argv, "-version") == 0) { fprintf(stderr, "%s: version %s patchlevel %d\n", Name, VERSION, PATCHLEVEL); exit(0); } else usage(); } } needMore(argc, argv) register int argc; register char ** argv; { if (argc == 1) { fprintf(stderr, "%s: option `%s' requires an argument\n", Name, *argv); usage(); } } getGeom(s, g) register char * s; register struct geom * g; { register int mask; unsigned int width; unsigned int height; mask = XParseGeometry(s, &g->x, &g->y, &width, &height); if (mask == 0) { fprintf(stderr, "%s: `%s' is a bad geometry specification\n", Name, s); exit(1); } if ((mask & WidthValue) || (mask & HeightValue)) fprintf(stderr, "%s: warning: width/height in geometry `%s' ignored\n", Name, s); g->mask = mask; } /* * Free resources. */ shutDown() { XFreeGC(dpy, GC_store); XFreeGC(dpy, GC_invert); XFreeGC(dpy, GC_bigf); XFreeGC(dpy, GC_smallf); XFreeFont(dpy, BigFont); XFreeFont(dpy, SmallFont); XFreePixmap(dpy, Mappix); XFreePixmap(dpy, Iconpix); XDestroyWindow(dpy, Clock); XDestroyWindow(dpy, Icon); XCloseDisplay(dpy); } usage() { fprintf(stderr, "usage: %s [-display dispname] [-geometry +x+y] [-icongeometry +x+y] [-iconic] [-version]\n", Name); exit(1); } /* * Set up stuff the window manager will want to know. Must be done * before mapping window, but after creating it. */ setAllHints(argc, argv) int argc; char ** argv; { XClassHint xch; XSizeHints xsh; XWMHints xwmh; xch.res_name = Name; xch.res_class = "Sunclock"; XSetClassHint(dpy, Clock, &xch); XStoreName(dpy, Clock, Name); XSetCommand(dpy, Clock, argv, argc); XSetIconName(dpy, Clock, Name); xsh.flags = PSize | PMinSize | PMaxSize; if (Geom.mask & (XValue | YValue)) { xsh.x = Geom.x; xsh.y = Geom.y; xsh.flags |= USPosition; } xsh.width = xsh.min_width = xsh.max_width = large_map_width; xsh.height = xsh.min_height = xsh.max_height = large_map_height; XSetNormalHints(dpy, Clock, &xsh); xwmh.flags = InputHint | StateHint | IconWindowHint; if (Icongeom.mask & (XValue | YValue)) { xwmh.icon_x = Icongeom.x; xwmh.icon_y = Icongeom.y; xwmh.flags |= IconPositionHint; } xwmh.input = False; xwmh.initial_state = Iconic ? IconicState : NormalState; xwmh.icon_window = Icon; XSetWMHints(dpy, Clock, &xwmh); } makeWindows() { register int ht; XSetWindowAttributes xswa; register int mask; ht = icon_map_height + SmallFont->max_bounds.ascent + SmallFont->max_bounds.descent + 2; xswa.background_pixel = White; xswa.border_pixel = Black; xswa.backing_store = WhenMapped; mask = CWBackPixel | CWBorderPixel | CWBackingStore; fixGeometry(&Geom, large_map_width, large_map_height); Clock = XCreateWindow(dpy, RootWindow(dpy, scr), Geom.x, Geom.y, large_map_width, large_map_height, 3, CopyFromParent, InputOutput, CopyFromParent, mask, &xswa); fixGeometry(&Icongeom, icon_map_width, ht); Icon = XCreateWindow(dpy, RootWindow(dpy, scr), Icongeom.x, Icongeom.y, icon_map_width, ht, 1, CopyFromParent, InputOutput, CopyFromParent, mask, &xswa); } fixGeometry(g, w, h) register struct geom * g; register int w; register int h; { if (g->mask & XNegative) g->x = DisplayWidth(dpy, scr) - w + g->x; if (g->mask & YNegative) g->y = DisplayHeight(dpy, scr) - h + g->y; } makeGCs(w, p) register Window w; register Pixmap p; { XGCValues gcv; gcv.foreground = Black; gcv.background = White; GC_store = XCreateGC(dpy, w, GCForeground | GCBackground, &gcv); gcv.function = GXinvert; GC_invert = XCreateGC(dpy, p, GCForeground | GCBackground | GCFunction, &gcv); gcv.font = BigFont->fid; GC_bigf = XCreateGC(dpy, w, GCForeground | GCBackground | GCFont, &gcv); gcv.font = SmallFont->fid; GC_smallf = XCreateGC(dpy, w, GCForeground | GCBackground | GCFont, &gcv); } getColors() { XColor c; XColor e; register Status s; s = XAllocNamedColor(dpy, DefaultColormap(dpy, scr), "Black", &c, &e); if (s != (Status)1) { fprintf(stderr, "%s: warning: can't allocate color `Black'\n"); Black = BlackPixel(dpy, scr); } else Black = c.pixel; s = XAllocNamedColor(dpy, DefaultColormap(dpy, scr), "White", &c, &e); if (s != (Status)1) { fprintf(stderr, "%s: can't allocate color `White'\n"); White = WhitePixel(dpy, scr); } else White = c.pixel; } getFonts() { BigFont = XLoadQueryFont(dpy, BIGFONT); if (BigFont == (XFontStruct *)NULL) { fprintf(stderr, "%s: can't open font `%s', using `%s'\n", Name, BIGFONT, FAILFONT); BigFont = XLoadQueryFont(dpy, FAILFONT); if (BigFont == (XFontStruct *)NULL) { fprintf(stderr, "%s: can't open font `%s', giving up\n", Name, FAILFONT); exit(1); } } SmallFont = XLoadQueryFont(dpy, SMALLFONT); if (SmallFont == (XFontStruct *)NULL) { fprintf(stderr, "%s: can't open font `%s', using `%s'\n", Name, SMALLFONT, FAILFONT); SmallFont = XLoadQueryFont(dpy, FAILFONT); if (SmallFont == (XFontStruct *)NULL) { fprintf(stderr, "%s: can't open font `%s', giving up\n", Name, FAILFONT); exit(1); } } } makePixmaps() { Mappix = XCreatePixmapFromBitmapData(dpy, RootWindow(dpy, scr), large_map_bits, large_map_width, large_map_height, 0, 1, 1); Iconpix = XCreatePixmapFromBitmapData(dpy, RootWindow(dpy, scr), icon_map_bits, icon_map_width, icon_map_height, 0, 1, 1); } makeClockContexts() { register struct sunclock * s; s = makeClockContext(large_map_width, large_map_height, Clock, Mappix, GC_bigf, bigtprint, 70, large_map_height - BigFont->max_bounds.descent - 1); Current = s; s = makeClockContext(icon_map_width, icon_map_height, Icon, Iconpix, GC_smallf, smalltprint, 6, icon_map_height + SmallFont->max_bounds.ascent + 1); Current->s_next = s; s->s_flags |= S_ICON; s->s_next = Current; } struct sunclock * makeClockContext(wid, ht, win, pix, gc, fun, txx, txy) int wid; int ht; Window win; Pixmap pix; GC gc; char * (*fun)(); int txx; int txy; { register struct sunclock * s; s = (struct sunclock *)salloc(sizeof (struct sunclock)); s->s_width = wid; s->s_height = ht; s->s_window = win; s->s_pixmap = pix; s->s_flags = S_DIRTY; s->s_noon = -1; s->s_wtab = (short *)salloc((int)(ht * sizeof (short *))); s->s_wtab1 = (short *)salloc((int)(ht * sizeof (short *))); s->s_increm = 0L; s->s_time = 0L; s->s_gc = gc; s->s_tfunc = fun; s->s_timeout = 0; s->s_projtime = -1L; s->s_text[0] = '\0'; s->s_textx = txx; s->s_texty = txy; return (s); } /* * Someone is sure to wonder why the event loop is coded this way, without * using select(). The answer is that this was developed on a System V * kernel, which has select() but the call has bugs; so, I was inspired * to make it portable to systems without select(). The slight delay in * expose event processing that results from using sleep(1) rather than * alarm() is a fine payoff for not having to worry about interrupted * system calls. * * I've got to use XCheckIfEvent with a degenerate predicate rather than * XCheckMaskEvent with a mask of -1L because the latter won't collect all * types of events, notably ClientMessage and Selection events. Sigh. */ eventLoop() { XEvent ev; for (;;) { if (XCheckIfEvent(dpy, &ev, evpred, (char *)0)) switch (ev.type) { case Expose: if (ev.xexpose.count == 0) doExpose(ev.xexpose.window); break; } else { sleep(1); doTimeout(); } } } Bool evpred(d, e, a) register Display * d; register XEvent * e; register char * a; { return (True); } /* * Got an expose event for window w. Do the right thing if it's not * currently the one we're displaying. */ doExpose(w) register Window w; { if (w != Current->s_window) { Current = Current->s_next; if (w != Current->s_window) { fprintf(stderr, "%s: expose event for unknown window, id = 0x%08lx\n", w); exit(1); } setTimeout(Current); } updimage(Current); Current->s_flags |= S_DIRTY; showImage(Current); } doTimeout() { if (QLength(dpy)) return; /* ensure events processed first */ if (--Current->s_timeout <= 0) { updimage(Current); showImage(Current); setTimeout(Current); } } setTimeout(s) register struct sunclock * s; { long t; if (s->s_flags & S_ICON) { time(&t); s->s_timeout = 60 - localtime(&t)->tm_sec; } else s->s_timeout = 1; } showImage(s) register struct sunclock * s; { register char * p; struct tm lt; register struct tm * gmtp; lt = *localtime(&s->s_time); gmtp = gmtime(&s->s_time); p = (*s->s_tfunc)(<, gmtp); if (s->s_flags & S_DIRTY) { XCopyPlane(dpy, s->s_pixmap, s->s_window, GC_store, 0, 0, s->s_width, s->s_height, 0, 0, 1); if (s->s_flags & S_ICON) XClearArea(dpy, s->s_window, 0, s->s_height + 1, 0, 0, False); s->s_flags &= ~S_DIRTY; } strcpy(s->s_text, p); showText(s); } showText(s) register struct sunclock * s; { XDrawImageString(dpy, s->s_window, s->s_gc, s->s_textx, s->s_texty, s->s_text, strlen(s->s_text)); } /* --- */ /* UPDIMAGE -- Update current displayed image. */ updimage(s) register struct sunclock * s; { register int i; int xl; struct tm * ct; double jt; double sunra; double sundec; double sunrv; double sunlong; double gt; struct tm lt; short * wtab_swap; /* If this is a full repaint of the window, force complete recalculation. */ if (s->s_noon < 0) { s->s_projtime = 0; for (i = 0; i < s->s_height; i++) { s->s_wtab1[i] = -1; } } if (s->s_flags & S_FAKE) { if (s->s_flags & S_ANIMATE) s->s_time += s->s_increm; if (s->s_time < 0) s->s_time = 0; } else time(&s->s_time); lt = *localtime(&s->s_time); ct = gmtime(&s->s_time); jt = jtime(ct); sunpos(jt, False, &sunra, &sundec, &sunrv, &sunlong); gt = gmst(jt); /* Projecting the illumination curve for the current seasonal instant is costly. If we're running in real time, only do it every PROJINT seconds. */ if ((s->s_flags & S_FAKE) || s->s_projtime < 0 || (s->s_time - s->s_projtime) > PROJINT) { projillum(s->s_wtab, s->s_width, s->s_height, sundec); wtab_swap = s->s_wtab; s->s_wtab = s->s_wtab1; s->s_wtab1 = wtab_swap; s->s_projtime = s->s_time; } sunlong = fixangle(180.0 + (sunra - (gt * 15))); xl = sunlong * (s->s_width / 360.0); /* If the subsolar point has moved at least one pixel, update the illuminated area on the screen. */ if ((s->s_flags & S_FAKE) || s->s_noon != xl) { moveterm(s->s_wtab1, xl, s->s_wtab, s->s_noon, s->s_width, s->s_height, s->s_pixmap); s->s_noon = xl; s->s_flags |= S_DIRTY; } } /* PROJILLUM -- Project illuminated area on the map. */ projillum(wtab, xdots, ydots, dec) short *wtab; int xdots, ydots; double dec; { int i, ftf = True, ilon, ilat, lilon, lilat, xt; double m, x, y, z, th, lon, lat, s, c; /* Clear unoccupied cells in width table */ for (i = 0; i < ydots; i++) wtab[i] = -1; /* Build transformation for declination */ s = sin(-dtr(dec)); c = cos(-dtr(dec)); /* Increment over a semicircle of illumination */ for (th = -(PI / 2); th <= PI / 2 + 0.001; th += PI / TERMINC) { /* Transform the point through the declination rotation. */ x = -s * sin(th); y = cos(th); z = c * sin(th); /* Transform the resulting co-ordinate through the map projection to obtain screen co-ordinates. */ lon = (y == 0 && x == 0) ? 0.0 : rtd(atan2(y, x)); lat = rtd(asin(z)); ilat = ydots - (lat + 90) * (ydots / 180.0); ilon = lon * (xdots / 360.0); if (ftf) { /* First time. Just save start co-ordinate. */ lilon = ilon; lilat = ilat; ftf = False; } else { /* Trace out the line and set the width table. */ if (lilat == ilat) { wtab[(ydots - 1) - ilat] = ilon == 0 ? 1 : ilon; } else { m = ((double) (ilon - lilon)) / (ilat - lilat); for (i = lilat; i != ilat; i += sgn(ilat - lilat)) { xt = lilon + floor((m * (i - lilat)) + 0.5); wtab[(ydots - 1) - i] = xt == 0 ? 1 : xt; } } lilon = ilon; lilat = ilat; } } /* Now tweak the widths to generate full illumination for the correct pole. */ if (dec < 0.0) { ilat = ydots - 1; lilat = -1; } else { ilat = 0; lilat = 1; } for (i = ilat; i != ydots / 2; i += lilat) { if (wtab[i] != -1) { while (True) { wtab[i] = xdots / 2; if (i == ilat) break; i -= lilat; } break; } } } /* XSPAN -- Complement a span of pixels. Called with line in which pixels are contained, leftmost pixel in the line, and the number of pixels to complement. Handles wrap-around at the right edge of the screen. */ xspan(pline, leftp, npix, xdots, p) register int pline; register int leftp; register int npix; register int xdots; register Pixmap p; { leftp = leftp % xdots; if (leftp + npix > xdots) { XDrawLine(dpy, p, GC_invert, leftp, pline, xdots - 1, pline); XDrawLine(dpy, p, GC_invert, 0, pline, (leftp + npix) - (xdots + 1), pline); } else XDrawLine(dpy, p, GC_invert, leftp, pline, leftp + (npix - 1), pline); } /* MOVETERM -- Update illuminated portion of the globe. */ moveterm(wtab, noon, otab, onoon, xdots, ydots, pixmap) short *wtab, *otab; int noon, onoon, xdots, ydots; Pixmap pixmap; { int i, ol, oh, nl, nh; for (i = 0; i < ydots; i++) { /* If line is off in new width table but is set in the old table, clear it. */ if (wtab[i] < 0) { if (otab[i] >= 0) { xspan(i, ((onoon - otab[i]) + xdots) % xdots, otab[i] * 2, xdots, pixmap); } } else { /* Line is on in new width table. If it was off in the old width table, just draw it. */ if (otab[i] < 0) { xspan(i, ((noon - wtab[i]) + xdots) % xdots, wtab[i] * 2, xdots, pixmap); } else { /* If both the old and new spans were the entire screen, they're equivalent. */ if (otab[i] == wtab[i] && wtab[i] == (xdots / 2)) continue; /* The line was on in both the old and new width tables. We must adjust the difference in the span. */ ol = ((onoon - otab[i]) + xdots) % xdots; oh = (ol + otab[i] * 2) - 1; nl = ((noon - wtab[i]) + xdots) % xdots; nh = (nl + wtab[i] * 2) - 1; /* If spans are disjoint, erase old span and set new span. */ if (oh < nl || nh < ol) { xspan(i, ol, (oh - ol) + 1, xdots, pixmap); xspan(i, nl, (nh - nl) + 1, xdots, pixmap); } else { /* Clear portion(s) of old span that extend beyond end of new span. */ if (ol < nl) { xspan(i, ol, nl - ol, xdots, pixmap); ol = nl; } if (oh > nh) { xspan(i, nh + 1, oh - nh, xdots, pixmap); oh = nh; } /* Extend existing (possibly trimmed) span to correct new length. */ if (nl < ol) { xspan(i, nl, ol - nl, xdots, pixmap); } if (nh > oh) { xspan(i, oh + 1, nh - oh, xdots, pixmap); } } } } otab[i] = wtab[i]; } } char * salloc(nbytes) register int nbytes; { register char * p; char * malloc(); p = malloc((unsigned)nbytes); if (p == (char *)NULL) { fprintf(stderr, "%s: out of memory\n"); exit(1); } return (p); } char * bigtprint(ltp, gmtp) register struct tm * ltp; register struct tm * gmtp; { static char s[80]; sprintf(s, "%02d:%02d:%02d %s %s %02d %s %02d %02d:%02d:%02d UTC %s %02d %s %02d", ltp->tm_hour, ltp->tm_min, ltp->tm_sec, #ifdef NEW_CTIME ltp->tm_zone, #else tzname[ltp->tm_isdst], #endif Wdayname[ltp->tm_wday], ltp->tm_mday, Monname[ltp->tm_mon], ltp->tm_year % 100, gmtp->tm_hour, gmtp->tm_min, gmtp->tm_sec, Wdayname[gmtp->tm_wday], gmtp->tm_mday, Monname[gmtp->tm_mon], gmtp->tm_year % 100); return (s); } char * smalltprint(ltp, gmtp) register struct tm * ltp; register struct tm * gmtp; { static char s[80]; sprintf(s, "%02d:%02d %s %02d %s %02d", ltp->tm_hour, ltp->tm_min, Wdayname[ltp->tm_wday], ltp->tm_mday, Monname[ltp->tm_mon], ltp->tm_year % 100); return (s); } @//E*O*F sunclock.c// chmod u=rw,g=,o= sunclock.c echo x - patchlevel.h sed 's/^@//' > "patchlevel.h" <<'@//E*O*F patchlevel.h//' /* * Track patch level of sunclock. */ #define PATCHLEVEL 0 @//E*O*F patchlevel.h// chmod u=rw,g=,o= patchlevel.h exit 0 dan ---------------------------------------------------- O'Reilly && Associates argv@sun.com / argv@ora.com Opinions expressed reflect those of the author only.