/ Usenet 1994 January / usenetsourcesnewsgroupsinfomagicjanuary1994.iso / sources / x / volume19 / xephem / part18

< prev

next >

view JSON data

view as text

open on a Mac

open on a PC

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 4e 65 77 73 67 72 6f 75 | 70 73 3a 20 63 6f 6d 70 |Newsgrou|ps: comp|
|00000010| 2e 73 6f 75 72 63 65 73 | 2e 78 0a 46 72 6f 6d 3a |.sources|.x.From:|
|00000020| 20 65 63 64 6f 77 6e 65 | 79 40 70 6f 62 6f 78 2e | ecdowne|y@pobox.|
|00000030| 63 63 61 2e 63 72 2e 72 | 6f 63 6b 77 65 6c 6c 2e |cca.cr.r|ockwell.|
|00000040| 63 6f 6d 20 28 45 6c 77 | 6f 6f 64 20 44 6f 77 6e |com (Elw|ood Down|
|00000050| 65 79 29 0a 53 75 62 6a | 65 63 74 3a 20 76 31 39 |ey).Subj|ect: v19|
|00000060| 69 31 30 36 3a 20 20 78 | 65 70 68 65 6d 20 2d 20 |i106:  x|ephem - |
|00000070| 61 73 74 72 6f 6e 6f 6d | 69 63 61 6c 20 65 70 68 |astronom|ical eph|
|00000080| 65 6d 65 72 69 73 20 70 | 72 6f 67 72 61 6d 2c 20 |emeris p|rogram, |
|00000090| 50 61 72 74 31 38 2f 32 | 31 0a 4d 65 73 73 61 67 |Part18/2|1.Messag|
|000000a0| 65 2d 49 44 3a 20 3c 31 | 39 39 33 4d 61 79 31 30 |e-ID: <1|993May10|
|000000b0| 2e 32 32 31 33 30 36 2e | 39 35 37 38 40 73 70 61 |.221306.|9578@spa|
|000000c0| 72 6b 79 2e 69 6d 64 2e | 73 74 65 72 6c 69 6e 67 |rky.imd.|sterling|
|000000d0| 2e 63 6f 6d 3e 0a 58 2d | 4d 64 34 2d 53 69 67 6e |.com>.X-|Md4-Sign|
|000000e0| 61 74 75 72 65 3a 20 32 | 62 62 38 62 39 62 61 30 |ature: 2|bb8b9ba0|
|000000f0| 63 61 36 35 64 36 65 35 | 62 33 33 31 66 36 64 36 |ca65d6e5|b331f6d6|
|00000100| 61 37 34 30 39 36 37 0a | 44 61 74 65 3a 20 4d 6f |a740967.|Date: Mo|
|00000110| 6e 2c 20 31 30 20 4d 61 | 79 20 31 39 39 33 20 32 |n, 10 Ma|y 1993 2|
|00000120| 32 3a 31 33 3a 30 36 20 | 47 4d 54 0a 41 70 70 72 |2:13:06 |GMT.Appr|
|00000130| 6f 76 65 64 3a 20 63 68 | 72 69 73 40 73 70 61 72 |oved: ch|ris@spar|
|00000140| 6b 79 2e 69 6d 64 2e 73 | 74 65 72 6c 69 6e 67 2e |ky.imd.s|terling.|
|00000150| 63 6f 6d 0a 0a 53 75 62 | 6d 69 74 74 65 64 2d 62 |com..Sub|mitted-b|
|00000160| 79 3a 20 65 63 64 6f 77 | 6e 65 79 40 70 6f 62 6f |y: ecdow|ney@pobo|
|00000170| 78 2e 63 63 61 2e 63 72 | 2e 72 6f 63 6b 77 65 6c |x.cca.cr|.rockwel|
|00000180| 6c 2e 63 6f 6d 20 28 45 | 6c 77 6f 6f 64 20 44 6f |l.com (E|lwood Do|
|00000190| 77 6e 65 79 29 0a 50 6f | 73 74 69 6e 67 2d 6e 75 |wney).Po|sting-nu|
|000001a0| 6d 62 65 72 3a 20 56 6f | 6c 75 6d 65 20 31 39 2c |mber: Vo|lume 19,|
|000001b0| 20 49 73 73 75 65 20 31 | 30 36 0a 41 72 63 68 69 | Issue 1|06.Archi|
|000001c0| 76 65 2d 6e 61 6d 65 3a | 20 78 65 70 68 65 6d 2f |ve-name:| xephem/|
|000001d0| 70 61 72 74 31 38 0a 45 | 6e 76 69 72 6f 6e 6d 65 |part18.E|nvironme|
|000001e0| 6e 74 3a 20 58 31 31 72 | 34 2c 20 4f 53 46 2f 4d |nt: X11r|4, OSF/M|
|000001f0| 6f 74 69 66 0a 53 75 70 | 65 72 73 65 64 65 73 3a |otif.Sup|ersedes:|
|00000200| 20 78 65 70 68 65 6d 3a | 20 56 6f 6c 75 6d 65 20 | xephem:| Volume |
|00000210| 31 36 2c 20 49 73 73 75 | 65 20 31 31 32 2d 31 33 |16, Issu|e 112-13|
|00000220| 34 0a 0a 23 21 20 2f 62 | 69 6e 2f 73 68 0a 23 20 |4..#! /b|in/sh.# |
|00000230| 54 68 69 73 20 69 73 20 | 61 20 73 68 65 6c 6c 20 |This is |a shell |
|00000240| 61 72 63 68 69 76 65 2e | 20 20 52 65 6d 6f 76 65 |archive.|  Remove|
|00000250| 20 61 6e 79 74 68 69 6e | 67 20 62 65 66 6f 72 65 | anythin|g before|
|00000260| 20 74 68 69 73 20 6c 69 | 6e 65 2c 20 74 68 65 6e | this li|ne, then|
|00000270| 20 66 65 65 64 20 69 74 | 0a 23 20 69 6e 74 6f 20 | feed it|.# into |
|00000280| 61 20 73 68 65 6c 6c 20 | 76 69 61 20 22 73 68 20 |a shell |via "sh |
|00000290| 66 69 6c 65 22 20 6f 72 | 20 73 69 6d 69 6c 61 72 |file" or| similar|
|000002a0| 2e 20 20 54 6f 20 6f 76 | 65 72 77 72 69 74 65 20 |.  To ov|erwrite |
|000002b0| 65 78 69 73 74 69 6e 67 | 20 66 69 6c 65 73 2c 0a |existing| files,.|
|000002c0| 23 20 74 79 70 65 20 22 | 73 68 20 66 69 6c 65 20 |# type "|sh file |
|000002d0| 2d 63 22 2e 0a 23 20 54 | 68 65 20 74 6f 6f 6c 20 |-c"..# T|he tool |
|000002e0| 74 68 61 74 20 67 65 6e | 65 72 61 74 65 64 20 74 |that gen|erated t|
|000002f0| 68 69 73 20 61 70 70 65 | 61 72 65 64 20 69 6e 20 |his appe|ared in |
|00000300| 74 68 65 20 63 6f 6d 70 | 2e 73 6f 75 72 63 65 73 |the comp|.sources|
|00000310| 2e 75 6e 69 78 20 6e 65 | 77 73 67 72 6f 75 70 3b |.unix ne|wsgroup;|
|00000320| 0a 23 20 73 65 6e 64 20 | 6d 61 69 6c 20 74 6f 20 |.# send |mail to |
|00000330| 63 6f 6d 70 2d 73 6f 75 | 72 63 65 73 2d 75 6e 69 |comp-sou|rces-uni|
|00000340| 78 40 75 75 6e 65 74 2e | 75 75 2e 6e 65 74 20 69 |x@uunet.|uu.net i|
|00000350| 66 20 79 6f 75 20 77 61 | 6e 74 20 74 68 61 74 20 |f you wa|nt that |
|00000360| 74 6f 6f 6c 2e 0a 23 20 | 43 6f 6e 74 65 6e 74 73 |tool..# |Contents|
|00000370| 3a 20 20 63 69 72 63 75 | 6d 2e 63 20 6d 69 73 63 |:  circu|m.c misc|
|00000380| 2e 63 20 70 72 65 63 65 | 73 73 2e 63 20 70 72 65 |.c prece|ss.c pre|
|00000390| 66 65 72 65 6e 63 65 73 | 2e 63 20 73 6b 79 76 69 |ferences|.c skyvi|
|000003a0| 65 77 6d 65 6e 75 2e 63 | 2e 33 0a 23 20 57 72 61 |ewmenu.c|.3.# Wra|
|000003b0| 70 70 65 64 20 62 79 20 | 63 68 72 69 73 40 6e 6f |pped by |chris@no|
|000003c0| 76 61 20 6f 6e 20 4d 6f | 6e 20 4d 61 79 20 31 30 |va on Mo|n May 10|
|000003d0| 20 31 36 3a 34 31 3a 35 | 33 20 31 39 39 33 0a 50 | 16:41:5|3 1993.P|
|000003e0| 41 54 48 3d 2f 62 69 6e | 3a 2f 75 73 72 2f 62 69 |ATH=/bin|:/usr/bi|
|000003f0| 6e 3a 2f 75 73 72 2f 75 | 63 62 20 3b 20 65 78 70 |n:/usr/u|cb ; exp|
|00000400| 6f 72 74 20 50 41 54 48 | 0a 65 63 68 6f 20 49 66 |ort PATH|.echo If|
|00000410| 20 74 68 69 73 20 61 72 | 63 68 69 76 65 20 69 73 | this ar|chive is|
|00000420| 20 63 6f 6d 70 6c 65 74 | 65 2c 20 79 6f 75 20 77 | complet|e, you w|
|00000430| 69 6c 6c 20 73 65 65 20 | 74 68 65 20 66 6f 6c 6c |ill see |the foll|
|00000440| 6f 77 69 6e 67 20 6d 65 | 73 73 61 67 65 3a 0a 65 |owing me|ssage:.e|
|00000450| 63 68 6f 20 27 20 20 20 | 20 20 20 20 20 20 20 22 |cho '   |       "|
|00000460| 73 68 61 72 3a 20 45 6e | 64 20 6f 66 20 61 72 63 |shar: En|d of arc|
|00000470| 68 69 76 65 20 31 38 20 | 28 6f 66 20 32 31 29 2e |hive 18 |(of 21).|
|00000480| 22 27 0a 69 66 20 74 65 | 73 74 20 2d 66 20 27 63 |"'.if te|st -f 'c|
|00000490| 69 72 63 75 6d 2e 63 27 | 20 2d 61 20 22 24 7b 31 |ircum.c'| -a "${1|
|000004a0| 7d 22 20 21 3d 20 22 2d | 63 22 20 3b 20 74 68 65 |}" != "-|c" ; the|
|000004b0| 6e 20 0a 20 20 65 63 68 | 6f 20 73 68 61 72 3a 20 |n .  ech|o shar: |
|000004c0| 57 69 6c 6c 20 6e 6f 74 | 20 63 6c 6f 62 62 65 72 |Will not| clobber|
|000004d0| 20 65 78 69 73 74 69 6e | 67 20 66 69 6c 65 20 5c | existin|g file \|
|000004e0| 22 27 63 69 72 63 75 6d | 2e 63 27 5c 22 0a 65 6c |"'circum|.c'\".el|
|000004f0| 73 65 0a 20 20 65 63 68 | 6f 20 73 68 61 72 3a 20 |se.  ech|o shar: |
|00000500| 45 78 74 72 61 63 74 69 | 6e 67 20 5c 22 27 63 69 |Extracti|ng \"'ci|
|00000510| 72 63 75 6d 2e 63 27 5c | 22 20 5c 28 32 30 30 30 |rcum.c'\|" \(2000|
|00000520| 38 20 63 68 61 72 61 63 | 74 65 72 73 5c 29 0a 20 |8 charac|ters\). |
|00000530| 20 73 65 64 20 22 73 2f | 5e 58 2f 2f 22 20 3e 27 | sed "s/|^X//" >'|
|00000540| 63 69 72 63 75 6d 2e 63 | 27 20 3c 3c 27 45 4e 44 |circum.c|' <<'END|
|00000550| 5f 4f 46 5f 46 49 4c 45 | 27 0a 58 2f 2a 20 67 69 |_OF_FILE|'.X/* gi|
|00000560| 76 65 6e 20 61 20 4e 6f | 77 20 61 6e 64 20 61 6e |ven a No|w and an|
|00000570| 20 4f 62 6a 20 77 69 74 | 68 20 74 68 65 20 6f 62 | Obj wit|h the ob|
|00000580| 6a 65 63 74 20 64 65 66 | 69 6e 69 74 69 6f 6e 20 |ject def|inition |
|00000590| 70 6f 72 74 69 6f 6e 20 | 66 69 6c 6c 65 64 20 69 |portion |filled i|
|000005a0| 6e 2c 0a 58 20 2a 20 66 | 69 6c 6c 20 69 6e 20 74 |n,.X * f|ill in t|
|000005b0| 68 65 20 73 6b 79 20 70 | 6f 73 69 74 69 6f 6e 20 |he sky p|osition |
|000005c0| 28 73 5f 2a 29 20 70 6f | 72 74 69 6f 6e 73 2e 0a |(s_*) po|rtions..|
|000005d0| 58 20 2a 2f 0a 58 0a 58 | 23 69 6e 63 6c 75 64 65 |X */.X.X|#include|
|000005e0| 20 3c 73 74 64 69 6f 2e | 68 3e 0a 58 23 69 6e 63 | <stdio.|h>.X#inc|
|000005f0| 6c 75 64 65 20 3c 6d 61 | 74 68 2e 68 3e 0a 58 23 |lude <ma|th.h>.X#|
|00000600| 69 66 20 64 65 66 69 6e | 65 64 28 5f 5f 53 54 44 |if defin|ed(__STD|
|00000610| 43 5f 5f 29 0a 58 23 69 | 6e 63 6c 75 64 65 20 3c |C__).X#i|nclude <|
|00000620| 73 74 64 6c 69 62 2e 68 | 3e 0a 58 23 65 6e 64 69 |stdlib.h|>.X#endi|
|00000630| 66 0a 58 23 69 6e 63 6c | 75 64 65 20 22 61 73 74 |f.X#incl|ude "ast|
|00000640| 72 6f 2e 68 22 0a 58 23 | 69 6e 63 6c 75 64 65 20 |ro.h".X#|include |
|00000650| 22 63 69 72 63 75 6d 2e | 68 22 0a 58 23 69 6e 63 |"circum.|h".X#inc|
|00000660| 6c 75 64 65 20 22 70 72 | 65 66 65 72 65 6e 63 65 |lude "pr|eference|
|00000670| 73 2e 68 22 0a 58 0a 58 | 23 69 66 20 64 65 66 69 |s.h".X.X|#if defi|
|00000680| 6e 65 64 28 5f 5f 53 54 | 44 43 5f 5f 29 20 7c 7c |ned(__ST|DC__) |||
|00000690| 20 64 65 66 69 6e 65 64 | 28 5f 5f 63 70 6c 75 73 | defined|(__cplus|
|000006a0| 70 6c 75 73 29 0a 58 23 | 64 65 66 69 6e 65 20 50 |plus).X#|define P|
|000006b0| 5f 28 73 29 20 73 0a 58 | 23 65 6c 73 65 0a 58 23 |_(s) s.X|#else.X#|
|000006c0| 64 65 66 69 6e 65 20 50 | 5f 28 73 29 20 28 29 0a |define P|_(s) ().|
|000006d0| 58 23 65 6e 64 69 66 0a | 58 0a 58 65 78 74 65 72 |X#endif.|X.Xexter|
|000006e0| 6e 20 76 6f 69 64 20 61 | 6e 6f 6d 61 6c 79 20 50 |n void a|nomaly P|
|000006f0| 5f 28 28 64 6f 75 62 6c | 65 20 6d 61 2c 20 64 6f |_((doubl|e ma, do|
|00000700| 75 62 6c 65 20 73 2c 20 | 64 6f 75 62 6c 65 20 2a |uble s, |double *|
|00000710| 6e 75 2c 20 64 6f 75 62 | 6c 65 20 2a 65 61 29 29 |nu, doub|le *ea))|
|00000720| 3b 0a 58 65 78 74 65 72 | 6e 20 76 6f 69 64 20 63 |;.Xexter|n void c|
|00000730| 6f 6d 65 74 20 50 5f 28 | 28 64 6f 75 62 6c 65 20 |omet P_(|(double |
|00000740| 4d 6a 64 2c 20 64 6f 75 | 62 6c 65 20 65 70 2c 20 |Mjd, dou|ble ep, |
|00000750| 64 6f 75 62 6c 65 20 69 | 6e 63 2c 20 64 6f 75 62 |double i|nc, doub|
|00000760| 6c 65 20 61 70 2c 20 64 | 6f 75 62 6c 65 20 71 70 |le ap, d|ouble qp|
|00000770| 2c 20 64 6f 75 62 6c 65 | 20 6f 6d 2c 20 64 6f 75 |, double| om, dou|
|00000780| 62 6c 65 20 2a 6c 70 64 | 2c 20 64 6f 75 62 6c 65 |ble *lpd|, double|
|00000790| 20 2a 70 73 69 2c 20 64 | 6f 75 62 6c 65 20 2a 72 | *psi, d|ouble *r|
|000007a0| 70 2c 20 64 6f 75 62 6c | 65 20 2a 72 68 6f 2c 20 |p, doubl|e *rho, |
|000007b0| 64 6f 75 62 6c 65 20 2a | 6c 61 6d 2c 20 64 6f 75 |double *|lam, dou|
|000007c0| 62 6c 65 20 2a 62 65 74 | 29 29 3b 0a 58 65 78 74 |ble *bet|));.Xext|
|000007d0| 65 72 6e 20 76 6f 69 64 | 20 65 63 6c 5f 65 71 20 |ern void| ecl_eq |
|000007e0| 50 5f 28 28 64 6f 75 62 | 6c 65 20 4d 6a 64 2c 20 |P_((doub|le Mjd, |
|000007f0| 64 6f 75 62 6c 65 20 6c | 61 74 69 2c 20 64 6f 75 |double l|ati, dou|
|00000800| 62 6c 65 20 6c 6e 67 69 | 2c 20 64 6f 75 62 6c 65 |ble lngi|, double|
|00000810| 20 2a 72 61 2c 20 64 6f | 75 62 6c 65 20 2a 64 65 | *ra, do|uble *de|
|00000820| 63 29 29 3b 0a 58 65 78 | 74 65 72 6e 20 76 6f 69 |c));.Xex|tern voi|
|00000830| 64 20 65 71 5f 65 63 6c | 20 50 5f 28 28 64 6f 75 |d eq_ecl| P_((dou|
|00000840| 62 6c 65 20 4d 6a 64 2c | 20 64 6f 75 62 6c 65 20 |ble Mjd,| double |
|00000850| 72 61 2c 20 64 6f 75 62 | 6c 65 20 64 65 63 2c 20 |ra, doub|le dec, |
|00000860| 64 6f 75 62 6c 65 20 2a | 6c 61 74 69 2c 20 64 6f |double *|lati, do|
|00000870| 75 62 6c 65 20 2a 6c 6e | 67 69 29 29 3b 0a 58 65 |uble *ln|gi));.Xe|
|00000880| 78 74 65 72 6e 20 76 6f | 69 64 20 67 6b 5f 6d 61 |xtern vo|id gk_ma|
|00000890| 67 20 50 5f 28 28 64 6f | 75 62 6c 65 20 67 2c 20 |g P_((do|uble g, |
|000008a0| 64 6f 75 62 6c 65 20 6b | 2c 20 64 6f 75 62 6c 65 |double k|, double|
|000008b0| 20 72 70 2c 20 64 6f 75 | 62 6c 65 20 72 68 6f 2c | rp, dou|ble rho,|
|000008c0| 20 64 6f 75 62 6c 65 20 | 2a 6d 70 29 29 3b 0a 58 | double |*mp));.X|
|000008d0| 65 78 74 65 72 6e 20 76 | 6f 69 64 20 68 61 64 65 |extern v|oid hade|
|000008e0| 63 5f 61 61 20 50 5f 28 | 28 64 6f 75 62 6c 65 20 |c_aa P_(|(double |
|000008f0| 6c 61 74 69 2c 20 64 6f | 75 62 6c 65 20 68 61 2c |lati, do|uble ha,|
|00000900| 20 64 6f 75 62 6c 65 20 | 64 65 63 2c 20 64 6f 75 | double |dec, dou|
|00000910| 62 6c 65 20 2a 61 6c 74 | 2c 20 64 6f 75 62 6c 65 |ble *alt|, double|
|00000920| 20 2a 61 7a 29 29 3b 0a | 58 65 78 74 65 72 6e 20 | *az));.|Xextern |
|00000930| 76 6f 69 64 20 68 67 5f | 6d 61 67 20 50 5f 28 28 |void hg_|mag P_((|
|00000940| 64 6f 75 62 6c 65 20 68 | 2c 20 64 6f 75 62 6c 65 |double h|, double|
|00000950| 20 67 2c 20 64 6f 75 62 | 6c 65 20 72 70 2c 20 64 | g, doub|le rp, d|
|00000960| 6f 75 62 6c 65 20 72 68 | 6f 2c 20 64 6f 75 62 6c |ouble rh|o, doubl|
|00000970| 65 20 72 73 6e 2c 20 64 | 6f 75 62 6c 65 20 2a 6d |e rsn, d|ouble *m|
|00000980| 70 29 29 3b 0a 58 65 78 | 74 65 72 6e 20 76 6f 69 |p));.Xex|tern voi|
|00000990| 64 20 6d 6f 6f 6e 20 50 | 5f 28 28 64 6f 75 62 6c |d moon P|_((doubl|
|000009a0| 65 20 4d 6a 64 2c 20 64 | 6f 75 62 6c 65 20 2a 6c |e Mjd, d|ouble *l|
|000009b0| 61 6d 2c 20 64 6f 75 62 | 6c 65 20 2a 62 65 74 2c |am, doub|le *bet,|
|000009c0| 20 64 6f 75 62 6c 65 20 | 2a 68 70 29 29 3b 0a 58 | double |*hp));.X|
|000009d0| 65 78 74 65 72 6e 20 76 | 6f 69 64 20 6e 6f 77 5f |extern v|oid now_|
|000009e0| 6c 73 74 20 50 5f 28 28 | 4e 6f 77 20 2a 6e 70 2c |lst P_((|Now *np,|
|000009f0| 20 64 6f 75 62 6c 65 20 | 2a 6c 73 74 29 29 3b 0a | double |*lst));.|
|00000a00| 58 65 78 74 65 72 6e 20 | 76 6f 69 64 20 6e 75 74 |Xextern |void nut|
|00000a10| 61 74 69 6f 6e 20 50 5f | 28 28 64 6f 75 62 6c 65 |ation P_|((double|
|00000a20| 20 4d 6a 64 2c 20 64 6f | 75 62 6c 65 20 2a 64 65 | Mjd, do|uble *de|
|00000a30| 70 73 2c 20 64 6f 75 62 | 6c 65 20 2a 64 70 73 69 |ps, doub|le *dpsi|
|00000a40| 29 29 3b 0a 58 65 78 74 | 65 72 6e 20 76 6f 69 64 |));.Xext|ern void|
|00000a50| 20 70 6c 61 6e 73 20 50 | 5f 28 28 64 6f 75 62 6c | plans P|_((doubl|
|00000a60| 65 20 4d 6a 64 2c 20 69 | 6e 74 20 70 2c 20 64 6f |e Mjd, i|nt p, do|
|00000a70| 75 62 6c 65 20 2a 6c 70 | 64 30 2c 20 64 6f 75 62 |uble *lp|d0, doub|
|00000a80| 6c 65 20 2a 70 73 69 30 | 2c 20 64 6f 75 62 6c 65 |le *psi0|, double|
|00000a90| 20 2a 72 70 30 2c 20 64 | 6f 75 62 6c 65 20 2a 72 | *rp0, d|ouble *r|
|00000aa0| 68 6f 30 2c 20 64 6f 75 | 62 6c 65 20 2a 6c 61 6d |ho0, dou|ble *lam|
|00000ab0| 2c 20 64 6f 75 62 6c 65 | 20 2a 62 65 74 2c 20 64 |, double| *bet, d|
|00000ac0| 6f 75 62 6c 65 20 2a 64 | 69 61 2c 20 64 6f 75 62 |ouble *d|ia, doub|
|00000ad0| 6c 65 20 2a 6d 61 67 29 | 29 3b 0a 58 65 78 74 65 |le *mag)|);.Xexte|
|00000ae0| 72 6e 20 76 6f 69 64 20 | 70 72 65 63 65 73 73 20 |rn void |precess |
|00000af0| 50 5f 28 28 64 6f 75 62 | 6c 65 20 6d 6a 64 31 2c |P_((doub|le mjd1,|
|00000b00| 20 64 6f 75 62 6c 65 20 | 6d 6a 64 32 2c 20 64 6f | double |mjd2, do|
|00000b10| 75 62 6c 65 20 2a 72 61 | 2c 20 64 6f 75 62 6c 65 |uble *ra|, double|
|00000b20| 20 2a 64 65 63 29 29 3b | 0a 58 65 78 74 65 72 6e | *dec));|.Xextern|
|00000b30| 20 76 6f 69 64 20 72 61 | 6e 67 65 20 50 5f 28 28 | void ra|nge P_((|
|00000b40| 64 6f 75 62 6c 65 20 2a | 76 2c 20 64 6f 75 62 6c |double *|v, doubl|
|00000b50| 65 20 72 29 29 3b 0a 58 | 65 78 74 65 72 6e 20 76 |e r));.X|extern v|
|00000b60| 6f 69 64 20 72 65 64 75 | 63 65 5f 65 6c 65 6d 65 |oid redu|ce_eleme|
|00000b70| 6e 74 73 20 50 5f 28 28 | 64 6f 75 62 6c 65 20 6d |nts P_((|double m|
|00000b80| 6a 64 30 2c 20 64 6f 75 | 62 6c 65 20 4d 6a 64 2c |jd0, dou|ble Mjd,|
|00000b90| 20 64 6f 75 62 6c 65 20 | 69 6e 63 30 2c 20 64 6f | double |inc0, do|
|00000ba0| 75 62 6c 65 20 61 70 30 | 2c 20 64 6f 75 62 6c 65 |uble ap0|, double|
|00000bb0| 20 6f 6d 30 2c 20 64 6f | 75 62 6c 65 20 2a 69 6e | om0, do|uble *in|
|00000bc0| 63 2c 20 64 6f 75 62 6c | 65 20 2a 61 70 2c 20 64 |c, doubl|e *ap, d|
|00000bd0| 6f 75 62 6c 65 20 2a 6f | 6d 29 29 3b 0a 58 65 78 |ouble *o|m));.Xex|
|00000be0| 74 65 72 6e 20 76 6f 69 | 64 20 72 65 66 72 61 63 |tern voi|d refrac|
|00000bf0| 74 20 50 5f 28 28 64 6f | 75 62 6c 65 20 70 72 2c |t P_((do|uble pr,|
|00000c00| 20 64 6f 75 62 6c 65 20 | 74 72 2c 20 64 6f 75 62 | double |tr, doub|
|00000c10| 6c 65 20 74 61 2c 20 64 | 6f 75 62 6c 65 20 2a 61 |le ta, d|ouble *a|
|00000c20| 61 29 29 3b 0a 58 65 78 | 74 65 72 6e 20 76 6f 69 |a));.Xex|tern voi|
|00000c30| 64 20 73 75 6e 70 6f 73 | 20 50 5f 28 28 64 6f 75 |d sunpos| P_((dou|
|00000c40| 62 6c 65 20 4d 6a 64 2c | 20 64 6f 75 62 6c 65 20 |ble Mjd,| double |
|00000c50| 2a 6c 73 6e 2c 20 64 6f | 75 62 6c 65 20 2a 72 73 |*lsn, do|uble *rs|
|00000c60| 6e 29 29 3b 0a 58 65 78 | 74 65 72 6e 20 76 6f 69 |n));.Xex|tern voi|
|00000c70| 64 20 74 61 5f 70 61 72 | 20 50 5f 28 28 64 6f 75 |d ta_par| P_((dou|
|00000c80| 62 6c 65 20 74 68 61 2c | 20 64 6f 75 62 6c 65 20 |ble tha,| double |
|00000c90| 74 64 65 63 2c 20 64 6f | 75 62 6c 65 20 70 68 69 |tdec, do|uble phi|
|00000ca0| 2c 20 64 6f 75 62 6c 65 | 20 68 74 2c 20 64 6f 75 |, double| ht, dou|
|00000cb0| 62 6c 65 20 65 68 70 2c | 20 64 6f 75 62 6c 65 20 |ble ehp,| double |
|00000cc0| 2a 61 68 61 2c 20 64 6f | 75 62 6c 65 20 2a 61 64 |*aha, do|uble *ad|
|00000cd0| 65 63 29 29 3b 0a 58 0a | 58 69 6e 74 20 6f 62 6a |ec));.X.|Xint obj|
|00000ce0| 5f 63 69 72 20 50 5f 28 | 28 4e 6f 77 20 2a 6e 70 |_cir P_(|(Now *np|
|00000cf0| 2c 20 4f 62 6a 20 2a 6f | 70 29 29 3b 0a 58 73 74 |, Obj *o|p));.Xst|
|00000d00| 61 74 69 63 20 6f 62 6a | 5f 70 6c 61 6e 65 74 20 |atic obj|_planet |
|00000d10| 50 5f 28 28 4e 6f 77 20 | 2a 6e 70 2c 20 4f 62 6a |P_((Now |*np, Obj|
|00000d20| 20 2a 6f 70 29 29 3b 0a | 58 73 74 61 74 69 63 20 | *op));.|Xstatic |
|00000d30| 6f 62 6a 5f 66 69 78 65 | 64 20 50 5f 28 28 4e 6f |obj_fixe|d P_((No|
|00000d40| 77 20 2a 6e 70 2c 20 4f | 62 6a 20 2a 6f 70 29 29 |w *np, O|bj *op))|
|00000d50| 3b 0a 58 73 74 61 74 69 | 63 20 6f 62 6a 5f 65 6c |;.Xstati|c obj_el|
|00000d60| 6c 69 70 74 69 63 61 6c | 20 50 5f 28 28 4e 6f 77 |liptical| P_((Now|
|00000d70| 20 2a 6e 70 2c 20 4f 62 | 6a 20 2a 6f 70 29 29 3b | *np, Ob|j *op));|
|00000d80| 0a 58 73 74 61 74 69 63 | 20 6f 62 6a 5f 68 79 70 |.Xstatic| obj_hyp|
|00000d90| 65 72 62 6f 6c 69 63 20 | 50 5f 28 28 4e 6f 77 20 |erbolic |P_((Now |
|00000da0| 2a 6e 70 2c 20 4f 62 6a | 20 2a 6f 70 29 29 3b 0a |*np, Obj| *op));.|
|00000db0| 58 73 74 61 74 69 63 20 | 6f 62 6a 5f 70 61 72 61 |Xstatic |obj_para|
|00000dc0| 62 6f 6c 69 63 20 50 5f | 28 28 4e 6f 77 20 2a 6e |bolic P_|((Now *n|
|00000dd0| 70 2c 20 4f 62 6a 20 2a | 6f 70 29 29 3b 0a 58 73 |p, Obj *|op));.Xs|
|00000de0| 74 61 74 69 63 20 73 75 | 6e 5f 63 69 72 20 50 5f |tatic su|n_cir P_|
|00000df0| 28 28 4e 6f 77 20 2a 6e | 70 2c 20 4f 62 6a 20 2a |((Now *n|p, Obj *|
|00000e00| 6f 70 29 29 3b 0a 58 73 | 74 61 74 69 63 20 6d 6f |op));.Xs|tatic mo|
|00000e10| 6f 6e 5f 63 69 72 20 50 | 5f 28 28 4e 6f 77 20 2a |on_cir P|_((Now *|
|00000e20| 6e 70 2c 20 4f 62 6a 20 | 2a 6f 70 29 29 3b 0a 58 |np, Obj |*op));.X|
|00000e30| 73 74 61 74 69 63 20 76 | 6f 69 64 20 63 69 72 5f |static v|oid cir_|
|00000e40| 73 6b 79 20 50 5f 28 28 | 4e 6f 77 20 2a 6e 70 2c |sky P_((|Now *np,|
|00000e50| 20 64 6f 75 62 6c 65 20 | 6c 70 64 2c 20 64 6f 75 | double |lpd, dou|
|00000e60| 62 6c 65 20 70 73 69 2c | 20 64 6f 75 62 6c 65 20 |ble psi,| double |
|00000e70| 72 70 2c 20 64 6f 75 62 | 6c 65 20 72 68 6f 2c 20 |rp, doub|le rho, |
|00000e80| 64 6f 75 62 6c 65 20 6c | 61 6d 2c 20 64 6f 75 62 |double l|am, doub|
|00000e90| 6c 65 20 62 65 74 2c 20 | 64 6f 75 62 6c 65 20 6c |le bet, |double l|
|00000ea0| 73 6e 2c 20 64 6f 75 62 | 6c 65 20 72 73 6e 2c 20 |sn, doub|le rsn, |
|00000eb0| 4f 62 6a 20 2a 6f 70 29 | 29 3b 0a 58 73 74 61 74 |Obj *op)|);.Xstat|
|00000ec0| 69 63 20 76 6f 69 64 20 | 65 6c 6f 6e 67 61 74 69 |ic void |elongati|
|00000ed0| 6f 6e 20 50 5f 28 28 64 | 6f 75 62 6c 65 20 6c 61 |on P_((d|ouble la|
|00000ee0| 6d 2c 20 64 6f 75 62 6c | 65 20 62 65 74 2c 20 64 |m, doubl|e bet, d|
|00000ef0| 6f 75 62 6c 65 20 6c 73 | 6e 2c 20 64 6f 75 62 6c |ouble ls|n, doubl|
|00000f00| 65 20 2a 65 6c 29 29 3b | 0a 58 0a 58 23 75 6e 64 |e *el));|.X.X#und|
|00000f10| 65 66 20 50 5f 0a 58 0a | 58 2f 2a 20 67 69 76 65 |ef P_.X.|X/* give|
|00000f20| 6e 20 61 20 4e 6f 77 20 | 61 6e 64 20 61 6e 20 4f |n a Now |and an O|
|00000f30| 62 6a 2c 20 66 69 6c 6c | 20 69 6e 20 74 68 65 20 |bj, fill| in the |
|00000f40| 61 70 70 72 6f 70 72 69 | 72 61 74 65 20 73 5f 2a |appropri|rate s_*|
|00000f50| 20 66 69 65 6c 64 73 20 | 77 69 74 68 69 6e 20 4f | fields |within O|
|00000f60| 62 6a 2e 0a 58 20 2a 20 | 72 65 74 75 72 6e 20 30 |bj..X * |return 0|
|00000f70| 20 69 66 20 61 6c 6c 20 | 6f 6b 2c 20 65 6c 73 65 | if all |ok, else|
|00000f80| 20 2d 31 2e 0a 58 20 2a | 2f 0a 58 6f 62 6a 5f 63 | -1..X *|/.Xobj_c|
|00000f90| 69 72 20 28 6e 70 2c 20 | 6f 70 29 0a 58 4e 6f 77 |ir (np, |op).XNow|
|00000fa0| 20 2a 6e 70 3b 0a 58 4f | 62 6a 20 2a 6f 70 3b 0a | *np;.XO|bj *op;.|
|00000fb0| 58 7b 0a 58 09 73 77 69 | 74 63 68 20 28 6f 70 2d |X{.X.swi|tch (op-|
|00000fc0| 3e 74 79 70 65 29 20 7b | 0a 58 09 63 61 73 65 20 |>type) {|.X.case |
|00000fd0| 46 49 58 45 44 3a 09 20 | 72 65 74 75 72 6e 20 28 |FIXED:. |return (|
|00000fe0| 6f 62 6a 5f 66 69 78 65 | 64 20 28 6e 70 2c 20 6f |obj_fixe|d (np, o|
|00000ff0| 70 29 29 3b 0a 58 09 63 | 61 73 65 20 45 4c 4c 49 |p));.X.c|ase ELLI|
|00001000| 50 54 49 43 41 4c 3a 20 | 72 65 74 75 72 6e 20 28 |PTICAL: |return (|
|00001010| 6f 62 6a 5f 65 6c 6c 69 | 70 74 69 63 61 6c 20 28 |obj_elli|ptical (|
|00001020| 6e 70 2c 20 6f 70 29 29 | 3b 0a 58 09 63 61 73 65 |np, op))|;.X.case|
|00001030| 20 48 59 50 45 52 42 4f | 4c 49 43 3a 20 72 65 74 | HYPERBO|LIC: ret|
|00001040| 75 72 6e 20 28 6f 62 6a | 5f 68 79 70 65 72 62 6f |urn (obj|_hyperbo|
|00001050| 6c 69 63 20 28 6e 70 2c | 20 6f 70 29 29 3b 0a 58 |lic (np,| op));.X|
|00001060| 09 63 61 73 65 20 50 41 | 52 41 42 4f 4c 49 43 3a |.case PA|RABOLIC:|
|00001070| 20 20 72 65 74 75 72 6e | 20 28 6f 62 6a 5f 70 61 |  return| (obj_pa|
|00001080| 72 61 62 6f 6c 69 63 20 | 28 6e 70 2c 20 6f 70 29 |rabolic |(np, op)|
|00001090| 29 3b 0a 58 09 63 61 73 | 65 20 50 4c 41 4e 45 54 |);.X.cas|e PLANET|
|000010a0| 3a 20 20 20 20 20 72 65 | 74 75 72 6e 20 28 6f 62 |:     re|turn (ob|
|000010b0| 6a 5f 70 6c 61 6e 65 74 | 20 28 6e 70 2c 20 6f 70 |j_planet| (np, op|
|000010c0| 29 29 3b 0a 58 09 64 65 | 66 61 75 6c 74 3a 0a 58 |));.X.de|fault:.X|
|000010d0| 09 20 20 20 20 70 72 69 | 6e 74 66 20 28 22 6f 62 |.    pri|ntf ("ob|
|000010e0| 6a 5f 63 69 72 28 29 20 | 63 61 6c 6c 65 64 20 77 |j_cir() |called w|
|000010f0| 69 74 68 20 74 79 70 65 | 20 25 64 5c 6e 22 2c 20 |ith type| %d\n", |
|00001100| 6f 70 2d 3e 74 79 70 65 | 29 3b 0a 58 09 20 20 20 |op->type|);.X.   |
|00001110| 20 65 78 69 74 28 31 29 | 3b 0a 58 09 20 20 20 20 | exit(1)|;.X.    |
|00001120| 72 65 74 75 72 6e 20 28 | 2d 31 29 3b 09 2f 2a 20 |return (|-1);./* |
|00001130| 6a 75 73 74 20 66 6f 72 | 20 6c 69 6e 74 20 2a 2f |just for| lint */|
|00001140| 0a 58 09 7d 0a 58 7d 0a | 58 0a 58 73 74 61 74 69 |.X.}.X}.|X.Xstati|
|00001150| 63 0a 58 6f 62 6a 5f 70 | 6c 61 6e 65 74 20 28 6e |c.Xobj_p|lanet (n|
|00001160| 70 2c 20 6f 70 29 0a 58 | 4e 6f 77 20 2a 6e 70 3b |p, op).X|Now *np;|
|00001170| 0a 58 4f 62 6a 20 2a 6f | 70 3b 0a 58 7b 0a 58 09 |.XObj *o|p;.X{.X.|
|00001180| 64 6f 75 62 6c 65 20 6c | 73 6e 2c 20 72 73 6e 3b |double l|sn, rsn;|
|00001190| 09 2f 2a 20 74 72 75 65 | 20 67 65 6f 63 20 6c 6e |./* true| geoc ln|
|000011a0| 67 20 6f 66 20 73 75 6e | 3b 20 64 69 73 74 20 66 |g of sun|; dist f|
|000011b0| 72 6f 6d 20 73 6e 20 74 | 6f 20 65 61 72 74 68 2a |rom sn t|o earth*|
|000011c0| 2f 0a 58 09 64 6f 75 62 | 6c 65 20 6c 70 64 2c 20 |/.X.doub|le lpd, |
|000011d0| 70 73 69 3b 09 2f 2a 20 | 68 65 6c 69 6f 63 65 6e |psi;./* |heliocen|
|000011e0| 74 72 69 63 20 65 63 6c | 69 70 74 69 63 20 6c 6f |tric ecl|iptic lo|
|000011f0| 6e 67 20 61 6e 64 20 6c | 61 74 20 2a 2f 0a 58 09 |ng and l|at */.X.|
|00001200| 64 6f 75 62 6c 65 20 72 | 70 3b 09 09 2f 2a 20 64 |double r|p;../* d|
|00001210| 69 73 74 20 66 72 6f 6d | 20 73 75 6e 20 2a 2f 0a |ist from| sun */.|
|00001220| 58 09 64 6f 75 62 6c 65 | 20 72 68 6f 3b 09 09 2f |X.double| rho;../|
|00001230| 2a 20 64 69 73 74 20 66 | 72 6f 6d 20 65 61 72 74 |* dist f|rom eart|
|00001240| 68 20 2a 2f 0a 58 09 64 | 6f 75 62 6c 65 20 6c 61 |h */.X.d|ouble la|
|00001250| 6d 2c 20 62 65 74 3b 09 | 2f 2a 20 67 65 6f 63 65 |m, bet;.|/* geoce|
|00001260| 6e 74 72 69 63 20 65 63 | 6c 69 70 74 69 63 20 6c |ntric ec|liptic l|
|00001270| 6f 6e 67 20 61 6e 64 20 | 6c 61 74 20 2a 2f 0a 58 |ong and |lat */.X|
|00001280| 09 64 6f 75 62 6c 65 20 | 64 69 61 2c 20 6d 61 67 |.double |dia, mag|
|00001290| 3b 09 2f 2a 20 61 6e 67 | 75 6c 61 72 20 64 69 61 |;./* ang|ular dia|
|000012a0| 6d 65 74 65 72 20 61 74 | 20 31 20 41 55 20 61 6e |meter at| 1 AU an|
|000012b0| 64 20 6d 61 67 6e 69 74 | 75 64 65 20 2a 2f 0a 58 |d magnit|ude */.X|
|000012c0| 09 64 6f 75 62 6c 65 20 | 66 3b 09 09 2f 2a 20 66 |.double |f;../* f|
|000012d0| 72 61 63 74 69 6f 6e 61 | 6c 20 70 68 61 73 65 20 |ractiona|l phase |
|000012e0| 66 72 6f 6d 20 65 61 72 | 74 68 20 2a 2f 0a 58 09 |from ear|th */.X.|
|000012f0| 69 6e 74 20 70 3b 0a 58 | 0a 58 09 2f 2a 20 76 61 |int p;.X|.X./* va|
|00001300| 6c 69 64 61 74 65 20 63 | 6f 64 65 20 61 6e 64 20 |lidate c|ode and |
|00001310| 63 68 65 63 6b 20 66 6f | 72 20 61 20 66 65 77 20 |check fo|r a few |
|00001320| 73 70 65 63 69 61 6c 20 | 63 61 73 65 73 20 2a 2f |special |cases */|
|00001330| 0a 58 09 70 20 3d 20 6f | 70 2d 3e 70 6c 2e 63 6f |.X.p = o|p->pl.co|
|00001340| 64 65 3b 0a 58 09 69 66 | 20 28 70 20 3c 20 30 20 |de;.X.if| (p < 0 |
|00001350| 7c 7c 20 70 20 3e 20 4d | 4f 4f 4e 29 20 7b 0a 58 ||| p > M|OON) {.X|
|00001360| 09 20 20 20 20 70 72 69 | 6e 74 66 20 28 22 75 6e |.    pri|ntf ("un|
|00001370| 6b 6e 6f 77 6e 20 70 6c | 61 6e 65 74 20 63 6f 64 |known pl|anet cod|
|00001380| 65 3a 20 25 64 20 5c 6e | 22 2c 20 70 29 3b 0a 58 |e: %d \n|", p);.X|
|00001390| 09 20 20 20 20 65 78 69 | 74 28 31 29 3b 0a 58 09 |.    exi|t(1);.X.|
|000013a0| 7d 0a 58 09 65 6c 73 65 | 20 69 66 20 28 70 20 3d |}.X.else| if (p =|
|000013b0| 3d 20 53 55 4e 29 0a 58 | 09 20 20 20 20 72 65 74 |= SUN).X|.    ret|
|000013c0| 75 72 6e 20 28 73 75 6e | 5f 63 69 72 20 28 6e 70 |urn (sun|_cir (np|
|000013d0| 2c 20 6f 70 29 29 3b 0a | 58 09 65 6c 73 65 20 69 |, op));.|X.else i|
|000013e0| 66 20 28 70 20 3d 3d 20 | 4d 4f 4f 4e 29 0a 58 09 |f (p == |MOON).X.|
|000013f0| 20 20 20 20 72 65 74 75 | 72 6e 20 28 6d 6f 6f 6e |    retu|rn (moon|
|00001400| 5f 63 69 72 20 28 6e 70 | 2c 20 6f 70 29 29 3b 0a |_cir (np|, op));.|
|00001410| 58 0a 58 0a 58 09 2f 2a | 20 66 69 6e 64 20 73 6f |X.X.X./*| find so|
|00001420| 6c 61 72 20 65 63 6c 69 | 70 74 69 63 61 6c 20 6c |lar ecli|ptical l|
|00001430| 6f 6e 67 69 74 75 64 65 | 20 61 6e 64 20 64 69 73 |ongitude| and dis|
|00001440| 74 61 6e 63 65 20 74 6f | 20 73 75 6e 20 66 72 6f |tance to| sun fro|
|00001450| 6d 20 65 61 72 74 68 20 | 2a 2f 0a 58 09 73 75 6e |m earth |*/.X.sun|
|00001460| 70 6f 73 20 28 6d 6a 64 | 2c 20 26 6c 73 6e 2c 20 |pos (mjd|, &lsn, |
|00001470| 26 72 73 6e 29 3b 0a 58 | 0a 58 09 2f 2a 20 66 69 |&rsn);.X|.X./* fi|
|00001480| 6e 64 20 68 65 6c 69 6f | 20 6c 6f 6e 67 2f 6c 61 |nd helio| long/la|
|00001490| 74 3b 20 73 75 6e 2f 70 | 6c 61 6e 65 74 20 61 6e |t; sun/p|lanet an|
|000014a0| 64 20 65 61 72 74 68 2f | 70 6c 61 6e 74 20 64 69 |d earth/|plant di|
|000014b0| 73 74 3b 20 65 63 6c 69 | 70 74 69 63 0a 58 09 20 |st; ecli|ptic.X. |
|000014c0| 2a 20 6c 6f 6e 67 2f 6c | 61 74 3b 20 64 69 61 6d |* long/l|at; diam|
|000014d0| 65 74 65 72 20 61 6e 64 | 20 6d 61 67 2e 0a 58 09 |eter and| mag..X.|
|000014e0| 20 2a 2f 0a 58 09 70 6c | 61 6e 73 28 6d 6a 64 2c | */.X.pl|ans(mjd,|
|000014f0| 20 70 2c 20 26 6c 70 64 | 2c 20 26 70 73 69 2c 20 | p, &lpd|, &psi, |
|00001500| 26 72 70 2c 20 26 72 68 | 6f 2c 20 26 6c 61 6d 2c |&rp, &rh|o, &lam,|
|00001510| 20 26 62 65 74 2c 20 26 | 64 69 61 2c 20 26 6d 61 | &bet, &|dia, &ma|
|00001520| 67 29 3b 0a 58 0a 58 09 | 2f 2a 20 66 69 6c 6c 20 |g);.X.X.|/* fill |
|00001530| 69 6e 20 61 6c 6c 20 6f | 66 20 6f 70 2d 3e 73 5f |in all o|f op->s_|
|00001540| 2a 20 73 74 75 66 66 20 | 65 78 63 65 70 74 20 73 |* stuff |except s|
|00001550| 5f 73 69 7a 65 20 61 6e | 64 20 73 5f 6d 61 67 20 |_size an|d s_mag |
|00001560| 2a 2f 0a 58 09 63 69 72 | 5f 73 6b 79 20 28 6e 70 |*/.X.cir|_sky (np|
|00001570| 2c 20 6c 70 64 2c 20 70 | 73 69 2c 20 72 70 2c 20 |, lpd, p|si, rp, |
|00001580| 72 68 6f 2c 20 6c 61 6d | 2c 20 62 65 74 2c 20 6c |rho, lam|, bet, l|
|00001590| 73 6e 2c 20 72 73 6e 2c | 20 6f 70 29 3b 0a 58 0a |sn, rsn,| op);.X.|
|000015a0| 58 09 2f 2a 20 63 6f 6d | 70 75 74 65 20 6d 61 67 |X./* com|pute mag|
|000015b0| 6e 69 74 75 64 65 20 61 | 6e 64 20 61 6e 67 75 6c |nitude a|nd angul|
|000015c0| 61 72 20 73 69 7a 65 20 | 2a 2f 0a 58 09 66 20 3d |ar size |*/.X.f =|
|000015d0| 20 30 2e 32 35 20 2a 20 | 28 28 28 72 70 2b 72 68 | 0.25 * |(((rp+rh|
|000015e0| 6f 29 2a 28 72 70 2b 72 | 68 6f 29 20 2d 20 72 73 |o)*(rp+r|ho) - rs|
|000015f0| 6e 2a 72 73 6e 29 2f 28 | 72 70 2a 72 68 6f 29 29 |n*rsn)/(|rp*rho))|
|00001600| 3b 0a 58 09 6f 70 2d 3e | 73 5f 6d 61 67 20 3d 20 |;.X.op->|s_mag = |
|00001610| 4d 41 47 53 43 41 4c 45 | 20 2a 20 28 6d 61 67 20 |MAGSCALE| * (mag |
|00001620| 2b 20 35 2e 30 2a 6c 6f | 67 28 72 70 2a 72 68 6f |+ 5.0*lo|g(rp*rho|
|00001630| 2f 73 71 72 74 28 66 29 | 29 2f 6c 6f 67 28 31 30 |/sqrt(f)|)/log(10|
|00001640| 2e 30 29 29 3b 0a 58 09 | 6f 70 2d 3e 73 5f 73 69 |.0));.X.|op->s_si|
|00001650| 7a 65 20 3d 20 64 69 61 | 2f 72 68 6f 20 2b 20 30 |ze = dia|/rho + 0|
|00001660| 2e 35 3b 0a 58 0a 58 09 | 72 65 74 75 72 6e 20 28 |.5;.X.X.|return (|
|00001670| 30 29 3b 0a 58 7d 0a 58 | 0a 58 73 74 61 74 69 63 |0);.X}.X|.Xstatic|
|00001680| 0a 58 6f 62 6a 5f 66 69 | 78 65 64 20 28 6e 70 2c |.Xobj_fi|xed (np,|
|00001690| 20 6f 70 29 0a 58 4e 6f | 77 20 2a 6e 70 3b 0a 58 | op).XNo|w *np;.X|
|000016a0| 4f 62 6a 20 2a 6f 70 3b | 0a 58 7b 0a 58 09 64 6f |Obj *op;|.X{.X.do|
|000016b0| 75 62 6c 65 20 6c 73 6e | 2c 20 72 73 6e 3b 09 2f |uble lsn|, rsn;./|
|000016c0| 2a 20 74 72 75 65 20 67 | 65 6f 63 20 6c 6e 67 20 |* true g|eoc lng |
|000016d0| 6f 66 20 73 75 6e 2c 20 | 64 69 73 74 20 66 72 6f |of sun, |dist fro|
|000016e0| 6d 20 73 6e 20 74 6f 20 | 65 61 72 74 68 2a 2f 0a |m sn to |earth*/.|
|000016f0| 58 09 64 6f 75 62 6c 65 | 20 6c 61 6d 2c 20 62 65 |X.double| lam, be|
|00001700| 74 3b 09 2f 2a 20 67 65 | 6f 63 65 6e 74 72 69 63 |t;./* ge|ocentric|
|00001710| 20 65 63 6c 69 70 74 69 | 63 20 6c 6f 6e 67 20 61 | eclipti|c long a|
|00001720| 6e 64 20 6c 61 74 20 2a | 2f 0a 58 09 64 6f 75 62 |nd lat *|/.X.doub|
|00001730| 6c 65 20 6c 73 74 3b 09 | 09 2f 2a 20 6c 6f 63 61 |le lst;.|./* loca|
|00001740| 6c 20 73 69 64 65 72 65 | 61 6c 20 74 69 6d 65 20 |l sidere|al time |
|00001750| 2a 2f 0a 58 09 64 6f 75 | 62 6c 65 20 68 61 3b 09 |*/.X.dou|ble ha;.|
|00001760| 09 2f 2a 20 6c 6f 63 61 | 6c 20 68 6f 75 72 20 61 |./* loca|l hour a|
|00001770| 6e 67 6c 65 20 2a 2f 0a | 58 09 64 6f 75 62 6c 65 |ngle */.|X.double|
|00001780| 20 65 6c 3b 09 09 2f 2a | 20 65 6c 6f 6e 67 61 74 | el;../*| elongat|
|00001790| 69 6f 6e 20 2a 2f 0a 58 | 09 64 6f 75 62 6c 65 20 |ion */.X|.double |
|000017a0| 61 6c 74 2c 20 61 7a 3b | 09 09 2f 2a 20 63 75 72 |alt, az;|../* cur|
|000017b0| 72 65 6e 74 20 61 6c 74 | 2c 20 61 7a 20 2a 2f 0a |rent alt|, az */.|
|000017c0| 58 09 64 6f 75 62 6c 65 | 20 72 61 2c 20 64 65 63 |X.double| ra, dec|
|000017d0| 3b 09 09 2f 2a 20 72 61 | 20 61 6e 64 20 64 65 63 |;../* ra| and dec|
|000017e0| 20 61 74 20 65 70 6f 63 | 68 20 6f 66 20 64 61 74 | at epoc|h of dat|
|000017f0| 65 20 2a 2f 0a 58 0a 58 | 09 2f 2a 20 73 65 74 20 |e */.X.X|./* set |
|00001800| 72 61 2f 64 65 63 20 74 | 6f 20 74 68 65 69 72 20 |ra/dec t|o their |
|00001810| 76 61 6c 75 65 73 20 61 | 74 20 65 70 6f 63 68 20 |values a|t epoch |
|00001820| 6f 66 20 64 61 74 65 20 | 2a 2f 0a 58 09 72 61 20 |of date |*/.X.ra |
|00001830| 3d 20 6f 70 2d 3e 66 5f | 52 41 3b 0a 58 09 64 65 |= op->f_|RA;.X.de|
|00001840| 63 20 3d 20 6f 70 2d 3e | 66 5f 64 65 63 3b 0a 58 |c = op->|f_dec;.X|
|00001850| 09 70 72 65 63 65 73 73 | 20 28 6f 70 2d 3e 66 5f |.precess| (op->f_|
|00001860| 65 70 6f 63 68 2c 20 6d | 6a 64 2c 20 26 72 61 2c |epoch, m|jd, &ra,|
|00001870| 20 26 64 65 63 29 3b 0a | 58 0a 58 09 2f 2a 20 73 | &dec);.|X.X./* s|
|00001880| 65 74 20 73 5f 72 61 2f | 64 65 63 20 61 74 20 64 |et s_ra/|dec at d|
|00001890| 65 73 69 72 65 64 20 65 | 70 6f 63 68 2e 20 20 2a |esired e|poch.  *|
|000018a0| 2f 0a 58 09 69 66 20 28 | 65 70 6f 63 68 20 3d 3d |/.X.if (|epoch ==|
|000018b0| 20 45 4f 44 29 20 7b 0a | 58 09 20 20 20 20 6f 70 | EOD) {.|X.    op|
|000018c0| 2d 3e 73 5f 72 61 20 3d | 20 72 61 3b 0a 58 09 20 |->s_ra =| ra;.X. |
|000018d0| 20 20 20 6f 70 2d 3e 73 | 5f 64 65 63 20 3d 20 64 |   op->s|_dec = d|
|000018e0| 65 63 3b 0a 58 09 7d 20 | 65 6c 73 65 20 7b 0a 58 |ec;.X.} |else {.X|
|000018f0| 09 20 20 20 20 2f 2a 20 | 77 61 6e 74 20 61 20 63 |.    /* |want a c|
|00001900| 65 72 74 61 69 6e 20 65 | 70 6f 63 68 20 2d 2d 20 |ertain e|poch -- |
|00001910| 69 66 20 69 74 27 73 20 | 6e 6f 74 20 77 68 61 74 |if it's |not what|
|00001920| 20 74 68 65 20 64 61 74 | 61 62 61 73 65 20 69 73 | the dat|abase is|
|00001930| 20 61 74 0a 58 09 20 20 | 20 20 20 2a 20 77 65 20 | at.X.  |   * we |
|00001940| 63 68 61 6e 67 65 20 74 | 68 65 20 6f 72 69 67 69 |change t|he origi|
|00001950| 6e 61 6c 20 74 6f 20 73 | 61 76 65 20 74 69 6d 65 |nal to s|ave time|
|00001960| 20 6e 65 78 74 20 74 69 | 6d 65 20 61 73 73 75 6d | next ti|me assum|
|00001970| 69 6e 67 20 74 68 65 0a | 58 09 20 20 20 20 20 2a |ing the.|X.     *|
|00001980| 20 75 73 65 72 20 69 73 | 20 6c 69 6b 65 6c 79 20 | user is| likely |
|00001990| 74 6f 20 73 74 69 63 6b | 20 77 69 74 68 20 74 68 |to stick| with th|
|000019a0| 69 73 20 66 6f 72 20 61 | 20 77 68 69 6c 65 2e 0a |is for a| while..|
|000019b0| 58 09 20 20 20 20 20 2a | 2f 0a 58 09 20 20 20 20 |X.     *|/.X.    |
|000019c0| 69 66 20 28 28 66 6c 6f | 61 74 29 65 70 6f 63 68 |if ((flo|at)epoch|
|000019d0| 20 21 3d 20 6f 70 2d 3e | 66 5f 65 70 6f 63 68 29 | != op->|f_epoch)|
|000019e0| 20 7b 0a 58 09 09 64 6f | 75 62 6c 65 20 74 72 61 | {.X..do|uble tra|
|000019f0| 20 3d 20 6f 70 2d 3e 66 | 5f 52 41 2c 20 74 64 65 | = op->f|_RA, tde|
|00001a00| 63 20 3d 20 6f 70 2d 3e | 66 5f 64 65 63 3b 0a 58 |c = op->|f_dec;.X|
|00001a10| 09 09 66 6c 6f 61 74 20 | 74 65 70 6f 63 68 20 3d |..float |tepoch =|
|00001a20| 20 65 70 6f 63 68 3b 09 | 2f 2a 20 63 6f 6d 70 61 | epoch;.|/* compa|
|00001a30| 72 65 20 74 6f 20 66 6c | 6f 61 74 20 70 72 65 63 |re to fl|oat prec|
|00001a40| 69 73 69 6f 6e 20 2a 2f | 0a 58 09 09 70 72 65 63 |ision */|.X..prec|
|00001a50| 65 73 73 20 28 6f 70 2d | 3e 66 5f 65 70 6f 63 68 |ess (op-|>f_epoch|
|00001a60| 2c 20 74 65 70 6f 63 68 | 2c 20 26 74 72 61 2c 20 |, tepoch|, &tra, |
|00001a70| 26 74 64 65 63 29 3b 0a | 58 09 09 6f 70 2d 3e 66 |&tdec);.|X..op->f|
|00001a80| 5f 65 70 6f 63 68 20 3d | 20 74 65 70 6f 63 68 3b |_epoch =| tepoch;|
|00001a90| 0a 58 09 09 6f 70 2d 3e | 66 5f 52 41 20 3d 20 74 |.X..op->|f_RA = t|
|00001aa0| 72 61 3b 0a 58 09 09 6f | 70 2d 3e 66 5f 64 65 63 |ra;.X..o|p->f_dec|
|00001ab0| 20 3d 20 74 64 65 63 3b | 0a 58 09 20 20 20 20 7d | = tdec;|.X.    }|
|00001ac0| 0a 58 09 20 20 20 20 6f | 70 2d 3e 73 5f 72 61 20 |.X.    o|p->s_ra |
|00001ad0| 3d 20 6f 70 2d 3e 66 5f | 52 41 3b 0a 58 09 20 20 |= op->f_|RA;.X.  |
|00001ae0| 20 20 6f 70 2d 3e 73 5f | 64 65 63 20 3d 20 6f 70 |  op->s_|dec = op|
|00001af0| 2d 3e 66 5f 64 65 63 3b | 0a 58 09 7d 0a 58 0a 58 |->f_dec;|.X.}.X.X|
|00001b00| 09 2f 2a 20 63 6f 6e 76 | 65 72 74 20 65 71 75 69 |./* conv|ert equi|
|00001b10| 74 6f 72 69 61 6c 20 72 | 61 2f 64 65 63 20 74 6f |torial r|a/dec to|
|00001b20| 20 67 65 6f 63 65 6e 74 | 72 69 63 20 65 63 6c 69 | geocent|ric ecli|
|00001b30| 70 74 69 63 20 6c 61 74 | 2f 6c 6f 6e 67 20 2a 2f |ptic lat|/long */|
|00001b40| 0a 58 09 65 71 5f 65 63 | 6c 20 28 6d 6a 64 2c 20 |.X.eq_ec|l (mjd, |
|00001b50| 72 61 2c 20 64 65 63 2c | 20 26 62 65 74 2c 20 26 |ra, dec,| &bet, &|
|00001b60| 6c 61 6d 29 3b 0a 58 0a | 58 09 2f 2a 20 66 69 6e |lam);.X.|X./* fin|
|00001b70| 64 20 73 6f 6c 61 72 20 | 65 63 6c 69 70 74 69 63 |d solar |ecliptic|
|00001b80| 61 6c 20 6c 6f 6e 67 69 | 74 75 64 65 20 61 6e 64 |al longi|tude and|
|00001b90| 20 64 69 73 74 61 6e 63 | 65 20 74 6f 20 73 75 6e | distanc|e to sun|
|00001ba0| 20 66 72 6f 6d 20 65 61 | 72 74 68 20 2a 2f 0a 58 | from ea|rth */.X|
|00001bb0| 09 73 75 6e 70 6f 73 20 | 28 6d 6a 64 2c 20 26 6c |.sunpos |(mjd, &l|
|00001bc0| 73 6e 2c 20 26 72 73 6e | 29 3b 0a 58 0a 58 09 2f |sn, &rsn|);.X.X./|
|00001bd0| 2a 20 63 6f 6d 70 75 74 | 65 20 65 6c 6f 6e 67 61 |* comput|e elonga|
|00001be0| 74 69 6f 6e 20 66 72 6f | 6d 20 65 63 6c 69 70 74 |tion fro|m eclipt|
|00001bf0| 69 63 20 6c 6f 6e 67 2f | 6c 61 74 20 61 6e 64 20 |ic long/|lat and |
|00001c00| 73 75 6e 20 67 65 6f 63 | 65 6e 74 72 69 63 20 6c |sun geoc|entric l|
|00001c10| 6f 6e 67 20 2a 2f 0a 58 | 09 65 6c 6f 6e 67 61 74 |ong */.X|.elongat|
|00001c20| 69 6f 6e 20 28 6c 61 6d | 2c 20 62 65 74 2c 20 6c |ion (lam|, bet, l|
|00001c30| 73 6e 2c 20 26 65 6c 29 | 3b 0a 58 09 65 6c 20 3d |sn, &el)|;.X.el =|
|00001c40| 20 72 61 64 64 65 67 28 | 65 6c 29 3b 0a 58 09 6f | raddeg(|el);.X.o|
|00001c50| 70 2d 3e 73 5f 65 6c 6f | 6e 67 20 3d 20 65 6c 3b |p->s_elo|ng = el;|
|00001c60| 0a 58 0a 58 09 2f 2a 20 | 74 68 65 73 65 20 61 72 |.X.X./* |these ar|
|00001c70| 65 20 72 65 61 6c 6c 79 | 20 74 68 65 20 73 61 6d |e really| the sam|
|00001c80| 65 20 66 69 65 6c 64 73 | 20 2e 2e 2e 0a 58 09 6f |e fields| ....X.o|
|00001c90| 70 2d 3e 73 5f 6d 61 67 | 20 3d 20 6f 70 2d 3e 66 |p->s_mag| = op->f|
|00001ca0| 5f 6d 61 67 3b 0a 58 09 | 6f 70 2d 3e 73 5f 73 69 |_mag;.X.|op->s_si|
|00001cb0| 7a 65 20 3d 20 6f 70 2d | 3e 66 5f 73 69 7a 65 3b |ze = op-|>f_size;|
|00001cc0| 0a 58 09 2a 2f 0a 58 0a | 58 09 2f 2a 20 61 6c 74 |.X.*/.X.|X./* alt|
|00001cd0| 2c 20 61 7a 3a 20 63 6f | 72 72 65 63 74 20 66 6f |, az: co|rrect fo|
|00001ce0| 72 20 72 65 66 72 61 63 | 74 69 6f 6e 3b 20 75 73 |r refrac|tion; us|
|00001cf0| 65 20 65 6f 64 20 72 61 | 2f 64 65 63 2e 20 2a 2f |e eod ra|/dec. */|
|00001d00| 0a 58 09 6e 6f 77 5f 6c | 73 74 20 28 6e 70 2c 20 |.X.now_l|st (np, |
|00001d10| 26 6c 73 74 29 3b 0a 58 | 09 68 61 20 3d 20 68 72 |&lst);.X|.ha = hr|
|00001d20| 72 61 64 28 6c 73 74 29 | 20 2d 20 72 61 3b 0a 58 |rad(lst)| - ra;.X|
|00001d30| 09 68 61 64 65 63 5f 61 | 61 20 28 6c 61 74 2c 20 |.hadec_a|a (lat, |
|00001d40| 68 61 2c 20 64 65 63 2c | 20 26 61 6c 74 2c 20 26 |ha, dec,| &alt, &|
|00001d50| 61 7a 29 3b 0a 58 09 72 | 65 66 72 61 63 74 20 28 |az);.X.r|efract (|
|00001d60| 70 72 65 73 73 75 72 65 | 2c 20 74 65 6d 70 2c 20 |pressure|, temp, |
|00001d70| 61 6c 74 2c 20 26 61 6c | 74 29 3b 0a 58 09 6f 70 |alt, &al|t);.X.op|
|00001d80| 2d 3e 73 5f 61 6c 74 20 | 3d 20 61 6c 74 3b 0a 58 |->s_alt |= alt;.X|
|00001d90| 09 6f 70 2d 3e 73 5f 61 | 7a 20 3d 20 61 7a 3b 0a |.op->s_a|z = az;.|
|00001da0| 58 0a 58 09 72 65 74 75 | 72 6e 20 28 30 29 3b 0a |X.X.retu|rn (0);.|
|00001db0| 58 7d 0a 58 0a 58 2f 2a | 20 63 6f 6d 70 75 74 65 |X}.X.X/*| compute|
|00001dc0| 20 73 6b 79 20 63 69 72 | 63 75 6d 73 74 61 6e 63 | sky cir|cumstanc|
|00001dd0| 65 73 20 6f 66 20 61 6e | 20 6f 62 6a 65 63 74 20 |es of an| object |
|00001de0| 69 6e 20 68 65 6c 69 6f | 63 65 6e 74 72 69 63 20 |in helio|centric |
|00001df0| 65 6c 6c 69 70 74 69 63 | 20 6f 72 62 69 74 20 61 |elliptic| orbit a|
|00001e00| 74 20 2a 6e 70 2e 0a 58 | 20 2a 2f 0a 58 73 74 61 |t *np..X| */.Xsta|
|00001e10| 74 69 63 0a 58 6f 62 6a | 5f 65 6c 6c 69 70 74 69 |tic.Xobj|_ellipti|
|00001e20| 63 61 6c 20 28 6e 70 2c | 20 6f 70 29 0a 58 4e 6f |cal (np,| op).XNo|
|00001e30| 77 20 2a 6e 70 3b 0a 58 | 4f 62 6a 20 2a 6f 70 3b |w *np;.X|Obj *op;|
|00001e40| 0a 58 7b 0a 58 09 64 6f | 75 62 6c 65 20 6c 73 6e |.X{.X.do|uble lsn|
|00001e50| 2c 20 72 73 6e 3b 09 2f | 2a 20 74 72 75 65 20 67 |, rsn;./|* true g|
|00001e60| 65 6f 63 20 6c 6e 67 20 | 6f 66 20 73 75 6e 3b 20 |eoc lng |of sun; |
|00001e70| 64 69 73 74 20 66 72 6f | 6d 20 73 6e 20 74 6f 20 |dist fro|m sn to |
|00001e80| 65 61 72 74 68 2a 2f 0a | 58 09 64 6f 75 62 6c 65 |earth*/.|X.double|
|00001e90| 20 64 74 3b 09 09 2f 2a | 20 6c 69 67 68 74 20 74 | dt;../*| light t|
|00001ea0| 72 61 76 65 6c 20 74 69 | 6d 65 20 74 6f 20 6f 62 |ravel ti|me to ob|
|00001eb0| 6a 65 63 74 20 2a 2f 0a | 58 09 64 6f 75 62 6c 65 |ject */.|X.double|
|00001ec0| 20 6c 67 3b 09 09 2f 2a | 20 68 65 6c 69 6f 20 6c | lg;../*| helio l|
|00001ed0| 6f 6e 67 20 6f 66 20 65 | 61 72 74 68 20 2a 2f 0a |ong of e|arth */.|
|00001ee0| 58 09 64 6f 75 62 6c 65 | 20 6e 75 2c 20 65 61 3b |X.double| nu, ea;|
|00001ef0| 09 09 2f 2a 20 74 72 75 | 65 20 61 6e 6f 6d 61 6c |../* tru|e anomal|
|00001f00| 79 20 61 6e 64 20 65 63 | 63 65 6e 74 72 69 63 20 |y and ec|centric |
|00001f10| 61 6e 6f 6d 61 6c 79 20 | 2a 2f 0a 58 09 64 6f 75 |anomaly |*/.X.dou|
|00001f20| 62 6c 65 20 6d 61 3b 09 | 09 2f 2a 20 6d 65 61 6e |ble ma;.|./* mean|
|00001f30| 20 61 6e 6f 6d 61 6c 79 | 20 2a 2f 0a 58 09 64 6f | anomaly| */.X.do|
|00001f40| 75 62 6c 65 20 72 70 3b | 09 09 2f 2a 20 64 69 73 |uble rp;|../* dis|
|00001f50| 74 61 6e 63 65 20 66 72 | 6f 6d 20 74 68 65 20 73 |tance fr|om the s|
|00001f60| 75 6e 20 2a 2f 0a 58 09 | 64 6f 75 62 6c 65 20 6c |un */.X.|double l|
|00001f70| 6f 2c 20 73 6c 6f 2c 20 | 63 6c 6f 3b 09 2f 2a 20 |o, slo, |clo;./* |
|00001f80| 61 6e 67 6c 65 20 66 72 | 6f 6d 20 61 73 63 65 6e |angle fr|om ascen|
|00001f90| 64 69 6e 67 20 6e 6f 64 | 65 20 2a 2f 0a 58 09 64 |ding nod|e */.X.d|
|00001fa0| 6f 75 62 6c 65 20 69 6e | 63 3b 09 09 2f 2a 20 69 |ouble in|c;../* i|
|00001fb0| 6e 63 6c 69 6e 61 74 69 | 6f 6e 20 2a 2f 0a 58 09 |nclinati|on */.X.|
|00001fc0| 64 6f 75 62 6c 65 20 70 | 73 69 2c 20 73 70 73 69 |double p|si, spsi|
|00001fd0| 2c 20 63 70 73 69 3b 09 | 2f 2a 20 68 65 6c 69 6f |, cpsi;.|/* helio|
|00001fe0| 63 65 6e 74 72 69 63 20 | 6c 61 74 69 74 75 64 65 |centric |latitude|
|00001ff0| 20 2a 2f 0a 58 09 64 6f | 75 62 6c 65 20 6c 70 64 | */.X.do|uble lpd|
|00002000| 3b 20 09 09 2f 2a 20 68 | 65 6c 69 6f 63 65 6e 74 |; ../* h|eliocent|
|00002010| 72 69 63 20 6c 6f 6e 67 | 69 74 75 64 65 20 2a 2f |ric long|itude */|
|00002020| 0a 58 09 64 6f 75 62 6c | 65 20 72 68 6f 3b 09 09 |.X.doubl|e rho;..|
|00002030| 2f 2a 20 64 69 73 74 61 | 6e 63 65 20 66 72 6f 6d |/* dista|nce from|
|00002040| 20 74 68 65 20 45 61 72 | 74 68 20 2a 2f 0a 58 09 | the Ear|th */.X.|
|00002050| 64 6f 75 62 6c 65 20 6f | 6d 3b 09 09 2f 2a 20 61 |double o|m;../* a|
|00002060| 72 67 20 6f 66 20 70 65 | 72 69 68 65 6c 69 6f 6e |rg of pe|rihelion|
|00002070| 20 2a 2f 0a 58 09 64 6f | 75 62 6c 65 20 4f 6d 3b | */.X.do|uble Om;|
|00002080| 09 09 2f 2a 20 6c 6f 6e | 67 20 6f 66 20 61 73 63 |../* lon|g of asc|
|00002090| 65 6e 64 69 6e 67 20 6e | 6f 64 65 2e 20 2a 2f 0a |ending n|ode. */.|
|000020a0| 58 09 64 6f 75 62 6c 65 | 20 6c 61 6d 3b 20 20 20 |X.double| lam;   |
|000020b0| 20 09 09 2f 2a 20 67 65 | 6f 63 65 6e 74 72 69 63 | ../* ge|ocentric|
|000020c0| 20 65 63 6c 69 70 74 69 | 63 20 6c 6f 6e 67 69 74 | eclipti|c longit|
|000020d0| 75 64 65 20 2a 2f 0a 58 | 09 64 6f 75 62 6c 65 20 |ude */.X|.double |
|000020e0| 62 65 74 3b 20 20 20 20 | 09 09 2f 2a 20 67 65 6f |bet;    |../* geo|
|000020f0| 63 65 6e 74 72 69 63 20 | 65 63 6c 69 70 74 69 63 |centric |ecliptic|
|00002100| 20 6c 61 74 69 74 75 64 | 65 20 2a 2f 0a 58 09 64 | latitud|e */.X.d|
|00002110| 6f 75 62 6c 65 20 65 3b | 09 09 2f 2a 20 66 61 73 |ouble e;|../* fas|
|00002120| 74 20 65 63 63 65 6e 74 | 72 69 63 69 74 79 20 2a |t eccent|ricity *|
|00002130| 2f 0a 58 09 64 6f 75 62 | 6c 65 20 6c 6c 2c 20 73 |/.X.doub|le ll, s|
|00002140| 6c 6c 2c 20 63 6c 6c 3b | 09 2f 2a 20 68 65 6c 69 |ll, cll;|./* heli|
|00002150| 6f 20 61 6e 67 6c 65 20 | 62 65 74 77 65 65 6e 20 |o angle |between |
|00002160| 6f 62 6a 65 63 74 20 61 | 6e 64 20 65 61 72 74 68 |object a|nd earth|
|00002170| 20 2a 2f 0a 58 09 64 6f | 75 62 6c 65 20 6d 61 67 | */.X.do|uble mag|
|00002180| 3b 09 09 2f 2a 20 6d 61 | 67 6e 69 74 75 64 65 20 |;../* ma|gnitude |
|00002190| 2a 2f 0a 58 09 64 6f 75 | 62 6c 65 20 72 70 64 3b |*/.X.dou|ble rpd;|
|000021a0| 0a 58 09 64 6f 75 62 6c | 65 20 79 3b 0a 58 09 69 |.X.doubl|e y;.X.i|
|000021b0| 6e 74 20 70 61 73 73 3b | 0a 58 0a 58 09 2f 2a 20 |nt pass;|.X.X./* |
|000021c0| 66 69 6e 64 20 6c 6f 63 | 61 74 69 6f 6e 20 6f 66 |find loc|ation of|
|000021d0| 20 65 61 72 74 68 20 66 | 72 6f 6d 20 73 75 6e 20 | earth f|rom sun |
|000021e0| 6e 6f 77 20 2a 2f 0a 58 | 09 73 75 6e 70 6f 73 20 |now */.X|.sunpos |
|000021f0| 28 6d 6a 64 2c 20 26 6c | 73 6e 2c 20 26 72 73 6e |(mjd, &l|sn, &rsn|
|00002200| 29 3b 0a 58 09 6c 67 20 | 3d 20 6c 73 6e 20 2b 20 |);.X.lg |= lsn + |
|00002210| 50 49 3b 0a 58 0a 58 09 | 2f 2a 20 66 61 73 74 65 |PI;.X.X.|/* faste|
|00002220| 72 20 61 63 63 65 73 73 | 20 74 6f 20 65 63 63 65 |r access| to ecce|
|00002230| 6e 74 72 69 63 74 79 20 | 2a 2f 0a 58 09 65 20 3d |ntricty |*/.X.e =|
|00002240| 20 6f 70 2d 3e 65 5f 65 | 3b 0a 58 0a 58 09 2f 2a | op->e_e|;.X.X./*|
|00002250| 20 6d 65 61 6e 20 64 61 | 69 6c 79 20 6d 6f 74 69 | mean da|ily moti|
|00002260| 6f 6e 20 69 73 20 6f 70 | 74 69 6f 6e 61 6c 20 2d |on is op|tional -|
|00002270| 2d 20 66 69 6c 6c 20 69 | 6e 20 66 72 6f 6d 20 70 |- fill i|n from p|
|00002280| 65 72 69 6f 64 20 69 66 | 20 30 20 2a 2f 0a 58 09 |eriod if| 0 */.X.|
|00002290| 69 66 20 28 6f 70 2d 3e | 65 5f 6e 20 3d 3d 20 30 |if (op->|e_n == 0|
|000022a0| 2e 30 29 20 7b 0a 58 09 | 20 20 20 20 64 6f 75 62 |.0) {.X.|    doub|
|000022b0| 6c 65 20 61 20 3d 20 6f | 70 2d 3e 65 5f 61 3b 0a |le a = o|p->e_a;.|
|000022c0| 58 09 20 20 20 20 6f 70 | 2d 3e 65 5f 6e 20 3d 20 |X.    op|->e_n = |
|000022d0| 30 2e 39 38 35 36 30 37 | 36 36 38 36 2f 73 71 72 |0.985607|6686/sqr|
|000022e0| 74 28 61 2a 61 2a 61 29 | 3b 0a 58 09 7d 0a 58 0a |t(a*a*a)|;.X.}.X.|
|000022f0| 58 09 2f 2a 20 63 6f 72 | 72 65 63 74 20 66 6f 72 |X./* cor|rect for|
|00002300| 20 6c 69 67 68 74 20 74 | 69 6d 65 20 62 79 20 63 | light t|ime by c|
|00002310| 6f 6d 70 75 74 69 6e 67 | 20 70 6f 73 69 74 69 6f |omputing| positio|
|00002320| 6e 20 61 74 20 74 69 6d | 65 20 6d 6a 64 2c 20 74 |n at tim|e mjd, t|
|00002330| 68 65 6e 0a 58 09 20 2a | 20 20 20 61 67 61 69 6e |hen.X. *|   again|
|00002340| 20 61 74 20 6d 6a 64 2d | 64 74 2c 20 77 68 65 72 | at mjd-|dt, wher|
|00002350| 65 0a 58 09 20 2a 20 20 | 20 64 74 20 3d 20 74 69 |e.X. *  | dt = ti|
|00002360| 6d 65 20 69 74 20 74 61 | 6b 65 73 20 6c 69 67 68 |me it ta|kes ligh|
|00002370| 74 20 74 6f 20 74 72 61 | 76 65 6c 20 65 61 72 74 |t to tra|vel eart|
|00002380| 68 2d 6f 62 6a 65 63 74 | 20 64 69 73 74 61 6e 63 |h-object| distanc|
|00002390| 65 2e 0a 58 09 20 2a 20 | 74 68 69 73 20 69 73 20 |e..X. * |this is |
|000023a0| 62 61 73 69 63 61 6c 6c | 79 20 74 68 65 20 73 61 |basicall|y the sa|
|000023b0| 6d 65 20 63 6f 64 65 20 | 61 73 20 70 65 6c 65 6d |me code |as pelem|
|000023c0| 65 6e 74 28 29 20 61 6e | 64 20 70 6c 61 6e 73 28 |ent() an|d plans(|
|000023d0| 29 0a 58 09 20 2a 20 20 | 20 63 6f 6d 62 69 6e 65 |).X. *  | combine|
|000023e0| 64 20 61 6e 64 20 73 69 | 6d 70 6c 69 66 69 65 64 |d and si|mplified|
|000023f0| 20 66 6f 72 20 74 68 65 | 20 73 70 65 63 69 61 6c | for the| special|
|00002400| 20 63 61 73 65 20 6f 66 | 20 6f 73 63 75 6c 61 74 | case of| osculat|
|00002410| 69 6e 67 0a 58 09 20 2a | 20 20 20 28 75 6e 70 65 |ing.X. *|   (unpe|
|00002420| 72 74 75 72 62 65 64 29 | 20 65 6c 65 6d 65 6e 74 |rturbed)| element|
|00002430| 73 2e 20 77 65 20 68 61 | 76 65 20 61 64 64 65 64 |s. we ha|ve added|
|00002440| 20 72 65 64 75 63 74 69 | 6f 6e 20 6f 66 20 65 6c | reducti|on of el|
|00002450| 65 6d 65 6e 74 73 20 75 | 73 69 6e 67 0a 58 09 20 |ements u|sing.X. |
|00002460| 2a 20 20 20 72 65 64 75 | 63 65 5f 65 6c 65 6d 65 |*   redu|ce_eleme|
|00002470| 6e 74 73 28 29 2e 0a 58 | 09 20 2a 2f 0a 58 09 64 |nts()..X|. */.X.d|
|00002480| 74 20 3d 20 30 3b 0a 58 | 09 66 6f 72 20 28 70 61 |t = 0;.X|.for (pa|
|00002490| 73 73 20 3d 20 70 72 65 | 66 5f 67 65 74 28 50 52 |ss = pre|f_get(PR|
|000024a0| 45 46 5f 41 4c 47 4f 29 | 3d 3d 50 52 45 46 5f 41 |EF_ALGO)|==PREF_A|
|000024b0| 43 43 55 52 41 54 45 20 | 3f 20 30 20 3a 20 31 3b |CCURATE |? 0 : 1;|
|000024c0| 20 70 61 73 73 3c 32 3b | 20 70 61 73 73 2b 2b 29 | pass<2;| pass++)|
|000024d0| 7b 0a 58 0a 58 09 20 20 | 20 20 72 65 64 75 63 65 |{.X.X.  |  reduce|
|000024e0| 5f 65 6c 65 6d 65 6e 74 | 73 20 28 6f 70 2d 3e 65 |_element|s (op->e|
|000024f0| 5f 65 70 6f 63 68 2c 20 | 6d 6a 64 2d 64 74 2c 20 |_epoch, |mjd-dt, |
|00002500| 64 65 67 72 61 64 28 6f | 70 2d 3e 65 5f 69 6e 63 |degrad(o|p->e_inc|
|00002510| 29 2c 0a 58 09 09 09 20 | 20 20 20 64 65 67 72 61 |),.X... |   degra|
|00002520| 64 20 28 6f 70 2d 3e 65 | 5f 6f 6d 29 2c 20 64 65 |d (op->e|_om), de|
|00002530| 67 72 61 64 20 28 6f 70 | 2d 3e 65 5f 4f 6d 29 2c |grad (op|->e_Om),|
|00002540| 0a 58 09 09 09 20 20 20 | 20 26 69 6e 63 2c 20 26 |.X...   | &inc, &|
|00002550| 6f 6d 2c 20 26 4f 6d 29 | 3b 0a 58 0a 58 09 20 20 |om, &Om)|;.X.X.  |
|00002560| 20 20 6d 61 20 3d 20 64 | 65 67 72 61 64 20 28 6f |  ma = d|egrad (o|
|00002570| 70 2d 3e 65 5f 4d 20 2b | 20 28 6d 6a 64 20 2d 20 |p->e_M +| (mjd - |
|00002580| 6f 70 2d 3e 65 5f 63 65 | 70 6f 63 68 20 2d 20 64 |op->e_ce|poch - d|
|00002590| 74 29 20 2a 20 6f 70 2d | 3e 65 5f 6e 29 3b 0a 58 |t) * op-|>e_n);.X|
|000025a0| 09 20 20 20 20 61 6e 6f | 6d 61 6c 79 20 28 6d 61 |.    ano|maly (ma|
|000025b0| 2c 20 65 2c 20 26 6e 75 | 2c 20 26 65 61 29 3b 0a |, e, &nu|, &ea);.|
|000025c0| 58 09 20 20 20 20 72 70 | 20 3d 20 6f 70 2d 3e 65 |X.    rp| = op->e|
|000025d0| 5f 61 20 2a 20 28 31 2d | 65 2a 65 29 20 2f 20 28 |_a * (1-|e*e) / (|
|000025e0| 31 2b 65 2a 63 6f 73 28 | 6e 75 29 29 3b 0a 58 09 |1+e*cos(|nu));.X.|
|000025f0| 20 20 20 20 6c 6f 20 3d | 20 6e 75 20 2b 20 6f 6d |    lo =| nu + om|
|00002600| 3b 0a 58 09 20 20 20 20 | 73 6c 6f 20 3d 20 73 69 |;.X.    |slo = si|
|00002610| 6e 28 6c 6f 29 3b 0a 58 | 09 20 20 20 20 63 6c 6f |n(lo);.X|.    clo|
|00002620| 20 3d 20 63 6f 73 28 6c | 6f 29 3b 0a 58 09 20 20 | = cos(l|o);.X.  |
|00002630| 20 20 73 70 73 69 20 3d | 20 73 6c 6f 2a 73 69 6e |  spsi =| slo*sin|
|00002640| 28 69 6e 63 29 3b 0a 58 | 09 20 20 20 20 79 20 3d |(inc);.X|.    y =|
|00002650| 20 73 6c 6f 2a 63 6f 73 | 28 69 6e 63 29 3b 0a 58 | slo*cos|(inc);.X|
|00002660| 09 20 20 20 20 70 73 69 | 20 3d 20 61 73 69 6e 28 |.    psi| = asin(|
|00002670| 73 70 73 69 29 3b 0a 58 | 09 20 20 20 20 6c 70 64 |spsi);.X|.    lpd|
|00002680| 20 3d 20 61 74 61 6e 28 | 79 2f 63 6c 6f 29 2b 4f | = atan(|y/clo)+O|
|00002690| 6d 3b 0a 58 09 20 20 20 | 20 69 66 20 28 63 6c 6f |m;.X.   | if (clo|
|000026a0| 3c 30 29 20 6c 70 64 20 | 2b 3d 20 50 49 3b 0a 58 |<0) lpd |+= PI;.X|
|000026b0| 09 20 20 20 20 72 61 6e | 67 65 20 28 26 6c 70 64 |.    ran|ge (&lpd|
|000026c0| 2c 20 32 2a 50 49 29 3b | 0a 58 09 20 20 20 20 63 |, 2*PI);|.X.    c|
|000026d0| 70 73 69 20 3d 20 63 6f | 73 28 70 73 69 29 3b 0a |psi = co|s(psi);.|
|000026e0| 58 09 20 20 20 20 72 70 | 64 20 3d 20 72 70 2a 63 |X.    rp|d = rp*c|
|000026f0| 70 73 69 3b 0a 58 09 20 | 20 20 20 6c 6c 20 3d 20 |psi;.X. |   ll = |
|00002700| 6c 70 64 2d 6c 67 3b 0a | 58 09 20 20 20 20 72 68 |lpd-lg;.|X.    rh|
|00002710| 6f 20 3d 20 73 71 72 74 | 28 72 73 6e 2a 72 73 6e |o = sqrt|(rsn*rsn|
|00002720| 2b 72 70 2a 72 70 2d 32 | 2a 72 73 6e 2a 72 70 2a |+rp*rp-2|*rsn*rp*|
|00002730| 63 70 73 69 2a 63 6f 73 | 28 6c 6c 29 29 3b 0a 58 |cpsi*cos|(ll));.X|
|00002740| 0a 58 09 20 20 20 20 64 | 74 20 3d 20 72 68 6f 2a |.X.    d|t = rho*|
|00002750| 4c 54 41 55 2f 33 36 30 | 30 2e 30 2f 32 34 2e 30 |LTAU/360|0.0/24.0|
|00002760| 3b 09 2f 2a 20 6c 69 67 | 68 74 20 74 72 61 76 65 |;./* lig|ht trave|
|00002770| 6c 20 74 69 6d 65 2c 20 | 69 6e 20 64 61 79 73 20 |l time, |in days |
|00002780| 2f 20 41 55 20 2a 2f 0a | 58 09 7d 0a 58 0a 58 09 |/ AU */.|X.}.X.X.|
|00002790| 2f 2a 20 63 6f 6d 70 75 | 74 65 20 73 69 6e 20 61 |/* compu|te sin a|
|000027a0| 6e 64 20 63 6f 73 20 6f | 66 20 6c 6c 20 2a 2f 0a |nd cos o|f ll */.|
|000027b0| 58 09 73 6c 6c 20 3d 20 | 73 69 6e 28 6c 6c 29 3b |X.sll = |sin(ll);|
|000027c0| 0a 58 09 63 6c 6c 20 3d | 20 63 6f 73 28 6c 6c 29 |.X.cll =| cos(ll)|
|000027d0| 3b 0a 58 0a 58 09 2f 2a | 20 66 69 6e 64 20 67 65 |;.X.X./*| find ge|
|000027e0| 6f 63 65 6e 74 72 69 63 | 20 65 63 6c 69 70 74 69 |ocentric| eclipti|
|000027f0| 63 20 6c 6f 6e 67 69 74 | 75 64 65 20 61 6e 64 20 |c longit|ude and |
|00002800| 6c 61 74 69 74 75 64 65 | 20 2a 2f 0a 58 09 69 66 |latitude| */.X.if|
|00002810| 20 28 72 70 64 20 3c 20 | 72 73 6e 29 0a 58 09 20 | (rpd < |rsn).X. |
|00002820| 20 20 20 6c 61 6d 20 3d | 20 61 74 61 6e 28 2d 31 |   lam =| atan(-1|
|00002830| 2a 72 70 64 2a 73 6c 6c | 2f 28 72 73 6e 2d 72 70 |*rpd*sll|/(rsn-rp|
|00002840| 64 2a 63 6c 6c 29 29 2b | 6c 67 2b 50 49 3b 0a 58 |d*cll))+|lg+PI;.X|
|00002850| 09 65 6c 73 65 0a 58 09 | 20 20 20 20 6c 61 6d 20 |.else.X.|    lam |
|00002860| 3d 20 61 74 61 6e 28 72 | 73 6e 2a 73 6c 6c 2f 28 |= atan(r|sn*sll/(|
|00002870| 72 70 64 2d 72 73 6e 2a | 63 6c 6c 29 29 2b 6c 70 |rpd-rsn*|cll))+lp|
|00002880| 64 3b 0a 58 09 72 61 6e | 67 65 20 28 26 6c 61 6d |d;.X.ran|ge (&lam|
|00002890| 2c 20 32 2a 50 49 29 3b | 0a 58 09 62 65 74 20 3d |, 2*PI);|.X.bet =|
|000028a0| 20 61 74 61 6e 28 72 70 | 64 2a 73 70 73 69 2a 73 | atan(rp|d*spsi*s|
|000028b0| 69 6e 28 6c 61 6d 2d 6c | 70 64 29 2f 28 63 70 73 |in(lam-l|pd)/(cps|
|000028c0| 69 2a 72 73 6e 2a 73 6c | 6c 29 29 3b 0a 58 0a 58 |i*rsn*sl|l));.X.X|
|000028d0| 09 2f 2a 20 66 69 6c 6c | 20 69 6e 20 61 6c 6c 20 |./* fill| in all |
|000028e0| 6f 66 20 6f 70 2d 3e 73 | 5f 2a 20 73 74 75 66 66 |of op->s|_* stuff|
|000028f0| 20 65 78 63 65 70 74 20 | 73 5f 73 69 7a 65 20 61 | except |s_size a|
|00002900| 6e 64 20 73 5f 6d 61 67 | 20 2a 2f 0a 58 09 63 69 |nd s_mag| */.X.ci|
|00002910| 72 5f 73 6b 79 20 28 6e | 70 2c 20 6c 70 64 2c 20 |r_sky (n|p, lpd, |
|00002920| 70 73 69 2c 20 72 70 2c | 20 72 68 6f 2c 20 6c 61 |psi, rp,| rho, la|
|00002930| 6d 2c 20 62 65 74 2c 20 | 6c 73 6e 2c 20 72 73 6e |m, bet, |lsn, rsn|
|00002940| 2c 20 6f 70 29 3b 0a 58 | 0a 58 09 2f 2a 20 63 6f |, op);.X|.X./* co|
|00002950| 6d 70 75 74 65 20 6d 61 | 67 6e 69 74 75 64 65 20 |mpute ma|gnitude |
|00002960| 61 6e 64 20 73 69 7a 65 | 20 2a 2f 0a 58 09 69 66 |and size| */.X.if|
|00002970| 20 28 6f 70 2d 3e 65 5f | 6d 61 67 2e 77 68 69 63 | (op->e_|mag.whic|
|00002980| 68 6d 20 3d 3d 20 4d 41 | 47 5f 48 47 29 20 7b 0a |hm == MA|G_HG) {.|
|00002990| 58 09 20 20 20 20 2f 2a | 20 74 68 65 20 48 20 61 |X.    /*| the H a|
|000029a0| 6e 64 20 47 20 70 61 72 | 61 6d 65 74 65 72 73 20 |nd G par|ameters |
|000029b0| 66 72 6f 6d 20 74 68 65 | 20 41 73 74 72 6f 2e 20 |from the| Astro. |
|000029c0| 41 6c 6d 61 6e 61 63 2e | 0a 58 09 20 20 20 20 20 |Almanac.|.X.     |
|000029d0| 2a 2f 0a 58 09 20 20 20 | 20 68 67 5f 6d 61 67 20 |*/.X.   | hg_mag |
|000029e0| 28 6f 70 2d 3e 65 5f 6d | 61 67 2e 6d 31 2c 20 6f |(op->e_m|ag.m1, o|
|000029f0| 70 2d 3e 65 5f 6d 61 67 | 2e 6d 32 2c 20 72 70 2c |p->e_mag|.m2, rp,|
|00002a00| 20 72 68 6f 2c 20 72 73 | 6e 2c 20 26 6d 61 67 29 | rho, rs|n, &mag)|
|00002a10| 3b 0a 58 09 7d 20 65 6c | 73 65 20 7b 0a 58 09 20 |;.X.} el|se {.X. |
|00002a20| 20 20 20 2f 2a 20 74 68 | 65 20 67 2f 6b 20 6d 6f |   /* th|e g/k mo|
|00002a30| 64 65 6c 20 6f 66 20 63 | 6f 6d 65 74 73 20 2a 2f |del of c|omets */|
|00002a40| 0a 58 09 20 20 20 20 67 | 6b 5f 6d 61 67 20 28 6f |.X.    g|k_mag (o|
|00002a50| 70 2d 3e 65 5f 6d 61 67 | 2e 6d 31 2c 20 6f 70 2d |p->e_mag|.m1, op-|
|00002a60| 3e 65 5f 6d 61 67 2e 6d | 32 2c 20 72 70 2c 20 72 |>e_mag.m|2, rp, r|
|00002a70| 68 6f 2c 20 26 6d 61 67 | 29 3b 0a 58 09 7d 0a 58 |ho, &mag|);.X.}.X|
|00002a80| 09 6f 70 2d 3e 73 5f 6d | 61 67 20 3d 20 6d 61 67 |.op->s_m|ag = mag|
|00002a90| 20 2a 20 4d 41 47 53 43 | 41 4c 45 3b 0a 58 09 6f | * MAGSC|ALE;.X.o|
|00002aa0| 70 2d 3e 73 5f 73 69 7a | 65 20 3d 20 6f 70 2d 3e |p->s_siz|e = op->|
|00002ab0| 65 5f 73 69 7a 65 20 2f | 20 72 68 6f 3b 0a 58 0a |e_size /| rho;.X.|
|00002ac0| 58 09 72 65 74 75 72 6e | 20 28 30 29 3b 0a 58 7d |X.return| (0);.X}|
|00002ad0| 0a 58 0a 58 2f 2a 20 63 | 6f 6d 70 75 74 65 20 73 |.X.X/* c|ompute s|
|00002ae0| 6b 79 20 63 69 72 63 75 | 6d 73 74 61 6e 63 65 73 |ky circu|mstances|
|00002af0| 20 6f 66 20 61 6e 20 6f | 62 6a 65 63 74 20 69 6e | of an o|bject in|
|00002b00| 20 68 65 6c 69 6f 63 65 | 6e 74 72 69 63 20 68 79 | helioce|ntric hy|
|00002b10| 70 65 72 62 6f 6c 69 63 | 20 6f 72 62 69 74 2e 0a |perbolic| orbit..|
|00002b20| 58 20 2a 2f 0a 58 73 74 | 61 74 69 63 0a 58 6f 62 |X */.Xst|atic.Xob|
|00002b30| 6a 5f 68 79 70 65 72 62 | 6f 6c 69 63 20 28 6e 70 |j_hyperb|olic (np|
|00002b40| 2c 20 6f 70 29 0a 58 4e | 6f 77 20 2a 6e 70 3b 0a |, op).XN|ow *np;.|
|00002b50| 58 4f 62 6a 20 2a 6f 70 | 3b 0a 58 7b 0a 58 09 64 |XObj *op|;.X{.X.d|
|00002b60| 6f 75 62 6c 65 20 6c 73 | 6e 2c 20 72 73 6e 3b 09 |ouble ls|n, rsn;.|
|00002b70| 2f 2a 20 74 72 75 65 20 | 67 65 6f 63 20 6c 6e 67 |/* true |geoc lng|
|00002b80| 20 6f 66 20 73 75 6e 3b | 20 64 69 73 74 20 66 72 | of sun;| dist fr|
|00002b90| 6f 6d 20 73 6e 20 74 6f | 20 65 61 72 74 68 2a 2f |om sn to| earth*/|
|00002ba0| 0a 58 09 64 6f 75 62 6c | 65 20 64 74 3b 09 09 2f |.X.doubl|e dt;../|
|00002bb0| 2a 20 6c 69 67 68 74 20 | 74 72 61 76 65 6c 20 74 |* light |travel t|
|00002bc0| 69 6d 65 20 74 6f 20 6f | 62 6a 65 63 74 20 2a 2f |ime to o|bject */|
|00002bd0| 0a 58 09 64 6f 75 62 6c | 65 20 6c 67 3b 09 09 2f |.X.doubl|e lg;../|
|00002be0| 2a 20 68 65 6c 69 6f 20 | 6c 6f 6e 67 20 6f 66 20 |* helio |long of |
|00002bf0| 65 61 72 74 68 20 2a 2f | 0a 58 09 64 6f 75 62 6c |earth */|.X.doubl|
|00002c00| 65 20 6e 75 2c 20 65 61 | 3b 09 09 2f 2a 20 74 72 |e nu, ea|;../* tr|
|00002c10| 75 65 20 61 6e 6f 6d 61 | 6c 79 20 61 6e 64 20 65 |ue anoma|ly and e|
|00002c20| 63 63 65 6e 74 72 69 63 | 20 61 6e 6f 6d 61 6c 79 |ccentric| anomaly|
|00002c30| 20 2a 2f 0a 58 09 64 6f | 75 62 6c 65 20 6d 61 3b | */.X.do|uble ma;|
|00002c40| 09 09 2f 2a 20 6d 65 61 | 6e 20 61 6e 6f 6d 61 6c |../* mea|n anomal|
|00002c50| 79 20 2a 2f 0a 58 09 64 | 6f 75 62 6c 65 20 72 70 |y */.X.d|ouble rp|
|00002c60| 3b 09 09 2f 2a 20 64 69 | 73 74 61 6e 63 65 20 66 |;../* di|stance f|
|00002c70| 72 6f 6d 20 74 68 65 20 | 73 75 6e 20 2a 2f 0a 58 |rom the |sun */.X|
|00002c80| 09 64 6f 75 62 6c 65 20 | 6c 6f 2c 20 73 6c 6f 2c |.double |lo, slo,|
|00002c90| 20 63 6c 6f 3b 09 2f 2a | 20 61 6e 67 6c 65 20 66 | clo;./*| angle f|
|00002ca0| 72 6f 6d 20 61 73 63 65 | 6e 64 69 6e 67 20 6e 6f |rom asce|nding no|
|00002cb0| 64 65 20 2a 2f 0a 58 09 | 64 6f 75 62 6c 65 20 69 |de */.X.|double i|
|00002cc0| 6e 63 3b 09 09 2f 2a 20 | 69 6e 63 6c 69 6e 61 74 |nc;../* |inclinat|
|00002cd0| 69 6f 6e 20 2a 2f 0a 58 | 09 64 6f 75 62 6c 65 20 |ion */.X|.double |
|00002ce0| 70 73 69 2c 20 73 70 73 | 69 2c 20 63 70 73 69 3b |psi, sps|i, cpsi;|
|00002cf0| 09 2f 2a 20 68 65 6c 69 | 6f 63 65 6e 74 72 69 63 |./* heli|ocentric|
|00002d00| 20 6c 61 74 69 74 75 64 | 65 20 2a 2f 0a 58 09 64 | latitud|e */.X.d|
|00002d10| 6f 75 62 6c 65 20 6c 70 | 64 3b 20 09 09 2f 2a 20 |ouble lp|d; ../* |
|00002d20| 68 65 6c 69 6f 63 65 6e | 74 72 69 63 20 6c 6f 6e |heliocen|tric lon|
|00002d30| 67 69 74 75 64 65 20 2a | 2f 0a 58 09 64 6f 75 62 |gitude *|/.X.doub|
|00002d40| 6c 65 20 72 68 6f 3b 09 | 09 2f 2a 20 64 69 73 74 |le rho;.|./* dist|
|00002d50| 61 6e 63 65 20 66 72 6f | 6d 20 74 68 65 20 45 61 |ance fro|m the Ea|
|00002d60| 72 74 68 20 2a 2f 0a 58 | 09 64 6f 75 62 6c 65 20 |rth */.X|.double |
|00002d70| 6f 6d 3b 09 09 2f 2a 20 | 61 72 67 20 6f 66 20 70 |om;../* |arg of p|
|00002d80| 65 72 69 68 65 6c 69 6f | 6e 20 2a 2f 0a 58 09 64 |erihelio|n */.X.d|
|00002d90| 6f 75 62 6c 65 20 4f 6d | 3b 09 09 2f 2a 20 6c 6f |ouble Om|;../* lo|
|00002da0| 6e 67 20 6f 66 20 61 73 | 63 65 6e 64 69 6e 67 20 |ng of as|cending |
|00002db0| 6e 6f 64 65 2e 20 2a 2f | 0a 58 09 64 6f 75 62 6c |node. */|.X.doubl|
|00002dc0| 65 20 6c 61 6d 3b 20 20 | 20 20 09 09 2f 2a 20 67 |e lam;  |  ../* g|
|00002dd0| 65 6f 63 65 6e 74 72 69 | 63 20 65 63 6c 69 70 74 |eocentri|c eclipt|
|00002de0| 69 63 20 6c 6f 6e 67 69 | 74 75 64 65 20 2a 2f 0a |ic longi|tude */.|
|00002df0| 58 09 64 6f 75 62 6c 65 | 20 62 65 74 3b 20 20 20 |X.double| bet;   |
|00002e00| 20 09 09 2f 2a 20 67 65 | 6f 63 65 6e 74 72 69 63 | ../* ge|ocentric|
|00002e10| 20 65 63 6c 69 70 74 69 | 63 20 6c 61 74 69 74 75 | eclipti|c latitu|
|00002e20| 64 65 20 2a 2f 0a 58 09 | 64 6f 75 62 6c 65 20 65 |de */.X.|double e|
|00002e30| 3b 09 09 2f 2a 20 66 61 | 73 74 20 65 63 63 65 6e |;../* fa|st eccen|
|00002e40| 74 72 69 63 69 74 79 20 | 2a 2f 0a 58 09 64 6f 75 |tricity |*/.X.dou|
|00002e50| 62 6c 65 20 6c 6c 2c 20 | 73 6c 6c 2c 20 63 6c 6c |ble ll, |sll, cll|
|00002e60| 3b 09 2f 2a 20 68 65 6c | 69 6f 20 61 6e 67 6c 65 |;./* hel|io angle|
|00002e70| 20 62 65 74 77 65 65 6e | 20 6f 62 6a 65 63 74 20 | between| object |
|00002e80| 61 6e 64 20 65 61 72 74 | 68 20 2a 2f 0a 58 09 64 |and eart|h */.X.d|
|00002e90| 6f 75 62 6c 65 20 6e 3b | 09 09 2f 2a 20 6d 65 61 |ouble n;|../* mea|
|00002ea0| 6e 20 64 61 69 6c 79 20 | 6d 6f 74 69 6f 6e 20 2a |n daily |motion *|
|00002eb0| 2f 0a 58 09 64 6f 75 62 | 6c 65 20 6d 61 67 3b 09 |/.X.doub|le mag;.|
|00002ec0| 09 2f 2a 20 6d 61 67 6e | 69 74 75 64 65 20 2a 2f |./* magn|itude */|
|00002ed0| 0a 58 09 64 6f 75 62 6c | 65 20 61 3b 09 09 2f 2a |.X.doubl|e a;../*|
|00002ee0| 20 6d 65 61 6e 20 64 69 | 73 74 61 6e 63 65 20 2a | mean di|stance *|
|00002ef0| 2f 0a 58 09 64 6f 75 62 | 6c 65 20 72 70 64 3b 0a |/.X.doub|le rpd;.|
|00002f00| 58 09 64 6f 75 62 6c 65 | 20 79 3b 0a 58 09 69 6e |X.double| y;.X.in|
|00002f10| 74 20 70 61 73 73 3b 0a | 58 0a 58 09 2f 2a 20 66 |t pass;.|X.X./* f|
|00002f20| 69 6e 64 20 73 6f 6c 61 | 72 20 65 63 6c 69 70 74 |ind sola|r eclipt|
|00002f30| 69 63 61 6c 20 6c 6f 6e | 67 69 74 75 64 65 20 61 |ical lon|gitude a|
|00002f40| 6e 64 20 64 69 73 74 61 | 6e 63 65 20 74 6f 20 73 |nd dista|nce to s|
|00002f50| 75 6e 20 66 72 6f 6d 20 | 65 61 72 74 68 20 2a 2f |un from |earth */|
|00002f60| 0a 58 09 73 75 6e 70 6f | 73 20 28 6d 6a 64 2c 20 |.X.sunpo|s (mjd, |
|00002f70| 26 6c 73 6e 2c 20 26 72 | 73 6e 29 3b 0a 58 0a 58 |&lsn, &r|sn);.X.X|
|00002f80| 09 6c 67 20 3d 20 6c 73 | 6e 20 2b 20 50 49 3b 0a |.lg = ls|n + PI;.|
|00002f90| 58 09 65 20 3d 20 6f 70 | 2d 3e 68 5f 65 3b 0a 58 |X.e = op|->h_e;.X|
|00002fa0| 09 61 20 3d 20 6f 70 2d | 3e 68 5f 71 70 2f 28 65 |.a = op-|>h_qp/(e|
|00002fb0| 20 2d 20 31 2e 30 29 3b | 0a 58 09 6e 20 3d 20 2e | - 1.0);|.X.n = .|
|00002fc0| 39 38 35 36 33 2f 73 71 | 72 74 28 61 2a 61 2a 61 |98563/sq|rt(a*a*a|
|00002fd0| 29 3b 0a 58 0a 58 09 2f | 2a 20 63 6f 72 72 65 63 |);.X.X./|* correc|
|00002fe0| 74 20 66 6f 72 20 6c 69 | 67 68 74 20 74 69 6d 65 |t for li|ght time|
|00002ff0| 20 62 79 20 63 6f 6d 70 | 75 74 69 6e 67 20 70 6f | by comp|uting po|
|00003000| 73 69 74 69 6f 6e 20 61 | 74 20 74 69 6d 65 20 6d |sition a|t time m|
|00003010| 6a 64 2c 20 74 68 65 6e | 0a 58 09 20 2a 20 20 20 |jd, then|.X. *   |
|00003020| 61 67 61 69 6e 20 61 74 | 20 6d 6a 64 2d 64 74 2c |again at| mjd-dt,|
|00003030| 20 77 68 65 72 65 0a 58 | 09 20 2a 20 20 20 64 74 | where.X|. *   dt|
|00003040| 20 3d 20 74 69 6d 65 20 | 69 74 20 74 61 6b 65 73 | = time |it takes|
|00003050| 20 6c 69 67 68 74 20 74 | 6f 20 74 72 61 76 65 6c | light t|o travel|
|00003060| 20 65 61 72 74 68 2d 6f | 62 6a 65 63 74 20 64 69 | earth-o|bject di|
|00003070| 73 74 61 6e 63 65 2e 0a | 58 09 20 2a 2f 0a 58 09 |stance..|X. */.X.|
|00003080| 64 74 20 3d 20 30 3b 0a | 58 09 66 6f 72 20 28 70 |dt = 0;.|X.for (p|
|00003090| 61 73 73 20 3d 20 70 72 | 65 66 5f 67 65 74 28 50 |ass = pr|ef_get(P|
|000030a0| 52 45 46 5f 41 4c 47 4f | 29 3d 3d 50 52 45 46 5f |REF_ALGO|)==PREF_|
|000030b0| 41 43 43 55 52 41 54 45 | 20 3f 20 30 20 3a 20 31 |ACCURATE| ? 0 : 1|
|000030c0| 3b 20 70 61 73 73 3c 32 | 3b 20 70 61 73 73 2b 2b |; pass<2|; pass++|
|000030d0| 29 7b 0a 58 0a 58 09 20 | 20 20 20 72 65 64 75 63 |){.X.X. |   reduc|
|000030e0| 65 5f 65 6c 65 6d 65 6e | 74 73 20 28 6f 70 2d 3e |e_elemen|ts (op->|
|000030f0| 68 5f 65 70 6f 63 68 2c | 20 6d 6a 64 2d 64 74 2c |h_epoch,| mjd-dt,|
|00003100| 20 64 65 67 72 61 64 28 | 6f 70 2d 3e 68 5f 69 6e | degrad(|op->h_in|
|00003110| 63 29 2c 0a 58 09 09 09 | 20 20 20 20 64 65 67 72 |c),.X...|    degr|
|00003120| 61 64 20 28 6f 70 2d 3e | 68 5f 6f 6d 29 2c 20 64 |ad (op->|h_om), d|
|00003130| 65 67 72 61 64 20 28 6f | 70 2d 3e 68 5f 4f 6d 29 |egrad (o|p->h_Om)|
|00003140| 2c 0a 58 09 09 09 20 20 | 20 20 26 69 6e 63 2c 20 |,.X...  |  &inc, |
|00003150| 26 6f 6d 2c 20 26 4f 6d | 29 3b 0a 58 0a 58 09 20 |&om, &Om|);.X.X. |
|00003160| 20 20 20 6d 61 20 3d 20 | 64 65 67 72 61 64 20 28 |   ma = |degrad (|
|00003170| 28 6d 6a 64 20 2d 20 6f | 70 2d 3e 68 5f 65 70 20 |(mjd - o|p->h_ep |
|00003180| 2d 20 64 74 29 20 2a 20 | 6e 29 3b 0a 58 09 20 20 |- dt) * |n);.X.  |
|00003190| 20 20 61 6e 6f 6d 61 6c | 79 20 28 6d 61 2c 20 65 |  anomal|y (ma, e|
|000031a0| 2c 20 26 6e 75 2c 20 26 | 65 61 29 3b 0a 58 09 20 |, &nu, &|ea);.X. |
|000031b0| 20 20 20 72 70 20 3d 20 | 61 20 2a 20 28 65 2a 65 |   rp = |a * (e*e|
|000031c0| 2d 31 2e 30 29 20 2f 20 | 28 31 2e 30 2b 65 2a 63 |-1.0) / |(1.0+e*c|
|000031d0| 6f 73 28 6e 75 29 29 3b | 0a 58 09 20 20 20 20 6c |os(nu));|.X.    l|
|000031e0| 6f 20 3d 20 6e 75 20 2b | 20 6f 6d 3b 0a 58 09 20 |o = nu +| om;.X. |
|000031f0| 20 20 20 73 6c 6f 20 3d | 20 73 69 6e 28 6c 6f 29 |   slo =| sin(lo)|
|00003200| 3b 0a 58 09 20 20 20 20 | 63 6c 6f 20 3d 20 63 6f |;.X.    |clo = co|
|00003210| 73 28 6c 6f 29 3b 0a 58 | 09 20 20 20 20 73 70 73 |s(lo);.X|.    sps|
|00003220| 69 20 3d 20 73 6c 6f 2a | 73 69 6e 28 69 6e 63 29 |i = slo*|sin(inc)|
|00003230| 3b 0a 58 09 20 20 20 20 | 79 20 3d 20 73 6c 6f 2a |;.X.    |y = slo*|
|00003240| 63 6f 73 28 69 6e 63 29 | 3b 0a 58 09 20 20 20 20 |cos(inc)|;.X.    |
|00003250| 70 73 69 20 3d 20 61 73 | 69 6e 28 73 70 73 69 29 |psi = as|in(spsi)|
|00003260| 3b 0a 58 09 20 20 20 20 | 6c 70 64 20 3d 20 61 74 |;.X.    |lpd = at|
|00003270| 61 6e 28 79 2f 63 6c 6f | 29 2b 4f 6d 3b 0a 58 09 |an(y/clo|)+Om;.X.|
|00003280| 20 20 20 20 69 66 20 28 | 63 6c 6f 3c 30 29 20 6c |    if (|clo<0) l|
|00003290| 70 64 20 2b 3d 20 50 49 | 3b 0a 58 09 20 20 20 20 |pd += PI|;.X.    |
|000032a0| 72 61 6e 67 65 20 28 26 | 6c 70 64 2c 20 32 2a 50 |range (&|lpd, 2*P|
|000032b0| 49 29 3b 0a 58 09 20 20 | 20 20 63 70 73 69 20 3d |I);.X.  |  cpsi =|
|000032c0| 20 63 6f 73 28 70 73 69 | 29 3b 0a 58 09 20 20 20 | cos(psi|);.X.   |
|000032d0| 20 72 70 64 20 3d 20 72 | 70 2a 63 70 73 69 3b 0a | rpd = r|p*cpsi;.|
|000032e0| 58 09 20 20 20 20 6c 6c | 20 3d 20 6c 70 64 2d 6c |X.    ll| = lpd-l|
|000032f0| 67 3b 0a 58 09 20 20 20 | 20 72 68 6f 20 3d 20 73 |g;.X.   | rho = s|
|00003300| 71 72 74 28 72 73 6e 2a | 72 73 6e 2b 72 70 2a 72 |qrt(rsn*|rsn+rp*r|
|00003310| 70 2d 32 2a 72 73 6e 2a | 72 70 2a 63 70 73 69 2a |p-2*rsn*|rp*cpsi*|
|00003320| 63 6f 73 28 6c 6c 29 29 | 3b 0a 58 0a 58 09 20 20 |cos(ll))|;.X.X.  |
|00003330| 20 20 64 74 20 3d 20 72 | 68 6f 2a 35 2e 37 37 35 |  dt = r|ho*5.775|
|00003340| 35 31 38 65 2d 33 3b 09 | 2f 2a 20 6c 69 67 68 74 |518e-3;.|/* light|
|00003350| 20 74 72 61 76 65 6c 20 | 74 69 6d 65 2c 20 69 6e | travel |time, in|
|00003360| 20 64 61 79 73 20 2a 2f | 0a 58 09 7d 0a 58 0a 58 | days */|.X.}.X.X|
|00003370| 09 2f 2a 20 63 6f 6d 70 | 75 74 65 20 73 69 6e 20 |./* comp|ute sin |
|00003380| 61 6e 64 20 63 6f 73 20 | 6f 66 20 6c 6c 20 2a 2f |and cos |of ll */|
|00003390| 0a 58 09 73 6c 6c 20 3d | 20 73 69 6e 28 6c 6c 29 |.X.sll =| sin(ll)|
|000033a0| 3b 0a 58 09 63 6c 6c 20 | 3d 20 63 6f 73 28 6c 6c |;.X.cll |= cos(ll|
|000033b0| 29 3b 0a 58 0a 58 09 2f | 2a 20 66 69 6e 64 20 67 |);.X.X./|* find g|
|000033c0| 65 6f 63 65 6e 74 72 69 | 63 20 65 63 6c 69 70 74 |eocentri|c eclipt|
|000033d0| 69 63 20 6c 6f 6e 67 69 | 74 75 64 65 20 61 6e 64 |ic longi|tude and|
|000033e0| 20 6c 61 74 69 74 75 64 | 65 20 2a 2f 0a 58 09 69 | latitud|e */.X.i|
|000033f0| 66 20 28 72 70 64 20 3c | 20 72 73 6e 29 0a 58 09 |f (rpd <| rsn).X.|
|00003400| 20 20 20 20 6c 61 6d 20 | 3d 20 61 74 61 6e 28 2d |    lam |= atan(-|
|00003410| 31 2a 72 70 64 2a 73 6c | 6c 2f 28 72 73 6e 2d 72 |1*rpd*sl|l/(rsn-r|
|00003420| 70 64 2a 63 6c 6c 29 29 | 2b 6c 67 2b 50 49 3b 0a |pd*cll))|+lg+PI;.|
|00003430| 58 09 65 6c 73 65 0a 58 | 09 20 20 20 20 6c 61 6d |X.else.X|.    lam|
|00003440| 20 3d 20 61 74 61 6e 28 | 72 73 6e 2a 73 6c 6c 2f | = atan(|rsn*sll/|
|00003450| 28 72 70 64 2d 72 73 6e | 2a 63 6c 6c 29 29 2b 6c |(rpd-rsn|*cll))+l|
|00003460| 70 64 3b 0a 58 09 72 61 | 6e 67 65 20 28 26 6c 61 |pd;.X.ra|nge (&la|
|00003470| 6d 2c 20 32 2a 50 49 29 | 3b 0a 58 09 62 65 74 20 |m, 2*PI)|;.X.bet |
|00003480| 3d 20 61 74 61 6e 28 72 | 70 64 2a 73 70 73 69 2a |= atan(r|pd*spsi*|
|00003490| 73 69 6e 28 6c 61 6d 2d | 6c 70 64 29 2f 28 63 70 |sin(lam-|lpd)/(cp|
|000034a0| 73 69 2a 72 73 6e 2a 73 | 6c 6c 29 29 3b 0a 58 0a |si*rsn*s|ll));.X.|
|000034b0| 58 09 2f 2a 20 66 69 6c | 6c 20 69 6e 20 61 6c 6c |X./* fil|l in all|
|000034c0| 20 6f 66 20 6f 70 2d 3e | 73 5f 2a 20 73 74 75 66 | of op->|s_* stuf|
|000034d0| 66 20 65 78 63 65 70 74 | 20 73 5f 73 69 7a 65 20 |f except| s_size |
|000034e0| 61 6e 64 20 73 5f 6d 61 | 67 20 2a 2f 0a 58 09 63 |and s_ma|g */.X.c|
|000034f0| 69 72 5f 73 6b 79 20 28 | 6e 70 2c 20 6c 70 64 2c |ir_sky (|np, lpd,|
|00003500| 20 70 73 69 2c 20 72 70 | 2c 20 72 68 6f 2c 20 6c | psi, rp|, rho, l|
|00003510| 61 6d 2c 20 62 65 74 2c | 20 6c 73 6e 2c 20 72 73 |am, bet,| lsn, rs|
|00003520| 6e 2c 20 6f 70 29 3b 0a | 58 0a 58 09 2f 2a 20 63 |n, op);.|X.X./* c|
|00003530| 6f 6d 70 75 74 65 20 6d | 61 67 6e 69 74 75 64 65 |ompute m|agnitude|
|00003540| 20 61 6e 64 20 73 69 7a | 65 20 2a 2f 0a 58 09 67 | and siz|e */.X.g|
|00003550| 6b 5f 6d 61 67 20 28 6f | 70 2d 3e 68 5f 67 2c 20 |k_mag (o|p->h_g, |
|00003560| 6f 70 2d 3e 68 5f 6b 2c | 20 72 70 2c 20 72 68 6f |op->h_k,| rp, rho|
|00003570| 2c 20 26 6d 61 67 29 3b | 0a 58 09 6f 70 2d 3e 73 |, &mag);|.X.op->s|
|00003580| 5f 6d 61 67 20 3d 20 6d | 61 67 20 2a 20 4d 41 47 |_mag = m|ag * MAG|
|00003590| 53 43 41 4c 45 3b 0a 58 | 09 6f 70 2d 3e 73 5f 73 |SCALE;.X|.op->s_s|
|000035a0| 69 7a 65 20 3d 20 6f 70 | 2d 3e 68 5f 73 69 7a 65 |ize = op|->h_size|
|000035b0| 20 2f 20 72 68 6f 3b 0a | 58 0a 58 09 72 65 74 75 | / rho;.|X.X.retu|
|000035c0| 72 6e 20 28 30 29 3b 0a | 58 7d 0a 58 0a 58 2f 2a |rn (0);.|X}.X.X/*|
|000035d0| 20 63 6f 6d 70 75 74 65 | 20 73 6b 79 20 63 69 72 | compute| sky cir|
|000035e0| 63 75 6d 73 74 61 6e 63 | 65 73 20 6f 66 20 61 6e |cumstanc|es of an|
|000035f0| 20 6f 62 6a 65 63 74 20 | 69 6e 20 68 65 6c 69 6f | object |in helio|
|00003600| 63 65 6e 74 72 69 63 20 | 68 79 70 65 72 62 6f 6c |centric |hyperbol|
|00003610| 69 63 20 6f 72 62 69 74 | 2e 0a 58 20 2a 2f 0a 58 |ic orbit|..X */.X|
|00003620| 73 74 61 74 69 63 0a 58 | 6f 62 6a 5f 70 61 72 61 |static.X|obj_para|
|00003630| 62 6f 6c 69 63 20 28 6e | 70 2c 20 6f 70 29 0a 58 |bolic (n|p, op).X|
|00003640| 4e 6f 77 20 2a 6e 70 3b | 0a 58 4f 62 6a 20 2a 6f |Now *np;|.XObj *o|
|00003650| 70 3b 0a 58 7b 0a 58 09 | 64 6f 75 62 6c 65 20 6c |p;.X{.X.|double l|
|00003660| 73 6e 2c 20 72 73 6e 3b | 09 2f 2a 20 74 72 75 65 |sn, rsn;|./* true|
|00003670| 20 67 65 6f 63 20 6c 6e | 67 20 6f 66 20 73 75 6e | geoc ln|g of sun|
|00003680| 3b 20 64 69 73 74 20 66 | 72 6f 6d 20 73 6e 20 74 |; dist f|rom sn t|
|00003690| 6f 20 65 61 72 74 68 2a | 2f 0a 58 09 64 6f 75 62 |o earth*|/.X.doub|
|000036a0| 6c 65 20 6c 61 6d 3b 20 | 20 20 20 09 09 2f 2a 20 |le lam; |   ../* |
|000036b0| 67 65 6f 63 65 6e 74 72 | 69 63 20 65 63 6c 69 70 |geocentr|ic eclip|
|000036c0| 74 69 63 20 6c 6f 6e 67 | 69 74 75 64 65 20 2a 2f |tic long|itude */|
|000036d0| 0a 58 09 64 6f 75 62 6c | 65 20 62 65 74 3b 20 20 |.X.doubl|e bet;  |
|000036e0| 20 20 09 09 2f 2a 20 67 | 65 6f 63 65 6e 74 72 69 |  ../* g|eocentri|
|000036f0| 63 20 65 63 6c 69 70 74 | 69 63 20 6c 61 74 69 74 |c eclipt|ic latit|
|00003700| 75 64 65 20 2a 2f 0a 58 | 09 64 6f 75 62 6c 65 20 |ude */.X|.double |
|00003710| 6d 61 67 3b 09 09 2f 2a | 20 6d 61 67 6e 69 74 75 |mag;../*| magnitu|
|00003720| 64 65 20 2a 2f 0a 58 09 | 64 6f 75 62 6c 65 20 69 |de */.X.|double i|
|00003730| 6e 63 2c 20 6f 6d 2c 20 | 4f 6d 3b 0a 58 09 64 6f |nc, om, |Om;.X.do|
|00003740| 75 62 6c 65 20 6c 70 64 | 2c 20 70 73 69 2c 20 72 |uble lpd|, psi, r|
|00003750| 70 2c 20 72 68 6f 3b 0a | 58 09 64 6f 75 62 6c 65 |p, rho;.|X.double|
|00003760| 20 64 74 3b 0a 58 09 69 | 6e 74 20 70 61 73 73 3b | dt;.X.i|nt pass;|
|00003770| 0a 58 0a 58 09 2f 2a 20 | 66 69 6e 64 20 73 6f 6c |.X.X./* |find sol|
|00003780| 61 72 20 65 63 6c 69 70 | 74 69 63 61 6c 20 6c 6f |ar eclip|tical lo|
|00003790| 6e 67 69 74 75 64 65 20 | 61 6e 64 20 64 69 73 74 |ngitude |and dist|
|000037a0| 61 6e 63 65 20 74 6f 20 | 73 75 6e 20 66 72 6f 6d |ance to |sun from|
|000037b0| 20 65 61 72 74 68 20 2a | 2f 0a 58 09 73 75 6e 70 | earth *|/.X.sunp|
|000037c0| 6f 73 20 28 6d 6a 64 2c | 20 26 6c 73 6e 2c 20 26 |os (mjd,| &lsn, &|
|000037d0| 72 73 6e 29 3b 0a 58 0a | 58 09 2f 2a 20 74 77 6f |rsn);.X.|X./* two|
|000037e0| 20 70 61 73 73 65 73 20 | 74 6f 20 63 6f 72 72 65 | passes |to corre|
|000037f0| 63 74 20 6c 61 6d 20 61 | 6e 64 20 62 65 74 20 66 |ct lam a|nd bet f|
|00003800| 6f 72 20 6c 69 67 68 74 | 20 74 72 61 76 65 6c 20 |or light| travel |
|00003810| 74 69 6d 65 2e 65 5f 20 | 2a 2f 0a 58 09 64 74 20 |time.e_ |*/.X.dt |
|00003820| 3d 20 30 2e 30 3b 0a 58 | 09 66 6f 72 20 28 70 61 |= 0.0;.X|.for (pa|
|00003830| 73 73 20 3d 20 70 72 65 | 66 5f 67 65 74 28 50 52 |ss = pre|f_get(PR|
|00003840| 45 46 5f 41 4c 47 4f 29 | 3d 3d 50 52 45 46 5f 41 |EF_ALGO)|==PREF_A|
|00003850| 43 43 55 52 41 54 45 20 | 3f 20 30 20 3a 20 31 3b |CCURATE |? 0 : 1;|
|00003860| 20 70 61 73 73 3c 32 3b | 20 70 61 73 73 2b 2b 29 | pass<2;| pass++)|
|00003870| 7b 0a 58 09 20 20 20 20 | 72 65 64 75 63 65 5f 65 |{.X.    |reduce_e|
|00003880| 6c 65 6d 65 6e 74 73 20 | 28 6f 70 2d 3e 70 5f 65 |lements |(op->p_e|
|00003890| 70 6f 63 68 2c 20 6d 6a | 64 2d 64 74 2c 20 64 65 |poch, mj|d-dt, de|
|000038a0| 67 72 61 64 28 6f 70 2d | 3e 70 5f 69 6e 63 29 2c |grad(op-|>p_inc),|
|000038b0| 0a 58 09 09 64 65 67 72 | 61 64 28 6f 70 2d 3e 70 |.X..degr|ad(op->p|
|000038c0| 5f 6f 6d 29 2c 20 64 65 | 67 72 61 64 28 6f 70 2d |_om), de|grad(op-|
|000038d0| 3e 70 5f 4f 6d 29 2c 20 | 26 69 6e 63 2c 20 26 6f |>p_Om), |&inc, &o|
|000038e0| 6d 2c 20 26 4f 6d 29 3b | 0a 58 09 20 20 20 20 63 |m, &Om);|.X.    c|
|000038f0| 6f 6d 65 74 20 28 6d 6a | 64 2d 64 74 2c 20 6f 70 |omet (mj|d-dt, op|
|00003900| 2d 3e 70 5f 65 70 2c 20 | 69 6e 63 2c 20 6f 6d 2c |->p_ep, |inc, om,|
|00003910| 20 6f 70 2d 3e 70 5f 71 | 70 2c 20 4f 6d 2c 0a 58 | op->p_q|p, Om,.X|
|00003920| 09 09 09 09 20 20 20 20 | 26 6c 70 64 2c 20 26 70 |....    |&lpd, &p|
|00003930| 73 69 2c 20 26 72 70 2c | 20 26 72 68 6f 2c 20 26 |si, &rp,| &rho, &|
|00003940| 6c 61 6d 2c 20 26 62 65 | 74 29 3b 0a 58 09 20 20 |lam, &be|t);.X.  |
|00003950| 20 20 64 74 20 3d 20 72 | 68 6f 2a 4c 54 41 55 2f |  dt = r|ho*LTAU/|
|00003960| 33 36 30 30 2e 30 2f 32 | 34 2e 30 3b 09 2f 2a 20 |3600.0/2|4.0;./* |
|00003970| 6c 69 67 68 74 20 74 72 | 61 76 65 6c 20 74 69 6d |light tr|avel tim|
|00003980| 65 2c 20 69 6e 20 64 61 | 79 73 20 2f 20 41 55 20 |e, in da|ys / AU |
|00003990| 2a 2f 0a 58 09 7d 0a 58 | 0a 58 09 2f 2a 20 66 69 |*/.X.}.X|.X./* fi|
|000039a0| 6c 6c 20 69 6e 20 61 6c | 6c 20 6f 66 20 6f 70 2d |ll in al|l of op-|
|000039b0| 3e 73 5f 2a 20 73 74 75 | 66 66 20 65 78 63 65 70 |>s_* stu|ff excep|
|000039c0| 74 20 73 5f 73 69 7a 65 | 20 61 6e 64 20 73 5f 6d |t s_size| and s_m|
|000039d0| 61 67 20 2a 2f 0a 58 09 | 63 69 72 5f 73 6b 79 20 |ag */.X.|cir_sky |
|000039e0| 28 6e 70 2c 20 6c 70 64 | 2c 20 70 73 69 2c 20 72 |(np, lpd|, psi, r|
|000039f0| 70 2c 20 72 68 6f 2c 20 | 6c 61 6d 2c 20 62 65 74 |p, rho, |lam, bet|
|00003a00| 2c 20 6c 73 6e 2c 20 72 | 73 6e 2c 20 6f 70 29 3b |, lsn, r|sn, op);|
|00003a10| 0a 58 0a 58 09 2f 2a 20 | 63 6f 6d 70 75 74 65 20 |.X.X./* |compute |
|00003a20| 6d 61 67 6e 69 74 75 64 | 65 20 61 6e 64 20 73 69 |magnitud|e and si|
|00003a30| 7a 65 20 2a 2f 0a 58 09 | 67 6b 5f 6d 61 67 20 28 |ze */.X.|gk_mag (|
|00003a40| 6f 70 2d 3e 70 5f 67 2c | 20 6f 70 2d 3e 70 5f 6b |op->p_g,| op->p_k|
|00003a50| 2c 20 72 70 2c 20 72 68 | 6f 2c 20 26 6d 61 67 29 |, rp, rh|o, &mag)|
|00003a60| 3b 0a 58 09 6f 70 2d 3e | 73 5f 6d 61 67 20 3d 20 |;.X.op->|s_mag = |
|00003a70| 6d 61 67 20 2a 20 4d 41 | 47 53 43 41 4c 45 3b 0a |mag * MA|GSCALE;.|
|00003a80| 58 09 6f 70 2d 3e 73 5f | 73 69 7a 65 20 3d 20 6f |X.op->s_|size = o|
|00003a90| 70 2d 3e 70 5f 73 69 7a | 65 20 2f 20 72 68 6f 3b |p->p_siz|e / rho;|
|00003aa0| 0a 58 0a 58 09 72 65 74 | 75 72 6e 20 28 30 29 3b |.X.X.ret|urn (0);|
|00003ab0| 0a 58 7d 0a 58 0a 58 2f | 2a 20 66 69 6e 64 20 73 |.X}.X.X/|* find s|
|00003ac0| 75 6e 27 73 20 63 69 72 | 63 75 6d 73 74 61 6e 63 |un's cir|cumstanc|
|00003ad0| 65 73 20 6e 6f 77 2e 0a | 58 20 2a 2f 0a 58 73 74 |es now..|X */.Xst|
|00003ae0| 61 74 69 63 0a 58 73 75 | 6e 5f 63 69 72 20 28 6e |atic.Xsu|n_cir (n|
|00003af0| 70 2c 20 6f 70 29 0a 58 | 4e 6f 77 20 2a 6e 70 3b |p, op).X|Now *np;|
|00003b00| 0a 58 4f 62 6a 20 2a 6f | 70 3b 0a 58 7b 0a 58 09 |.XObj *o|p;.X{.X.|
|00003b10| 64 6f 75 62 6c 65 20 6c | 73 6e 2c 20 72 73 6e 3b |double l|sn, rsn;|
|00003b20| 09 2f 2a 20 74 72 75 65 | 20 67 65 6f 63 20 6c 6e |./* true| geoc ln|
|00003b30| 67 20 6f 66 20 73 75 6e | 3b 20 64 69 73 74 20 66 |g of sun|; dist f|
|00003b40| 72 6f 6d 20 73 6e 20 74 | 6f 20 65 61 72 74 68 2a |rom sn t|o earth*|
|00003b50| 2f 0a 58 09 64 6f 75 62 | 6c 65 20 64 65 70 73 2c |/.X.doub|le deps,|
|00003b60| 20 64 70 73 69 3b 09 2f | 2a 20 6e 75 74 61 74 69 | dpsi;./|* nutati|
|00003b70| 6f 6e 20 6f 6e 20 6f 62 | 6c 69 71 75 69 74 79 20 |on on ob|liquity |
|00003b80| 61 6e 64 20 6c 6f 6e 67 | 69 74 75 64 65 20 2a 2f |and long|itude */|
|00003b90| 0a 58 09 64 6f 75 62 6c | 65 20 6c 73 74 3b 09 09 |.X.doubl|e lst;..|
|00003ba0| 2f 2a 20 6c 6f 63 61 6c | 20 73 69 64 65 72 65 61 |/* local| siderea|
|00003bb0| 6c 20 74 69 6d 65 20 2a | 2f 0a 58 09 64 6f 75 62 |l time *|/.X.doub|
|00003bc0| 6c 65 20 65 68 70 3b 09 | 09 2f 2a 20 61 6e 67 75 |le ehp;.|./* angu|
|00003bd0| 6c 61 72 20 64 69 61 6d | 74 65 72 20 6f 66 20 65 |lar diam|ter of e|
|00003be0| 61 72 74 68 20 66 72 6f | 6d 20 6f 62 6a 65 63 74 |arth fro|m object|
|00003bf0| 20 2a 2f 0a 58 09 64 6f | 75 62 6c 65 20 68 61 3b | */.X.do|uble ha;|
|00003c00| 09 09 2f 2a 20 68 6f 75 | 72 20 61 6e 67 6c 65 20 |../* hou|r angle |
|00003c10| 2a 2f 0a 58 09 64 6f 75 | 62 6c 65 20 72 61 2c 20 |*/.X.dou|ble ra, |
|00003c20| 64 65 63 3b 09 09 2f 2a | 20 72 61 20 61 6e 64 20 |dec;../*| ra and |
|00003c30| 64 65 63 20 6e 6f 77 20 | 2a 2f 0a 58 09 64 6f 75 |dec now |*/.X.dou|
|00003c40| 62 6c 65 20 61 6c 74 2c | 20 61 7a 3b 09 09 2f 2a |ble alt,| az;../*|
|00003c50| 20 61 6c 74 20 61 6e 64 | 20 61 7a 20 2a 2f 0a 58 | alt and| az */.X|
|00003c60| 09 64 6f 75 62 6c 65 20 | 64 68 6c 6f 6e 67 3b 0a |.double |dhlong;.|
|00003c70| 58 0a 58 09 73 75 6e 70 | 6f 73 20 28 6d 6a 64 2c |X.X.sunp|os (mjd,|
|00003c80| 20 26 6c 73 6e 2c 20 26 | 72 73 6e 29 3b 09 2f 2a | &lsn, &|rsn);./*|
|00003c90| 20 73 75 6e 27 73 20 74 | 72 75 65 20 65 63 6c 69 | sun's t|rue ecli|
|00003ca0| 70 74 69 63 20 6c 6f 6e | 67 20 61 6e 64 20 64 69 |ptic lon|g and di|
|00003cb0| 73 74 20 2a 2f 0a 58 09 | 6e 75 74 61 74 69 6f 6e |st */.X.|nutation|
|00003cc0| 20 28 6d 6a 64 2c 20 26 | 64 65 70 73 2c 20 26 64 | (mjd, &|deps, &d|
|00003cd0| 70 73 69 29 3b 09 2f 2a | 20 63 6f 72 72 65 63 74 |psi);./*| correct|
|00003ce0| 20 66 6f 72 20 6e 75 74 | 61 74 69 6f 6e 20 2e 2e | for nut|ation ..|
|00003cf0| 2e 20 2a 2f 0a 58 09 6c | 73 6e 20 2b 3d 20 64 70 |. */.X.l|sn += dp|
|00003d00| 73 69 3b 0a 58 09 6c 73 | 6e 20 2d 3d 20 64 65 67 |si;.X.ls|n -= deg|
|00003d10| 72 61 64 28 32 30 2e 34 | 2f 33 36 30 30 29 3b 09 |rad(20.4|/3600);.|
|00003d20| 2f 2a 20 61 6e 64 20 61 | 62 65 72 61 74 69 6f 6e |/* and a|beration|
|00003d30| 20 2a 2f 0a 58 0a 58 09 | 6f 70 2d 3e 73 5f 65 64 | */.X.X.|op->s_ed|
|00003d40| 69 73 74 20 3d 20 72 73 | 6e 3b 0a 58 09 6f 70 2d |ist = rs|n;.X.op-|
|00003d50| 3e 73 5f 73 64 69 73 74 | 20 3d 20 30 2e 30 3b 0a |>s_sdist| = 0.0;.|
|00003d60| 58 09 6f 70 2d 3e 73 5f | 65 6c 6f 6e 67 20 3d 20 |X.op->s_|elong = |
|00003d70| 30 2e 30 3b 0a 58 09 6f | 70 2d 3e 73 5f 70 68 61 |0.0;.X.o|p->s_pha|
|00003d80| 73 65 20 3d 20 31 30 30 | 2e 30 3b 0a 58 09 6f 70 |se = 100|.0;.X.op|
|00003d90| 2d 3e 73 5f 73 69 7a 65 | 20 3d 20 72 61 64 64 65 |->s_size| = radde|
|00003da0| 67 28 34 2e 36 35 32 34 | 32 65 2d 33 2f 72 73 6e |g(4.6524|2e-3/rsn|
|00003db0| 29 2a 33 36 30 30 2a 32 | 3b 0a 58 09 6f 70 2d 3e |)*3600*2|;.X.op->|
|00003dc0| 73 5f 6d 61 67 20 3d 20 | 4d 41 47 53 43 41 4c 45 |s_mag = |MAGSCALE|
|00003dd0| 20 2a 20 2d 32 36 2e 38 | 3b 09 2f 2a 20 54 4f 44 | * -26.8|;./* TOD|
|00003de0| 4f 20 2a 2f 0a 58 09 64 | 68 6c 6f 6e 67 20 3d 20 |O */.X.d|hlong = |
|00003df0| 6c 73 6e 2d 50 49 3b 09 | 2f 2a 20 67 65 6f 2d 20 |lsn-PI;.|/* geo- |
|00003e00| 74 6f 20 68 65 6c 69 6f | 2d 20 63 65 6e 74 72 69 |to helio|- centri|
|00003e10| 63 20 2a 2f 0a 58 09 72 | 61 6e 67 65 20 28 26 64 |c */.X.r|ange (&d|
|00003e20| 68 6c 6f 6e 67 2c 20 32 | 2a 50 49 29 3b 0a 58 09 |hlong, 2|*PI);.X.|
|00003e30| 6f 70 2d 3e 73 5f 68 6c | 6f 6e 67 20 3d 20 64 68 |op->s_hl|ong = dh|
|00003e40| 6c 6f 6e 67 3b 0a 58 09 | 6f 70 2d 3e 73 5f 68 6c |long;.X.|op->s_hl|
|00003e50| 61 74 20 3d 20 30 2e 30 | 3b 0a 58 0a 58 0a 58 09 |at = 0.0|;.X.X.X.|
|00003e60| 2f 2a 20 63 6f 6e 76 65 | 72 74 20 67 65 6f 63 65 |/* conve|rt geoce|
|00003e70| 6e 74 72 69 63 20 65 63 | 6c 69 70 74 69 63 20 6c |ntric ec|liptic l|
|00003e80| 61 74 2f 6c 6f 6e 67 20 | 74 6f 20 65 71 75 69 74 |at/long |to equit|
|00003e90| 6f 72 69 61 6c 20 72 61 | 2f 64 65 63 20 6f 66 20 |orial ra|/dec of |
|00003ea0| 64 61 74 65 20 2a 2f 0a | 58 09 65 63 6c 5f 65 71 |date */.|X.ecl_eq|
|00003eb0| 20 28 6d 6a 64 2c 20 30 | 2e 30 2c 20 6c 73 6e 2c | (mjd, 0|.0, lsn,|
|00003ec0| 20 26 72 61 2c 20 26 64 | 65 63 29 3b 0a 58 0a 58 | &ra, &d|ec);.X.X|
|00003ed0| 09 2f 2a 20 66 69 6e 64 | 20 73 5f 72 61 2f 64 65 |./* find| s_ra/de|
|00003ee0| 63 20 61 74 20 64 65 73 | 69 72 65 64 20 65 70 6f |c at des|ired epo|
|00003ef0| 63 68 20 2a 2f 0a 58 09 | 69 66 20 28 65 70 6f 63 |ch */.X.|if (epoc|
|00003f00| 68 20 3d 3d 20 45 4f 44 | 29 20 7b 0a 58 09 20 20 |h == EOD|) {.X.  |
|00003f10| 20 20 6f 70 2d 3e 73 5f | 72 61 20 3d 20 72 61 3b |  op->s_|ra = ra;|
|00003f20| 0a 58 09 20 20 20 20 6f | 70 2d 3e 73 5f 64 65 63 |.X.    o|p->s_dec|
|00003f30| 20 3d 20 64 65 63 3b 0a | 58 09 7d 20 65 6c 73 65 | = dec;.|X.} else|
|00003f40| 20 7b 0a 58 09 20 20 20 | 20 64 6f 75 62 6c 65 20 | {.X.   | double |
|00003f50| 74 72 61 20 3d 20 72 61 | 2c 20 74 64 65 63 20 3d |tra = ra|, tdec =|
|00003f60| 20 64 65 63 3b 0a 58 09 | 20 20 20 20 70 72 65 63 | dec;.X.|    prec|
|00003f70| 65 73 73 20 28 6d 6a 64 | 2c 20 65 70 6f 63 68 2c |ess (mjd|, epoch,|
|00003f80| 20 26 74 72 61 2c 20 26 | 74 64 65 63 29 3b 0a 58 | &tra, &|tdec);.X|
|00003f90| 09 20 20 20 20 6f 70 2d | 3e 73 5f 72 61 20 3d 20 |.    op-|>s_ra = |
|00003fa0| 74 72 61 3b 0a 58 09 20 | 20 20 20 6f 70 2d 3e 73 |tra;.X. |   op->s|
|00003fb0| 5f 64 65 63 20 3d 20 74 | 64 65 63 3b 0a 58 09 7d |_dec = t|dec;.X.}|
|00003fc0| 0a 58 0a 58 09 2f 2a 20 | 66 69 6e 64 20 61 6c 74 |.X.X./* |find alt|
|00003fd0| 2f 61 7a 20 62 61 73 65 | 64 20 6f 6e 20 75 6e 70 |/az base|d on unp|
|00003fe0| 72 65 63 65 73 73 65 64 | 20 72 61 2f 64 65 63 20 |recessed| ra/dec |
|00003ff0| 2a 2f 0a 58 09 6e 6f 77 | 5f 6c 73 74 20 28 6e 70 |*/.X.now|_lst (np|
|00004000| 2c 20 26 6c 73 74 29 3b | 0a 58 09 68 61 20 3d 20 |, &lst);|.X.ha = |
|00004010| 68 72 72 61 64 28 6c 73 | 74 29 20 2d 20 72 61 3b |hrrad(ls|t) - ra;|
|00004020| 0a 58 09 65 68 70 20 3d | 20 28 32 2e 30 20 2a 20 |.X.ehp =| (2.0 * |
|00004030| 36 33 37 38 2e 30 20 2f | 20 31 34 36 2e 30 65 36 |6378.0 /| 146.0e6|
|00004040| 29 20 2f 20 6f 70 2d 3e | 73 5f 65 64 69 73 74 3b |) / op->|s_edist;|
|00004050| 0a 58 09 74 61 5f 70 61 | 72 20 28 68 61 2c 20 64 |.X.ta_pa|r (ha, d|
|00004060| 65 63 2c 20 6c 61 74 2c | 20 65 6c 65 76 2c 20 65 |ec, lat,| elev, e|
|00004070| 68 70 2c 20 26 68 61 2c | 20 26 64 65 63 29 3b 0a |hp, &ha,| &dec);.|
|00004080| 58 09 68 61 64 65 63 5f | 61 61 20 28 6c 61 74 2c |X.hadec_|aa (lat,|
|00004090| 20 68 61 2c 20 64 65 63 | 2c 20 26 61 6c 74 2c 20 | ha, dec|, &alt, |
|000040a0| 26 61 7a 29 3b 0a 58 09 | 72 65 66 72 61 63 74 20 |&az);.X.|refract |
|000040b0| 28 70 72 65 73 73 75 72 | 65 2c 20 74 65 6d 70 2c |(pressur|e, temp,|
|000040c0| 20 61 6c 74 2c 20 26 61 | 6c 74 29 3b 0a 58 09 6f | alt, &a|lt);.X.o|
|000040d0| 70 2d 3e 73 5f 61 6c 74 | 20 3d 20 61 6c 74 3b 0a |p->s_alt| = alt;.|
|000040e0| 58 09 6f 70 2d 3e 73 5f | 61 7a 20 3d 20 61 7a 3b |X.op->s_|az = az;|
|000040f0| 0a 58 09 72 65 74 75 72 | 6e 20 28 30 29 3b 0a 58 |.X.retur|n (0);.X|
|00004100| 7d 0a 58 0a 58 2f 2a 20 | 66 69 6e 64 20 6d 6f 6f |}.X.X/* |find moo|
|00004110| 6e 27 73 20 63 69 72 63 | 75 6d 73 74 61 6e 63 65 |n's circ|umstance|
|00004120| 73 20 6e 6f 77 2e 0a 58 | 20 2a 2f 0a 58 73 74 61 |s now..X| */.Xsta|
|00004130| 74 69 63 0a 58 6d 6f 6f | 6e 5f 63 69 72 20 28 6e |tic.Xmoo|n_cir (n|
|00004140| 70 2c 20 6f 70 29 0a 58 | 4e 6f 77 20 2a 6e 70 3b |p, op).X|Now *np;|
|00004150| 0a 58 4f 62 6a 20 2a 6f | 70 3b 0a 58 7b 0a 58 09 |.XObj *o|p;.X{.X.|
|00004160| 64 6f 75 62 6c 65 20 6c | 73 6e 2c 20 72 73 6e 3b |double l|sn, rsn;|
|00004170| 09 2f 2a 20 74 72 75 65 | 20 67 65 6f 63 20 6c 6e |./* true| geoc ln|
|00004180| 67 20 6f 66 20 73 75 6e | 3b 20 64 69 73 74 20 66 |g of sun|; dist f|
|00004190| 72 6f 6d 20 73 6e 20 74 | 6f 20 65 61 72 74 68 2a |rom sn t|o earth*|
|000041a0| 2f 0a 58 09 64 6f 75 62 | 6c 65 20 64 65 70 73 2c |/.X.doub|le deps,|
|000041b0| 20 64 70 73 69 3b 09 2f | 2a 20 6e 75 74 61 74 69 | dpsi;./|* nutati|
|000041c0| 6f 6e 20 6f 6e 20 6f 62 | 6c 69 71 75 69 74 79 20 |on on ob|liquity |
|000041d0| 61 6e 64 20 6c 6f 6e 67 | 69 74 75 64 65 20 2a 2f |and long|itude */|
|000041e0| 0a 58 09 64 6f 75 62 6c | 65 20 6c 73 74 3b 09 09 |.X.doubl|e lst;..|
|000041f0| 2f 2a 20 6c 6f 63 61 6c | 20 73 69 64 65 72 65 61 |/* local| siderea|
|00004200| 6c 20 74 69 6d 65 20 2a | 2f 0a 58 09 64 6f 75 62 |l time *|/.X.doub|
|00004210| 6c 65 20 65 68 70 3b 09 | 09 2f 2a 20 61 6e 67 75 |le ehp;.|./* angu|
|00004220| 6c 61 72 20 64 69 61 6d | 74 65 72 20 6f 66 20 65 |lar diam|ter of e|
|00004230| 61 72 74 68 20 66 72 6f | 6d 20 6f 62 6a 65 63 74 |arth fro|m object|
|00004240| 20 2a 2f 0a 58 09 64 6f | 75 62 6c 65 20 68 61 3b | */.X.do|uble ha;|
|00004250| 09 09 2f 2a 20 68 6f 75 | 72 20 61 6e 67 6c 65 20 |../* hou|r angle |
|00004260| 2a 2f 0a 58 09 64 6f 75 | 62 6c 65 20 72 61 2c 20 |*/.X.dou|ble ra, |
|00004270| 64 65 63 3b 09 09 2f 2a | 20 72 61 20 61 6e 64 20 |dec;../*| ra and |
|00004280| 64 65 63 20 6e 6f 77 20 | 2a 2f 0a 58 09 64 6f 75 |dec now |*/.X.dou|
|00004290| 62 6c 65 20 61 6c 74 2c | 20 61 7a 3b 09 09 2f 2a |ble alt,| az;../*|
|000042a0| 20 61 6c 74 20 61 6e 64 | 20 61 7a 20 2a 2f 0a 58 | alt and| az */.X|
|000042b0| 09 64 6f 75 62 6c 65 20 | 6c 61 6d 3b 20 20 20 20 |.double |lam;    |
|000042c0| 09 09 2f 2a 20 67 65 6f | 63 65 6e 74 72 69 63 20 |../* geo|centric |
|000042d0| 65 63 6c 69 70 74 69 63 | 20 6c 6f 6e 67 69 74 75 |ecliptic| longitu|
|000042e0| 64 65 20 2a 2f 0a 58 09 | 64 6f 75 62 6c 65 20 62 |de */.X.|double b|
|000042f0| 65 74 3b 20 20 20 20 09 | 09 2f 2a 20 67 65 6f 63 |et;    .|./* geoc|
|00004300| 65 6e 74 72 69 63 20 65 | 63 6c 69 70 74 69 63 20 |entric e|cliptic |
|00004310| 6c 61 74 69 74 75 64 65 | 20 2a 2f 0a 58 09 64 6f |latitude| */.X.do|
|00004320| 75 62 6c 65 20 65 64 69 | 73 74 61 75 3b 09 09 2f |uble edi|stau;../|
|00004330| 2a 20 65 61 72 74 68 2d | 6d 6f 6f 6e 20 64 69 73 |* earth-|moon dis|
|00004340| 74 2c 20 69 6e 20 61 75 | 20 2a 2f 0a 58 09 64 6f |t, in au| */.X.do|
|00004350| 75 62 6c 65 20 65 6c 3b | 09 09 2f 2a 20 65 6c 6f |uble el;|../* elo|
|00004360| 6e 67 61 74 69 6f 6e 2c | 20 72 61 64 73 20 65 61 |ngation,| rads ea|
|00004370| 73 74 20 2a 2f 0a 58 0a | 58 09 6d 6f 6f 6e 20 28 |st */.X.|X.moon (|
|00004380| 6d 6a 64 2c 20 26 6c 61 | 6d 2c 20 26 62 65 74 2c |mjd, &la|m, &bet,|
|00004390| 20 26 65 68 70 29 3b 09 | 2f 2a 20 6d 6f 6f 6e 27 | &ehp);.|/* moon'|
|000043a0| 73 20 74 72 75 65 20 65 | 63 6c 69 70 74 69 63 20 |s true e|cliptic |
|000043b0| 6c 6f 63 20 2a 2f 0a 58 | 09 6e 75 74 61 74 69 6f |loc */.X|.nutatio|
|000043c0| 6e 20 28 6d 6a 64 2c 20 | 26 64 65 70 73 2c 20 26 |n (mjd, |&deps, &|
|000043d0| 64 70 73 69 29 3b 09 2f | 2a 20 63 6f 72 72 65 63 |dpsi);./|* correc|
|000043e0| 74 20 66 6f 72 20 6e 75 | 74 61 74 69 6f 6e 20 2a |t for nu|tation *|
|000043f0| 2f 0a 58 09 6c 61 6d 20 | 2b 3d 20 64 70 73 69 3b |/.X.lam |+= dpsi;|
|00004400| 0a 58 09 72 61 6e 67 65 | 20 28 26 6c 61 6d 2c 20 |.X.range| (&lam, |
|00004410| 32 2a 50 49 29 3b 0a 58 | 0a 58 09 6f 70 2d 3e 73 |2*PI);.X|.X.op->s|
|00004420| 5f 65 64 69 73 74 20 3d | 20 36 33 37 38 2e 31 34 |_edist =| 6378.14|
|00004430| 2f 73 69 6e 28 65 68 70 | 29 3b 09 2f 2a 20 65 61 |/sin(ehp|);./* ea|
|00004440| 72 74 68 2d 6d 6f 6f 6e | 20 64 69 73 74 2c 20 77 |rth-moon| dist, w|
|00004450| 61 6e 74 20 6b 6d 20 2a | 2f 0a 58 09 6f 70 2d 3e |ant km *|/.X.op->|
|00004460| 73 5f 73 69 7a 65 20 3d | 20 33 36 30 30 2a 33 31 |s_size =| 3600*31|
|00004470| 2e 32 32 35 31 32 2a 73 | 69 6e 28 65 68 70 29 3b |.22512*s|in(ehp);|
|00004480| 2f 2a 20 6d 6f 6f 6e 20 | 61 6e 67 75 6c 61 72 20 |/* moon |angular |
|00004490| 64 69 61 2c 20 73 65 63 | 6f 6e 64 73 20 2a 2f 0a |dia, sec|onds */.|
|000044a0| 58 09 6f 70 2d 3e 73 5f | 68 6c 6f 6e 67 20 3d 20 |X.op->s_|hlong = |
|000044b0| 6c 61 6d 3b 09 09 2f 2a | 20 73 61 76 65 20 67 65 |lam;../*| save ge|
|000044c0| 6f 20 69 6e 20 68 65 6c | 69 6f 20 66 69 65 6c 64 |o in hel|io field|
|000044d0| 73 20 2a 2f 0a 58 09 6f | 70 2d 3e 73 5f 68 6c 61 |s */.X.o|p->s_hla|
|000044e0| 74 20 3d 20 62 65 74 3b | 0a 58 0a 58 09 2f 2a 20 |t = bet;|.X.X./* |
|000044f0| 63 6f 6e 76 65 72 74 20 | 67 65 6f 63 65 6e 74 72 |convert |geocentr|
|00004500| 69 63 20 65 63 6c 69 70 | 74 69 63 20 6c 61 74 2f |ic eclip|tic lat/|
|00004510| 6c 6f 6e 67 20 74 6f 20 | 65 71 75 69 74 6f 72 69 |long to |equitori|
|00004520| 61 6c 20 72 61 2f 64 65 | 63 20 6f 66 20 64 61 74 |al ra/de|c of dat|
|00004530| 65 20 2a 2f 0a 58 09 65 | 63 6c 5f 65 71 20 28 6d |e */.X.e|cl_eq (m|
|00004540| 6a 64 2c 20 62 65 74 2c | 20 6c 61 6d 2c 20 26 72 |jd, bet,| lam, &r|
|00004550| 61 2c 20 26 64 65 63 29 | 3b 0a 58 0a 58 09 2f 2a |a, &dec)|;.X.X./*|
|00004560| 20 66 69 6e 64 20 73 5f | 72 61 2f 64 65 63 20 61 | find s_|ra/dec a|
|00004570| 74 20 64 65 73 69 72 65 | 64 20 65 70 6f 63 68 20 |t desire|d epoch |
|00004580| 2a 2f 0a 58 09 69 66 20 | 28 65 70 6f 63 68 20 3d |*/.X.if |(epoch =|
|00004590| 3d 20 45 4f 44 29 20 7b | 0a 58 09 20 20 20 20 6f |= EOD) {|.X.    o|
|000045a0| 70 2d 3e 73 5f 72 61 20 | 3d 20 72 61 3b 0a 58 09 |p->s_ra |= ra;.X.|
|000045b0| 20 20 20 20 6f 70 2d 3e | 73 5f 64 65 63 20 3d 20 |    op->|s_dec = |
|000045c0| 64 65 63 3b 0a 58 09 7d | 20 65 6c 73 65 20 7b 0a |dec;.X.}| else {.|
|000045d0| 58 09 20 20 20 20 64 6f | 75 62 6c 65 20 74 72 61 |X.    do|uble tra|
|000045e0| 20 3d 20 72 61 2c 20 74 | 64 65 63 20 3d 20 64 65 | = ra, t|dec = de|
|000045f0| 63 3b 0a 58 09 20 20 20 | 20 70 72 65 63 65 73 73 |c;.X.   | precess|
|00004600| 20 28 6d 6a 64 2c 20 65 | 70 6f 63 68 2c 20 26 74 | (mjd, e|poch, &t|
|00004610| 72 61 2c 20 26 74 64 65 | 63 29 3b 0a 58 09 20 20 |ra, &tde|c);.X.  |
|00004620| 20 20 6f 70 2d 3e 73 5f | 72 61 20 3d 20 74 72 61 |  op->s_|ra = tra|
|00004630| 3b 0a 58 09 20 20 20 20 | 6f 70 2d 3e 73 5f 64 65 |;.X.    |op->s_de|
|00004640| 63 20 3d 20 74 64 65 63 | 3b 0a 58 09 7d 0a 58 0a |c = tdec|;.X.}.X.|
|00004650| 58 09 73 75 6e 70 6f 73 | 20 28 6d 6a 64 2c 20 26 |X.sunpos| (mjd, &|
|00004660| 6c 73 6e 2c 20 26 72 73 | 6e 29 3b 0a 58 09 72 61 |lsn, &rs|n);.X.ra|
|00004670| 6e 67 65 20 28 26 6c 73 | 6e 2c 20 32 2a 50 49 29 |nge (&ls|n, 2*PI)|
|00004680| 3b 0a 58 09 65 6c 6f 6e | 67 61 74 69 6f 6e 20 28 |;.X.elon|gation (|
|00004690| 6c 61 6d 2c 20 62 65 74 | 2c 20 6c 73 6e 2c 20 26 |lam, bet|, lsn, &|
|000046a0| 65 6c 29 3b 0a 58 0a 58 | 09 2f 2a 20 73 6f 6c 76 |el);.X.X|./* solv|
|000046b0| 65 20 74 72 69 61 6e 67 | 6c 65 20 6f 66 20 65 61 |e triang|le of ea|
|000046c0| 72 74 68 2c 20 73 75 6e | 2c 20 61 6e 64 20 65 6c |rth, sun|, and el|
|000046d0| 6f 6e 67 61 74 69 6f 6e | 20 66 6f 72 20 6d 6f 6f |ongation| for moo|
|000046e0| 6e 2d 73 75 6e 20 64 69 | 73 74 20 2a 2f 0a 58 09 |n-sun di|st */.X.|
|000046f0| 65 64 69 73 74 61 75 20 | 3d 20 6f 70 2d 3e 73 5f |edistau |= op->s_|
|00004700| 65 64 69 73 74 2f 31 2e | 34 39 35 39 37 39 65 38 |edist/1.|495979e8|
|00004710| 3b 09 2f 2a 20 6b 6d 20 | 2d 3e 20 61 75 20 2a 2f |;./* km |-> au */|
|00004720| 0a 58 09 6f 70 2d 3e 73 | 5f 73 64 69 73 74 3d 20 |.X.op->s|_sdist= |
|00004730| 73 71 72 74 20 28 65 64 | 69 73 74 61 75 2a 65 64 |sqrt (ed|istau*ed|
|00004740| 69 73 74 61 75 20 2b 20 | 72 73 6e 2a 72 73 6e 20 |istau + |rsn*rsn |
|00004750| 2d 20 32 2e 30 2a 65 64 | 69 73 74 61 75 2a 72 73 |- 2.0*ed|istau*rs|
|00004760| 6e 2a 63 6f 73 28 65 6c | 29 29 3b 0a 58 0a 58 09 |n*cos(el|));.X.X.|
|00004770| 2f 2a 20 54 4f 44 4f 3a | 20 69 6d 70 72 6f 76 65 |/* TODO:| improve|
|00004780| 20 6d 61 67 3b 20 74 68 | 69 73 20 69 73 20 62 61 | mag; th|is is ba|
|00004790| 73 65 64 20 6f 6e 20 61 | 20 66 6c 61 74 20 6d 6f |sed on a| flat mo|
|000047a0| 6f 6e 20 6d 6f 64 65 6c | 2e 20 2a 2f 0a 58 09 6f |on model|. */.X.o|
|000047b0| 70 2d 3e 73 5f 6d 61 67 | 20 3d 20 4d 41 47 53 43 |p->s_mag| = MAGSC|
|000047c0| 41 4c 45 20 2a 20 0a 58 | 09 20 20 20 20 28 2d 31 |ALE * .X|.    (-1|
|000047d0| 32 2e 37 20 2b 20 32 2e | 35 2a 28 6c 6f 67 31 30 |2.7 + 2.|5*(log10|
|000047e0| 28 50 49 29 20 2d 20 6c | 6f 67 31 30 28 50 49 2f |(PI) - l|og10(PI/|
|000047f0| 32 2a 28 31 2b 31 2e 65 | 2d 36 2d 63 6f 73 28 65 |2*(1+1.e|-6-cos(e|
|00004800| 6c 29 29 29 29 29 3b 0a | 58 0a 58 09 6f 70 2d 3e |l)))));.|X.X.op->|
|00004810| 73 5f 65 6c 6f 6e 67 20 | 3d 20 72 61 64 64 65 67 |s_elong |= raddeg|
|00004820| 28 65 6c 29 3b 09 09 2f | 2a 20 77 61 6e 74 20 64 |(el);../|* want d|
|00004830| 65 67 72 65 65 73 20 2a | 2f 0a 58 09 6f 70 2d 3e |egrees *|/.X.op->|
|00004840| 73 5f 70 68 61 73 65 20 | 3d 20 66 61 62 73 28 65 |s_phase |= fabs(e|
|00004850| 6c 29 2f 50 49 2a 31 30 | 30 2e 30 3b 09 2f 2a 20 |l)/PI*10|0.0;./* |
|00004860| 77 61 6e 74 20 6e 6f 6e | 2d 6e 65 67 61 74 69 76 |want non|-negativ|
|00004870| 65 20 25 20 2a 2f 0a 58 | 0a 58 09 2f 2a 20 73 68 |e % */.X|.X./* sh|
|00004880| 6f 77 20 74 6f 70 6f 63 | 65 6e 74 72 69 63 20 61 |ow topoc|entric a|
|00004890| 6c 74 2f 61 7a 20 62 79 | 20 63 6f 72 72 65 63 74 |lt/az by| correct|
|000048a0| 69 6e 67 20 72 61 2f 64 | 65 63 20 66 6f 72 20 70 |ing ra/d|ec for p|
|000048b0| 61 72 61 6c 6c 61 78 20 | 0a 58 09 20 2a 20 61 73 |arallax |.X. * as|
|000048c0| 20 77 65 6c 6c 20 61 73 | 20 72 65 66 72 61 63 74 | well as| refract|
|000048d0| 69 6f 6e 2e 0a 58 09 20 | 2a 20 75 73 65 20 75 6e |ion..X. |* use un|
|000048e0| 70 72 65 63 65 73 73 65 | 64 20 72 61 2f 64 65 63 |precesse|d ra/dec|
|000048f0| 2e 0a 58 09 20 2a 2f 0a | 58 09 6e 6f 77 5f 6c 73 |..X. */.|X.now_ls|
|00004900| 74 20 28 6e 70 2c 20 26 | 6c 73 74 29 3b 0a 58 09 |t (np, &|lst);.X.|
|00004910| 68 61 20 3d 20 68 72 72 | 61 64 28 6c 73 74 29 20 |ha = hrr|ad(lst) |
|00004920| 2d 20 72 61 3b 0a 58 09 | 74 61 5f 70 61 72 20 28 |- ra;.X.|ta_par (|
|00004930| 68 61 2c 20 64 65 63 2c | 20 6c 61 74 2c 20 65 6c |ha, dec,| lat, el|
|00004940| 65 76 2c 20 65 68 70 2c | 20 26 68 61 2c 20 26 64 |ev, ehp,| &ha, &d|
|00004950| 65 63 29 3b 0a 58 09 68 | 61 64 65 63 5f 61 61 20 |ec);.X.h|adec_aa |
|00004960| 28 6c 61 74 2c 20 68 61 | 2c 20 64 65 63 2c 20 26 |(lat, ha|, dec, &|
|00004970| 61 6c 74 2c 20 26 61 7a | 29 3b 0a 58 09 72 65 66 |alt, &az|);.X.ref|
|00004980| 72 61 63 74 20 28 70 72 | 65 73 73 75 72 65 2c 20 |ract (pr|essure, |
|00004990| 74 65 6d 70 2c 20 61 6c | 74 2c 20 26 61 6c 74 29 |temp, al|t, &alt)|
|000049a0| 3b 0a 58 09 6f 70 2d 3e | 73 5f 61 6c 74 20 3d 20 |;.X.op->|s_alt = |
|000049b0| 61 6c 74 3b 0a 58 09 6f | 70 2d 3e 73 5f 61 7a 20 |alt;.X.o|p->s_az |
|000049c0| 3d 20 61 7a 3b 0a 58 0a | 58 09 72 65 74 75 72 6e |= az;.X.|X.return|
|000049d0| 20 28 30 29 3b 0a 58 7d | 0a 58 0a 58 2f 2a 20 66 | (0);.X}|.X.X/* f|
|000049e0| 69 6c 6c 20 69 6e 20 61 | 6c 6c 20 6f 66 20 6f 70 |ill in a|ll of op|
|000049f0| 2d 3e 73 5f 2a 20 73 74 | 75 66 66 20 65 78 63 65 |->s_* st|uff exce|
|00004a00| 70 74 20 73 5f 73 69 7a | 65 20 61 6e 64 20 73 5f |pt s_siz|e and s_|
|00004a10| 6d 61 67 20 2a 2f 0a 58 | 73 74 61 74 69 63 20 76 |mag */.X|static v|
|00004a20| 6f 69 64 0a 58 63 69 72 | 5f 73 6b 79 20 28 6e 70 |oid.Xcir|_sky (np|
|00004a30| 2c 20 6c 70 64 2c 20 70 | 73 69 2c 20 72 70 2c 20 |, lpd, p|si, rp, |
|00004a40| 72 68 6f 2c 20 6c 61 6d | 2c 20 62 65 74 2c 20 6c |rho, lam|, bet, l|
|00004a50| 73 6e 2c 20 72 73 6e 2c | 20 6f 70 29 0a 58 4e 6f |sn, rsn,| op).XNo|
|00004a60| 77 20 2a 6e 70 3b 0a 58 | 64 6f 75 62 6c 65 20 6c |w *np;.X|double l|
|00004a70| 70 64 2c 20 70 73 69 3b | 09 2f 2a 20 68 65 6c 69 |pd, psi;|./* heli|
|00004a80| 6f 63 65 6e 74 72 69 63 | 20 65 63 6c 69 70 74 69 |ocentric| eclipti|
|00004a90| 63 20 6c 6f 6e 67 20 61 | 6e 64 20 6c 61 74 20 2a |c long a|nd lat *|
|00004aa0| 2f 0a 58 64 6f 75 62 6c | 65 20 72 70 3b 09 09 2f |/.Xdoubl|e rp;../|
|00004ab0| 2a 20 64 69 73 74 20 66 | 72 6f 6d 20 73 75 6e 20 |* dist f|rom sun |
|00004ac0| 2a 2f 0a 58 64 6f 75 62 | 6c 65 20 72 68 6f 3b 09 |*/.Xdoub|le rho;.|
|00004ad0| 09 2f 2a 20 64 69 73 74 | 20 66 72 6f 6d 20 65 61 |./* dist| from ea|
|00004ae0| 72 74 68 20 2a 2f 0a 58 | 64 6f 75 62 6c 65 20 6c |rth */.X|double l|
|00004af0| 61 6d 2c 20 62 65 74 3b | 09 2f 2a 20 74 72 75 65 |am, bet;|./* true|
|00004b00| 20 67 65 6f 63 65 6e 74 | 72 69 63 20 65 63 6c 69 | geocent|ric ecli|
|00004b10| 70 74 69 63 20 6c 6f 6e | 67 20 61 6e 64 20 6c 61 |ptic lon|g and la|
|00004b20| 74 20 2a 2f 0a 58 64 6f | 75 62 6c 65 20 6c 73 6e |t */.Xdo|uble lsn|
|00004b30| 2c 20 72 73 6e 3b 09 2f | 2a 20 74 72 75 65 20 67 |, rsn;./|* true g|
|00004b40| 65 6f 63 20 6c 6e 67 20 | 6f 66 20 73 75 6e 3b 20 |eoc lng |of sun; |
|00004b50| 64 69 73 74 20 66 72 6f | 6d 20 73 6e 20 74 6f 20 |dist fro|m sn to |
|00004b60| 65 61 72 74 68 2a 2f 0a | 58 4f 62 6a 20 2a 6f 70 |earth*/.|XObj *op|
|00004b70| 3b 0a 58 7b 0a 58 09 64 | 6f 75 62 6c 65 20 61 3b |;.X{.X.d|ouble a;|
|00004b80| 09 09 2f 2a 20 67 65 6f | 63 20 61 6e 67 6c 65 20 |../* geo|c angle |
|00004b90| 62 65 74 77 65 65 6e 20 | 6f 62 6a 65 63 74 20 61 |between |object a|
|00004ba0| 6e 64 20 73 75 6e 20 2a | 2f 0a 58 09 64 6f 75 62 |nd sun *|/.X.doub|
|00004bb0| 6c 65 20 72 61 2c 20 64 | 65 63 3b 09 09 2f 2a 20 |le ra, d|ec;../* |
|00004bc0| 72 61 20 61 6e 64 20 64 | 65 63 20 61 74 20 65 70 |ra and d|ec at ep|
|00004bd0| 6f 63 68 20 6f 66 20 64 | 61 74 65 20 2a 2f 0a 58 |och of d|ate */.X|
|00004be0| 09 64 6f 75 62 6c 65 20 | 65 6c 3b 09 09 2f 2a 20 |.double |el;../* |
|00004bf0| 65 6c 6f 6e 67 61 74 69 | 6f 6e 20 2a 2f 0a 58 09 |elongati|on */.X.|
|00004c00| 64 6f 75 62 6c 65 20 6c | 73 74 3b 09 09 2f 2a 20 |double l|st;../* |
|00004c10| 6c 6f 63 61 6c 20 73 69 | 64 65 72 65 61 6c 20 74 |local si|dereal t|
|00004c20| 69 6d 65 20 2a 2f 0a 58 | 09 64 6f 75 62 6c 65 20 |ime */.X|.double |
|00004c30| 65 68 70 3b 09 09 2f 2a | 20 65 68 70 3a 20 61 6e |ehp;../*| ehp: an|
|00004c40| 67 75 6c 61 72 20 64 69 | 61 20 6f 66 20 65 61 72 |gular di|a of ear|
|00004c50| 74 68 20 66 72 6f 6d 20 | 62 6f 64 79 20 2a 2f 0a |th from |body */.|
|00004c60| 58 09 64 6f 75 62 6c 65 | 20 68 61 3b 09 09 2f 2a |X.double| ha;../*|
|00004c70| 20 6c 6f 63 61 6c 20 68 | 6f 75 72 20 61 6e 67 6c | local h|our angl|
|00004c80| 65 20 2a 2f 0a 58 09 64 | 6f 75 62 6c 65 20 66 3b |e */.X.d|ouble f;|
|00004c90| 09 09 2f 2a 20 66 72 61 | 63 74 69 6f 6e 61 6c 20 |../* fra|ctional |
|00004ca0| 70 68 61 73 65 20 66 72 | 6f 6d 20 65 61 72 74 68 |phase fr|om earth|
|00004cb0| 20 2a 2f 0a 58 09 64 6f | 75 62 6c 65 20 61 6c 74 | */.X.do|uble alt|
|00004cc0| 2c 20 61 7a 3b 09 09 2f | 2a 20 63 75 72 72 65 6e |, az;../|* curren|
|00004cd0| 74 20 61 6c 74 2c 20 61 | 7a 20 2a 2f 0a 58 09 64 |t alt, a|z */.X.d|
|00004ce0| 6f 75 62 6c 65 20 64 65 | 70 73 2c 20 64 70 73 69 |ouble de|ps, dpsi|
|00004cf0| 3b 09 2f 2a 20 6e 75 74 | 61 74 69 6f 6e 20 6f 6e |;./* nut|ation on|
|00004d00| 20 6f 62 6c 69 71 75 69 | 74 79 20 61 6e 64 20 6c | obliqui|ty and l|
|00004d10| 6f 6e 67 69 74 75 64 65 | 20 2a 2f 0a 58 0a 58 09 |ongitude| */.X.X.|
|00004d20| 2f 2a 20 63 6f 72 72 65 | 63 74 20 66 6f 72 20 6e |/* corre|ct for n|
|00004d30| 75 74 61 74 69 6f 6e 20 | 2e 2e 2e 20 2a 2f 0a 58 |utation |... */.X|
|00004d40| 09 6e 75 74 61 74 69 6f | 6e 20 28 6d 6a 64 2c 20 |.nutatio|n (mjd, |
|00004d50| 26 64 65 70 73 2c 20 26 | 64 70 73 69 29 3b 0a 58 |&deps, &|dpsi);.X|
|00004d60| 09 6c 61 6d 20 2b 3d 20 | 64 70 73 69 3b 0a 58 0a |.lam += |dpsi;.X.|
|00004d70| 58 09 2f 2a 20 61 6e 64 | 20 63 6f 72 72 65 63 74 |X./* and| correct|
|00004d80| 20 66 6f 72 20 32 30 2e | 34 22 20 61 62 65 72 61 | for 20.|4" abera|
|00004d90| 74 69 6f 6e 20 2a 2f 0a | 58 09 61 20 3d 20 6c 73 |tion */.|X.a = ls|
|00004da0| 6e 20 2d 20 6c 61 6d 3b | 0a 58 09 6c 61 6d 20 2d |n - lam;|.X.lam -|
|00004db0| 3d 20 64 65 67 72 61 64 | 28 32 30 2e 34 2f 33 36 |= degrad|(20.4/36|
|00004dc0| 30 30 29 2a 63 6f 73 28 | 61 29 2f 63 6f 73 28 62 |00)*cos(|a)/cos(b|
|00004dd0| 65 74 29 3b 0a 58 09 62 | 65 74 20 2d 3d 20 64 65 |et);.X.b|et -= de|
|00004de0| 67 72 61 64 28 32 30 2e | 34 2f 33 36 30 30 29 2a |grad(20.|4/3600)*|
|00004df0| 73 69 6e 28 61 29 2a 73 | 69 6e 28 62 65 74 29 3b |sin(a)*s|in(bet);|
|00004e00| 0a 58 0a 58 09 2f 2a 20 | 63 6f 6e 76 65 72 74 20 |.X.X./* |convert |
|00004e10| 67 65 6f 63 65 6e 74 72 | 69 63 20 65 63 6c 69 70 |geocentr|ic eclip|
|00004e20| 74 69 63 20 6c 61 74 2f | 6c 6f 6e 67 20 74 6f 20 |tic lat/|long to |
|00004e30| 65 71 75 69 74 6f 72 69 | 61 6c 20 72 61 2f 64 65 |equitori|al ra/de|
|00004e40| 63 20 6f 66 20 64 61 74 | 65 20 2a 2f 0a 58 09 65 |c of dat|e */.X.e|
|00004e50| 63 6c 5f 65 71 20 28 6d | 6a 64 2c 20 62 65 74 2c |cl_eq (m|jd, bet,|
|00004e60| 20 6c 61 6d 2c 20 26 72 | 61 2c 20 26 64 65 63 29 | lam, &r|a, &dec)|
|00004e70| 3b 0a 58 0a 58 09 2f 2a | 20 73 65 74 20 73 5f 72 |;.X.X./*| set s_r|
|00004e80| 61 2f 73 5f 64 65 63 20 | 74 6f 20 74 68 61 74 20 |a/s_dec |to that |
|00004e90| 6f 66 20 64 65 73 69 72 | 65 64 20 64 69 73 70 6c |of desir|ed displ|
|00004ea0| 61 79 20 65 70 6f 63 68 | 2e 20 2a 2f 0a 58 09 69 |ay epoch|. */.X.i|
|00004eb0| 66 20 28 65 70 6f 63 68 | 20 3d 3d 20 45 4f 44 29 |f (epoch| == EOD)|
|00004ec0| 20 7b 0a 58 09 20 20 20 | 20 6f 70 2d 3e 73 5f 72 | {.X.   | op->s_r|
|00004ed0| 61 20 3d 20 72 61 3b 0a | 58 09 20 20 20 20 6f 70 |a = ra;.|X.    op|
|00004ee0| 2d 3e 73 5f 64 65 63 20 | 3d 20 64 65 63 3b 0a 58 |->s_dec |= dec;.X|
|00004ef0| 09 7d 20 65 6c 73 65 20 | 7b 0a 58 09 20 20 20 20 |.} else |{.X.    |
|00004f00| 64 6f 75 62 6c 65 20 74 | 72 61 20 3d 20 72 61 2c |double t|ra = ra,|
|00004f10| 20 74 64 65 63 20 3d 20 | 64 65 63 3b 0a 58 09 20 | tdec = |dec;.X. |
|00004f20| 20 20 20 70 72 65 63 65 | 73 73 20 28 6d 6a 64 2c |   prece|ss (mjd,|
|00004f30| 20 65 70 6f 63 68 2c 20 | 26 74 72 61 2c 20 26 74 | epoch, |&tra, &t|
|00004f40| 64 65 63 29 3b 0a 58 09 | 20 20 20 20 6f 70 2d 3e |dec);.X.|    op->|
|00004f50| 73 5f 72 61 20 3d 20 74 | 72 61 3b 0a 58 09 20 20 |s_ra = t|ra;.X.  |
|00004f60| 20 20 6f 70 2d 3e 73 5f | 64 65 63 20 3d 20 74 64 |  op->s_|dec = td|
|00004f70| 65 63 3b 0a 58 09 7d 0a | 58 0a 58 09 2f 2a 20 73 |ec;.X.}.|X.X./* s|
|00004f80| 65 74 20 65 61 72 74 68 | 2f 70 6c 61 6e 65 74 20 |et earth|/planet |
|00004f90| 61 6e 64 20 73 75 6e 2f | 70 6c 61 6e 65 74 20 64 |and sun/|planet d|
|00004fa0| 69 73 74 61 6e 63 65 20 | 2a 2f 0a 58 09 6f 70 2d |istance |*/.X.op-|
|00004fb0| 3e 73 5f 65 64 69 73 74 | 20 3d 20 72 68 6f 3b 0a |>s_edist| = rho;.|
|00004fc0| 58 09 6f 70 2d 3e 73 5f | 73 64 69 73 74 20 3d 20 |X.op->s_|sdist = |
|00004fd0| 72 70 3b 0a 58 0a 58 09 | 2f 2a 20 63 6f 6d 70 75 |rp;.X.X.|/* compu|
|00004fe0| 74 65 20 65 6c 6f 6e 67 | 61 74 69 6f 6e 20 61 6e |te elong|ation an|
|00004ff0| 64 20 70 68 61 73 65 20 | 2a 2f 0a 58 09 65 6c 6f |d phase |*/.X.elo|
|00005000| 6e 67 61 74 69 6f 6e 20 | 28 6c 61 6d 2c 20 62 65 |ngation |(lam, be|
|00005010| 74 2c 20 6c 73 6e 2c 20 | 26 65 6c 29 3b 0a 58 09 |t, lsn, |&el);.X.|
|00005020| 65 6c 20 3d 20 72 61 64 | 64 65 67 28 65 6c 29 3b |el = rad|deg(el);|
|00005030| 0a 58 09 6f 70 2d 3e 73 | 5f 65 6c 6f 6e 67 20 3d |.X.op->s|_elong =|
|00005040| 20 65 6c 3b 0a 58 09 66 | 20 3d 20 30 2e 32 35 20 | el;.X.f| = 0.25 |
|00005050| 2a 20 28 28 72 70 2b 72 | 68 6f 29 2a 28 72 70 2b |* ((rp+r|ho)*(rp+|
|00005060| 72 68 6f 29 20 2d 20 72 | 73 6e 2a 72 73 6e 29 2f |rho) - r|sn*rsn)/|
|00005070| 28 72 70 2a 72 68 6f 29 | 3b 0a 58 09 6f 70 2d 3e |(rp*rho)|;.X.op->|
|00005080| 73 5f 70 68 61 73 65 20 | 3d 20 66 2a 31 30 30 2e |s_phase |= f*100.|
|00005090| 30 3b 20 2f 2a 20 70 65 | 72 63 65 6e 74 20 2a 2f |0; /* pe|rcent */|
|000050a0| 0a 58 0a 58 09 2f 2a 20 | 73 65 74 20 68 65 6c 69 |.X.X./* |set heli|
|000050b0| 6f 63 65 6e 74 72 69 63 | 20 6c 6f 6e 67 2f 6c 61 |ocentric| long/la|
|000050c0| 74 20 2a 2f 0a 58 09 6f | 70 2d 3e 73 5f 68 6c 6f |t */.X.o|p->s_hlo|
|000050d0| 6e 67 20 3d 20 6c 70 64 | 3b 0a 58 09 6f 70 2d 3e |ng = lpd|;.X.op->|
|000050e0| 73 5f 68 6c 61 74 20 3d | 20 70 73 69 3b 0a 58 0a |s_hlat =| psi;.X.|
|000050f0| 58 09 2f 2a 20 61 6c 74 | 2c 20 61 7a 3a 20 63 6f |X./* alt|, az: co|
|00005100| 72 72 65 63 74 20 66 6f | 72 20 70 61 72 61 6c 6c |rrect fo|r parall|
|00005110| 61 78 20 61 6e 64 20 72 | 65 66 72 61 63 74 69 6f |ax and r|efractio|
|00005120| 6e 3b 20 75 73 65 20 65 | 6f 64 20 72 61 2f 64 65 |n; use e|od ra/de|
|00005130| 63 2e 20 2a 2f 0a 58 09 | 6e 6f 77 5f 6c 73 74 20 |c. */.X.|now_lst |
|00005140| 28 6e 70 2c 20 26 6c 73 | 74 29 3b 0a 58 09 68 61 |(np, &ls|t);.X.ha|
|00005150| 20 3d 20 68 72 72 61 64 | 28 6c 73 74 29 20 2d 20 | = hrrad|(lst) - |
|00005160| 72 61 3b 0a 58 09 65 68 | 70 20 3d 20 28 32 2e 30 |ra;.X.eh|p = (2.0|
|00005170| 2a 36 33 37 38 2e 30 2f | 31 34 36 2e 30 65 36 29 |*6378.0/|146.0e6)|
|00005180| 20 2f 20 72 68 6f 3b 0a | 58 09 74 61 5f 70 61 72 | / rho;.|X.ta_par|
|00005190| 20 28 68 61 2c 20 64 65 | 63 2c 20 6c 61 74 2c 20 | (ha, de|c, lat, |
|000051a0| 65 6c 65 76 2c 20 65 68 | 70 2c 20 26 68 61 2c 20 |elev, eh|p, &ha, |
|000051b0| 26 64 65 63 29 3b 0a 58 | 09 68 61 64 65 63 5f 61 |&dec);.X|.hadec_a|
|000051c0| 61 20 28 6c 61 74 2c 20 | 68 61 2c 20 64 65 63 2c |a (lat, |ha, dec,|
|000051d0| 20 26 61 6c 74 2c 20 26 | 61 7a 29 3b 0a 58 09 72 | &alt, &|az);.X.r|
|000051e0| 65 66 72 61 63 74 20 28 | 70 72 65 73 73 75 72 65 |efract (|pressure|
|000051f0| 2c 20 74 65 6d 70 2c 20 | 61 6c 74 2c 20 26 61 6c |, temp, |alt, &al|
|00005200| 74 29 3b 0a 58 09 6f 70 | 2d 3e 73 5f 61 6c 74 20 |t);.X.op|->s_alt |
|00005210| 3d 20 61 6c 74 3b 0a 58 | 09 6f 70 2d 3e 73 5f 61 |= alt;.X|.op->s_a|
|00005220| 7a 20 3d 20 61 7a 3b 0a | 58 7d 0a 58 0a 58 2f 2a |z = az;.|X}.X.X/*|
|00005230| 20 67 69 76 65 6e 20 67 | 65 6f 63 65 6e 74 72 69 | given g|eocentri|
|00005240| 63 20 65 63 6c 69 70 74 | 69 63 20 6c 6f 6e 67 69 |c eclipt|ic longi|
|00005250| 74 75 64 65 20 61 6e 64 | 20 6c 61 74 69 74 75 64 |tude and| latitud|
|00005260| 65 2c 20 6c 61 6d 20 61 | 6e 64 20 62 65 74 2c 20 |e, lam a|nd bet, |
|00005270| 6f 66 20 73 6f 6d 65 20 | 6f 62 6a 65 63 74 0a 58 |of some |object.X|
|00005280| 20 2a 20 61 6e 64 20 74 | 68 65 20 6c 6f 6e 67 69 | * and t|he longi|
|00005290| 74 75 64 65 20 6f 66 20 | 74 68 65 20 73 75 6e 2c |tude of |the sun,|
|000052a0| 20 6c 73 6e 2c 20 66 69 | 6e 64 20 74 68 65 20 65 | lsn, fi|nd the e|
|000052b0| 6c 6f 6e 67 61 74 69 6f | 6e 2c 20 65 6c 2e 20 74 |longatio|n, el. t|
|000052c0| 68 69 73 20 69 73 20 74 | 68 65 0a 58 20 2a 20 61 |his is t|he.X * a|
|000052d0| 63 74 75 61 6c 20 61 6e | 67 75 6c 61 72 20 73 65 |ctual an|gular se|
|000052e0| 70 61 72 61 74 69 6f 6e | 20 6f 66 20 74 68 65 20 |paration| of the |
|000052f0| 6f 62 6a 65 63 74 20 66 | 72 6f 6d 20 74 68 65 20 |object f|rom the |
|00005300| 73 75 6e 2c 20 6e 6f 74 | 20 6a 75 73 74 20 74 68 |sun, not| just th|
|00005310| 65 20 64 69 66 66 65 72 | 65 6e 63 65 0a 58 20 2a |e differ|ence.X *|
|00005320| 20 69 6e 20 74 68 65 20 | 6c 6f 6e 67 69 74 75 64 | in the |longitud|
|00005330| 65 2e 20 74 68 65 20 73 | 69 67 6e 2c 20 68 6f 77 |e. the s|ign, how|
|00005340| 65 76 65 72 2c 20 49 53 | 20 73 65 74 20 73 69 6d |ever, IS| set sim|
|00005350| 70 6c 79 20 61 73 20 61 | 20 74 65 73 74 20 6f 6e |ply as a| test on|
|00005360| 20 6c 6f 6e 67 69 74 75 | 64 65 0a 58 20 2a 20 73 | longitu|de.X * s|
|00005370| 75 63 68 20 74 68 61 74 | 20 65 6c 20 77 69 6c 6c |uch that| el will|
|00005380| 20 62 65 20 3e 30 20 66 | 6f 72 20 61 6e 20 65 76 | be >0 f|or an ev|
|00005390| 65 6e 69 6e 67 20 6f 62 | 6a 65 63 74 20 3c 30 20 |ening ob|ject <0 |
|000053a0| 66 6f 72 20 61 20 6d 6f | 72 6e 69 6e 67 20 6f 62 |for a mo|rning ob|
|000053b0| 6a 65 63 74 2e 0a 58 20 | 2a 20 74 6f 20 75 6e 64 |ject..X |* to und|
|000053c0| 65 72 73 74 61 6e 64 20 | 74 68 65 20 74 65 73 74 |erstand |the test|
|000053d0| 20 66 6f 72 20 65 6c 20 | 73 69 67 6e 2c 20 64 72 | for el |sign, dr|
|000053e0| 61 77 20 61 20 67 72 61 | 70 68 20 77 69 74 68 20 |aw a gra|ph with |
|000053f0| 6c 61 6d 20 67 6f 69 6e | 67 20 66 72 6f 6d 20 30 |lam goin|g from 0|
|00005400| 2d 32 2a 50 49 0a 58 20 | 2a 20 20 20 64 6f 77 6e |-2*PI.X |*   down|
|00005410| 20 74 68 65 20 76 65 72 | 74 69 63 61 6c 20 61 78 | the ver|tical ax|
|00005420| 69 73 2c 20 6c 73 6e 20 | 67 6f 69 6e 67 20 66 72 |is, lsn |going fr|
|00005430| 6f 6d 20 30 2d 32 2a 50 | 49 20 61 63 72 6f 73 73 |om 0-2*P|I across|
|00005440| 20 74 68 65 20 68 6f 72 | 20 61 78 69 73 2e 20 74 | the hor| axis. t|
|00005450| 68 65 6e 0a 58 20 2a 20 | 20 20 64 65 66 69 6e 65 |hen.X * |  define|
|00005460| 20 74 68 65 20 64 69 61 | 67 6f 6e 61 6c 20 72 65 | the dia|gonal re|
|00005470| 67 69 6f 6e 73 20 62 6f | 75 6e 64 65 64 20 62 79 |gions bo|unded by|
|00005480| 20 74 68 65 20 6c 69 6e | 65 73 20 6c 61 6d 3d 6c | the lin|es lam=l|
|00005490| 73 6e 2b 50 49 2c 20 6c | 61 6d 3d 6c 73 6e 20 61 |sn+PI, l|am=lsn a|
|000054a0| 6e 64 0a 58 20 2a 20 20 | 20 6c 61 6d 3d 6c 73 6e |nd.X *  | lam=lsn|
|000054b0| 2d 50 49 2e 20 74 68 65 | 20 22 6d 6f 72 6e 69 6e |-PI. the| "mornin|
|000054c0| 67 22 20 72 65 67 69 6f | 6e 73 20 61 72 65 20 61 |g" regio|ns are a|
|000054d0| 6e 79 20 76 61 6c 75 65 | 73 20 74 6f 20 74 68 65 |ny value|s to the|
|000054e0| 20 6c 6f 77 65 72 20 6c | 65 66 74 20 6f 66 20 74 | lower l|eft of t|
|000054f0| 68 65 0a 58 20 2a 20 20 | 20 66 69 72 73 74 20 6c |he.X *  | first l|
|00005500| 69 6e 65 20 61 6e 64 20 | 62 6f 75 6e 64 65 64 20 |ine and |bounded |
|00005510| 77 69 74 68 69 6e 20 74 | 68 65 20 73 65 63 6f 6e |within t|he secon|
|00005520| 64 20 70 61 69 72 20 6f | 66 20 6c 69 6e 65 73 2e |d pair o|f lines.|
|00005530| 0a 58 20 2a 20 61 6c 6c | 20 61 6e 67 6c 65 73 20 |.X * all| angles |
|00005540| 69 6e 20 72 61 64 69 61 | 6e 73 2e 0a 58 20 2a 2f |in radia|ns..X */|
|00005550| 0a 58 73 74 61 74 69 63 | 20 76 6f 69 64 0a 58 65 |.Xstatic| void.Xe|
|00005560| 6c 6f 6e 67 61 74 69 6f | 6e 20 28 6c 61 6d 2c 20 |longatio|n (lam, |
|00005570| 62 65 74 2c 20 6c 73 6e | 2c 20 65 6c 29 0a 58 64 |bet, lsn|, el).Xd|
|00005580| 6f 75 62 6c 65 20 6c 61 | 6d 2c 20 62 65 74 2c 20 |ouble la|m, bet, |
|00005590| 6c 73 6e 3b 0a 58 64 6f | 75 62 6c 65 20 2a 65 6c |lsn;.Xdo|uble *el|
|000055a0| 3b 0a 58 7b 0a 58 09 2a | 65 6c 20 3d 20 61 63 6f |;.X{.X.*|el = aco|
|000055b0| 73 28 63 6f 73 28 62 65 | 74 29 2a 63 6f 73 28 6c |s(cos(be|t)*cos(l|
|000055c0| 61 6d 2d 6c 73 6e 29 29 | 3b 0a 58 09 69 66 20 28 |am-lsn))|;.X.if (|
|000055d0| 6c 61 6d 3e 6c 73 6e 2b | 50 49 20 7c 7c 20 6c 61 |lam>lsn+|PI || la|
|000055e0| 6d 3e 6c 73 6e 2d 50 49 | 20 26 26 20 6c 61 6d 3c |m>lsn-PI| && lam<|
|000055f0| 6c 73 6e 29 20 2a 65 6c | 20 3d 20 2d 20 2a 65 6c |lsn) *el| = - *el|
|00005600| 3b 0a 58 7d 0a 45 4e 44 | 5f 4f 46 5f 46 49 4c 45 |;.X}.END|_OF_FILE|
|00005610| 0a 20 20 69 66 20 74 65 | 73 74 20 32 30 30 30 38 |.  if te|st 20008|
|00005620| 20 2d 6e 65 20 60 77 63 | 20 2d 63 20 3c 27 63 69 | -ne `wc| -c <'ci|
|00005630| 72 63 75 6d 2e 63 27 60 | 3b 20 74 68 65 6e 0a 20 |rcum.c'`|; then. |
|00005640| 20 20 20 65 63 68 6f 20 | 73 68 61 72 3a 20 5c 22 |   echo |shar: \"|
|00005650| 27 63 69 72 63 75 6d 2e | 63 27 5c 22 20 75 6e 70 |'circum.|c'\" unp|
|00005660| 61 63 6b 65 64 20 77 69 | 74 68 20 77 72 6f 6e 67 |acked wi|th wrong|
|00005670| 20 73 69 7a 65 21 0a 20 | 20 66 69 0a 20 20 23 20 | size!. | fi.  # |
|00005680| 65 6e 64 20 6f 66 20 27 | 63 69 72 63 75 6d 2e 63 |end of '|circum.c|
|00005690| 27 0a 66 69 0a 69 66 20 | 74 65 73 74 20 2d 66 20 |'.fi.if |test -f |
|000056a0| 27 6d 69 73 63 2e 63 27 | 20 2d 61 20 22 24 7b 31 |'misc.c'| -a "${1|
|000056b0| 7d 22 20 21 3d 20 22 2d | 63 22 20 3b 20 74 68 65 |}" != "-|c" ; the|
|000056c0| 6e 20 0a 20 20 65 63 68 | 6f 20 73 68 61 72 3a 20 |n .  ech|o shar: |
|000056d0| 57 69 6c 6c 20 6e 6f 74 | 20 63 6c 6f 62 62 65 72 |Will not| clobber|
|000056e0| 20 65 78 69 73 74 69 6e | 67 20 66 69 6c 65 20 5c | existin|g file \|
|000056f0| 22 27 6d 69 73 63 2e 63 | 27 5c 22 0a 65 6c 73 65 |"'misc.c|'\".else|
|00005700| 0a 20 20 65 63 68 6f 20 | 73 68 61 72 3a 20 45 78 |.  echo |shar: Ex|
|00005710| 74 72 61 63 74 69 6e 67 | 20 5c 22 27 6d 69 73 63 |tracting| \"'misc|
|00005720| 2e 63 27 5c 22 20 5c 28 | 31 39 35 36 37 20 63 68 |.c'\" \(|19567 ch|
|00005730| 61 72 61 63 74 65 72 73 | 5c 29 0a 20 20 73 65 64 |aracters|\).  sed|
|00005740| 20 22 73 2f 5e 58 2f 2f | 22 20 3e 27 6d 69 73 63 | "s/^X//|" >'misc|
|00005750| 2e 63 27 20 3c 3c 27 45 | 4e 44 5f 4f 46 5f 46 49 |.c' <<'E|ND_OF_FI|
|00005760| 4c 45 27 0a 58 2f 2a 20 | 6d 69 73 63 20 68 61 6e |LE'.X/* |misc han|
|00005770| 64 79 20 66 75 6e 63 74 | 69 6f 6e 73 2e 0a 58 20 |dy funct|ions..X |
|00005780| 2a 20 65 76 65 72 79 20 | 73 79 73 74 65 6d 20 68 |* every |system h|
|00005790| 61 73 20 73 75 63 68 2c | 20 6e 6f 3f 0a 58 20 2a |as such,| no?.X *|
|000057a0| 2f 0a 58 0a 58 23 69 6e | 63 6c 75 64 65 20 3c 73 |/.X.X#in|clude <s|
|000057b0| 74 64 69 6f 2e 68 3e 0a | 58 23 69 6e 63 6c 75 64 |tdio.h>.|X#includ|
|000057c0| 65 20 3c 63 74 79 70 65 | 2e 68 3e 0a 58 23 69 6e |e <ctype|.h>.X#in|
|000057d0| 63 6c 75 64 65 20 3c 6d | 61 74 68 2e 68 3e 0a 58 |clude <m|ath.h>.X|
|000057e0| 23 69 66 20 64 65 66 69 | 6e 65 64 28 5f 5f 53 54 |#if defi|ned(__ST|
|000057f0| 44 43 5f 5f 29 0a 58 23 | 69 6e 63 6c 75 64 65 20 |DC__).X#|include |
|00005800| 3c 73 74 64 6c 69 62 2e | 68 3e 0a 58 23 65 6e 64 |<stdlib.|h>.X#end|
|00005810| 69 66 0a 58 23 69 6e 63 | 6c 75 64 65 20 3c 58 31 |if.X#inc|lude <X1|
|00005820| 31 2f 58 6c 69 62 2e 68 | 3e 0a 58 23 69 6e 63 6c |1/Xlib.h|>.X#incl|
|00005830| 75 64 65 20 3c 58 31 31 | 2f 63 75 72 73 6f 72 66 |ude <X11|/cursorf|
|00005840| 6f 6e 74 2e 68 3e 0a 58 | 23 69 6e 63 6c 75 64 65 |ont.h>.X|#include|
|00005850| 20 3c 58 6d 2f 58 6d 2e | 68 3e 0a 58 23 69 6e 63 | <Xm/Xm.|h>.X#inc|
|00005860| 6c 75 64 65 20 3c 58 6d | 2f 50 75 73 68 42 2e 68 |lude <Xm|/PushB.h|
|00005870| 3e 0a 58 23 69 6e 63 6c | 75 64 65 20 22 61 73 74 |>.X#incl|ude "ast|
|00005880| 72 6f 2e 68 22 0a 58 23 | 69 6e 63 6c 75 64 65 20 |ro.h".X#|include |
|00005890| 22 63 69 72 63 75 6d 2e | 68 22 0a 58 0a 58 23 69 |"circum.|h".X.X#i|
|000058a0| 66 20 64 65 66 69 6e 65 | 64 28 5f 5f 53 54 44 43 |f define|d(__STDC|
|000058b0| 5f 5f 29 20 7c 7c 20 64 | 65 66 69 6e 65 64 28 5f |__) || d|efined(_|
|000058c0| 5f 63 70 6c 75 73 70 6c | 75 73 29 0a 58 23 64 65 |_cpluspl|us).X#de|
|000058d0| 66 69 6e 65 20 50 5f 28 | 73 29 20 73 0a 58 23 65 |fine P_(|s) s.X#e|
|000058e0| 6c 73 65 0a 58 23 64 65 | 66 69 6e 65 20 50 5f 28 |lse.X#de|fine P_(|
|000058f0| 73 29 20 28 29 0a 58 23 | 65 6e 64 69 66 0a 58 0a |s) ().X#|endif.X.|
|00005900| 58 65 78 74 65 72 6e 20 | 4e 6f 77 20 2a 6d 6d 5f |Xextern |Now *mm_|
|00005910| 67 65 74 5f 6e 6f 77 20 | 50 5f 28 28 76 6f 69 64 |get_now |P_((void|
|00005920| 29 29 3b 0a 58 65 78 74 | 65 72 6e 20 69 6e 74 20 |));.Xext|ern int |
|00005930| 6c 69 73 74 69 6e 67 5f | 69 73 6f 6e 20 50 5f 28 |listing_|ison P_(|
|00005940| 28 76 6f 69 64 29 29 3b | 0a 58 65 78 74 65 72 6e |(void));|.Xextern|
|00005950| 20 69 6e 74 20 70 6c 6f | 74 5f 69 73 6f 6e 20 50 | int plo|t_ison P|
|00005960| 5f 28 28 76 6f 69 64 29 | 29 3b 0a 58 65 78 74 65 |_((void)|);.Xexte|
|00005970| 72 6e 20 69 6e 74 20 73 | 72 63 68 5f 69 73 6f 6e |rn int s|rch_ison|
|00005980| 20 50 5f 28 28 76 6f 69 | 64 29 29 3b 0a 58 65 78 | P_((voi|d));.Xex|
|00005990| 74 65 72 6e 20 76 6f 69 | 64 20 64 6d 5f 6e 65 77 |tern voi|d dm_new|
|000059a0| 6f 62 6a 20 50 5f 28 28 | 69 6e 74 20 64 62 69 64 |obj P_((|int dbid|
|000059b0| 78 29 29 3b 0a 58 65 78 | 74 65 72 6e 20 76 6f 69 |x));.Xex|tern voi|
|000059c0| 64 20 64 6d 5f 73 65 6c | 65 63 74 69 6f 6e 5f 6d |d dm_sel|ection_m|
|000059d0| 6f 64 65 20 50 5f 28 28 | 69 6e 74 20 77 68 65 74 |ode P_((|int whet|
|000059e0| 68 65 72 29 29 3b 0a 58 | 65 78 74 65 72 6e 20 76 |her));.X|extern v|
|000059f0| 6f 69 64 20 64 6d 5f 75 | 70 64 61 74 65 20 50 5f |oid dm_u|pdate P_|
|00005a00| 28 28 4e 6f 77 20 2a 6e | 70 2c 20 69 6e 74 20 68 |((Now *n|p, int h|
|00005a10| 6f 77 5f 6d 75 63 68 29 | 29 3b 0a 58 65 78 74 65 |ow_much)|);.Xexte|
|00005a20| 72 6e 20 76 6f 69 64 20 | 65 5f 63 75 72 73 6f 72 |rn void |e_cursor|
|00005a30| 20 50 5f 28 28 43 75 72 | 73 6f 72 20 63 29 29 3b | P_((Cur|sor c));|
|00005a40| 0a 58 65 78 74 65 72 6e | 20 76 6f 69 64 20 65 5f |.Xextern| void e_|
|00005a50| 73 65 6c 65 63 74 69 6f | 6e 5f 6d 6f 64 65 20 50 |selectio|n_mode P|
|00005a60| 5f 28 28 69 6e 74 20 77 | 68 65 74 68 65 72 29 29 |_((int w|hether))|
|00005a70| 3b 0a 58 65 78 74 65 72 | 6e 20 76 6f 69 64 20 65 |;.Xexter|n void e|
|00005a80| 5f 75 70 64 61 74 65 20 | 50 5f 28 28 4e 6f 77 20 |_update |P_((Now |
|00005a90| 2a 6e 70 2c 20 69 6e 74 | 20 66 6f 72 63 65 29 29 |*np, int| force))|
|00005aa0| 3b 0a 58 65 78 74 65 72 | 6e 20 76 6f 69 64 20 66 |;.Xexter|n void f|
|00005ab0| 5f 73 68 6f 77 69 74 20 | 50 5f 28 28 57 69 64 67 |_showit |P_((Widg|
|00005ac0| 65 74 20 77 2c 20 63 68 | 61 72 20 2a 73 29 29 3b |et w, ch|ar *s));|
|00005ad0| 0a 58 65 78 74 65 72 6e | 20 76 6f 69 64 20 66 73 |.Xextern| void fs|
|00005ae0| 5f 64 61 74 65 20 50 5f | 28 28 63 68 61 72 20 6f |_date P_|((char o|
|00005af0| 75 74 5b 5d 2c 20 64 6f | 75 62 6c 65 20 6a 64 29 |ut[], do|uble jd)|
|00005b00| 29 3b 0a 58 65 78 74 65 | 72 6e 20 76 6f 69 64 20 |);.Xexte|rn void |
|00005b10| 66 73 5f 74 69 6d 65 20 | 50 5f 28 28 63 68 61 72 |fs_time |P_((char|
|00005b20| 20 6f 75 74 5b 5d 2c 20 | 64 6f 75 62 6c 65 20 74 | out[], |double t|
|00005b30| 29 29 3b 0a 58 65 78 74 | 65 72 6e 20 76 6f 69 64 |));.Xext|ern void|
|00005b40| 20 6a 6d 5f 73 65 6c 65 | 63 74 69 6f 6e 5f 6d 6f | jm_sele|ction_mo|
|00005b50| 64 65 20 50 5f 28 28 69 | 6e 74 20 77 68 65 74 68 |de P_((i|nt wheth|
|00005b60| 65 72 29 29 3b 0a 58 65 | 78 74 65 72 6e 20 76 6f |er));.Xe|xtern vo|
|00005b70| 69 64 20 6a 6d 5f 75 70 | 64 61 74 65 20 50 5f 28 |id jm_up|date P_(|
|00005b80| 28 4e 6f 77 20 2a 6e 70 | 2c 20 69 6e 74 20 68 6f |(Now *np|, int ho|
|00005b90| 77 5f 6d 75 63 68 29 29 | 3b 0a 58 65 78 74 65 72 |w_much))|;.Xexter|
|00005ba0| 6e 20 76 6f 69 64 20 6c | 73 74 5f 6c 6f 67 20 50 |n void l|st_log P|
|00005bb0| 5f 28 28 63 68 61 72 20 | 2a 6e 61 6d 65 2c 20 63 |_((char |*name, c|
|00005bc0| 68 61 72 20 2a 73 74 72 | 29 29 3b 0a 58 65 78 74 |har *str|));.Xext|
|00005bd0| 65 72 6e 20 76 6f 69 64 | 20 6c 73 74 5f 73 65 6c |ern void| lst_sel|
|00005be0| 65 63 74 69 6f 6e 20 50 | 5f 28 28 63 68 61 72 20 |ection P|_((char |
|00005bf0| 2a 6e 61 6d 65 29 29 3b | 0a 58 65 78 74 65 72 6e |*name));|.Xextern|
|00005c00| 20 76 6f 69 64 20 6d 5f | 75 70 64 61 74 65 20 50 | void m_|update P|
|00005c10| 5f 28 28 4e 6f 77 20 2a | 6e 70 2c 20 69 6e 74 20 |_((Now *|np, int |
|00005c20| 68 6f 77 5f 6d 75 63 68 | 29 29 3b 0a 58 65 78 74 |how_much|));.Xext|
|00005c30| 65 72 6e 20 76 6f 69 64 | 20 6d 61 72 73 5f 63 75 |ern void| mars_cu|
|00005c40| 72 73 6f 72 20 50 5f 28 | 28 43 75 72 73 6f 72 20 |rsor P_(|(Cursor |
|00005c50| 63 29 29 3b 0a 58 65 78 | 74 65 72 6e 20 76 6f 69 |c));.Xex|tern voi|
|00005c60| 64 20 6d 61 72 73 5f 73 | 65 6c 65 63 74 69 6f 6e |d mars_s|election|
|00005c70| 5f 6d 6f 64 65 20 50 5f | 28 28 69 6e 74 20 77 68 |_mode P_|((int wh|
|00005c80| 65 74 68 65 72 29 29 3b | 0a 58 65 78 74 65 72 6e |ether));|.Xextern|
|00005c90| 20 76 6f 69 64 20 6d 61 | 72 73 5f 75 70 64 61 74 | void ma|rs_updat|
|00005ca0| 65 20 50 5f 28 28 4e 6f | 77 20 2a 6e 70 2c 20 69 |e P_((No|w *np, i|
|00005cb0| 6e 74 20 66 6f 72 63 65 | 29 29 3b 0a 58 65 78 74 |nt force|));.Xext|
|00005cc0| 65 72 6e 20 76 6f 69 64 | 20 6d 6d 5f 73 65 6c 65 |ern void| mm_sele|
|00005cd0| 63 74 69 6f 6e 5f 6d 6f | 64 65 20 50 5f 28 28 69 |ction_mo|de P_((i|
|00005ce0| 6e 74 20 77 68 65 74 68 | 65 72 29 29 3b 0a 58 65 |nt wheth|er));.Xe|
|00005cf0| 78 74 65 72 6e 20 76 6f | 69 64 20 6f 62 6a 5f 6e |xtern vo|id obj_n|
|00005d00| 65 77 64 62 20 50 5f 28 | 28 69 6e 74 20 61 70 70 |ewdb P_(|(int app|
|00005d10| 65 6e 64 65 64 29 29 3b | 0a 58 65 78 74 65 72 6e |ended));|.Xextern|
|00005d20| 20 76 6f 69 64 20 6f 62 | 6a 5f 75 70 64 61 74 65 | void ob|j_update|
|00005d30| 20 50 5f 28 28 4e 6f 77 | 20 2a 6e 70 2c 20 69 6e | P_((Now| *np, in|
|00005d40| 74 20 68 6f 77 6d 75 63 | 68 29 29 3b 0a 58 65 78 |t howmuc|h));.Xex|
|00005d50| 74 65 72 6e 20 76 6f 69 | 64 20 70 6c 74 5f 6c 6f |tern voi|d plt_lo|
|00005d60| 67 20 50 5f 28 28 63 68 | 61 72 20 2a 6e 61 6d 65 |g P_((ch|ar *name|
|00005d70| 2c 20 64 6f 75 62 6c 65 | 20 76 61 6c 75 65 29 29 |, double| value))|
|00005d80| 3b 0a 58 65 78 74 65 72 | 6e 20 76 6f 69 64 20 70 |;.Xexter|n void p|
|00005d90| 6c 74 5f 73 65 6c 65 63 | 74 69 6f 6e 20 50 5f 28 |lt_selec|tion P_(|
|00005da0| 28 63 68 61 72 20 2a 6e | 61 6d 65 29 29 3b 0a 58 |(char *n|ame));.X|
|00005db0| 65 78 74 65 72 6e 20 76 | 6f 69 64 20 73 6d 5f 73 |extern v|oid sm_s|
|00005dc0| 65 6c 65 63 74 69 6f 6e | 5f 6d 6f 64 65 20 50 5f |election|_mode P_|
|00005dd0| 28 28 69 6e 74 20 77 68 | 65 74 68 65 72 29 29 3b |((int wh|ether));|
|00005de0| 0a 58 65 78 74 65 72 6e | 20 76 6f 69 64 20 73 6d |.Xextern| void sm|
|00005df0| 5f 75 70 64 61 74 65 20 | 50 5f 28 28 4e 6f 77 20 |_update |P_((Now |
|00005e00| 2a 6e 70 2c 20 69 6e 74 | 20 68 6f 77 5f 6d 75 63 |*np, int| how_muc|
|00005e10| 68 29 29 3b 0a 58 65 78 | 74 65 72 6e 20 76 6f 69 |h));.Xex|tern voi|
|00005e20| 64 20 73 72 63 68 5f 6c | 6f 67 20 50 5f 28 28 63 |d srch_l|og P_((c|
|00005e30| 68 61 72 20 2a 6e 61 6d | 65 2c 20 64 6f 75 62 6c |har *nam|e, doubl|
|00005e40| 65 20 76 61 6c 75 65 29 | 29 3b 0a 58 65 78 74 65 |e value)|);.Xexte|
|00005e50| 72 6e 20 76 6f 69 64 20 | 73 72 63 68 5f 73 65 6c |rn void |srch_sel|
|00005e60| 65 63 74 69 6f 6e 20 50 | 5f 28 28 63 68 61 72 20 |ection P|_((char |
|00005e70| 2a 6e 61 6d 65 29 29 3b | 0a 58 65 78 74 65 72 6e |*name));|.Xextern|
|00005e80| 20 76 6f 69 64 20 73 72 | 63 68 5f 73 65 6c 65 63 | void sr|ch_selec|
|00005e90| 74 69 6f 6e 5f 6d 6f 64 | 65 20 50 5f 28 28 69 6e |tion_mod|e P_((in|
|00005ea0| 74 20 77 68 65 74 68 65 | 72 29 29 3b 0a 58 65 78 |t whethe|r));.Xex|
|00005eb0| 74 65 72 6e 20 76 6f 69 | 64 20 73 73 5f 6e 65 77 |tern voi|d ss_new|
|00005ec0| 6f 62 6a 20 50 5f 28 28 | 69 6e 74 20 64 62 69 64 |obj P_((|int dbid|
|00005ed0| 78 29 29 3b 0a 58 65 78 | 74 65 72 6e 20 76 6f 69 |x));.Xex|tern voi|
|00005ee0| 64 20 73 73 5f 75 70 64 | 61 74 65 20 50 5f 28 28 |d ss_upd|ate P_((|
|00005ef0| 4e 6f 77 20 2a 6e 70 2c | 20 69 6e 74 20 68 6f 77 |Now *np,| int how|
|00005f00| 5f 6d 75 63 68 29 29 3b | 0a 58 65 78 74 65 72 6e |_much));|.Xextern|
|00005f10| 20 76 6f 69 64 20 73 76 | 5f 6e 65 77 64 62 20 50 | void sv|_newdb P|
|00005f20| 5f 28 28 69 6e 74 20 61 | 70 70 65 6e 64 65 64 29 |_((int a|ppended)|
|00005f30| 29 3b 0a 58 65 78 74 65 | 72 6e 20 76 6f 69 64 20 |);.Xexte|rn void |
|00005f40| 73 76 5f 6e 65 77 6f 62 | 6a 20 50 5f 28 28 69 6e |sv_newob|j P_((in|
|00005f50| 74 20 64 62 69 64 78 29 | 29 3b 0a 58 65 78 74 65 |t dbidx)|);.Xexte|
|00005f60| 72 6e 20 76 6f 69 64 20 | 73 76 5f 75 70 64 61 74 |rn void |sv_updat|
|00005f70| 65 20 50 5f 28 28 4e 6f | 77 20 2a 6e 70 2c 20 69 |e P_((No|w *np, i|
|00005f80| 6e 74 20 68 6f 77 5f 6d | 75 63 68 29 29 3b 0a 58 |nt how_m|uch));.X|
|00005f90| 65 78 74 65 72 6e 20 76 | 6f 69 64 20 75 74 63 5f |extern v|oid utc_|
|00005fa0| 67 73 74 20 50 5f 28 28 | 64 6f 75 62 6c 65 20 4d |gst P_((|double M|
|00005fb0| 6a 64 2c 20 64 6f 75 62 | 6c 65 20 75 74 63 2c 20 |jd, doub|le utc, |
|00005fc0| 64 6f 75 62 6c 65 20 2a | 67 73 74 29 29 3b 0a 58 |double *|gst));.X|
|00005fd0| 65 78 74 65 72 6e 20 76 | 6f 69 64 20 78 65 5f 6d |extern v|oid xe_m|
|00005fe0| 73 67 20 50 5f 28 28 63 | 68 61 72 20 2a 6d 73 67 |sg P_((c|har *msg|
|00005ff0| 2c 20 69 6e 74 20 61 70 | 70 5f 6d 6f 64 61 6c 29 |, int ap|p_modal)|
|00006000| 29 3b 0a 58 0a 58 76 6f | 69 64 20 73 65 74 5f 73 |);.X.Xvo|id set_s|
|00006010| 6f 6d 65 74 68 69 6e 67 | 20 50 5f 28 28 57 69 64 |omething| P_((Wid|
|00006020| 67 65 74 20 77 2c 20 63 | 68 61 72 20 2a 72 65 73 |get w, c|har *res|
|00006030| 6f 75 72 63 65 2c 20 63 | 68 61 72 20 2a 76 61 6c |ource, c|har *val|
|00006040| 75 65 29 29 3b 0a 58 76 | 6f 69 64 20 67 65 74 5f |ue));.Xv|oid get_|
|00006050| 73 6f 6d 65 74 68 69 6e | 67 20 50 5f 28 28 57 69 |somethin|g P_((Wi|
|00006060| 64 67 65 74 20 77 2c 20 | 63 68 61 72 20 2a 72 65 |dget w, |char *re|
|00006070| 73 6f 75 72 63 65 2c 20 | 63 68 61 72 20 2a 76 61 |source, |char *va|
|00006080| 6c 75 65 29 29 3b 0a 58 | 76 6f 69 64 20 67 65 74 |lue));.X|void get|
|00006090| 5f 78 6d 73 74 72 69 6e | 67 20 50 5f 28 28 57 69 |_xmstrin|g P_((Wi|
|000060a0| 64 67 65 74 20 77 2c 20 | 63 68 61 72 20 2a 72 65 |dget w, |char *re|
|000060b0| 73 6f 75 72 63 65 2c 20 | 63 68 61 72 20 2a 2a 74 |source, |char **t|
|000060c0| 78 74 70 29 29 3b 0a 58 | 76 6f 69 64 20 73 65 74 |xtp));.X|void set|
|000060d0| 5f 78 6d 73 74 72 69 6e | 67 20 50 5f 28 28 57 69 |_xmstrin|g P_((Wi|
|000060e0| 64 67 65 74 20 77 2c 20 | 63 68 61 72 20 2a 72 65 |dget w, |char *re|
|000060f0| 73 6f 75 72 63 65 2c 20 | 63 68 61 72 20 2a 74 78 |source, |char *tx|
|00006100| 74 29 29 3b 0a 58 76 6f | 69 64 20 72 61 6e 67 65 |t));.Xvo|id range|
|00006110| 20 50 5f 28 28 64 6f 75 | 62 6c 65 20 2a 76 2c 20 | P_((dou|ble *v, |
|00006120| 64 6f 75 62 6c 65 20 72 | 29 29 3b 0a 58 76 6f 69 |double r|));.Xvoi|
|00006130| 64 20 6e 6f 77 5f 6c 73 | 74 20 50 5f 28 28 4e 6f |d now_ls|t P_((No|
|00006140| 77 20 2a 6e 70 2c 20 64 | 6f 75 62 6c 65 20 2a 6c |w *np, d|ouble *l|
|00006150| 73 74 29 29 3b 0a 58 76 | 6f 69 64 20 72 6e 64 5f |st));.Xv|oid rnd_|
|00006160| 73 65 63 6f 6e 64 20 50 | 5f 28 28 64 6f 75 62 6c |second P|_((doubl|
|00006170| 65 20 2a 74 29 29 3b 0a | 58 64 6f 75 62 6c 65 20 |e *t));.|Xdouble |
|00006180| 6d 6a 64 5f 64 61 79 20 | 50 5f 28 28 64 6f 75 62 |mjd_day |P_((doub|
|00006190| 6c 65 20 6a 64 29 29 3b | 0a 58 64 6f 75 62 6c 65 |le jd));|.Xdouble|
|000061a0| 20 6d 6a 64 5f 68 72 20 | 50 5f 28 28 64 6f 75 62 | mjd_hr |P_((doub|
|000061b0| 6c 65 20 6a 64 29 29 3b | 0a 58 76 6f 69 64 20 7a |le jd));|.Xvoid z|
|000061c0| 65 72 6f 5f 6d 65 6d 20 | 50 5f 28 28 63 68 61 72 |ero_mem |P_((char|
|000061d0| 20 2a 6c 6f 63 2c 20 75 | 6e 73 69 67 6e 65 64 20 | *loc, u|nsigned |
|000061e0| 6c 65 6e 29 29 3b 0a 58 | 76 6f 69 64 20 77 61 74 |len));.X|void wat|
|000061f0| 63 68 5f 63 75 72 73 6f | 72 20 50 5f 28 28 69 6e |ch_curso|r P_((in|
|00006200| 74 20 77 61 6e 74 29 29 | 3b 0a 58 76 6f 69 64 20 |t want))|;.Xvoid |
|00006210| 61 6c 6c 5f 75 70 64 61 | 74 65 20 50 5f 28 28 4e |all_upda|te P_((N|
|00006220| 6f 77 20 2a 6e 70 2c 20 | 69 6e 74 20 68 6f 77 5f |ow *np, |int how_|
|00006230| 6d 75 63 68 29 29 3b 0a | 58 76 6f 69 64 20 61 6c |much));.|Xvoid al|
|00006240| 6c 5f 6e 65 77 6f 62 6a | 20 50 5f 28 28 69 6e 74 |l_newobj| P_((int|
|00006250| 20 64 62 69 64 78 29 29 | 3b 0a 58 76 6f 69 64 20 | dbidx))|;.Xvoid |
|00006260| 61 6c 6c 5f 6e 65 77 64 | 62 20 50 5f 28 28 69 6e |all_newd|b P_((in|
|00006270| 74 20 61 70 70 65 6e 64 | 65 64 29 29 3b 0a 58 76 |t append|ed));.Xv|
|00006280| 6f 69 64 20 61 6c 6c 5f | 73 65 6c 65 63 74 69 6f |oid all_|selectio|
|00006290| 6e 5f 6d 6f 64 65 20 50 | 5f 28 28 69 6e 74 20 77 |n_mode P|_((int w|
|000062a0| 68 65 74 68 65 72 29 29 | 3b 0a 58 76 6f 69 64 20 |hether))|;.Xvoid |
|000062b0| 72 65 67 69 73 74 65 72 | 5f 73 65 6c 65 63 74 69 |register|_selecti|
|000062c0| 6f 6e 20 50 5f 28 28 63 | 68 61 72 20 2a 6e 61 6d |on P_((c|har *nam|
|000062d0| 65 29 29 3b 0a 58 76 6f | 69 64 20 66 69 65 6c 64 |e));.Xvo|id field|
|000062e0| 5f 6c 6f 67 20 50 5f 28 | 28 57 69 64 67 65 74 20 |_log P_(|(Widget |
|000062f0| 77 2c 20 64 6f 75 62 6c | 65 20 76 61 6c 75 65 2c |w, doubl|e value,|
|00006300| 20 69 6e 74 20 6c 6f 67 | 76 2c 20 63 68 61 72 20 | int log|v, char |
|00006310| 2a 73 74 72 29 29 3b 0a | 58 76 6f 69 64 20 70 72 |*str));.|Xvoid pr|
|00006320| 6f 6d 70 74 5f 6d 61 70 | 5f 63 62 20 50 5f 28 28 |ompt_map|_cb P_((|
|00006330| 57 69 64 67 65 74 20 77 | 2c 20 58 74 50 6f 69 6e |Widget w|, XtPoin|
|00006340| 74 65 72 20 63 6c 69 65 | 6e 74 2c 20 58 74 50 6f |ter clie|nt, XtPo|
|00006350| 69 6e 74 65 72 20 63 61 | 6c 6c 29 29 3b 0a 58 73 |inter ca|ll));.Xs|
|00006360| 74 61 74 69 63 20 76 6f | 69 64 20 67 65 74 5f 63 |tatic vo|id get_c|
|00006370| 6f 6c 6f 72 5f 72 65 73 | 6f 75 72 63 65 20 50 5f |olor_res|ource P_|
|00006380| 28 28 57 69 64 67 65 74 | 20 77 2c 20 63 68 61 72 |((Widget| w, char|
|00006390| 20 2a 63 6e 61 6d 65 2c | 20 50 69 78 65 6c 20 2a | *cname,| Pixel *|
|000063a0| 70 29 29 3b 0a 58 73 74 | 61 74 69 63 20 67 65 74 |p));.Xst|atic get|
|000063b0| 5f 76 69 65 77 73 5f 66 | 6f 6e 74 20 50 5f 28 28 |_views_f|ont P_((|
|000063c0| 44 69 73 70 6c 61 79 20 | 2a 64 73 70 2c 20 46 6f |Display |*dsp, Fo|
|000063d0| 6e 74 20 2a 66 70 29 29 | 3b 0a 58 76 6f 69 64 20 |nt *fp))|;.Xvoid |
|000063e0| 6f 62 6a 5f 70 69 63 6b | 67 63 20 50 5f 28 28 4f |obj_pick|gc P_((O|
|000063f0| 62 6a 20 2a 6f 70 2c 20 | 57 69 64 67 65 74 20 77 |bj *op, |Widget w|
+--------+-------------------------+-------------------------+--------+--------+