home *** CD-ROM | disk | FTP | other *** search
open in: 
MacOS 8.1
|
Win98
|
DOS
browse contents |
view JSON data
|
view as text
This file was processed as: SHell self-extracting ARchive
(archive/shar).
| Confidence | Program | Detection | Match Type | Support
|
|---|
| 100%
| dexvert
| SHell self-extracting ARchive (archive/shar)
| magic
| Supported |
| 1%
| dexvert
| Text File (text/txt)
| fallback
| Supported |
| 100%
| file
| ASCII text
| default
| |
| 100%
| checkBytes
| Printable ASCII
| default
| |
| 100%
| perlTextCheck
| Likely Text (Perl)
| default
| |
| 100%
| siegfried
| fmt/329 Shell Archive Format
| default
| |
| 100%
| detectItEasy
| Format: plain text[LF]
| default (weak)
|
|
hex view+--------+-------------------------+-------------------------+--------+--------+
|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 6d 69 73 63 0a 46 72 |.sources|.misc.Fr|
|00000020| 6f 6d 3a 20 64 66 73 40 | 64 6f 65 2e 63 61 72 6c |om: dfs@|doe.carl|
|00000030| 65 74 6f 6e 2e 63 61 20 | 28 44 61 76 69 64 20 46 |eton.ca |(David F|
|00000040| 2e 20 53 6b 6f 6c 6c 29 | 0a 53 75 62 6a 65 63 74 |. Skoll)|.Subject|
|00000050| 3a 20 76 33 39 69 31 31 | 38 3a 20 20 72 65 6d 69 |: v39i11|8: remi|
|00000060| 6e 64 20 2d 20 41 20 72 | 65 70 6c 61 63 65 6d 65 |nd - A r|eplaceme|
|00000070| 6e 74 20 66 6f 72 20 63 | 61 6c 65 6e 64 61 72 2c |nt for c|alendar,|
|00000080| 20 50 61 74 63 68 30 39 | 64 2f 34 0a 4d 65 73 73 | Patch09|d/4.Mess|
|00000090| 61 67 65 2d 49 44 3a 20 | 3c 31 39 39 33 4f 63 74 |age-ID: |<1993Oct|
|000000a0| 34 2e 32 30 34 30 31 33 | 2e 38 32 32 30 40 73 70 |4.204013|.8220@sp|
|000000b0| 61 72 6b 79 2e 73 74 65 | 72 6c 69 6e 67 2e 63 6f |arky.ste|rling.co|
|000000c0| 6d 3e 0a 58 2d 4d 64 34 | 2d 53 69 67 6e 61 74 75 |m>.X-Md4|-Signatu|
|000000d0| 72 65 3a 20 62 63 32 30 | 37 31 36 37 32 63 34 62 |re: bc20|71672c4b|
|000000e0| 62 65 32 32 34 31 32 39 | 38 62 65 61 33 64 63 30 |be224129|8bea3dc0|
|000000f0| 66 35 63 33 0a 53 65 6e | 64 65 72 3a 20 6b 65 6e |f5c3.Sen|der: ken|
|00000100| 74 40 73 70 61 72 6b 79 | 2e 73 74 65 72 6c 69 6e |t@sparky|.sterlin|
|00000110| 67 2e 63 6f 6d 20 28 4b | 65 6e 74 20 4c 61 6e 64 |g.com (K|ent Land|
|00000120| 66 69 65 6c 64 29 0a 4f | 72 67 61 6e 69 7a 61 74 |field).O|rganizat|
|00000130| 69 6f 6e 3a 20 44 65 70 | 74 2e 20 6f 66 20 45 6c |ion: Dep|t. of El|
|00000140| 65 63 74 72 6f 6e 69 63 | 73 2c 20 43 61 72 6c 65 |ectronic|s, Carle|
|00000150| 74 6f 6e 20 55 6e 69 76 | 65 72 73 69 74 79 0a 44 |ton Univ|ersity.D|
|00000160| 61 74 65 3a 20 4d 6f 6e | 2c 20 34 20 4f 63 74 20 |ate: Mon|, 4 Oct |
|00000170| 31 39 39 33 20 32 30 3a | 34 30 3a 31 33 20 47 4d |1993 20:|40:13 GM|
|00000180| 54 0a 41 70 70 72 6f 76 | 65 64 3a 20 6b 65 6e 74 |T.Approv|ed: kent|
|00000190| 40 73 70 61 72 6b 79 2e | 73 74 65 72 6c 69 6e 67 |@sparky.|sterling|
|000001a0| 2e 63 6f 6d 0a 0a 53 75 | 62 6d 69 74 74 65 64 2d |.com..Su|bmitted-|
|000001b0| 62 79 3a 20 64 66 73 40 | 64 6f 65 2e 63 61 72 6c |by: dfs@|doe.carl|
|000001c0| 65 74 6f 6e 2e 63 61 20 | 28 44 61 76 69 64 20 46 |eton.ca |(David F|
|000001d0| 2e 20 53 6b 6f 6c 6c 29 | 0a 50 6f 73 74 69 6e 67 |. Skoll)|.Posting|
|000001e0| 2d 6e 75 6d 62 65 72 3a | 20 56 6f 6c 75 6d 65 20 |-number:| Volume |
|000001f0| 33 39 2c 20 49 73 73 75 | 65 20 31 31 38 0a 41 72 |39, Issu|e 118.Ar|
|00000200| 63 68 69 76 65 2d 6e 61 | 6d 65 3a 20 72 65 6d 69 |chive-na|me: remi|
|00000210| 6e 64 2f 70 61 74 63 68 | 30 39 64 0a 45 6e 76 69 |nd/patch|09d.Envi|
|00000220| 72 6f 6e 6d 65 6e 74 3a | 20 55 4e 49 58 2c 20 4d |ronment:| UNIX, M|
|00000230| 53 2d 44 4f 53 2c 20 4f | 53 2f 32 0a 50 61 74 63 |S-DOS, O|S/2.Patc|
|00000240| 68 2d 54 6f 3a 20 72 65 | 6d 69 6e 64 3a 20 56 6f |h-To: re|mind: Vo|
|00000250| 6c 75 6d 65 20 33 33 2c | 20 49 73 73 75 65 20 35 |lume 33,| Issue 5|
|00000260| 38 2d 36 39 0a 0a 23 21 | 20 2f 62 69 6e 2f 73 68 |8-69..#!| /bin/sh|
|00000270| 0a 23 20 54 68 69 73 20 | 69 73 20 61 20 73 68 65 |.# This |is a she|
|00000280| 6c 6c 20 61 72 63 68 69 | 76 65 2e 20 20 52 65 6d |ll archi|ve. Rem|
|00000290| 6f 76 65 20 61 6e 79 74 | 68 69 6e 67 20 62 65 66 |ove anyt|hing bef|
|000002a0| 6f 72 65 20 74 68 69 73 | 20 6c 69 6e 65 2c 20 74 |ore this| line, t|
|000002b0| 68 65 6e 20 66 65 65 64 | 20 69 74 0a 23 20 69 6e |hen feed| it.# in|
|000002c0| 74 6f 20 61 20 73 68 65 | 6c 6c 20 76 69 61 20 22 |to a she|ll via "|
|000002d0| 73 68 20 66 69 6c 65 22 | 20 6f 72 20 73 69 6d 69 |sh file"| or simi|
|000002e0| 6c 61 72 2e 20 20 54 6f | 20 6f 76 65 72 77 72 69 |lar. To| overwri|
|000002f0| 74 65 20 65 78 69 73 74 | 69 6e 67 20 66 69 6c 65 |te exist|ing file|
|00000300| 73 2c 0a 23 20 74 79 70 | 65 20 22 73 68 20 66 69 |s,.# typ|e "sh fi|
|00000310| 6c 65 20 2d 63 22 2e 0a | 23 20 43 6f 6e 74 65 6e |le -c"..|# Conten|
|00000320| 74 73 3a 20 20 6d 6f 6f | 6e 2e 63 0a 23 20 57 72 |ts: moo|n.c.# Wr|
|00000330| 61 70 70 65 64 20 62 79 | 20 6b 65 6e 74 40 73 70 |apped by| kent@sp|
|00000340| 61 72 6b 79 20 6f 6e 20 | 4d 6f 6e 20 4f 63 74 20 |arky on |Mon Oct |
|00000350| 20 34 20 31 35 3a 30 33 | 3a 33 31 20 31 39 39 33 | 4 15:03|:31 1993|
|00000360| 0a 50 41 54 48 3d 2f 62 | 69 6e 3a 2f 75 73 72 2f |.PATH=/b|in:/usr/|
|00000370| 62 69 6e 3a 2f 75 73 72 | 2f 75 63 62 3a 2f 75 73 |bin:/usr|/ucb:/us|
|00000380| 72 2f 6c 6f 63 61 6c 2f | 62 69 6e 3a 2f 75 73 72 |r/local/|bin:/usr|
|00000390| 2f 6c 62 69 6e 20 3b 20 | 65 78 70 6f 72 74 20 50 |/lbin ; |export P|
|000003a0| 41 54 48 0a 65 63 68 6f | 20 49 66 20 74 68 69 73 |ATH.echo| If this|
|000003b0| 20 61 72 63 68 69 76 65 | 20 69 73 20 63 6f 6d 70 | archive| is comp|
|000003c0| 6c 65 74 65 2c 20 79 6f | 75 20 77 69 6c 6c 20 73 |lete, yo|u will s|
|000003d0| 65 65 20 74 68 65 20 66 | 6f 6c 6c 6f 77 69 6e 67 |ee the f|ollowing|
|000003e0| 20 6d 65 73 73 61 67 65 | 3a 0a 65 63 68 6f 20 27 | message|:.echo '|
|000003f0| 20 20 20 20 20 20 20 20 | 20 20 22 73 68 61 72 3a | | "shar:|
|00000400| 20 45 6e 64 20 6f 66 20 | 61 72 63 68 69 76 65 20 | End of |archive |
|00000410| 34 20 28 6f 66 20 34 29 | 2e 22 27 0a 69 66 20 74 |4 (of 4)|."'.if t|
|00000420| 65 73 74 20 2d 66 20 27 | 6d 6f 6f 6e 2e 63 27 20 |est -f '|moon.c' |
|00000430| 2d 61 20 22 24 7b 31 7d | 22 20 21 3d 20 22 2d 63 |-a "${1}|" != "-c|
|00000440| 22 20 3b 20 74 68 65 6e | 20 0a 20 20 65 63 68 6f |" ; then| . echo|
|00000450| 20 73 68 61 72 3a 20 57 | 69 6c 6c 20 6e 6f 74 20 | shar: W|ill not |
|00000460| 63 6c 6f 62 62 65 72 20 | 65 78 69 73 74 69 6e 67 |clobber |existing|
|00000470| 20 66 69 6c 65 20 5c 22 | 27 6d 6f 6f 6e 2e 63 27 | file \"|'moon.c'|
|00000480| 5c 22 0a 65 6c 73 65 0a | 20 20 65 63 68 6f 20 73 |\".else.| echo s|
|00000490| 68 61 72 3a 20 45 78 74 | 72 61 63 74 69 6e 67 20 |har: Ext|racting |
|000004a0| 5c 22 27 6d 6f 6f 6e 2e | 63 27 5c 22 20 5c 28 32 |\"'moon.|c'\" \(2|
|000004b0| 32 34 36 31 20 63 68 61 | 72 61 63 74 65 72 73 5c |2461 cha|racters\|
|000004c0| 29 0a 20 20 73 65 64 20 | 22 73 2f 5e 58 2f 2f 22 |). sed |"s/^X//"|
|000004d0| 20 3e 27 6d 6f 6f 6e 2e | 63 27 20 3c 3c 27 45 4e | >'moon.|c' <<'EN|
|000004e0| 44 5f 4f 46 5f 46 49 4c | 45 27 0a 58 2f 2a 2a 2a |D_OF_FIL|E'.X/***|
|000004f0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000500| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000510| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000520| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2f 0a 58 2f |********|****/.X/|
|00000530| 2a 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |* | |
|00000540| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00000550| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00000560| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 2a 2f | | */|
|00000570| 0a 58 2f 2a 20 20 4d 4f | 4f 4e 2e 43 20 20 20 20 |.X/* MO|ON.C |
|00000580| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00000590| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000005a0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000005b0| 20 2a 2f 0a 58 2f 2a 20 | 20 20 20 20 20 20 20 20 | */.X/* | |
|000005c0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000005d0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000005e0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000005f0| 20 20 20 20 2a 2f 0a 58 | 2f 2a 20 20 43 61 6c 63 | */.X|/* Calc|
|00000600| 75 6c 61 74 69 6f 6e 73 | 20 66 6f 72 20 66 69 67 |ulations| for fig|
|00000610| 75 72 69 6e 67 20 6f 75 | 74 20 6d 6f 6f 6e 20 70 |uring ou|t moon p|
|00000620| 68 61 73 65 73 2e 20 20 | 20 20 20 20 20 20 20 20 |hases. | |
|00000630| 20 20 20 20 20 20 20 2a | 2f 0a 58 2f 2a 20 20 20 | *|/.X/* |
|00000640| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00000650| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00000660| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00000670| 20 20 20 20 20 20 20 20 | 20 20 2a 2f 0a 58 2f 2a | | */.X/*|
|00000680| 20 20 54 68 69 73 20 66 | 69 6c 65 20 69 73 20 70 | This f|ile is p|
|00000690| 61 72 74 20 6f 66 20 52 | 45 4d 49 4e 44 2e 20 20 |art of R|EMIND. |
|000006a0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000006b0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 2a 2f 0a | | */.|
|000006c0| 58 2f 2a 20 20 43 6f 70 | 79 72 69 67 68 74 20 28 |X/* Cop|yright (|
|000006d0| 43 29 20 31 39 39 32 2c | 20 31 39 39 33 20 62 79 |C) 1992,| 1993 by|
|000006e0| 20 44 61 76 69 64 20 46 | 2e 20 53 6b 6f 6c 6c 2e | David F|. Skoll.|
|000006f0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00000700| 2a 2f 0a 58 2f 2a 20 20 | 20 20 20 20 20 20 20 20 |*/.X/* | |
|00000710| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00000720| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00000730| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00000740| 20 20 20 2a 2f 0a 58 2f | 2a 2a 2a 2a 2a 2a 2a 2a | */.X/|********|
|00000750| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000760| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000770| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00000780| 2a 2a 2a 2a 2a 2a 2a 2f | 0a 58 0a 58 2f 2a 20 41 |*******/|.X.X/* A|
|00000790| 6c 6c 20 6f 66 20 74 68 | 65 73 65 20 72 6f 75 74 |ll of th|ese rout|
|000007a0| 69 6e 65 73 20 77 65 72 | 65 20 61 64 61 70 74 65 |ines wer|e adapte|
|000007b0| 64 20 66 72 6f 6d 20 74 | 68 65 20 70 72 6f 67 72 |d from t|he progr|
|000007c0| 61 6d 20 22 6d 6f 6f 6e | 74 6f 6f 6c 22 0a 58 20 |am "moon|tool".X |
|000007d0| 20 20 62 79 20 4a 6f 68 | 6e 20 57 61 6c 6b 65 72 | by Joh|n Walker|
|000007e0| 2c 20 46 65 62 72 75 61 | 72 79 20 31 39 38 38 2e |, Februa|ry 1988.|
|000007f0| 20 20 48 65 72 65 27 73 | 20 74 68 65 20 62 6c 75 | Here's| the blu|
|00000800| 72 62 20 66 72 6f 6d 20 | 6d 6f 6f 6e 74 6f 6f 6c |rb from |moontool|
|00000810| 3a 0a 58 0a 58 20 20 20 | 2e 2e 2e 20 54 68 65 20 |:.X.X |... The |
|00000820| 69 6e 66 6f 72 6d 61 74 | 69 6f 6e 20 69 73 20 67 |informat|ion is g|
|00000830| 65 6e 65 72 61 6c 6c 79 | 20 61 63 63 75 72 61 74 |enerally| accurat|
|00000840| 65 20 74 6f 20 77 69 74 | 68 69 6e 20 74 65 6e 0a |e to wit|hin ten.|
|00000850| 58 20 20 20 6d 69 6e 75 | 74 65 73 2e 0a 58 0a 58 |X minu|tes..X.X|
|00000860| 20 20 20 54 68 65 20 61 | 6c 67 6f 72 69 74 68 6d | The a|lgorithm|
|00000870| 73 20 75 73 65 64 20 69 | 6e 20 74 68 69 73 20 70 |s used i|n this p|
|00000880| 72 6f 67 72 61 6d 20 74 | 6f 20 63 61 6c 63 75 6c |rogram t|o calcul|
|00000890| 61 74 65 20 74 68 65 20 | 70 6f 73 69 74 69 6f 6e |ate the |position|
|000008a0| 73 20 53 75 6e 20 61 6e | 64 0a 58 20 20 20 4d 6f |s Sun an|d.X Mo|
|000008b0| 6f 6e 20 61 73 20 73 65 | 65 6e 20 66 72 6f 6d 20 |on as se|en from |
|000008c0| 74 68 65 20 45 61 72 74 | 68 20 61 72 65 20 67 69 |the Eart|h are gi|
|000008d0| 76 65 6e 20 69 6e 20 74 | 68 65 20 62 6f 6f 6b 20 |ven in t|he book |
|000008e0| 22 50 72 61 63 74 69 63 | 61 6c 20 41 73 74 72 6f |"Practic|al Astro|
|000008f0| 6e 6f 6d 79 0a 58 20 20 | 20 57 69 74 68 20 59 6f |nomy.X | With Yo|
|00000900| 75 72 20 43 61 6c 63 75 | 6c 61 74 6f 72 22 20 62 |ur Calcu|lator" b|
|00000910| 79 20 50 65 74 65 72 20 | 44 75 66 66 65 74 74 2d |y Peter |Duffett-|
|00000920| 53 6d 69 74 68 2c 20 53 | 65 63 6f 6e 64 20 45 64 |Smith, S|econd Ed|
|00000930| 69 74 69 6f 6e 2c 0a 58 | 20 20 20 43 61 6d 62 72 |ition,.X| Cambr|
|00000940| 69 64 67 65 20 55 6e 69 | 76 65 72 73 69 74 79 20 |idge Uni|versity |
|00000950| 50 72 65 73 73 2c 20 31 | 39 38 31 2e 20 49 67 6e |Press, 1|981. Ign|
|00000960| 6f 72 65 20 74 68 65 20 | 77 6f 72 64 20 22 43 61 |ore the |word "Ca|
|00000970| 6c 63 75 6c 61 74 6f 72 | 22 20 69 6e 20 74 68 65 |lculator|" in the|
|00000980| 0a 58 20 20 20 74 69 74 | 6c 65 3b 20 74 68 69 73 |.X tit|le; this|
|00000990| 20 69 73 20 61 6e 20 65 | 73 73 65 6e 74 69 61 6c | is an e|ssential|
|000009a0| 20 72 65 66 65 72 65 6e | 63 65 20 69 66 20 79 6f | referen|ce if yo|
|000009b0| 75 27 72 65 20 69 6e 74 | 65 72 65 73 74 65 64 20 |u're int|erested |
|000009c0| 69 6e 0a 58 20 20 20 64 | 65 76 65 6c 6f 70 69 6e |in.X d|evelopin|
|000009d0| 67 20 73 6f 66 74 77 61 | 72 65 20 77 68 69 63 68 |g softwa|re which|
|000009e0| 20 63 61 6c 63 75 6c 61 | 74 65 73 20 70 6c 61 6e | calcula|tes plan|
|000009f0| 65 74 61 72 79 20 70 6f | 73 69 74 69 6f 6e 73 2c |etary po|sitions,|
|00000a00| 20 6f 72 62 69 74 73 2c | 0a 58 20 20 20 65 63 6c | orbits,|.X ecl|
|00000a10| 69 70 73 65 73 2c 20 61 | 6e 64 20 74 68 65 20 6c |ipses, a|nd the l|
|00000a20| 69 6b 65 2e 20 49 66 20 | 79 6f 75 27 72 65 20 69 |ike. If |you're i|
|00000a30| 6e 74 65 72 65 73 74 65 | 64 20 69 6e 20 70 75 72 |ntereste|d in pur|
|00000a40| 73 75 69 6e 67 20 73 75 | 63 68 0a 58 20 20 20 70 |suing su|ch.X p|
|00000a50| 72 6f 67 72 61 6d 6d 69 | 6e 67 2c 20 79 6f 75 20 |rogrammi|ng, you |
|00000a60| 73 68 6f 75 6c 64 20 61 | 6c 73 6f 20 6f 62 74 61 |should a|lso obta|
|00000a70| 69 6e 3a 0a 58 0a 58 20 | 20 20 22 41 73 74 72 6f |in:.X.X | "Astro|
|00000a80| 6e 6f 6d 69 63 61 6c 20 | 46 6f 72 6d 75 6c 61 65 |nomical |Formulae|
|00000a90| 20 66 6f 72 20 43 61 6c | 63 75 6c 61 74 6f 72 73 | for Cal|culators|
|00000aa0| 22 20 62 79 20 4a 65 61 | 6e 20 4d 65 65 75 73 2c |" by Jea|n Meeus,|
|00000ab0| 20 54 68 69 72 64 20 45 | 64 69 74 69 6f 6e 2c 0a | Third E|dition,.|
|00000ac0| 58 20 20 20 57 69 6c 6c | 6d 61 6e 6e 2d 42 65 6c |X Will|mann-Bel|
|00000ad0| 6c 2c 20 31 39 38 35 2e | 20 41 20 6d 75 73 74 2d |l, 1985.| A must-|
|00000ae0| 68 61 76 65 2e 0a 58 0a | 58 20 20 20 22 50 6c 61 |have..X.|X "Pla|
|00000af0| 6e 65 74 61 72 79 20 50 | 72 6f 67 72 61 6d 73 20 |netary P|rograms |
|00000b00| 61 6e 64 20 54 61 62 6c | 65 73 20 66 72 6f 6d 20 |and Tabl|es from |
|00000b10| 2d 34 30 30 30 20 74 6f | 20 2b 32 38 30 30 22 20 |-4000 to| +2800" |
|00000b20| 62 79 20 50 69 65 72 72 | 65 0a 58 20 20 20 42 72 |by Pierr|e.X Br|
|00000b30| 65 74 61 67 6e 6f 6e 20 | 61 6e 64 20 4a 65 61 6e |etagnon |and Jean|
|00000b40| 2d 4c 6f 75 69 73 20 53 | 69 6d 6f 6e 2c 20 57 69 |-Louis S|imon, Wi|
|00000b50| 6c 6c 6d 61 6e 6e 2d 42 | 65 6c 6c 2c 20 31 39 38 |llmann-B|ell, 198|
|00000b60| 36 2e 20 49 66 20 79 6f | 75 20 77 61 6e 74 20 74 |6. If yo|u want t|
|00000b70| 68 65 0a 58 20 20 20 75 | 74 6d 6f 73 74 20 28 6f |he.X u|tmost (o|
|00000b80| 75 74 73 69 64 65 20 6f | 66 20 4a 50 4c 29 20 61 |utside o|f JPL) a|
|00000b90| 63 63 75 72 61 63 79 20 | 66 6f 72 20 74 68 65 20 |ccuracy |for the |
|00000ba0| 70 6c 61 6e 65 74 73 2c | 20 69 74 27 73 20 68 65 |planets,| it's he|
|00000bb0| 72 65 2e 0a 58 0a 58 20 | 20 20 22 43 65 6c 65 73 |re..X.X | "Celes|
|00000bc0| 74 69 61 6c 20 42 41 53 | 49 43 22 20 62 79 20 45 |tial BAS|IC" by E|
|00000bd0| 72 69 63 20 42 75 72 67 | 65 73 73 2c 20 52 65 76 |ric Burg|ess, Rev|
|00000be0| 69 73 65 64 20 45 64 69 | 74 69 6f 6e 2c 20 53 79 |ised Edi|tion, Sy|
|00000bf0| 62 65 78 2c 20 31 39 38 | 35 2e 20 56 65 72 79 0a |bex, 198|5. Very.|
|00000c00| 58 20 20 20 63 6f 6f 6b | 62 6f 6f 6b 20 6f 72 69 |X cook|book ori|
|00000c10| 65 6e 74 65 64 2c 20 61 | 6e 64 20 6d 61 6e 79 20 |ented, a|nd many |
|00000c20| 6f 66 20 74 68 65 20 61 | 6c 67 6f 72 69 74 68 6d |of the a|lgorithm|
|00000c30| 73 20 61 72 65 20 68 61 | 72 64 20 74 6f 20 64 69 |s are ha|rd to di|
|00000c40| 67 20 6f 75 74 20 6f 66 | 0a 58 20 20 20 74 68 65 |g out of|.X the|
|00000c50| 20 74 75 72 67 69 64 20 | 42 41 53 49 43 20 63 6f | turgid |BASIC co|
|00000c60| 64 65 2c 20 62 75 74 20 | 79 6f 75 27 6c 6c 20 70 |de, but |you'll p|
|00000c70| 72 6f 62 61 62 6c 79 20 | 77 61 6e 74 20 69 74 20 |robably |want it |
|00000c80| 61 6e 79 77 61 79 2e 0a | 58 0a 58 20 20 20 4d 61 |anyway..|X.X Ma|
|00000c90| 6e 79 20 6f 66 20 74 68 | 65 73 65 20 72 65 66 65 |ny of th|ese refe|
|00000ca0| 72 65 6e 63 65 73 20 63 | 61 6e 20 62 65 20 6f 62 |rences c|an be ob|
|00000cb0| 74 61 69 6e 65 64 20 66 | 72 6f 6d 20 57 69 6c 6c |tained f|rom Will|
|00000cc0| 6d 61 6e 6e 2d 42 65 6c | 6c 2c 20 50 2e 4f 2e 20 |mann-Bel|l, P.O. |
|00000cd0| 42 6f 78 0a 58 20 20 20 | 33 35 30 32 35 2c 20 52 |Box.X |35025, R|
|00000ce0| 69 63 68 6d 6f 6e 64 2c | 20 56 41 20 32 33 32 33 |ichmond,| VA 2323|
|00000cf0| 35 2c 20 55 53 41 2e 20 | 50 68 6f 6e 65 3a 20 28 |5, USA. |Phone: (|
|00000d00| 38 30 34 29 20 33 32 30 | 2d 37 30 31 36 2e 20 49 |804) 320|-7016. I|
|00000d10| 6e 20 61 64 64 69 74 69 | 6f 6e 0a 58 20 20 20 74 |n additi|on.X t|
|00000d20| 6f 20 74 68 65 69 72 20 | 6f 77 6e 20 70 75 62 6c |o their |own publ|
|00000d30| 69 63 61 74 69 6f 6e 73 | 2c 20 74 68 65 79 20 73 |ications|, they s|
|00000d40| 74 6f 63 6b 20 6d 6f 73 | 74 20 6f 66 20 74 68 65 |tock mos|t of the|
|00000d50| 20 73 74 61 6e 64 61 72 | 64 20 72 65 66 65 72 65 | standar|d refere|
|00000d60| 6e 63 65 73 0a 58 20 20 | 20 66 6f 72 20 6d 61 74 |nces.X | for mat|
|00000d70| 68 65 6d 61 74 69 63 61 | 6c 20 61 6e 64 20 70 6f |hematica|l and po|
|00000d80| 73 69 74 69 6f 6e 61 6c | 20 61 73 74 72 6f 6e 6f |sitional| astrono|
|00000d90| 6d 79 2e 0a 58 0a 58 20 | 20 20 54 68 69 73 20 70 |my..X.X | This p|
|00000da0| 72 6f 67 72 61 6d 20 77 | 61 73 20 77 72 69 74 74 |rogram w|as writt|
|00000db0| 65 6e 20 62 79 3a 0a 58 | 0a 58 20 20 20 20 20 20 |en by:.X|.X |
|00000dc0| 4a 6f 68 6e 20 57 61 6c | 6b 65 72 0a 58 20 20 20 |John Wal|ker.X |
|00000dd0| 20 20 20 41 75 74 6f 64 | 65 73 6b 2c 20 49 6e 63 | Autod|esk, Inc|
|00000de0| 2e 0a 58 20 20 20 20 20 | 20 32 33 32 30 20 4d 61 |..X | 2320 Ma|
|00000df0| 72 69 6e 73 68 69 70 20 | 57 61 79 0a 58 20 20 20 |rinship |Way.X |
|00000e00| 20 20 20 53 61 75 73 61 | 6c 69 74 6f 2c 20 43 41 | Sausa|lito, CA|
|00000e10| 20 39 34 39 36 35 0a 58 | 20 20 20 20 20 20 28 34 | 94965.X| (4|
|00000e20| 31 35 29 20 33 33 32 2d | 32 33 34 34 20 45 78 74 |15) 332-|2344 Ext|
|00000e30| 2e 20 38 32 39 0a 58 0a | 58 20 20 20 20 20 20 55 |. 829.X.|X U|
|00000e40| 73 65 6e 65 74 3a 20 7b | 73 75 6e 21 77 65 6c 6c |senet: {|sun!well|
|00000e50| 7d 21 61 63 61 64 21 6b | 65 6c 76 69 6e 0a 58 0a |}!acad!k|elvin.X.|
|00000e60| 58 20 20 20 54 68 69 73 | 20 70 72 6f 67 72 61 6d |X This| program|
|00000e70| 20 69 73 20 69 6e 20 74 | 68 65 20 70 75 62 6c 69 | is in t|he publi|
|00000e80| 63 20 64 6f 6d 61 69 6e | 3a 20 22 44 6f 20 77 68 |c domain|: "Do wh|
|00000e90| 61 74 20 74 68 6f 75 20 | 77 69 6c 74 20 73 68 61 |at thou |wilt sha|
|00000ea0| 6c 6c 20 62 65 20 74 68 | 65 0a 58 20 20 20 77 68 |ll be th|e.X wh|
|00000eb0| 6f 6c 65 20 6f 66 20 74 | 68 65 20 6c 61 77 22 2e |ole of t|he law".|
|00000ec0| 20 49 27 64 20 61 70 70 | 72 65 63 69 61 74 65 20 | I'd app|reciate |
|00000ed0| 72 65 63 65 69 76 69 6e | 67 20 61 6e 79 20 62 75 |receivin|g any bu|
|00000ee0| 67 20 66 69 78 65 73 20 | 61 6e 64 2f 6f 72 0a 58 |g fixes |and/or.X|
|00000ef0| 20 20 20 65 6e 68 61 6e | 63 65 6d 65 6e 74 73 2c | enhan|cements,|
|00000f00| 20 77 68 69 63 68 20 49 | 27 6c 6c 20 69 6e 63 6f | which I|'ll inco|
|00000f10| 72 70 6f 72 61 74 65 20 | 69 6e 20 66 75 74 75 72 |rporate |in futur|
|00000f20| 65 20 76 65 72 73 69 6f | 6e 73 20 6f 66 20 74 68 |e versio|ns of th|
|00000f30| 65 0a 58 20 20 20 70 72 | 6f 67 72 61 6d 2e 20 50 |e.X pr|ogram. P|
|00000f40| 6c 65 61 73 65 20 6c 65 | 61 76 65 20 74 68 65 20 |lease le|ave the |
|00000f50| 6f 72 69 67 69 6e 61 6c | 20 61 74 74 72 69 62 75 |original| attribu|
|00000f60| 74 69 6f 6e 20 69 6e 66 | 6f 72 6d 61 74 69 6f 6e |tion inf|ormation|
|00000f70| 20 69 6e 74 61 63 74 20 | 73 6f 0a 58 20 20 20 74 | intact |so.X t|
|00000f80| 68 61 74 20 63 72 65 64 | 69 74 20 61 6e 64 20 62 |hat cred|it and b|
|00000f90| 6c 61 6d 65 20 6d 61 79 | 20 62 65 20 70 72 6f 70 |lame may| be prop|
|00000fa0| 65 72 6c 79 20 61 70 70 | 6f 72 74 69 6f 6e 65 64 |erly app|ortioned|
|00000fb0| 2e 0a 58 0a 58 2a 2f 0a | 58 23 69 6e 63 6c 75 64 |..X.X*/.|X#includ|
|00000fc0| 65 20 22 63 6f 6e 66 69 | 67 2e 68 22 0a 58 23 69 |e "confi|g.h".X#i|
|00000fd0| 66 64 65 66 20 48 41 56 | 45 5f 53 54 44 4c 49 42 |fdef HAV|E_STDLIB|
|00000fe0| 5f 48 0a 58 23 69 6e 63 | 6c 75 64 65 20 3c 73 74 |_H.X#inc|lude <st|
|00000ff0| 64 6c 69 62 2e 68 3e 0a | 58 23 65 6e 64 69 66 0a |dlib.h>.|X#endif.|
|00001000| 58 23 69 6e 63 6c 75 64 | 65 20 3c 73 74 64 69 6f |X#includ|e <stdio|
|00001010| 2e 68 3e 0a 58 23 69 6e | 63 6c 75 64 65 20 3c 6d |.h>.X#in|clude <m|
|00001020| 61 74 68 2e 68 3e 0a 58 | 23 69 6e 63 6c 75 64 65 |ath.h>.X|#include|
|00001030| 20 3c 74 69 6d 65 2e 68 | 3e 0a 58 23 69 6e 63 6c | <time.h|>.X#incl|
|00001040| 75 64 65 20 22 74 79 70 | 65 73 2e 68 22 0a 58 23 |ude "typ|es.h".X#|
|00001050| 69 6e 63 6c 75 64 65 20 | 22 70 72 6f 74 6f 73 2e |include |"protos.|
|00001060| 68 22 0a 58 23 69 6e 63 | 6c 75 64 65 20 22 65 78 |h".X#inc|lude "ex|
|00001070| 70 72 2e 68 22 0a 58 23 | 69 6e 63 6c 75 64 65 20 |pr.h".X#|include |
|00001080| 22 67 6c 6f 62 61 6c 73 | 2e 68 22 0a 58 23 69 6e |"globals|.h".X#in|
|00001090| 63 6c 75 64 65 20 22 65 | 72 72 2e 68 22 0a 58 0a |clude "e|rr.h".X.|
|000010a0| 58 2f 2a 20 46 75 6e 63 | 74 69 6f 6e 20 70 72 6f |X/* Func|tion pro|
|000010b0| 74 6f 74 79 70 65 73 20 | 2a 2f 0a 58 50 52 49 56 |totypes |*/.XPRIV|
|000010c0| 41 54 45 20 6c 6f 6e 67 | 20 6a 64 61 74 65 20 41 |ATE long| jdate A|
|000010d0| 52 47 53 28 28 69 6e 74 | 20 79 2c 20 69 6e 74 20 |RGS((int| y, int |
|000010e0| 6d 6f 6e 2c 20 69 6e 74 | 20 64 61 79 29 29 3b 0a |mon, int| day));.|
|000010f0| 58 50 52 49 56 41 54 45 | 20 64 6f 75 62 6c 65 20 |XPRIVATE| double |
|00001100| 6a 74 69 6d 65 20 41 52 | 47 53 28 28 69 6e 74 20 |jtime AR|GS((int |
|00001110| 79 2c 20 69 6e 74 20 6d | 6f 6e 2c 20 69 6e 74 20 |y, int m|on, int |
|00001120| 64 61 79 2c 20 69 6e 74 | 20 68 6f 75 72 2c 20 69 |day, int| hour, i|
|00001130| 6e 74 20 6d 69 6e 2c 20 | 69 6e 74 20 73 65 63 29 |nt min, |int sec)|
|00001140| 29 3b 0a 58 50 52 49 56 | 41 54 45 20 76 6f 69 64 |);.XPRIV|ATE void|
|00001150| 20 6a 79 65 61 72 20 41 | 52 47 53 28 28 64 6f 75 | jyear A|RGS((dou|
|00001160| 62 6c 65 20 74 64 2c 20 | 69 6e 74 20 2a 79 79 2c |ble td, |int *yy,|
|00001170| 20 69 6e 74 20 2a 6d 6d | 2c 20 69 6e 74 20 2a 64 | int *mm|, int *d|
|00001180| 64 29 29 3b 0a 58 50 52 | 49 56 41 54 45 20 76 6f |d));.XPR|IVATE vo|
|00001190| 69 64 20 6a 68 6d 73 20 | 41 52 47 53 28 28 64 6f |id jhms |ARGS((do|
|000011a0| 75 62 6c 65 20 6a 2c 20 | 69 6e 74 20 2a 68 2c 20 |uble j, |int *h, |
|000011b0| 69 6e 74 20 2a 6d 2c 20 | 69 6e 74 20 2a 73 29 29 |int *m, |int *s))|
|000011c0| 3b 0a 58 50 52 49 56 41 | 54 45 20 64 6f 75 62 6c |;.XPRIVA|TE doubl|
|000011d0| 65 20 6d 65 61 6e 70 68 | 61 73 65 20 41 52 47 53 |e meanph|ase ARGS|
|000011e0| 28 28 64 6f 75 62 6c 65 | 20 73 64 61 74 65 2c 20 |((double| sdate, |
|000011f0| 64 6f 75 62 6c 65 20 70 | 68 61 73 65 2c 20 64 6f |double p|hase, do|
|00001200| 75 62 6c 65 20 2a 75 73 | 65 6b 29 29 3b 0a 58 50 |uble *us|ek));.XP|
|00001210| 52 49 56 41 54 45 20 64 | 6f 75 62 6c 65 20 74 72 |RIVATE d|ouble tr|
|00001220| 75 65 70 68 61 73 65 20 | 41 52 47 53 28 28 64 6f |uephase |ARGS((do|
|00001230| 75 62 6c 65 20 6b 2c 20 | 64 6f 75 62 6c 65 20 70 |uble k, |double p|
|00001240| 68 61 73 65 29 29 3b 0a | 58 50 52 49 56 41 54 45 |hase));.|XPRIVATE|
|00001250| 20 64 6f 75 62 6c 65 20 | 6b 65 70 6c 65 72 20 41 | double |kepler A|
|00001260| 52 47 53 28 28 64 6f 75 | 62 6c 65 20 6d 2c 20 64 |RGS((dou|ble m, d|
|00001270| 6f 75 62 6c 65 20 65 63 | 63 29 29 3b 0a 58 50 52 |ouble ec|c));.XPR|
|00001280| 49 56 41 54 45 20 64 6f | 75 62 6c 65 20 70 68 61 |IVATE do|uble pha|
|00001290| 73 65 20 41 52 47 53 28 | 28 64 6f 75 62 6c 65 2c |se ARGS(|(double,|
|000012a0| 20 64 6f 75 62 6c 65 20 | 2a 2c 20 64 6f 75 62 6c | double |*, doubl|
|000012b0| 65 20 2a 2c 20 64 6f 75 | 62 6c 65 20 2a 2c 20 64 |e *, dou|ble *, d|
|000012c0| 6f 75 62 6c 65 20 2a 2c | 20 64 6f 75 62 6c 65 20 |ouble *,| double |
|000012d0| 2a 2c 20 64 6f 75 62 6c | 65 20 2a 29 29 3b 0a 58 |*, doubl|e *));.X|
|000012e0| 0a 58 0a 58 2f 2a 20 20 | 41 73 74 72 6f 6e 6f 6d |.X.X/* |Astronom|
|000012f0| 69 63 61 6c 20 63 6f 6e | 73 74 61 6e 74 73 20 20 |ical con|stants |
|00001300| 2a 2f 0a 58 0a 58 23 64 | 65 66 69 6e 65 20 65 70 |*/.X.X#d|efine ep|
|00001310| 6f 63 68 09 20 20 20 20 | 32 34 34 34 32 33 38 2e |och. |2444238.|
|00001320| 35 09 20 20 20 2f 2a 20 | 31 39 38 30 20 4a 61 6e |5. /* |1980 Jan|
|00001330| 75 61 72 79 20 30 2e 30 | 20 2a 2f 0a 58 0a 58 2f |uary 0.0| */.X.X/|
|00001340| 2a 20 20 43 6f 6e 73 74 | 61 6e 74 73 20 64 65 66 |* Const|ants def|
|00001350| 69 6e 69 6e 67 20 74 68 | 65 20 53 75 6e 27 73 20 |ining th|e Sun's |
|00001360| 61 70 70 61 72 65 6e 74 | 20 6f 72 62 69 74 20 20 |apparent| orbit |
|00001370| 2a 2f 0a 58 0a 58 23 64 | 65 66 69 6e 65 20 65 6c |*/.X.X#d|efine el|
|00001380| 6f 6e 67 65 09 20 20 20 | 20 32 37 38 2e 38 33 33 |onge. | 278.833|
|00001390| 35 34 30 09 20 20 20 2f | 2a 20 45 63 6c 69 70 74 |540. /|* Eclipt|
|000013a0| 69 63 20 6c 6f 6e 67 69 | 74 75 64 65 20 6f 66 20 |ic longi|tude of |
|000013b0| 74 68 65 20 53 75 6e 0a | 58 09 09 09 09 20 20 20 |the Sun.|X.... |
|000013c0| 20 20 20 61 74 20 65 70 | 6f 63 68 20 31 39 38 30 | at ep|och 1980|
|000013d0| 2e 30 20 2a 2f 0a 58 23 | 64 65 66 69 6e 65 20 65 |.0 */.X#|define e|
|000013e0| 6c 6f 6e 67 70 09 20 20 | 20 20 32 38 32 2e 35 39 |longp. | 282.59|
|000013f0| 36 34 30 33 09 20 20 20 | 2f 2a 20 45 63 6c 69 70 |6403. |/* Eclip|
|00001400| 74 69 63 20 6c 6f 6e 67 | 69 74 75 64 65 20 6f 66 |tic long|itude of|
|00001410| 20 74 68 65 20 53 75 6e | 20 61 74 0a 58 09 09 09 | the Sun| at.X...|
|00001420| 09 20 20 20 20 20 20 70 | 65 72 69 67 65 65 20 2a |. p|erigee *|
|00001430| 2f 0a 58 23 64 65 66 69 | 6e 65 20 65 63 63 65 6e |/.X#defi|ne eccen|
|00001440| 74 20 20 20 20 20 20 30 | 2e 30 31 36 37 31 38 20 |t 0|.016718 |
|00001450| 20 20 20 20 20 20 2f 2a | 20 45 63 63 65 6e 74 72 | /*| Eccentr|
|00001460| 69 63 69 74 79 20 6f 66 | 20 45 61 72 74 68 27 73 |icity of| Earth's|
|00001470| 20 6f 72 62 69 74 20 2a | 2f 0a 58 23 64 65 66 69 | orbit *|/.X#defi|
|00001480| 6e 65 20 73 75 6e 73 6d | 61 78 20 20 20 20 20 31 |ne sunsm|ax 1|
|00001490| 2e 34 39 35 39 38 35 65 | 38 20 20 20 20 20 2f 2a |.495985e|8 /*|
|000014a0| 20 53 65 6d 69 2d 6d 61 | 6a 6f 72 20 61 78 69 73 | Semi-ma|jor axis|
|000014b0| 20 6f 66 20 45 61 72 74 | 68 27 73 20 6f 72 62 69 | of Eart|h's orbi|
|000014c0| 74 2c 20 6b 6d 20 2a 2f | 0a 58 23 64 65 66 69 6e |t, km */|.X#defin|
|000014d0| 65 20 73 75 6e 61 6e 67 | 73 69 7a 20 20 20 30 2e |e sunang|siz 0.|
|000014e0| 35 33 33 31 32 38 20 20 | 20 20 20 20 20 2f 2a 20 |533128 | /* |
|000014f0| 53 75 6e 27 73 20 61 6e | 67 75 6c 61 72 20 73 69 |Sun's an|gular si|
|00001500| 7a 65 2c 20 64 65 67 72 | 65 65 73 2c 20 61 74 0a |ze, degr|ees, at.|
|00001510| 58 09 09 09 09 20 20 20 | 20 20 20 73 65 6d 69 2d |X.... | semi-|
|00001520| 6d 61 6a 6f 72 20 61 78 | 69 73 20 64 69 73 74 61 |major ax|is dista|
|00001530| 6e 63 65 20 2a 2f 0a 58 | 0a 58 2f 2a 20 20 45 6c |nce */.X|.X/* El|
|00001540| 65 6d 65 6e 74 73 20 6f | 66 20 74 68 65 20 4d 6f |ements o|f the Mo|
|00001550| 6f 6e 27 73 20 6f 72 62 | 69 74 2c 20 65 70 6f 63 |on's orb|it, epoc|
|00001560| 68 20 31 39 38 30 2e 30 | 20 20 2a 2f 0a 58 0a 58 |h 1980.0| */.X.X|
|00001570| 23 64 65 66 69 6e 65 20 | 6d 6d 6c 6f 6e 67 20 20 |#define |mmlong |
|00001580| 20 20 20 20 36 34 2e 39 | 37 35 34 36 34 20 20 20 | 64.9|75464 |
|00001590| 20 20 20 2f 2a 20 4d 6f | 6f 6e 27 73 20 6d 65 61 | /* Mo|on's mea|
|000015a0| 6e 20 6c 6f 6e 69 67 69 | 74 75 64 65 20 61 74 20 |n lonigi|tude at |
|000015b0| 74 68 65 20 65 70 6f 63 | 68 20 2a 2f 0a 58 23 64 |the epoc|h */.X#d|
|000015c0| 65 66 69 6e 65 20 6d 6d | 6c 6f 6e 67 70 20 20 20 |efine mm|longp |
|000015d0| 20 20 33 34 39 2e 33 38 | 33 30 36 33 09 20 20 20 | 349.38|3063. |
|000015e0| 2f 2a 20 4d 65 61 6e 20 | 6c 6f 6e 67 69 74 75 64 |/* Mean |longitud|
|000015f0| 65 20 6f 66 20 74 68 65 | 20 70 65 72 69 67 65 65 |e of the| perigee|
|00001600| 20 61 74 20 74 68 65 0a | 58 09 09 09 09 20 20 20 | at the.|X.... |
|00001610| 20 20 20 65 70 6f 63 68 | 20 2a 2f 0a 58 23 64 65 | epoch| */.X#de|
|00001620| 66 69 6e 65 20 6d 6c 6e | 6f 64 65 09 20 20 20 20 |fine mln|ode. |
|00001630| 31 35 31 2e 39 35 30 34 | 32 39 09 20 20 20 2f 2a |151.9504|29. /*|
|00001640| 20 4d 65 61 6e 20 6c 6f | 6e 67 69 74 75 64 65 20 | Mean lo|ngitude |
|00001650| 6f 66 20 74 68 65 20 6e | 6f 64 65 20 61 74 20 74 |of the n|ode at t|
|00001660| 68 65 0a 58 09 09 09 09 | 20 20 20 20 20 20 65 70 |he.X....| ep|
|00001670| 6f 63 68 20 2a 2f 0a 58 | 23 64 65 66 69 6e 65 20 |och */.X|#define |
|00001680| 6d 69 6e 63 20 20 20 20 | 20 20 20 20 35 2e 31 34 |minc | 5.14|
|00001690| 35 33 39 36 20 20 20 20 | 20 20 20 2f 2a 20 49 6e |5396 | /* In|
|000016a0| 63 6c 69 6e 61 74 69 6f | 6e 20 6f 66 20 74 68 65 |clinatio|n of the|
|000016b0| 20 4d 6f 6f 6e 27 73 20 | 6f 72 62 69 74 20 2a 2f | Moon's |orbit */|
|000016c0| 0a 58 23 64 65 66 69 6e | 65 20 6d 65 63 63 20 20 |.X#defin|e mecc |
|000016d0| 20 20 20 20 20 20 30 2e | 30 35 34 39 30 30 20 20 | 0.|054900 |
|000016e0| 20 20 20 20 20 2f 2a 20 | 45 63 63 65 6e 74 72 69 | /* |Eccentri|
|000016f0| 63 69 74 79 20 6f 66 20 | 74 68 65 20 4d 6f 6f 6e |city of |the Moon|
|00001700| 27 73 20 6f 72 62 69 74 | 20 2a 2f 0a 58 23 64 65 |'s orbit| */.X#de|
|00001710| 66 69 6e 65 20 6d 61 6e | 67 73 69 7a 20 20 20 20 |fine man|gsiz |
|00001720| 20 30 2e 35 31 38 31 20 | 20 20 20 20 20 20 20 20 | 0.5181 | |
|00001730| 2f 2a 20 4d 6f 6f 6e 27 | 73 20 61 6e 67 75 6c 61 |/* Moon'|s angula|
|00001740| 72 20 73 69 7a 65 20 61 | 74 20 64 69 73 74 61 6e |r size a|t distan|
|00001750| 63 65 20 61 0a 58 09 09 | 09 09 20 20 20 20 20 20 |ce a.X..|.. |
|00001760| 66 72 6f 6d 20 45 61 72 | 74 68 20 2a 2f 0a 58 23 |from Ear|th */.X#|
|00001770| 64 65 66 69 6e 65 20 6d | 73 6d 61 78 20 20 20 20 |define m|smax |
|00001780| 20 20 20 33 38 34 34 30 | 31 2e 30 20 20 20 20 20 | 38440|1.0 |
|00001790| 20 20 2f 2a 20 53 65 6d | 69 2d 6d 61 6a 6f 72 20 | /* Sem|i-major |
|000017a0| 61 78 69 73 20 6f 66 20 | 4d 6f 6f 6e 27 73 20 6f |axis of |Moon's o|
|000017b0| 72 62 69 74 20 69 6e 20 | 6b 6d 20 2a 2f 0a 58 23 |rbit in |km */.X#|
|000017c0| 64 65 66 69 6e 65 20 6d | 70 61 72 61 6c 6c 61 78 |define m|parallax|
|000017d0| 20 20 20 30 2e 39 35 30 | 37 09 20 20 20 2f 2a 20 | 0.950|7. /* |
|000017e0| 50 61 72 61 6c 6c 61 78 | 20 61 74 20 64 69 73 74 |Parallax| at dist|
|000017f0| 61 6e 63 65 20 61 20 66 | 72 6f 6d 20 45 61 72 74 |ance a f|rom Eart|
|00001800| 68 20 2a 2f 0a 58 23 64 | 65 66 69 6e 65 20 73 79 |h */.X#d|efine sy|
|00001810| 6e 6d 6f 6e 74 68 20 20 | 20 20 32 39 2e 35 33 30 |nmonth | 29.530|
|00001820| 35 38 38 36 38 20 20 20 | 20 2f 2a 20 53 79 6e 6f |58868 | /* Syno|
|00001830| 64 69 63 20 6d 6f 6e 74 | 68 20 28 6e 65 77 20 4d |dic mont|h (new M|
|00001840| 6f 6f 6e 20 74 6f 20 6e | 65 77 20 4d 6f 6f 6e 29 |oon to n|ew Moon)|
|00001850| 20 2a 2f 0a 58 23 64 65 | 66 69 6e 65 20 6c 75 6e | */.X#de|fine lun|
|00001860| 61 74 62 61 73 65 20 20 | 20 32 34 32 33 34 33 36 |atbase | 2423436|
|00001870| 2e 30 20 20 20 20 20 20 | 2f 2a 20 42 61 73 65 20 |.0 |/* Base |
|00001880| 64 61 74 65 20 66 6f 72 | 20 45 2e 20 57 2e 20 42 |date for| E. W. B|
|00001890| 72 6f 77 6e 27 73 20 6e | 75 6d 62 65 72 65 64 0a |rown's n|umbered.|
|000018a0| 58 09 09 09 09 20 20 20 | 20 20 20 73 65 72 69 65 |X.... | serie|
|000018b0| 73 20 6f 66 20 6c 75 6e | 61 74 69 6f 6e 73 20 28 |s of lun|ations (|
|000018c0| 31 39 32 33 20 4a 61 6e | 75 61 72 79 20 31 36 29 |1923 Jan|uary 16)|
|000018d0| 20 2a 2f 0a 58 0a 58 2f | 2a 20 20 50 72 6f 70 65 | */.X.X/|* Prope|
|000018e0| 72 74 69 65 73 20 6f 66 | 20 74 68 65 20 45 61 72 |rties of| the Ear|
|000018f0| 74 68 20 20 2a 2f 0a 58 | 0a 58 23 64 65 66 69 6e |th */.X|.X#defin|
|00001900| 65 20 65 61 72 74 68 72 | 61 64 20 20 20 20 36 33 |e earthr|ad 63|
|00001910| 37 38 2e 31 36 09 20 20 | 20 2f 2a 20 52 61 64 69 |78.16. | /* Radi|
|00001920| 75 73 20 6f 66 20 45 61 | 72 74 68 20 69 6e 20 6b |us of Ea|rth in k|
|00001930| 69 6c 6f 6d 65 74 72 65 | 73 20 2a 2f 0a 58 0a 58 |ilometre|s */.X.X|
|00001940| 23 64 65 66 69 6e 65 20 | 50 49 20 33 2e 31 34 31 |#define |PI 3.141|
|00001950| 35 39 32 36 35 33 35 38 | 39 37 39 33 32 33 38 34 |59265358|97932384|
|00001960| 36 0a 58 0a 58 2f 2a 20 | 20 48 61 6e 64 79 20 6d |6.X.X/* | Handy m|
|00001970| 61 74 68 65 6d 61 74 69 | 63 61 6c 20 66 75 6e 63 |athemati|cal func|
|00001980| 74 69 6f 6e 73 20 20 2a | 2f 0a 58 0a 58 23 69 66 |tions *|/.X.X#if|
|00001990| 64 65 66 20 73 67 6e 0a | 58 23 75 6e 64 65 66 20 |def sgn.|X#undef |
|000019a0| 73 67 6e 0a 58 23 65 6e | 64 69 66 0a 58 23 64 65 |sgn.X#en|dif.X#de|
|000019b0| 66 69 6e 65 20 73 67 6e | 28 78 29 20 28 28 28 78 |fine sgn|(x) (((x|
|000019c0| 29 20 3c 20 30 29 20 3f | 20 2d 31 20 3a 20 28 28 |) < 0) ?| -1 : ((|
|000019d0| 78 29 20 3e 20 30 20 3f | 20 31 20 3a 20 30 29 29 |x) > 0 ?| 1 : 0))|
|000019e0| 09 20 20 2f 2a 20 45 78 | 74 72 61 63 74 20 73 69 |. /* Ex|tract si|
|000019f0| 67 6e 20 2a 2f 0a 58 0a | 58 23 69 66 64 65 66 20 |gn */.X.|X#ifdef |
|00001a00| 61 62 73 0a 58 23 75 6e | 64 65 66 20 61 62 73 0a |abs.X#un|def abs.|
|00001a10| 58 23 65 6e 64 69 66 0a | 58 23 64 65 66 69 6e 65 |X#endif.|X#define|
|00001a20| 20 61 62 73 28 78 29 20 | 28 28 78 29 20 3c 20 30 | abs(x) |((x) < 0|
|00001a30| 20 3f 20 28 2d 28 78 29 | 29 20 3a 20 28 78 29 29 | ? (-(x)|) : (x))|
|00001a40| 20 09 09 20 20 2f 2a 20 | 41 62 73 6f 6c 75 74 65 | .. /* |Absolute|
|00001a50| 20 76 61 6c 20 2a 2f 0a | 58 0a 58 23 64 65 66 69 | val */.|X.X#defi|
|00001a60| 6e 65 20 66 69 78 61 6e | 67 6c 65 28 61 29 20 28 |ne fixan|gle(a) (|
|00001a70| 28 61 29 20 2d 20 33 36 | 30 2e 30 20 2a 20 28 66 |(a) - 36|0.0 * (f|
|00001a80| 6c 6f 6f 72 28 28 61 29 | 20 2f 20 33 36 30 2e 30 |loor((a)| / 360.0|
|00001a90| 29 29 29 20 20 2f 2a 20 | 46 69 78 20 61 6e 67 6c |))) /* |Fix angl|
|00001aa0| 65 09 20 20 2a 2f 0a 58 | 23 64 65 66 69 6e 65 20 |e. */.X|#define |
|00001ab0| 74 6f 72 61 64 28 64 29 | 20 28 28 64 29 20 2a 20 |torad(d)| ((d) * |
|00001ac0| 28 50 49 20 2f 20 31 38 | 30 2e 30 29 29 09 09 09 |(PI / 18|0.0))...|
|00001ad0| 20 20 2f 2a 20 44 65 67 | 2d 3e 52 61 64 09 20 20 | /* Deg|->Rad. |
|00001ae0| 2a 2f 0a 58 23 64 65 66 | 69 6e 65 20 74 6f 64 65 |*/.X#def|ine tode|
|00001af0| 67 28 64 29 20 28 28 64 | 29 20 2a 20 28 31 38 30 |g(d) ((d|) * (180|
|00001b00| 2e 30 20 2f 20 50 49 29 | 29 09 09 09 20 20 2f 2a |.0 / PI)|)... /*|
|00001b10| 20 52 61 64 2d 3e 44 65 | 67 09 20 20 2a 2f 0a 58 | Rad->De|g. */.X|
|00001b20| 23 64 65 66 69 6e 65 20 | 64 73 69 6e 28 78 29 20 |#define |dsin(x) |
|00001b30| 28 73 69 6e 28 74 6f 72 | 61 64 28 28 78 29 29 29 |(sin(tor|ad((x)))|
|00001b40| 29 09 09 09 20 20 2f 2a | 20 53 69 6e 20 66 72 6f |)... /*| Sin fro|
|00001b50| 6d 20 64 65 67 20 2a 2f | 0a 58 23 64 65 66 69 6e |m deg */|.X#defin|
|00001b60| 65 20 64 63 6f 73 28 78 | 29 20 28 63 6f 73 28 74 |e dcos(x|) (cos(t|
|00001b70| 6f 72 61 64 28 28 78 29 | 29 29 29 09 09 09 20 20 |orad((x)|)))... |
|00001b80| 2f 2a 20 43 6f 73 20 66 | 72 6f 6d 20 64 65 67 20 |/* Cos f|rom deg |
|00001b90| 2a 2f 0a 58 0a 58 2f 2a | 2a 2a 2a 2a 2a 2a 2a 2a |*/.X.X/*|********|
|00001ba0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001bb0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001bc0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001bd0| 2a 2a 2a 2a 2a 2a 2f 0a | 58 2f 2a 20 20 20 20 20 |******/.|X/* |
|00001be0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001bf0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001c00| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001c10| 20 20 20 20 20 20 20 20 | 2a 2f 0a 58 2f 2a 20 20 | |*/.X/* |
|00001c20| 6a 64 61 74 65 20 20 20 | 20 20 20 20 20 20 20 20 |jdate | |
|00001c30| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001c40| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001c50| 20 20 20 20 20 20 20 20 | 20 20 20 2a 2f 0a 58 2f | | */.X/|
|00001c60| 2a 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |* | |
|00001c70| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001c80| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001c90| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 2a 2f | | */|
|00001ca0| 0a 58 2f 2a 20 20 43 6f | 6e 76 65 72 74 20 61 20 |.X/* Co|nvert a |
|00001cb0| 64 61 74 65 20 61 6e 64 | 20 74 69 6d 65 20 74 6f |date and| time to|
|00001cc0| 20 4a 75 6c 69 61 6e 20 | 64 61 79 20 61 6e 64 20 | Julian |day and |
|00001cd0| 66 72 61 63 74 69 6f 6e | 2e 20 20 20 20 20 20 20 |fraction|. |
|00001ce0| 20 2a 2f 0a 58 2f 2a 20 | 20 20 20 20 20 20 20 20 | */.X/* | |
|00001cf0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001d00| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001d10| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001d20| 20 20 20 20 2a 2f 0a 58 | 2f 2a 2a 2a 2a 2a 2a 2a | */.X|/*******|
|00001d30| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001d40| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001d50| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001d60| 2a 2a 2a 2a 2a 2a 2a 2a | 2f 0a 58 23 69 66 64 65 |********|/.X#ifde|
|00001d70| 66 20 48 41 56 45 5f 50 | 52 4f 54 4f 53 0a 58 50 |f HAVE_P|ROTOS.XP|
|00001d80| 52 49 56 41 54 45 20 6c | 6f 6e 67 20 6a 64 61 74 |RIVATE l|ong jdat|
|00001d90| 65 28 69 6e 74 20 79 2c | 20 69 6e 74 20 6d 6f 6e |e(int y,| int mon|
|00001da0| 2c 20 69 6e 74 20 64 61 | 79 29 0a 58 23 65 6c 73 |, int da|y).X#els|
|00001db0| 65 0a 58 73 74 61 74 69 | 63 20 6c 6f 6e 67 20 6a |e.Xstati|c long j|
|00001dc0| 64 61 74 65 28 79 2c 20 | 6d 6f 6e 2c 20 64 61 79 |date(y, |mon, day|
|00001dd0| 29 0a 58 69 6e 74 20 79 | 2c 20 6d 6f 6e 2c 20 64 |).Xint y|, mon, d|
|00001de0| 61 79 3b 0a 58 23 65 6e | 64 69 66 0a 58 7b 0a 58 |ay;.X#en|dif.X{.X|
|00001df0| 20 20 20 6c 6f 6e 67 20 | 63 2c 20 6d 3b 0a 58 0a | long |c, m;.X.|
|00001e00| 58 20 20 20 6d 20 3d 20 | 6d 6f 6e 2b 31 3b 0a 58 |X m = |mon+1;.X|
|00001e10| 20 20 20 69 66 20 28 6d | 3e 32 29 20 7b 0a 58 20 | if (m|>2) {.X |
|00001e20| 20 20 20 20 20 6d 20 2d | 3d 20 33 3b 0a 58 20 20 | m -|= 3;.X |
|00001e30| 20 7d 20 65 6c 73 65 20 | 7b 0a 58 20 20 20 20 20 | } else |{.X |
|00001e40| 20 6d 20 2b 3d 20 39 3b | 0a 58 20 20 20 20 20 20 | m += 9;|.X |
|00001e50| 79 2d 2d 3b 0a 58 20 20 | 20 7d 0a 58 20 20 20 63 |y--;.X | }.X c|
|00001e60| 20 3d 20 79 2f 31 30 30 | 4c 3b 20 20 20 2f 2a 20 | = y/100|L; /* |
|00001e70| 43 65 6e 74 75 72 79 20 | 2a 2f 0a 58 20 20 20 79 |Century |*/.X y|
|00001e80| 20 2d 3d 20 31 30 30 4c | 20 2a 20 63 3b 0a 58 20 | -= 100L| * c;.X |
|00001e90| 20 20 72 65 74 75 72 6e | 20 64 61 79 20 2b 20 28 | return| day + (|
|00001ea0| 63 2a 31 34 36 30 39 37 | 4c 29 2f 34 20 2b 20 28 |c*146097|L)/4 + (|
|00001eb0| 79 2a 31 34 36 31 4c 29 | 2f 34 20 2b 20 28 6d 2a |y*1461L)|/4 + (m*|
|00001ec0| 31 35 33 4c 2b 32 29 2f | 35 20 2b 20 31 37 32 31 |153L+2)/|5 + 1721|
|00001ed0| 31 31 39 4c 3b 0a 58 7d | 0a 58 0a 58 2f 2a 2a 2a |119L;.X}|.X.X/***|
|00001ee0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001ef0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001f00| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00001f10| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2f 0a 58 2f |********|****/.X/|
|00001f20| 2a 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |* | |
|00001f30| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001f40| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001f50| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 2a 2f | | */|
|00001f60| 0a 58 2f 2a 20 20 6a 74 | 69 6d 65 20 20 20 20 20 |.X/* jt|ime |
|00001f70| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001f80| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001f90| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001fa0| 20 2a 2f 0a 58 2f 2a 20 | 20 20 20 20 20 20 20 20 | */.X/* | |
|00001fb0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001fc0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001fd0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001fe0| 20 20 20 20 2a 2f 0a 58 | 2f 2a 20 20 43 6f 6e 76 | */.X|/* Conv|
|00001ff0| 65 72 74 20 61 20 47 4d | 54 20 64 61 74 65 20 61 |ert a GM|T date a|
|00002000| 6e 64 20 74 69 6d 65 20 | 74 6f 20 61 73 74 72 6f |nd time |to astro|
|00002010| 6e 6f 6d 69 63 61 6c 20 | 4a 75 6c 69 61 6e 20 74 |nomical |Julian t|
|00002020| 69 6d 65 2c 20 20 20 2a | 2f 0a 58 2f 2a 20 20 69 |ime, *|/.X/* i|
|00002030| 2e 65 2e 20 4a 75 6c 69 | 61 6e 20 64 61 74 65 20 |.e. Juli|an date |
|00002040| 70 6c 75 73 20 64 61 79 | 20 66 72 61 63 74 69 6f |plus day| fractio|
|00002050| 6e 2c 20 65 78 70 72 65 | 73 73 65 64 20 61 73 20 |n, expre|ssed as |
|00002060| 61 20 64 6f 75 62 6c 65 | 20 20 2a 2f 0a 58 2f 2a |a double| */.X/*|
|00002070| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002080| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002090| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000020a0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 2a 2f 0a | | */.|
|000020b0| 58 2f 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |X/******|********|
|000020c0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000020d0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000020e0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000020f0| 2a 2f 0a 58 23 69 66 64 | 65 66 20 48 41 56 45 5f |*/.X#ifd|ef HAVE_|
|00002100| 50 52 4f 54 4f 53 0a 58 | 50 52 49 56 41 54 45 20 |PROTOS.X|PRIVATE |
|00002110| 64 6f 75 62 6c 65 20 6a | 74 69 6d 65 28 69 6e 74 |double j|time(int|
|00002120| 20 79 2c 20 69 6e 74 20 | 6d 6f 6e 2c 20 69 6e 74 | y, int |mon, int|
|00002130| 20 64 61 79 2c 20 69 6e | 74 20 68 6f 75 72 2c 20 | day, in|t hour, |
|00002140| 69 6e 74 20 6d 69 6e 2c | 20 69 6e 74 20 73 65 63 |int min,| int sec|
|00002150| 29 0a 58 23 65 6c 73 65 | 0a 58 73 74 61 74 69 63 |).X#else|.Xstatic|
|00002160| 20 64 6f 75 62 6c 65 20 | 6a 74 69 6d 65 28 79 2c | double |jtime(y,|
|00002170| 20 6d 6f 6e 2c 20 64 61 | 79 2c 20 68 6f 75 72 2c | mon, da|y, hour,|
|00002180| 20 6d 69 6e 2c 20 73 65 | 63 29 0a 58 69 6e 74 20 | min, se|c).Xint |
|00002190| 79 2c 20 6d 6f 6e 2c 20 | 64 61 79 2c 20 68 6f 75 |y, mon, |day, hou|
|000021a0| 72 2c 20 6d 69 6e 2c 20 | 73 65 63 3b 0a 58 23 65 |r, min, |sec;.X#e|
|000021b0| 6e 64 69 66 0a 58 7b 0a | 58 20 20 20 72 65 74 75 |ndif.X{.|X retu|
|000021c0| 72 6e 20 28 6a 64 61 74 | 65 28 79 2c 20 6d 6f 6e |rn (jdat|e(y, mon|
|000021d0| 2c 20 64 61 79 29 2d 30 | 2e 35 29 20 2b 0a 58 20 |, day)-0|.5) +.X |
|000021e0| 20 20 20 20 20 20 20 20 | 20 20 20 20 28 73 65 63 | | (sec|
|000021f0| 20 2b 20 36 30 4c 20 2a | 20 28 6c 6f 6e 67 29 20 | + 60L *| (long) |
|00002200| 6d 69 6e 20 2b 20 33 36 | 30 30 4c 20 2a 20 28 6c |min + 36|00L * (l|
|00002210| 6f 6e 67 29 20 68 6f 75 | 72 29 20 2f 20 38 36 34 |ong) hou|r) / 864|
|00002220| 30 30 2e 30 3b 0a 58 7d | 0a 58 0a 58 2f 2a 2a 2a |00.0;.X}|.X.X/***|
|00002230| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002240| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002250| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002260| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2f 0a 58 2f |********|****/.X/|
|00002270| 2a 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |* | |
|00002280| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002290| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000022a0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 2a 2f | | */|
|000022b0| 0a 58 2f 2a 20 20 6a 79 | 65 61 72 20 20 20 20 20 |.X/* jy|ear |
|000022c0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000022d0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000022e0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000022f0| 20 2a 2f 0a 58 2f 2a 20 | 20 20 20 20 20 20 20 20 | */.X/* | |
|00002300| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002310| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002320| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002330| 20 20 20 20 2a 2f 0a 58 | 2f 2a 20 20 43 6f 6e 76 | */.X|/* Conv|
|00002340| 65 72 74 20 61 20 4a 75 | 6c 69 61 6e 20 64 61 74 |ert a Ju|lian dat|
|00002350| 65 20 74 6f 20 79 65 61 | 72 2c 20 6d 6f 6e 74 68 |e to yea|r, month|
|00002360| 2c 20 64 61 79 2e 20 20 | 20 20 20 20 20 20 20 20 |, day. | |
|00002370| 20 20 20 20 20 20 20 2a | 2f 0a 58 2f 2a 20 20 20 | *|/.X/* |
|00002380| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002390| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000023a0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000023b0| 20 20 20 20 20 20 20 20 | 20 20 2a 2f 0a 58 2f 2a | | */.X/*|
|000023c0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000023d0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000023e0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000023f0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2f 0a |********|******/.|
|00002400| 58 23 69 66 64 65 66 20 | 48 41 56 45 5f 50 52 4f |X#ifdef |HAVE_PRO|
|00002410| 54 4f 53 0a 58 50 52 49 | 56 41 54 45 20 76 6f 69 |TOS.XPRI|VATE voi|
|00002420| 64 20 6a 79 65 61 72 28 | 64 6f 75 62 6c 65 20 74 |d jyear(|double t|
|00002430| 64 2c 20 69 6e 74 20 2a | 79 79 2c 20 69 6e 74 20 |d, int *|yy, int |
|00002440| 2a 6d 6d 2c 20 69 6e 74 | 20 2a 64 64 29 0a 58 23 |*mm, int| *dd).X#|
|00002450| 65 6c 73 65 0a 58 73 74 | 61 74 69 63 20 76 6f 69 |else.Xst|atic voi|
|00002460| 64 20 6a 79 65 61 72 28 | 74 64 2c 20 79 79 2c 20 |d jyear(|td, yy, |
|00002470| 6d 6d 2c 20 64 64 29 0a | 58 64 6f 75 62 6c 65 20 |mm, dd).|Xdouble |
|00002480| 74 64 3b 0a 58 69 6e 74 | 20 2a 79 79 2c 20 2a 6d |td;.Xint| *yy, *m|
|00002490| 6d 2c 20 2a 64 64 3b 0a | 58 23 65 6e 64 69 66 0a |m, *dd;.|X#endif.|
|000024a0| 58 7b 0a 58 20 20 20 64 | 6f 75 62 6c 65 20 6a 2c |X{.X d|ouble j,|
|000024b0| 20 64 2c 20 79 2c 20 6d | 3b 0a 58 0a 58 20 20 20 | d, y, m|;.X.X |
|000024c0| 74 64 20 2b 3d 20 30 2e | 35 3b 20 20 20 20 20 20 |td += 0.|5; |
|000024d0| 20 20 20 2f 2a 20 41 73 | 74 72 6f 6e 6f 6d 69 63 | /* As|tronomic|
|000024e0| 61 6c 20 74 6f 20 63 69 | 76 69 6c 20 2a 2f 0a 58 |al to ci|vil */.X|
|000024f0| 20 20 20 6a 20 3d 20 66 | 6c 6f 6f 72 28 74 64 29 | j = f|loor(td)|
|00002500| 3b 0a 58 20 20 20 6a 20 | 3d 20 6a 20 2d 20 31 37 |;.X j |= j - 17|
|00002510| 32 31 31 31 39 2e 30 3b | 0a 58 20 20 20 79 20 3d |21119.0;|.X y =|
|00002520| 20 66 6c 6f 6f 72 28 28 | 28 34 20 2a 20 6a 29 20 | floor((|(4 * j) |
|00002530| 2d 20 31 29 20 2f 20 31 | 34 36 30 39 37 2e 30 29 |- 1) / 1|46097.0)|
|00002540| 3b 0a 58 20 20 20 6a 20 | 3d 20 28 6a 20 2a 20 34 |;.X j |= (j * 4|
|00002550| 2e 30 29 20 2d 20 28 31 | 2e 30 20 2b 20 28 31 34 |.0) - (1|.0 + (14|
|00002560| 36 30 39 37 2e 30 20 2a | 20 79 29 29 3b 0a 58 20 |6097.0 *| y));.X |
|00002570| 20 20 64 20 3d 20 66 6c | 6f 6f 72 28 6a 20 2f 20 | d = fl|oor(j / |
|00002580| 34 2e 30 29 3b 0a 58 20 | 20 20 6a 20 3d 20 66 6c |4.0);.X | j = fl|
|00002590| 6f 6f 72 28 28 28 34 2e | 30 20 2a 20 64 29 20 2b |oor(((4.|0 * d) +|
|000025a0| 20 33 2e 30 29 20 2f 20 | 31 34 36 31 2e 30 29 3b | 3.0) / |1461.0);|
|000025b0| 0a 58 20 20 20 64 20 3d | 20 28 28 34 2e 30 20 2a |.X d =| ((4.0 *|
|000025c0| 20 64 29 20 2b 20 33 2e | 30 29 20 2d 20 28 31 34 | d) + 3.|0) - (14|
|000025d0| 36 31 2e 30 20 2a 20 6a | 29 3b 0a 58 20 20 20 64 |61.0 * j|);.X d|
|000025e0| 20 3d 20 66 6c 6f 6f 72 | 28 28 64 20 2b 20 34 2e | = floor|((d + 4.|
|000025f0| 30 29 20 2f 20 34 2e 30 | 29 3b 0a 58 20 20 20 6d |0) / 4.0|);.X m|
|00002600| 20 3d 20 66 6c 6f 6f 72 | 28 28 28 35 2e 30 20 2a | = floor|(((5.0 *|
|00002610| 20 64 29 20 2d 20 33 29 | 20 2f 20 31 35 33 2e 30 | d) - 3)| / 153.0|
|00002620| 29 3b 0a 58 20 20 20 64 | 20 3d 20 28 35 2e 30 20 |);.X d| = (5.0 |
|00002630| 2a 20 64 29 20 2d 20 28 | 33 2e 30 20 2b 20 28 31 |* d) - (|3.0 + (1|
|00002640| 35 33 2e 30 20 2a 20 6d | 29 29 3b 0a 58 20 20 20 |53.0 * m|));.X |
|00002650| 64 20 3d 20 66 6c 6f 6f | 72 28 28 64 20 2b 20 35 |d = floo|r((d + 5|
|00002660| 2e 30 29 20 2f 20 35 2e | 30 29 3b 0a 58 20 20 20 |.0) / 5.|0);.X |
|00002670| 79 20 3d 20 28 31 30 30 | 2e 30 20 2a 20 79 29 20 |y = (100|.0 * y) |
|00002680| 2b 20 6a 3b 0a 58 20 20 | 20 69 66 20 28 6d 20 3c |+ j;.X | if (m <|
|00002690| 20 31 30 2e 30 29 0a 58 | 20 20 20 20 20 20 6d 20 | 10.0).X| m |
|000026a0| 3d 20 6d 20 2b 20 32 3b | 0a 58 20 20 20 65 6c 73 |= m + 2;|.X els|
|000026b0| 65 20 7b 0a 58 20 20 20 | 20 20 20 6d 20 3d 20 6d |e {.X | m = m|
|000026c0| 20 2d 20 31 30 3b 0a 58 | 20 20 20 20 20 20 79 20 | - 10;.X| y |
|000026d0| 3d 20 79 20 2b 20 31 3b | 0a 58 20 20 20 7d 0a 58 |= y + 1;|.X }.X|
|000026e0| 20 20 20 2a 79 79 20 3d | 20 79 3b 0a 58 20 20 20 | *yy =| y;.X |
|000026f0| 2a 6d 6d 20 3d 20 6d 3b | 0a 58 20 20 20 2a 64 64 |*mm = m;|.X *dd|
|00002700| 20 3d 20 64 3b 0a 58 7d | 0a 58 0a 58 2f 2a 2a 2a | = d;.X}|.X.X/***|
|00002710| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002720| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002730| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002740| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2f 0a 58 2f |********|****/.X/|
|00002750| 2a 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |* | |
|00002760| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002770| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002780| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 2a 2f | | */|
|00002790| 0a 58 2f 2a 20 20 6a 68 | 6d 73 20 20 20 20 20 20 |.X/* jh|ms |
|000027a0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000027b0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000027c0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000027d0| 20 2a 2f 0a 58 2f 2a 20 | 20 20 20 20 20 20 20 20 | */.X/* | |
|000027e0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000027f0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002800| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002810| 20 20 20 20 2a 2f 0a 58 | 2f 2a 20 20 43 6f 6e 76 | */.X|/* Conv|
|00002820| 65 72 74 20 61 20 4a 75 | 6c 69 61 6e 20 74 69 6d |ert a Ju|lian tim|
|00002830| 65 20 74 6f 20 68 6f 75 | 72 2c 20 6d 69 6e 75 74 |e to hou|r, minut|
|00002840| 65 73 20 61 6e 64 20 73 | 65 63 6f 6e 64 73 2e 20 |es and s|econds. |
|00002850| 20 20 20 20 20 20 20 2a | 2f 0a 58 2f 2a 20 20 20 | *|/.X/* |
|00002860| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002870| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002880| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002890| 20 20 20 20 20 20 20 20 | 20 20 2a 2f 0a 58 2f 2a | | */.X/*|
|000028a0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000028b0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000028c0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000028d0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2f 0a |********|******/.|
|000028e0| 58 23 69 66 64 65 66 20 | 48 41 56 45 5f 50 52 4f |X#ifdef |HAVE_PRO|
|000028f0| 54 4f 53 0a 58 50 52 49 | 56 41 54 45 20 76 6f 69 |TOS.XPRI|VATE voi|
|00002900| 64 20 6a 68 6d 73 28 64 | 6f 75 62 6c 65 20 6a 2c |d jhms(d|ouble j,|
|00002910| 20 69 6e 74 20 2a 68 2c | 20 69 6e 74 20 2a 6d 2c | int *h,| int *m,|
|00002920| 20 69 6e 74 20 2a 73 29 | 0a 58 23 65 6c 73 65 0a | int *s)|.X#else.|
|00002930| 58 73 74 61 74 69 63 20 | 76 6f 69 64 20 6a 68 6d |Xstatic |void jhm|
|00002940| 73 28 6a 2c 20 68 2c 20 | 6d 2c 20 73 29 0a 58 64 |s(j, h, |m, s).Xd|
|00002950| 6f 75 62 6c 65 20 6a 3b | 0a 58 69 6e 74 20 2a 68 |ouble j;|.Xint *h|
|00002960| 2c 20 2a 6d 2c 20 2a 73 | 3b 0a 58 23 65 6e 64 69 |, *m, *s|;.X#endi|
|00002970| 66 0a 58 7b 0a 58 20 20 | 20 6c 6f 6e 67 20 69 6a |f.X{.X | long ij|
|00002980| 3b 0a 58 0a 58 20 20 20 | 6a 20 2b 3d 20 30 2e 35 |;.X.X |j += 0.5|
|00002990| 3b 20 20 20 20 20 20 20 | 20 20 2f 2a 20 41 73 74 |; | /* Ast|
|000029a0| 72 6f 6e 6f 6d 69 63 61 | 6c 20 74 6f 20 63 69 76 |ronomica|l to civ|
|000029b0| 69 6c 20 2a 2f 0a 58 20 | 20 20 69 6a 20 3d 20 28 |il */.X | ij = (|
|000029c0| 6a 20 2d 20 66 6c 6f 6f | 72 28 6a 29 29 20 2a 20 |j - floo|r(j)) * |
|000029d0| 38 36 34 30 30 2e 30 3b | 0a 58 20 20 20 2a 68 20 |86400.0;|.X *h |
|000029e0| 3d 20 69 6a 20 2f 20 33 | 36 30 30 4c 3b 0a 58 20 |= ij / 3|600L;.X |
|000029f0| 20 20 2a 6d 20 3d 20 28 | 69 6a 20 2f 20 36 30 4c | *m = (|ij / 60L|
|00002a00| 29 20 25 20 36 30 4c 3b | 0a 58 20 20 20 2a 73 20 |) % 60L;|.X *s |
|00002a10| 3d 20 69 6a 20 25 20 36 | 30 4c 3b 0a 58 7d 0a 58 |= ij % 6|0L;.X}.X|
|00002a20| 0a 58 2f 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |.X/*****|********|
|00002a30| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002a40| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002a50| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002a60| 2a 2a 2f 0a 58 2f 2a 20 | 20 20 20 20 20 20 20 20 |**/.X/* | |
|00002a70| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002a80| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002a90| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002aa0| 20 20 20 20 2a 2f 0a 58 | 2f 2a 20 20 6d 65 61 6e | */.X|/* mean|
|00002ab0| 70 68 61 73 65 20 20 20 | 20 20 20 20 20 20 20 20 |phase | |
|00002ac0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002ad0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002ae0| 20 20 20 20 20 20 20 2a | 2f 0a 58 2f 2a 20 20 20 | *|/.X/* |
|00002af0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002b00| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002b10| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002b20| 20 20 20 20 20 20 20 20 | 20 20 2a 2f 0a 58 2f 2a | | */.X/*|
|00002b30| 20 20 43 61 6c 63 75 6c | 61 74 65 73 20 6d 65 61 | Calcul|ates mea|
|00002b40| 6e 20 70 68 61 73 65 20 | 6f 66 20 74 68 65 20 4d |n phase |of the M|
|00002b50| 6f 6f 6e 20 66 6f 72 20 | 61 20 20 20 20 20 20 20 |oon for |a |
|00002b60| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 2a 2f 0a | | */.|
|00002b70| 58 2f 2a 20 20 67 69 76 | 65 6e 20 62 61 73 65 20 |X/* giv|en base |
|00002b80| 64 61 74 65 20 61 6e 64 | 20 64 65 73 69 72 65 64 |date and| desired|
|00002b90| 20 70 68 61 73 65 3a 09 | 09 09 20 20 20 20 20 20 | phase:.|.. |
|00002ba0| 20 2a 2f 0a 58 2f 2a 20 | 20 20 20 20 30 2e 30 20 | */.X/* | 0.0 |
|00002bb0| 20 20 4e 65 77 20 4d 6f | 6f 6e 09 09 09 09 09 20 | New Mo|on..... |
|00002bc0| 20 20 20 20 20 20 2a 2f | 0a 58 2f 2a 20 20 20 20 | */|.X/* |
|00002bd0| 20 30 2e 32 35 20 20 46 | 69 72 73 74 20 71 75 61 | 0.25 F|irst qua|
|00002be0| 72 74 65 72 09 09 09 09 | 20 20 20 20 20 20 20 2a |rter....| *|
|00002bf0| 2f 0a 58 2f 2a 20 20 20 | 20 20 30 2e 35 20 20 20 |/.X/* | 0.5 |
|00002c00| 46 75 6c 6c 20 6d 6f 6f | 6e 09 09 09 09 09 20 20 |Full moo|n..... |
|00002c10| 20 20 20 20 20 2a 2f 0a | 58 2f 2a 20 20 20 20 20 | */.|X/* |
|00002c20| 30 2e 37 35 20 20 4c 61 | 73 74 20 71 75 61 72 74 |0.75 La|st quart|
|00002c30| 65 72 09 09 09 09 20 20 | 20 20 20 20 20 2a 2f 0a |er.... | */.|
|00002c40| 58 2f 2a 20 20 42 65 77 | 61 72 65 21 21 21 20 20 |X/* Bew|are!!! |
|00002c50| 54 68 69 73 20 72 6f 75 | 74 69 6e 65 20 72 65 74 |This rou|tine ret|
|00002c60| 75 72 6e 73 20 6d 65 61 | 6e 69 6e 67 6c 65 73 73 |urns mea|ningless|
|00002c70| 09 09 20 20 20 20 20 20 | 20 2a 2f 0a 58 2f 2a 20 |.. | */.X/* |
|00002c80| 20 72 65 73 75 6c 74 73 | 20 66 6f 72 20 61 6e 79 | results| for any|
|00002c90| 20 6f 74 68 65 72 20 70 | 68 61 73 65 20 61 72 67 | other p|hase arg|
|00002ca0| 75 6d 65 6e 74 73 2e 20 | 20 44 6f 6e 27 74 09 20 |uments. | Don't. |
|00002cb0| 20 20 20 20 20 20 2a 2f | 0a 58 2f 2a 20 20 61 74 | */|.X/* at|
|00002cc0| 74 65 6d 70 74 20 74 6f | 20 67 65 6e 65 72 61 6c |tempt to| general|
|00002cd0| 69 73 65 20 69 74 20 77 | 69 74 68 6f 75 74 20 75 |ise it w|ithout u|
|00002ce0| 6e 64 65 72 73 74 61 6e | 64 69 6e 67 09 20 20 20 |nderstan|ding. |
|00002cf0| 20 20 20 20 2a 2f 0a 58 | 2f 2a 20 20 74 68 61 74 | */.X|/* that|
|00002d00| 20 74 68 65 20 6d 6f 74 | 69 6f 6e 20 6f 66 20 74 | the mot|ion of t|
|00002d10| 68 65 20 6d 6f 6f 6e 20 | 69 73 20 66 61 72 20 6d |he moon |is far m|
|00002d20| 6f 72 65 20 63 6f 6d 70 | 6c 69 63 61 74 65 64 09 |ore comp|licated.|
|00002d30| 20 20 20 20 20 20 20 2a | 2f 0a 58 2f 2a 20 20 74 | *|/.X/* t|
|00002d40| 68 61 6e 20 74 68 69 73 | 20 63 61 6c 63 75 6c 61 |han this| calcula|
|00002d50| 74 69 6f 6e 20 72 65 76 | 65 61 6c 73 2e 09 09 09 |tion rev|eals....|
|00002d60| 20 20 20 20 20 20 20 2a | 2f 0a 58 2f 2a 20 20 20 | *|/.X/* |
|00002d70| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002d80| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002d90| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002da0| 20 20 20 20 20 20 20 20 | 20 20 2a 2f 0a 58 2f 2a | | */.X/*|
|00002db0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002dc0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002dd0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00002de0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2f 0a |********|******/.|
|00002df0| 58 23 69 66 64 65 66 20 | 48 41 56 45 5f 50 52 4f |X#ifdef |HAVE_PRO|
|00002e00| 54 4f 53 0a 58 50 52 49 | 56 41 54 45 20 64 6f 75 |TOS.XPRI|VATE dou|
|00002e10| 62 6c 65 20 6d 65 61 6e | 70 68 61 73 65 28 64 6f |ble mean|phase(do|
|00002e20| 75 62 6c 65 20 73 64 61 | 74 65 2c 20 64 6f 75 62 |uble sda|te, doub|
|00002e30| 6c 65 20 70 68 61 73 65 | 2c 20 64 6f 75 62 6c 65 |le phase|, double|
|00002e40| 20 2a 75 73 65 6b 29 0a | 58 23 65 6c 73 65 0a 58 | *usek).|X#else.X|
|00002e50| 73 74 61 74 69 63 20 64 | 6f 75 62 6c 65 20 6d 65 |static d|ouble me|
|00002e60| 61 6e 70 68 61 73 65 28 | 73 64 61 74 65 2c 20 70 |anphase(|sdate, p|
|00002e70| 68 61 73 65 2c 20 75 73 | 65 6b 29 0a 58 64 6f 75 |hase, us|ek).Xdou|
|00002e80| 62 6c 65 20 73 64 61 74 | 65 2c 20 70 68 61 73 65 |ble sdat|e, phase|
|00002e90| 3b 0a 58 64 6f 75 62 6c | 65 20 2a 75 73 65 6b 3b |;.Xdoubl|e *usek;|
|00002ea0| 0a 58 23 65 6e 64 69 66 | 0a 58 7b 0a 58 20 20 20 |.X#endif|.X{.X |
|00002eb0| 69 6e 74 20 79 79 2c 20 | 6d 6d 2c 20 64 64 3b 0a |int yy, |mm, dd;.|
|00002ec0| 58 20 20 20 64 6f 75 62 | 6c 65 20 6b 2c 20 74 2c |X doub|le k, t,|
|00002ed0| 20 74 32 2c 20 74 33 2c | 20 6e 74 31 3b 0a 58 0a | t2, t3,| nt1;.X.|
|00002ee0| 58 20 20 20 6a 79 65 61 | 72 28 73 64 61 74 65 2c |X jyea|r(sdate,|
|00002ef0| 20 26 79 79 2c 20 26 6d | 6d 2c 20 26 64 64 29 3b | &yy, &m|m, &dd);|
|00002f00| 0a 58 0a 58 20 20 20 6b | 20 3d 20 28 79 79 20 2b |.X.X k| = (yy +|
|00002f10| 20 28 28 6d 6d 20 2d 20 | 31 29 20 2a 20 28 31 2e | ((mm - |1) * (1.|
|00002f20| 30 20 2f 20 31 32 2e 30 | 29 29 20 2d 20 31 39 30 |0 / 12.0|)) - 190|
|00002f30| 30 29 20 2a 20 31 32 2e | 33 36 38 35 3b 0a 58 0a |0) * 12.|3685;.X.|
|00002f40| 58 20 20 20 2f 2a 20 54 | 69 6d 65 20 69 6e 20 4a |X /* T|ime in J|
|00002f50| 75 6c 69 61 6e 20 63 65 | 6e 74 75 72 69 65 73 20 |ulian ce|nturies |
|00002f60| 66 72 6f 6d 20 31 39 30 | 30 20 4a 61 6e 75 61 72 |from 190|0 Januar|
|00002f70| 79 20 30 2e 35 20 2a 2f | 0a 58 20 20 20 74 20 3d |y 0.5 */|.X t =|
|00002f80| 20 28 73 64 61 74 65 20 | 2d 20 32 34 31 35 30 32 | (sdate |- 241502|
|00002f90| 30 2e 30 29 20 2f 20 33 | 36 35 32 35 3b 0a 58 20 |0.0) / 3|6525;.X |
|00002fa0| 20 20 74 32 20 3d 20 74 | 20 2a 20 74 3b 09 09 20 | t2 = t| * t;.. |
|00002fb0| 20 20 2f 2a 20 53 71 75 | 61 72 65 20 66 6f 72 20 | /* Squ|are for |
|00002fc0| 66 72 65 71 75 65 6e 74 | 20 75 73 65 20 2a 2f 0a |frequent| use */.|
|00002fd0| 58 20 20 20 74 33 20 3d | 20 74 32 20 2a 20 74 3b |X t3 =| t2 * t;|
|00002fe0| 09 09 20 20 20 2f 2a 20 | 43 75 62 65 20 66 6f 72 |.. /* |Cube for|
|00002ff0| 20 66 72 65 71 75 65 6e | 74 20 75 73 65 20 2a 2f | frequen|t use */|
|00003000| 0a 58 0a 58 20 20 20 2a | 75 73 65 6b 20 3d 20 6b |.X.X *|usek = k|
|00003010| 20 3d 20 66 6c 6f 6f 72 | 28 6b 29 20 2b 20 70 68 | = floor|(k) + ph|
|00003020| 61 73 65 3b 0a 58 20 20 | 20 6e 74 31 20 3d 20 32 |ase;.X | nt1 = 2|
|00003030| 34 31 35 30 32 30 2e 37 | 35 39 33 33 20 2b 20 73 |415020.7|5933 + s|
|00003040| 79 6e 6d 6f 6e 74 68 20 | 2a 20 6b 0a 58 20 20 20 |ynmonth |* k.X |
|00003050| 20 20 20 20 20 20 2b 20 | 30 2e 30 30 30 31 31 37 | + |0.000117|
|00003060| 38 20 2a 20 74 32 0a 58 | 20 20 20 20 20 20 20 20 |8 * t2.X| |
|00003070| 20 2d 20 30 2e 30 30 30 | 30 30 30 31 35 35 20 2a | - 0.000|000155 *|
|00003080| 20 74 33 0a 58 20 20 20 | 20 20 20 20 20 20 2b 20 | t3.X | + |
|00003090| 30 2e 30 30 30 33 33 20 | 2a 20 64 73 69 6e 28 31 |0.00033 |* dsin(1|
|000030a0| 36 36 2e 35 36 20 2b 20 | 31 33 32 2e 38 37 20 2a |66.56 + |132.87 *|
|000030b0| 20 74 20 2d 20 30 2e 30 | 30 39 31 37 33 20 2a 20 | t - 0.0|09173 * |
|000030c0| 74 32 29 3b 0a 58 0a 58 | 20 20 20 72 65 74 75 72 |t2);.X.X| retur|
|000030d0| 6e 20 6e 74 31 3b 0a 58 | 7d 0a 58 0a 58 2f 2a 2a |n nt1;.X|}.X.X/**|
|000030e0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000030f0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003100| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003110| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2f 0a 58 |********|*****/.X|
|00003120| 2f 2a 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |/* | |
|00003130| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003140| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003150| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 2a | | *|
|00003160| 2f 0a 58 2f 2a 20 20 74 | 72 75 65 70 68 61 73 65 |/.X/* t|ruephase|
|00003170| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003180| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003190| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000031a0| 20 20 2a 2f 0a 58 2f 2a | 20 20 20 20 20 20 20 20 | */.X/*| |
|000031b0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000031c0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000031d0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000031e0| 20 20 20 20 20 2a 2f 0a | 58 2f 2a 20 20 47 69 76 | */.|X/* Giv|
|000031f0| 65 6e 20 61 20 4b 20 76 | 61 6c 75 65 20 75 73 65 |en a K v|alue use|
|00003200| 64 20 74 6f 20 64 65 74 | 65 72 6d 69 6e 65 20 74 |d to det|ermine t|
|00003210| 68 65 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |he | |
|00003220| 20 20 20 20 20 20 20 20 | 2a 2f 0a 58 2f 2a 20 20 | |*/.X/* |
|00003230| 6d 65 61 6e 20 70 68 61 | 73 65 20 6f 66 20 74 68 |mean pha|se of th|
|00003240| 65 20 6e 65 77 20 6d 6f | 6f 6e 2c 20 61 6e 64 20 |e new mo|on, and |
|00003250| 61 20 70 68 61 73 65 20 | 20 20 20 20 20 20 20 20 |a phase | |
|00003260| 20 20 20 20 20 20 20 20 | 20 20 20 2a 2f 0a 58 2f | | */.X/|
|00003270| 2a 20 20 73 65 6c 65 63 | 74 6f 72 20 28 30 2e 30 |* selec|tor (0.0|
|00003280| 2c 20 30 2e 32 35 2c 20 | 30 2e 35 2c 20 30 2e 37 |, 0.25, |0.5, 0.7|
|00003290| 35 29 2c 20 6f 62 74 61 | 69 6e 20 20 20 20 20 20 |5), obta|in |
|000032a0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 2a 2f | | */|
|000032b0| 0a 58 2f 2a 20 20 74 68 | 65 20 74 72 75 65 2c 20 |.X/* th|e true, |
|000032c0| 63 6f 72 72 65 63 74 65 | 64 20 70 68 61 73 65 20 |correcte|d phase |
|000032d0| 74 69 6d 65 2e 20 20 20 | 20 20 20 20 20 20 20 20 |time. | |
|000032e0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000032f0| 20 2a 2f 0a 58 2f 2a 20 | 20 20 20 20 20 20 20 20 | */.X/* | |
|00003300| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003310| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003320| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003330| 20 20 20 20 2a 2f 0a 58 | 2f 2a 2a 2a 2a 2a 2a 2a | */.X|/*******|
|00003340| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003350| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003360| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003370| 2a 2a 2a 2a 2a 2a 2a 2a | 2f 0a 58 23 69 66 64 65 |********|/.X#ifde|
|00003380| 66 20 48 41 56 45 5f 50 | 52 4f 54 4f 53 0a 58 50 |f HAVE_P|ROTOS.XP|
|00003390| 52 49 56 41 54 45 20 64 | 6f 75 62 6c 65 20 74 72 |RIVATE d|ouble tr|
|000033a0| 75 65 70 68 61 73 65 28 | 64 6f 75 62 6c 65 20 6b |uephase(|double k|
|000033b0| 2c 20 64 6f 75 62 6c 65 | 20 70 68 61 73 65 29 0a |, double| phase).|
|000033c0| 58 23 65 6c 73 65 0a 58 | 73 74 61 74 69 63 20 64 |X#else.X|static d|
|000033d0| 6f 75 62 6c 65 20 74 72 | 75 65 70 68 61 73 65 28 |ouble tr|uephase(|
|000033e0| 6b 2c 20 70 68 61 73 65 | 29 0a 58 64 6f 75 62 6c |k, phase|).Xdoubl|
|000033f0| 65 20 6b 2c 20 70 68 61 | 73 65 3b 0a 58 23 65 6e |e k, pha|se;.X#en|
|00003400| 64 69 66 0a 58 7b 0a 58 | 20 20 20 64 6f 75 62 6c |dif.X{.X| doubl|
|00003410| 65 20 74 2c 20 74 32 2c | 20 74 33 2c 20 70 74 2c |e t, t2,| t3, pt,|
|00003420| 20 6d 2c 20 6d 70 72 69 | 6d 65 2c 20 66 3b 0a 58 | m, mpri|me, f;.X|
|00003430| 20 20 20 69 6e 74 20 61 | 70 63 6f 72 20 3d 20 30 | int a|pcor = 0|
|00003440| 3b 0a 58 0a 58 20 20 20 | 6b 20 2b 3d 20 70 68 61 |;.X.X |k += pha|
|00003450| 73 65 3b 09 09 20 20 20 | 2f 2a 20 41 64 64 20 70 |se;.. |/* Add p|
|00003460| 68 61 73 65 20 74 6f 20 | 6e 65 77 20 6d 6f 6f 6e |hase to |new moon|
|00003470| 20 74 69 6d 65 20 2a 2f | 0a 58 20 20 20 74 20 3d | time */|.X t =|
|00003480| 20 6b 20 2f 20 31 32 33 | 36 2e 38 35 3b 09 20 20 | k / 123|6.85;. |
|00003490| 20 2f 2a 20 54 69 6d 65 | 20 69 6e 20 4a 75 6c 69 | /* Time| in Juli|
|000034a0| 61 6e 20 63 65 6e 74 75 | 72 69 65 73 20 66 72 6f |an centu|ries fro|
|000034b0| 6d 0a 58 20 20 20 09 09 | 09 20 20 20 20 20 20 31 |m.X ..|. 1|
|000034c0| 39 30 30 20 4a 61 6e 75 | 61 72 79 20 30 2e 35 20 |900 Janu|ary 0.5 |
|000034d0| 2a 2f 0a 58 20 20 20 74 | 32 20 3d 20 74 20 2a 20 |*/.X t|2 = t * |
|000034e0| 74 3b 09 09 20 20 20 2f | 2a 20 53 71 75 61 72 65 |t;.. /|* Square|
|000034f0| 20 66 6f 72 20 66 72 65 | 71 75 65 6e 74 20 75 73 | for fre|quent us|
|00003500| 65 20 2a 2f 0a 58 20 20 | 20 74 33 20 3d 20 74 32 |e */.X | t3 = t2|
|00003510| 20 2a 20 74 3b 09 09 20 | 20 20 2f 2a 20 43 75 62 | * t;.. | /* Cub|
|00003520| 65 20 66 6f 72 20 66 72 | 65 71 75 65 6e 74 20 75 |e for fr|equent u|
|00003530| 73 65 20 2a 2f 0a 58 20 | 20 20 70 74 20 3d 20 32 |se */.X | pt = 2|
|00003540| 34 31 35 30 32 30 2e 37 | 35 39 33 33 09 20 20 20 |415020.7|5933. |
|00003550| 2f 2a 20 4d 65 61 6e 20 | 74 69 6d 65 20 6f 66 20 |/* Mean |time of |
|00003560| 70 68 61 73 65 20 2a 2f | 0a 58 20 20 20 20 20 20 |phase */|.X |
|00003570| 20 20 2b 20 73 79 6e 6d | 6f 6e 74 68 20 2a 20 6b | + synm|onth * k|
|00003580| 0a 58 20 20 20 20 20 20 | 20 20 2b 20 30 2e 30 30 |.X | + 0.00|
|00003590| 30 31 31 37 38 20 2a 20 | 74 32 0a 58 20 20 20 20 |01178 * |t2.X |
|000035a0| 20 20 20 20 2d 20 30 2e | 30 30 30 30 30 30 31 35 | - 0.|00000015|
|000035b0| 35 20 2a 20 74 33 0a 58 | 20 20 20 20 20 20 20 20 |5 * t3.X| |
|000035c0| 2b 20 30 2e 30 30 30 33 | 33 20 2a 20 64 73 69 6e |+ 0.0003|3 * dsin|
|000035d0| 28 31 36 36 2e 35 36 20 | 2b 20 31 33 32 2e 38 37 |(166.56 |+ 132.87|
|000035e0| 20 2a 20 74 20 2d 20 30 | 2e 30 30 39 31 37 33 20 | * t - 0|.009173 |
|000035f0| 2a 20 74 32 29 3b 0a 58 | 0a 58 20 20 20 20 20 20 |* t2);.X|.X |
|00003600| 20 20 6d 20 3d 20 33 35 | 39 2e 32 32 34 32 20 20 | m = 35|9.2242 |
|00003610| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 2f 2a 20 | | /* |
|00003620| 53 75 6e 27 73 20 6d 65 | 61 6e 20 61 6e 6f 6d 61 |Sun's me|an anoma|
|00003630| 6c 79 20 2a 2f 0a 58 20 | 20 20 20 20 20 20 2b 20 |ly */.X | + |
|00003640| 32 39 2e 31 30 35 33 35 | 36 30 38 20 2a 20 6b 0a |29.10535|608 * k.|
|00003650| 58 20 20 20 20 20 20 20 | 2d 20 30 2e 30 30 30 30 |X |- 0.0000|
|00003660| 33 33 33 20 2a 20 74 32 | 0a 58 20 20 20 20 20 20 |333 * t2|.X |
|00003670| 20 2d 20 30 2e 30 30 30 | 30 30 33 34 37 20 2a 20 | - 0.000|00347 * |
|00003680| 74 33 3b 0a 58 20 20 20 | 20 20 20 20 20 6d 70 72 |t3;.X | mpr|
|00003690| 69 6d 65 20 3d 20 33 30 | 36 2e 30 32 35 33 20 20 |ime = 30|6.0253 |
|000036a0| 20 20 20 20 20 20 20 20 | 2f 2a 20 4d 6f 6f 6e 27 | |/* Moon'|
|000036b0| 73 20 6d 65 61 6e 20 61 | 6e 6f 6d 61 6c 79 20 2a |s mean a|nomaly *|
|000036c0| 2f 0a 58 20 20 20 20 20 | 20 20 2b 20 33 38 35 2e |/.X | + 385.|
|000036d0| 38 31 36 39 31 38 30 36 | 20 2a 20 6b 0a 58 20 20 |81691806| * k.X |
|000036e0| 20 20 20 20 20 2b 20 30 | 2e 30 31 30 37 33 30 36 | + 0|.0107306|
|000036f0| 20 2a 20 74 32 0a 58 20 | 20 20 20 20 20 20 2b 20 | * t2.X | + |
|00003700| 30 2e 30 30 30 30 31 32 | 33 36 20 2a 20 74 33 3b |0.000012|36 * t3;|
|00003710| 0a 58 20 20 20 20 20 20 | 20 20 66 20 3d 20 32 31 |.X | f = 21|
|00003720| 2e 32 39 36 34 20 20 20 | 20 20 20 20 20 20 20 20 |.2964 | |
|00003730| 20 20 20 20 20 2f 2a 20 | 4d 6f 6f 6e 27 73 20 61 | /* |Moon's a|
|00003740| 72 67 75 6d 65 6e 74 20 | 6f 66 20 6c 61 74 69 74 |rgument |of latit|
|00003750| 75 64 65 20 2a 2f 0a 58 | 20 20 20 20 20 20 20 2b |ude */.X| +|
|00003760| 20 33 39 30 2e 36 37 30 | 35 30 36 34 36 20 2a 20 | 390.670|50646 * |
|00003770| 6b 0a 58 20 20 20 20 20 | 20 20 2d 20 30 2e 30 30 |k.X | - 0.00|
|00003780| 31 36 35 32 38 20 2a 20 | 74 32 0a 58 20 20 20 20 |16528 * |t2.X |
|00003790| 20 20 20 2d 20 30 2e 30 | 30 30 30 30 32 33 39 20 | - 0.0|0000239 |
|000037a0| 2a 20 74 33 3b 0a 58 20 | 20 20 69 66 20 28 28 70 |* t3;.X | if ((p|
|000037b0| 68 61 73 65 20 3c 20 30 | 2e 30 31 29 20 7c 7c 20 |hase < 0|.01) || |
|000037c0| 28 61 62 73 28 70 68 61 | 73 65 20 2d 20 30 2e 35 |(abs(pha|se - 0.5|
|000037d0| 29 20 3c 20 30 2e 30 31 | 29 29 20 7b 0a 58 0a 58 |) < 0.01|)) {.X.X|
|000037e0| 20 20 20 20 20 20 2f 2a | 20 43 6f 72 72 65 63 74 | /*| Correct|
|000037f0| 69 6f 6e 73 20 66 6f 72 | 20 4e 65 77 20 61 6e 64 |ions for| New and|
|00003800| 20 46 75 6c 6c 20 4d 6f | 6f 6e 20 2a 2f 0a 58 0a | Full Mo|on */.X.|
|00003810| 58 20 20 20 20 20 20 70 | 74 20 2b 3d 20 20 20 20 |X p|t += |
|00003820| 20 28 30 2e 31 37 33 34 | 20 2d 20 30 2e 30 30 30 | (0.1734| - 0.000|
|00003830| 33 39 33 20 2a 20 74 29 | 20 2a 20 64 73 69 6e 28 |393 * t)| * dsin(|
|00003840| 6d 29 0a 58 20 20 20 09 | 20 20 20 20 2b 20 30 2e |m).X .| + 0.|
|00003850| 30 30 32 31 20 2a 20 64 | 73 69 6e 28 32 20 2a 20 |0021 * d|sin(2 * |
|00003860| 6d 29 0a 58 20 20 20 09 | 20 20 20 20 2d 20 30 2e |m).X .| - 0.|
|00003870| 34 30 36 38 20 2a 20 64 | 73 69 6e 28 6d 70 72 69 |4068 * d|sin(mpri|
|00003880| 6d 65 29 0a 58 20 20 20 | 09 20 20 20 20 2b 20 30 |me).X |. + 0|
|00003890| 2e 30 31 36 31 20 2a 20 | 64 73 69 6e 28 32 20 2a |.0161 * |dsin(2 *|
|000038a0| 20 6d 70 72 69 6d 65 29 | 0a 58 20 20 20 09 20 20 | mprime)|.X . |
|000038b0| 20 20 2d 20 30 2e 30 30 | 30 34 20 2a 20 64 73 69 | - 0.00|04 * dsi|
|000038c0| 6e 28 33 20 2a 20 6d 70 | 72 69 6d 65 29 0a 58 20 |n(3 * mp|rime).X |
|000038d0| 20 20 09 20 20 20 20 2b | 20 30 2e 30 31 30 34 20 | . +| 0.0104 |
|000038e0| 2a 20 64 73 69 6e 28 32 | 20 2a 20 66 29 0a 58 20 |* dsin(2| * f).X |
|000038f0| 20 20 09 20 20 20 20 2d | 20 30 2e 30 30 35 31 20 | . -| 0.0051 |
|00003900| 2a 20 64 73 69 6e 28 6d | 20 2b 20 6d 70 72 69 6d |* dsin(m| + mprim|
|00003910| 65 29 0a 58 20 20 20 09 | 20 20 20 20 2d 20 30 2e |e).X .| - 0.|
|00003920| 30 30 37 34 20 2a 20 64 | 73 69 6e 28 6d 20 2d 20 |0074 * d|sin(m - |
|00003930| 6d 70 72 69 6d 65 29 0a | 58 20 20 20 09 20 20 20 |mprime).|X . |
|00003940| 20 2b 20 30 2e 30 30 30 | 34 20 2a 20 64 73 69 6e | + 0.000|4 * dsin|
|00003950| 28 32 20 2a 20 66 20 2b | 20 6d 29 0a 58 20 20 20 |(2 * f +| m).X |
|00003960| 09 20 20 20 20 2d 20 30 | 2e 30 30 30 34 20 2a 20 |. - 0|.0004 * |
|00003970| 64 73 69 6e 28 32 20 2a | 20 66 20 2d 20 6d 29 0a |dsin(2 *| f - m).|
|00003980| 58 20 20 20 09 20 20 20 | 20 2d 20 30 2e 30 30 30 |X . | - 0.000|
|00003990| 36 20 2a 20 64 73 69 6e | 28 32 20 2a 20 66 20 2b |6 * dsin|(2 * f +|
|000039a0| 20 6d 70 72 69 6d 65 29 | 0a 58 20 20 20 09 20 20 | mprime)|.X . |
|000039b0| 20 20 2b 20 30 2e 30 30 | 31 30 20 2a 20 64 73 69 | + 0.00|10 * dsi|
|000039c0| 6e 28 32 20 2a 20 66 20 | 2d 20 6d 70 72 69 6d 65 |n(2 * f |- mprime|
|000039d0| 29 0a 58 20 20 20 09 20 | 20 20 20 2b 20 30 2e 30 |).X . | + 0.0|
|000039e0| 30 30 35 20 2a 20 64 73 | 69 6e 28 6d 20 2b 20 32 |005 * ds|in(m + 2|
|000039f0| 20 2a 20 6d 70 72 69 6d | 65 29 3b 0a 58 20 20 20 | * mprim|e);.X |
|00003a00| 20 20 20 61 70 63 6f 72 | 20 3d 20 31 3b 0a 58 20 | apcor| = 1;.X |
|00003a10| 20 20 7d 20 65 6c 73 65 | 20 69 66 20 28 28 61 62 | } else| if ((ab|
|00003a20| 73 28 70 68 61 73 65 20 | 2d 20 30 2e 32 35 29 20 |s(phase |- 0.25) |
|00003a30| 3c 20 30 2e 30 31 20 7c | 7c 20 28 61 62 73 28 70 |< 0.01 ||| (abs(p|
|00003a40| 68 61 73 65 20 2d 20 30 | 2e 37 35 29 20 3c 20 30 |hase - 0|.75) < 0|
|00003a50| 2e 30 31 29 29 29 20 7b | 0a 58 20 20 20 20 20 20 |.01))) {|.X |
|00003a60| 70 74 20 2b 3d 20 20 20 | 20 20 28 30 2e 31 37 32 |pt += | (0.172|
|00003a70| 31 20 2d 20 30 2e 30 30 | 30 34 20 2a 20 74 29 20 |1 - 0.00|04 * t) |
|00003a80| 2a 20 64 73 69 6e 28 6d | 29 0a 58 20 20 20 09 20 |* dsin(m|).X . |
|00003a90| 20 20 20 2b 20 30 2e 30 | 30 32 31 20 2a 20 64 73 | + 0.0|021 * ds|
|00003aa0| 69 6e 28 32 20 2a 20 6d | 29 0a 58 20 20 20 09 20 |in(2 * m|).X . |
|00003ab0| 20 20 20 2d 20 30 2e 36 | 32 38 30 20 2a 20 64 73 | - 0.6|280 * ds|
|00003ac0| 69 6e 28 6d 70 72 69 6d | 65 29 0a 58 20 20 20 09 |in(mprim|e).X .|
|00003ad0| 20 20 20 20 2b 20 30 2e | 30 30 38 39 20 2a 20 64 | + 0.|0089 * d|
|00003ae0| 73 69 6e 28 32 20 2a 20 | 6d 70 72 69 6d 65 29 0a |sin(2 * |mprime).|
|00003af0| 58 20 20 20 09 20 20 20 | 20 2d 20 30 2e 30 30 30 |X . | - 0.000|
|00003b00| 34 20 2a 20 64 73 69 6e | 28 33 20 2a 20 6d 70 72 |4 * dsin|(3 * mpr|
|00003b10| 69 6d 65 29 0a 58 20 20 | 20 09 20 20 20 20 2b 20 |ime).X | . + |
|00003b20| 30 2e 30 30 37 39 20 2a | 20 64 73 69 6e 28 32 20 |0.0079 *| dsin(2 |
|00003b30| 2a 20 66 29 0a 58 20 20 | 20 09 20 20 20 20 2d 20 |* f).X | . - |
|00003b40| 30 2e 30 31 31 39 20 2a | 20 64 73 69 6e 28 6d 20 |0.0119 *| dsin(m |
|00003b50| 2b 20 6d 70 72 69 6d 65 | 29 0a 58 20 20 20 09 20 |+ mprime|).X . |
|00003b60| 20 20 20 2d 20 30 2e 30 | 30 34 37 20 2a 20 64 73 | - 0.0|047 * ds|
|00003b70| 69 6e 28 6d 20 2d 20 6d | 70 72 69 6d 65 29 0a 58 |in(m - m|prime).X|
|00003b80| 20 20 20 09 20 20 20 20 | 2b 20 30 2e 30 30 30 33 | . |+ 0.0003|
|00003b90| 20 2a 20 64 73 69 6e 28 | 32 20 2a 20 66 20 2b 20 | * dsin(|2 * f + |
|00003ba0| 6d 29 0a 58 20 20 20 09 | 20 20 20 20 2d 20 30 2e |m).X .| - 0.|
|00003bb0| 30 30 30 34 20 2a 20 64 | 73 69 6e 28 32 20 2a 20 |0004 * d|sin(2 * |
|00003bc0| 66 20 2d 20 6d 29 0a 58 | 20 20 20 09 20 20 20 20 |f - m).X| . |
|00003bd0| 2d 20 30 2e 30 30 30 36 | 20 2a 20 64 73 69 6e 28 |- 0.0006| * dsin(|
|00003be0| 32 20 2a 20 66 20 2b 20 | 6d 70 72 69 6d 65 29 0a |2 * f + |mprime).|
|00003bf0| 58 20 20 20 09 20 20 20 | 20 2b 20 30 2e 30 30 32 |X . | + 0.002|
|00003c00| 31 20 2a 20 64 73 69 6e | 28 32 20 2a 20 66 20 2d |1 * dsin|(2 * f -|
|00003c10| 20 6d 70 72 69 6d 65 29 | 0a 58 20 20 20 09 20 20 | mprime)|.X . |
|00003c20| 20 20 2b 20 30 2e 30 30 | 30 33 20 2a 20 64 73 69 | + 0.00|03 * dsi|
|00003c30| 6e 28 6d 20 2b 20 32 20 | 2a 20 6d 70 72 69 6d 65 |n(m + 2 |* mprime|
|00003c40| 29 0a 58 20 20 20 09 20 | 20 20 20 2b 20 30 2e 30 |).X . | + 0.0|
|00003c50| 30 30 34 20 2a 20 64 73 | 69 6e 28 6d 20 2d 20 32 |004 * ds|in(m - 2|
|00003c60| 20 2a 20 6d 70 72 69 6d | 65 29 0a 58 20 20 20 09 | * mprim|e).X .|
|00003c70| 20 20 20 20 2d 20 30 2e | 30 30 30 33 20 2a 20 64 | - 0.|0003 * d|
|00003c80| 73 69 6e 28 32 20 2a 20 | 6d 20 2b 20 6d 70 72 69 |sin(2 * |m + mpri|
|00003c90| 6d 65 29 3b 0a 58 20 20 | 20 20 20 20 69 66 20 28 |me);.X | if (|
|00003ca0| 70 68 61 73 65 20 3c 20 | 30 2e 35 29 0a 58 20 20 |phase < |0.5).X |
|00003cb0| 20 20 20 20 20 20 20 2f | 2a 20 46 69 72 73 74 20 | /|* First |
|00003cc0| 71 75 61 72 74 65 72 20 | 63 6f 72 72 65 63 74 69 |quarter |correcti|
|00003cd0| 6f 6e 20 2a 2f 0a 58 20 | 20 20 20 20 20 20 20 20 |on */.X | |
|00003ce0| 70 74 20 2b 3d 20 30 2e | 30 30 32 38 20 2d 20 30 |pt += 0.|0028 - 0|
|00003cf0| 2e 30 30 30 34 20 2a 20 | 64 63 6f 73 28 6d 29 20 |.0004 * |dcos(m) |
|00003d00| 2b 20 30 2e 30 30 30 33 | 20 2a 20 64 63 6f 73 28 |+ 0.0003| * dcos(|
|00003d10| 6d 70 72 69 6d 65 29 3b | 0a 58 20 20 20 20 20 20 |mprime);|.X |
|00003d20| 65 6c 73 65 0a 58 20 20 | 20 20 20 20 20 20 20 2f |else.X | /|
|00003d30| 2a 20 4c 61 73 74 20 71 | 75 61 72 74 65 72 20 63 |* Last q|uarter c|
|00003d40| 6f 72 72 65 63 74 69 6f | 6e 20 2a 2f 0a 58 20 20 |orrectio|n */.X |
|00003d50| 20 20 20 20 20 20 20 70 | 74 20 2b 3d 20 2d 30 2e | p|t += -0.|
|00003d60| 30 30 32 38 20 2b 20 30 | 2e 30 30 30 34 20 2a 20 |0028 + 0|.0004 * |
|00003d70| 64 63 6f 73 28 6d 29 20 | 2d 20 30 2e 30 30 30 33 |dcos(m) |- 0.0003|
|00003d80| 20 2a 20 64 63 6f 73 28 | 6d 70 72 69 6d 65 29 3b | * dcos(|mprime);|
|00003d90| 0a 58 20 20 20 20 20 20 | 61 70 63 6f 72 20 3d 20 |.X |apcor = |
|00003da0| 31 3b 0a 58 20 20 20 7d | 0a 58 20 20 20 69 66 20 |1;.X }|.X if |
|00003db0| 28 21 61 70 63 6f 72 29 | 20 72 65 74 75 72 6e 20 |(!apcor)| return |
|00003dc0| 30 2e 30 3b 0a 58 20 20 | 20 72 65 74 75 72 6e 20 |0.0;.X | return |
|00003dd0| 70 74 3b 0a 58 7d 0a 58 | 0a 58 2f 2a 2a 2a 2a 2a |pt;.X}.X|.X/*****|
|00003de0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003df0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003e00| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003e10| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2f 0a 58 2f 2a 20 |********|**/.X/* |
|00003e20| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003e30| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003e40| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003e50| 20 20 20 20 20 20 20 20 | 20 20 20 20 2a 2f 0a 58 | | */.X|
|00003e60| 2f 2a 20 20 6b 65 70 6c | 65 72 20 20 20 20 20 20 |/* kepl|er |
|00003e70| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003e80| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003e90| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 2a | | *|
|00003ea0| 2f 0a 58 2f 2a 20 20 20 | 20 20 20 20 20 20 20 20 |/.X/* | |
|00003eb0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003ec0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003ed0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003ee0| 20 20 2a 2f 0a 58 2f 2a | 20 20 53 6f 6c 76 65 20 | */.X/*| Solve |
|00003ef0| 74 68 65 20 65 71 75 61 | 74 69 6f 6e 20 6f 66 20 |the equa|tion of |
|00003f00| 4b 65 70 6c 65 72 2e 20 | 20 20 20 20 20 20 20 20 |Kepler. | |
|00003f10| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003f20| 20 20 20 20 20 2a 2f 0a | 58 2f 2a 20 20 20 20 20 | */.|X/* |
|00003f30| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003f40| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003f50| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003f60| 20 20 20 20 20 20 20 20 | 2a 2f 0a 58 2f 2a 2a 2a | |*/.X/***|
|00003f70| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003f80| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003f90| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00003fa0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2f 0a 58 23 |********|****/.X#|
|00003fb0| 69 66 64 65 66 20 48 41 | 56 45 5f 50 52 4f 54 4f |ifdef HA|VE_PROTO|
|00003fc0| 53 0a 58 50 52 49 56 41 | 54 45 20 64 6f 75 62 6c |S.XPRIVA|TE doubl|
|00003fd0| 65 20 6b 65 70 6c 65 72 | 28 64 6f 75 62 6c 65 20 |e kepler|(double |
|00003fe0| 6d 2c 20 64 6f 75 62 6c | 65 20 65 63 63 29 0a 58 |m, doubl|e ecc).X|
|00003ff0| 23 65 6c 73 65 0a 58 73 | 74 61 74 69 63 20 64 6f |#else.Xs|tatic do|
|00004000| 75 62 6c 65 20 6b 65 70 | 6c 65 72 28 6d 2c 20 65 |uble kep|ler(m, e|
|00004010| 63 63 29 0a 58 64 6f 75 | 62 6c 65 20 6d 2c 20 65 |cc).Xdou|ble m, e|
|00004020| 63 63 3b 0a 58 23 65 6e | 64 69 66 0a 58 7b 0a 58 |cc;.X#en|dif.X{.X|
|00004030| 20 20 20 64 6f 75 62 6c | 65 20 65 2c 20 64 65 6c | doubl|e e, del|
|00004040| 74 61 3b 0a 58 23 64 65 | 66 69 6e 65 20 45 50 53 |ta;.X#de|fine EPS|
|00004050| 49 4c 4f 4e 20 31 45 2d | 36 0a 58 0a 58 20 20 20 |ILON 1E-|6.X.X |
|00004060| 65 20 3d 20 6d 20 3d 20 | 74 6f 72 61 64 28 6d 29 |e = m = |torad(m)|
|00004070| 3b 0a 58 20 20 20 64 6f | 20 7b 0a 58 20 20 20 20 |;.X do| {.X |
|00004080| 20 20 64 65 6c 74 61 20 | 3d 20 65 20 2d 20 65 63 | delta |= e - ec|
|00004090| 63 20 2a 20 73 69 6e 28 | 65 29 20 2d 20 6d 3b 0a |c * sin(|e) - m;.|
|000040a0| 58 20 20 20 20 20 20 65 | 20 2d 3d 20 64 65 6c 74 |X e| -= delt|
|000040b0| 61 20 2f 20 28 31 20 2d | 20 65 63 63 20 2a 20 63 |a / (1 -| ecc * c|
|000040c0| 6f 73 28 65 29 29 3b 0a | 58 20 20 20 7d 20 77 68 |os(e));.|X } wh|
|000040d0| 69 6c 65 20 28 61 62 73 | 28 64 65 6c 74 61 29 20 |ile (abs|(delta) |
|000040e0| 3e 20 45 50 53 49 4c 4f | 4e 29 3b 0a 58 20 20 20 |> EPSILO|N);.X |
|000040f0| 72 65 74 75 72 6e 20 65 | 3b 0a 58 7d 0a 58 2f 2a |return e|;.X}.X/*|
|00004100| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00004110| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00004120| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00004130| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2f 0a |********|******/.|
|00004140| 58 2f 2a 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |X/* | |
|00004150| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00004160| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00004170| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00004180| 2a 2f 0a 58 2f 2a 20 20 | 50 48 41 53 45 20 20 2d |*/.X/* |PHASE -|
|00004190| 2d 20 20 43 61 6c 63 75 | 6c 61 74 65 20 70 68 61 |- Calcu|late pha|
|000041a0| 73 65 20 6f 66 20 6d 6f | 6f 6e 20 61 73 20 61 20 |se of mo|on as a |
|000041b0| 66 72 61 63 74 69 6f 6e | 3a 20 20 20 20 20 20 20 |fraction|: |
|000041c0| 20 20 20 2a 2f 0a 58 2f | 2a 20 20 20 20 20 20 20 | */.X/|* |
|000041d0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000041e0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000041f0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00004200| 20 20 20 20 20 20 2a 2f | 0a 58 2f 2a 20 20 20 54 | */|.X/* T|
|00004210| 68 65 20 61 72 67 75 6d | 65 6e 74 20 69 73 20 74 |he argum|ent is t|
|00004220| 68 65 20 74 69 6d 65 20 | 66 6f 72 20 77 68 69 63 |he time |for whic|
|00004230| 68 20 74 68 65 20 70 68 | 61 73 65 20 69 73 20 20 |h the ph|ase is |
|00004240| 20 09 20 20 20 20 20 20 | 20 2a 2f 0a 58 2f 2a 20 | . | */.X/* |
|00004250| 20 20 52 65 71 75 65 73 | 74 65 64 2c 20 65 78 70 | Reques|ted, exp|
|00004260| 72 65 73 73 65 64 20 61 | 73 20 61 20 4a 75 6c 69 |ressed a|s a Juli|
|00004270| 61 6e 20 64 61 74 65 20 | 61 6e 64 09 09 20 20 20 |an date |and.. |
|00004280| 20 20 20 20 2a 2f 0a 58 | 2f 2a 20 20 20 66 72 61 | */.X|/* fra|
|00004290| 63 74 69 6f 6e 2e 20 20 | 52 65 74 75 72 6e 73 20 |ction. |Returns |
|000042a0| 74 68 65 20 74 65 72 6d | 69 6e 61 74 6f 72 20 70 |the term|inator p|
|000042b0| 68 61 73 65 20 61 6e 67 | 6c 65 20 61 73 20 61 09 |hase ang|le as a.|
|000042c0| 20 20 20 20 20 20 20 2a | 2f 0a 58 2f 2a 20 20 20 | *|/.X/* |
|000042d0| 70 65 72 63 65 6e 74 61 | 67 65 20 6f 66 20 61 20 |percenta|ge of a |
|000042e0| 66 75 6c 6c 20 63 69 72 | 63 6c 65 20 28 69 2e 65 |full cir|cle (i.e|
|000042f0| 2e 2c 20 30 20 74 6f 20 | 31 29 2c 20 61 6e 64 09 |., 0 to |1), and.|
|00004300| 20 20 20 20 20 20 20 2a | 2f 0a 58 2f 2a 20 20 20 | *|/.X/* |
|00004310| 73 74 6f 72 65 73 20 69 | 6e 74 6f 20 70 6f 69 6e |stores i|nto poin|
|00004320| 74 65 72 20 61 72 67 75 | 6d 65 6e 74 73 20 74 68 |ter argu|ments th|
|00004330| 65 20 69 6c 6c 75 6d 69 | 6e 61 74 65 64 09 20 20 |e illumi|nated. |
|00004340| 20 20 20 20 20 2a 2f 0a | 58 2f 2a 20 20 20 66 72 | */.|X/* fr|
|00004350| 61 63 74 69 6f 6e 20 6f | 66 20 74 68 65 20 4d 6f |action o|f the Mo|
|00004360| 6f 6e 27 73 20 64 69 73 | 63 2c 20 74 68 65 20 4d |on's dis|c, the M|
|00004370| 6f 6f 6e 27 73 20 61 67 | 65 20 69 6e 09 20 20 20 |oon's ag|e in. |
|00004380| 20 20 20 20 2a 2f 0a 58 | 2f 2a 20 20 20 64 61 79 | */.X|/* day|
|00004390| 73 20 61 6e 64 20 66 72 | 61 63 74 69 6f 6e 2c 20 |s and fr|action, |
|000043a0| 74 68 65 20 64 69 73 74 | 61 6e 63 65 20 6f 66 20 |the dist|ance of |
|000043b0| 74 68 65 20 4d 6f 6f 6e | 20 66 72 6f 6d 09 20 20 |the Moon| from. |
|000043c0| 20 20 20 20 20 2a 2f 0a | 58 2f 2a 20 20 20 74 68 | */.|X/* th|
|000043d0| 65 20 63 65 6e 74 72 65 | 20 6f 66 20 74 68 65 20 |e centre| of the |
|000043e0| 45 61 72 74 68 2c 20 61 | 6e 64 20 74 68 65 20 61 |Earth, a|nd the a|
|000043f0| 6e 67 75 6c 61 72 20 64 | 69 61 6d 65 74 65 72 09 |ngular d|iameter.|
|00004400| 20 20 20 20 20 20 20 2a | 2f 0a 58 2f 2a 20 20 20 | *|/.X/* |
|00004410| 73 75 62 74 65 6e 64 65 | 64 20 62 79 20 74 68 65 |subtende|d by the|
|00004420| 20 4d 6f 6f 6e 20 61 73 | 20 73 65 65 6e 20 62 79 | Moon as| seen by|
|00004430| 20 61 6e 20 6f 62 73 65 | 72 76 65 72 20 61 74 09 | an obse|rver at.|
|00004440| 20 20 20 20 20 20 20 2a | 2f 0a 58 2f 2a 20 20 20 | *|/.X/* |
|00004450| 74 68 65 20 63 65 6e 74 | 72 65 20 6f 66 20 74 68 |the cent|re of th|
|00004460| 65 20 45 61 72 74 68 2e | 09 09 09 09 20 20 20 20 |e Earth.|.... |
|00004470| 20 20 20 2a 2f 0a 58 2f | 2a 20 20 20 20 20 20 20 | */.X/|* |
|00004480| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00004490| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000044a0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000044b0| 20 20 20 20 20 20 2a 2f | 0a 58 2f 2a 2a 2a 2a 2a | */|.X/*****|
|000044c0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000044d0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000044e0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|000044f0| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2f 0a 58 23 69 66 |********|**/.X#if|
|00004500| 64 65 66 20 48 41 56 45 | 5f 50 52 4f 54 4f 53 0a |def HAVE|_PROTOS.|
|00004510| 58 50 52 49 56 41 54 45 | 20 64 6f 75 62 6c 65 20 |XPRIVATE| double |
|00004520| 70 68 61 73 65 28 64 6f | 75 62 6c 65 20 70 64 61 |phase(do|uble pda|
|00004530| 74 65 2c 0a 58 09 20 20 | 20 20 20 20 20 20 20 20 |te,.X. | |
|00004540| 20 20 20 64 6f 75 62 6c | 65 20 2a 70 70 68 61 73 | doubl|e *pphas|
|00004550| 65 2c 0a 58 09 09 20 20 | 20 20 20 64 6f 75 62 6c |e,.X.. | doubl|
|00004560| 65 20 2a 6d 61 67 65 2c | 0a 58 09 09 20 20 20 20 |e *mage,|.X.. |
|00004570| 20 64 6f 75 62 6c 65 20 | 2a 64 69 73 74 2c 0a 58 | double |*dist,.X|
|00004580| 09 09 20 20 20 20 20 64 | 6f 75 62 6c 65 20 2a 61 |.. d|ouble *a|
|00004590| 6e 67 64 69 61 2c 0a 58 | 09 09 20 20 20 20 20 64 |ngdia,.X|.. d|
|000045a0| 6f 75 62 6c 65 20 2a 73 | 75 64 69 73 74 2c 0a 58 |ouble *s|udist,.X|
|000045b0| 09 09 20 20 20 20 20 64 | 6f 75 62 6c 65 20 2a 73 |.. d|ouble *s|
|000045c0| 75 61 6e 67 64 69 61 29 | 0a 58 23 65 6c 73 65 0a |uangdia)|.X#else.|
|000045d0| 58 73 74 61 74 69 63 20 | 64 6f 75 62 6c 65 20 70 |Xstatic |double p|
|000045e0| 68 61 73 65 28 70 64 61 | 74 65 2c 20 70 70 68 61 |hase(pda|te, ppha|
|000045f0| 73 65 2c 20 6d 61 67 65 | 2c 20 64 69 73 74 2c 20 |se, mage|, dist, |
|00004600| 61 6e 67 64 69 61 2c 20 | 73 75 64 69 73 74 2c 20 |angdia, |sudist, |
|00004610| 73 75 61 6e 67 64 69 61 | 29 0a 58 64 6f 75 62 6c |suangdia|).Xdoubl|
|00004620| 65 20 70 64 61 74 65 3b | 0a 58 64 6f 75 62 6c 65 |e pdate;|.Xdouble|
|00004630| 20 2a 70 70 68 61 73 65 | 3b 20 09 09 20 20 20 2f | *pphase|; .. /|
|00004640| 2a 20 49 6c 6c 75 6d 69 | 6e 61 74 65 64 20 66 72 |* Illumi|nated fr|
|00004650| 61 63 74 69 6f 6e 20 2a | 2f 0a 58 64 6f 75 62 6c |action *|/.Xdoubl|
|00004660| 65 20 2a 6d 61 67 65 3b | 09 09 09 20 20 20 2f 2a |e *mage;|... /*|
|00004670| 20 41 67 65 20 6f 66 20 | 6d 6f 6f 6e 20 69 6e 20 | Age of |moon in |
|00004680| 64 61 79 73 20 2a 2f 0a | 58 64 6f 75 62 6c 65 20 |days */.|Xdouble |
|00004690| 2a 64 69 73 74 3b 09 09 | 09 20 20 20 2f 2a 20 44 |*dist;..|. /* D|
|000046a0| 69 73 74 61 6e 63 65 20 | 69 6e 20 6b 69 6c 6f 6d |istance |in kilom|
|000046b0| 65 74 72 65 73 20 2a 2f | 0a 58 64 6f 75 62 6c 65 |etres */|.Xdouble|
|000046c0| 20 2a 61 6e 67 64 69 61 | 3b 20 09 09 20 20 20 2f | *angdia|; .. /|
|000046d0| 2a 20 41 6e 67 75 6c 61 | 72 20 64 69 61 6d 65 74 |* Angula|r diamet|
|000046e0| 65 72 20 69 6e 20 64 65 | 67 72 65 65 73 20 2a 2f |er in de|grees */|
|000046f0| 0a 58 64 6f 75 62 6c 65 | 20 2a 73 75 64 69 73 74 |.Xdouble| *sudist|
|00004700| 3b 20 09 09 20 20 20 2f | 2a 20 44 69 73 74 61 6e |; .. /|* Distan|
|00004710| 63 65 20 74 6f 20 53 75 | 6e 20 2a 2f 0a 58 64 6f |ce to Su|n */.Xdo|
|00004720| 75 62 6c 65 20 2a 73 75 | 61 6e 67 64 69 61 3b 20 |uble *su|angdia; |
|00004730| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00004740| 20 2f 2a 20 53 75 6e 27 | 73 20 61 6e 67 75 6c 61 | /* Sun'|s angula|
|00004750| 72 20 64 69 61 6d 65 74 | 65 72 20 2a 2f 0a 58 23 |r diamet|er */.X#|
|00004760| 65 6e 64 69 66 0a 58 7b | 0a 58 0a 58 20 20 20 64 |endif.X{|.X.X d|
|00004770| 6f 75 62 6c 65 20 44 61 | 79 2c 20 4e 2c 20 4d 2c |ouble Da|y, N, M,|
|00004780| 20 45 63 2c 20 4c 61 6d | 62 64 61 73 75 6e 2c 20 | Ec, Lam|bdasun, |
|00004790| 6d 6c 2c 20 4d 4d 2c 20 | 4d 4e 2c 20 45 76 2c 20 |ml, MM, |MN, Ev, |
|000047a0| 41 65 2c 20 41 33 2c 20 | 4d 6d 50 2c 0a 58 20 20 |Ae, A3, |MmP,.X |
|000047b0| 20 20 20 20 20 20 20 20 | 6d 45 63 2c 20 41 34 2c | |mEc, A4,|
|000047c0| 20 6c 50 2c 20 56 2c 20 | 6c 50 50 2c 20 4e 50 2c | lP, V, |lPP, NP,|
|000047d0| 20 79 2c 20 78 2c 20 4c | 61 6d 62 64 61 6d 6f 6f | y, x, L|ambdamoo|
|000047e0| 6e 2c 0a 58 20 20 20 20 | 20 20 20 20 20 20 4d 6f |n,.X | Mo|
|000047f0| 6f 6e 41 67 65 2c 20 4d | 6f 6f 6e 50 68 61 73 65 |onAge, M|oonPhase|
|00004800| 2c 0a 58 20 20 20 20 20 | 20 20 20 20 20 4d 6f 6f |,.X | Moo|
|00004810| 6e 44 69 73 74 2c 20 4d | 6f 6f 6e 44 46 72 61 63 |nDist, M|oonDFrac|
|00004820| 2c 20 4d 6f 6f 6e 41 6e | 67 2c 0a 58 20 20 20 20 |, MoonAn|g,.X |
|00004830| 20 20 20 20 20 20 46 2c | 20 53 75 6e 44 69 73 74 | F,| SunDist|
|00004840| 2c 20 53 75 6e 41 6e 67 | 3b 0a 58 0a 58 20 20 20 |, SunAng|;.X.X |
|00004850| 20 20 20 20 20 2f 2a 20 | 43 61 6c 63 75 6c 61 74 | /* |Calculat|
|00004860| 69 6f 6e 20 6f 66 20 74 | 68 65 20 53 75 6e 27 73 |ion of t|he Sun's|
|00004870| 20 70 6f 73 69 74 69 6f | 6e 20 2a 2f 0a 58 0a 58 | positio|n */.X.X|
|00004880| 20 20 20 44 61 79 20 3d | 20 70 64 61 74 65 20 2d | Day =| pdate -|
|00004890| 20 65 70 6f 63 68 3b 09 | 20 20 20 20 2f 2a 20 44 | epoch;.| /* D|
|000048a0| 61 74 65 20 77 69 74 68 | 69 6e 20 65 70 6f 63 68 |ate with|in epoch|
|000048b0| 20 2a 2f 0a 58 20 20 20 | 4e 20 3d 20 66 69 78 61 | */.X |N = fixa|
|000048c0| 6e 67 6c 65 28 28 33 36 | 30 20 2f 20 33 36 35 2e |ngle((36|0 / 365.|
|000048d0| 32 34 32 32 29 20 2a 20 | 44 61 79 29 3b 20 2f 2a |2422) * |Day); /*|
|000048e0| 20 4d 65 61 6e 20 61 6e | 6f 6d 61 6c 79 20 6f 66 | Mean an|omaly of|
|000048f0| 20 74 68 65 20 53 75 6e | 20 2a 2f 0a 58 20 20 20 | the Sun| */.X |
|00004900| 4d 20 3d 20 66 69 78 61 | 6e 67 6c 65 28 4e 20 2b |M = fixa|ngle(N +|
|00004910| 20 65 6c 6f 6e 67 65 20 | 2d 20 65 6c 6f 6e 67 70 | elonge |- elongp|
|00004920| 29 3b 20 20 20 20 2f 2a | 20 43 6f 6e 76 65 72 74 |); /*| Convert|
|00004930| 20 66 72 6f 6d 20 70 65 | 72 69 67 65 65 0a 58 20 | from pe|rigee.X |
|00004940| 20 20 09 09 09 20 20 20 | 20 20 20 20 63 6f 2d 6f | ... | co-o|
|00004950| 72 64 69 6e 61 74 65 73 | 20 74 6f 20 65 70 6f 63 |rdinates| to epoc|
|00004960| 68 20 31 39 38 30 2e 30 | 20 2a 2f 0a 58 20 20 20 |h 1980.0| */.X |
|00004970| 45 63 20 3d 20 6b 65 70 | 6c 65 72 28 4d 2c 20 65 |Ec = kep|ler(M, e|
|00004980| 63 63 65 6e 74 29 3b 20 | 20 20 20 20 2f 2a 20 53 |ccent); | /* S|
|00004990| 6f 6c 76 65 20 65 71 75 | 61 74 69 6f 6e 20 6f 66 |olve equ|ation of|
|000049a0| 20 4b 65 70 6c 65 72 20 | 2a 2f 0a 58 20 20 20 45 | Kepler |*/.X E|
|000049b0| 63 20 3d 20 73 71 72 74 | 28 28 31 20 2b 20 65 63 |c = sqrt|((1 + ec|
|000049c0| 63 65 6e 74 29 20 2f 20 | 28 31 20 2d 20 65 63 63 |cent) / |(1 - ecc|
|000049d0| 65 6e 74 29 29 20 2a 20 | 74 61 6e 28 45 63 20 2f |ent)) * |tan(Ec /|
|000049e0| 20 32 29 3b 0a 58 20 20 | 20 45 63 20 3d 20 32 20 | 2);.X | Ec = 2 |
|000049f0| 2a 20 74 6f 64 65 67 28 | 61 74 61 6e 28 45 63 29 |* todeg(|atan(Ec)|
|00004a00| 29 3b 20 20 20 2f 2a 20 | 31 20 61 6e 6f 6d 61 6c |); /* |1 anomal|
|00004a10| 79 20 2a 2f 0a 58 20 20 | 20 4c 61 6d 62 64 61 73 |y */.X | Lambdas|
|00004a20| 75 6e 20 3d 20 66 69 78 | 61 6e 67 6c 65 28 45 63 |un = fix|angle(Ec|
|00004a30| 20 2b 20 65 6c 6f 6e 67 | 70 29 3b 20 20 2f 2a 20 | + elong|p); /* |
|00004a40| 53 75 6e 27 73 20 67 65 | 6f 63 65 6e 74 72 69 63 |Sun's ge|ocentric|
|00004a50| 20 65 63 6c 69 70 74 69 | 63 0a 58 20 20 20 09 09 | eclipti|c.X ..|
|00004a60| 09 09 20 20 20 20 20 20 | 20 6c 6f 6e 67 69 74 75 |.. | longitu|
|00004a70| 64 65 20 2a 2f 0a 58 20 | 20 20 2f 2a 20 4f 72 62 |de */.X | /* Orb|
|00004a80| 69 74 61 6c 20 64 69 73 | 74 61 6e 63 65 20 66 61 |ital dis|tance fa|
|00004a90| 63 74 6f 72 20 2a 2f 0a | 58 20 20 20 46 20 3d 20 |ctor */.|X F = |
|00004aa0| 28 28 31 20 2b 20 65 63 | 63 65 6e 74 20 2a 20 63 |((1 + ec|cent * c|
|00004ab0| 6f 73 28 74 6f 72 61 64 | 28 45 63 29 29 29 20 2f |os(torad|(Ec))) /|
|00004ac0| 20 28 31 20 2d 20 65 63 | 63 65 6e 74 20 2a 20 65 | (1 - ec|cent * e|
|00004ad0| 63 63 65 6e 74 29 29 3b | 0a 58 20 20 20 53 75 6e |ccent));|.X Sun|
|00004ae0| 44 69 73 74 20 3d 20 73 | 75 6e 73 6d 61 78 20 2f |Dist = s|unsmax /|
|00004af0| 20 46 3b 09 20 20 20 20 | 2f 2a 20 44 69 73 74 61 | F;. |/* Dista|
|00004b00| 6e 63 65 20 74 6f 20 53 | 75 6e 20 69 6e 20 6b 6d |nce to S|un in km|
|00004b10| 20 2a 2f 0a 58 20 20 20 | 20 20 20 20 20 53 75 6e | */.X | Sun|
|00004b20| 41 6e 67 20 3d 20 46 20 | 2a 20 73 75 6e 61 6e 67 |Ang = F |* sunang|
|00004b30| 73 69 7a 3b 20 20 20 20 | 20 2f 2a 20 53 75 6e 27 |siz; | /* Sun'|
|00004b40| 73 20 61 6e 67 75 6c 61 | 72 20 73 69 7a 65 20 69 |s angula|r size i|
|00004b50| 6e 20 64 65 67 72 65 65 | 73 20 2a 2f 0a 58 0a 58 |n degree|s */.X.X|
|00004b60| 0a 58 20 20 20 20 20 20 | 20 20 2f 2a 20 43 61 6c |.X | /* Cal|
|00004b70| 63 75 6c 61 74 69 6f 6e | 20 6f 66 20 74 68 65 20 |culation| of the |
|00004b80| 4d 6f 6f 6e 27 73 20 70 | 6f 73 69 74 69 6f 6e 20 |Moon's p|osition |
|00004b90| 2a 2f 0a 58 0a 58 20 20 | 20 20 20 20 20 20 2f 2a |*/.X.X | /*|
|00004ba0| 20 4d 6f 6f 6e 27 73 20 | 6d 65 61 6e 20 6c 6f 6e | Moon's |mean lon|
|00004bb0| 67 69 74 75 64 65 20 2a | 2f 0a 58 20 20 20 6d 6c |gitude *|/.X ml|
|00004bc0| 20 3d 20 66 69 78 61 6e | 67 6c 65 28 31 33 2e 31 | = fixan|gle(13.1|
|00004bd0| 37 36 33 39 36 36 20 2a | 20 44 61 79 20 2b 20 6d |763966 *| Day + m|
|00004be0| 6d 6c 6f 6e 67 29 3b 0a | 58 0a 58 20 20 20 20 20 |mlong);.|X.X |
|00004bf0| 20 20 20 2f 2a 20 4d 6f | 6f 6e 27 73 20 6d 65 61 | /* Mo|on's mea|
|00004c00| 6e 20 61 6e 6f 6d 61 6c | 79 20 2a 2f 0a 58 20 20 |n anomal|y */.X |
|00004c10| 20 4d 4d 20 3d 20 66 69 | 78 61 6e 67 6c 65 28 6d | MM = fi|xangle(m|
|00004c20| 6c 20 2d 20 30 2e 31 31 | 31 34 30 34 31 20 2a 20 |l - 0.11|14041 * |
|00004c30| 44 61 79 20 2d 20 6d 6d | 6c 6f 6e 67 70 29 3b 0a |Day - mm|longp);.|
|00004c40| 58 0a 58 20 20 20 20 20 | 20 20 20 2f 2a 20 4d 6f |X.X | /* Mo|
|00004c50| 6f 6e 27 73 20 61 73 63 | 65 6e 64 69 6e 67 20 6e |on's asc|ending n|
|00004c60| 6f 64 65 20 6d 65 61 6e | 20 6c 6f 6e 67 69 74 75 |ode mean| longitu|
|00004c70| 64 65 20 2a 2f 0a 58 20 | 20 20 4d 4e 20 3d 20 66 |de */.X | MN = f|
|00004c80| 69 78 61 6e 67 6c 65 28 | 6d 6c 6e 6f 64 65 20 2d |ixangle(|mlnode -|
|00004c90| 20 30 2e 30 35 32 39 35 | 33 39 20 2a 20 44 61 79 | 0.05295|39 * Day|
|00004ca0| 29 3b 0a 58 0a 58 20 20 | 20 2f 2a 20 45 76 65 63 |);.X.X | /* Evec|
|00004cb0| 74 69 6f 6e 20 2a 2f 0a | 58 20 20 20 45 76 20 3d |tion */.|X Ev =|
|00004cc0| 20 31 2e 32 37 33 39 20 | 2a 20 73 69 6e 28 74 6f | 1.2739 |* sin(to|
|00004cd0| 72 61 64 28 32 20 2a 20 | 28 6d 6c 20 2d 20 4c 61 |rad(2 * |(ml - La|
|00004ce0| 6d 62 64 61 73 75 6e 29 | 20 2d 20 4d 4d 29 29 3b |mbdasun)| - MM));|
|00004cf0| 0a 58 0a 58 20 20 20 2f | 2a 20 41 6e 6e 75 61 6c |.X.X /|* Annual|
|00004d00| 20 65 71 75 61 74 69 6f | 6e 20 2a 2f 0a 58 20 20 | equatio|n */.X |
|00004d10| 20 41 65 20 3d 20 30 2e | 31 38 35 38 20 2a 20 73 | Ae = 0.|1858 * s|
|00004d20| 69 6e 28 74 6f 72 61 64 | 28 4d 29 29 3b 0a 58 0a |in(torad|(M));.X.|
|00004d30| 58 20 20 20 2f 2a 20 43 | 6f 72 72 65 63 74 69 6f |X /* C|orrectio|
|00004d40| 6e 20 74 65 72 6d 20 2a | 2f 0a 58 20 20 20 41 33 |n term *|/.X A3|
|00004d50| 20 3d 20 30 2e 33 37 20 | 2a 20 73 69 6e 28 74 6f | = 0.37 |* sin(to|
|00004d60| 72 61 64 28 4d 29 29 3b | 0a 58 0a 58 20 20 20 2f |rad(M));|.X.X /|
|00004d70| 2a 20 43 6f 72 72 65 63 | 74 65 64 20 61 6e 6f 6d |* Correc|ted anom|
|00004d80| 61 6c 79 20 2a 2f 0a 58 | 20 20 20 4d 6d 50 20 3d |aly */.X| MmP =|
|00004d90| 20 4d 4d 20 2b 20 45 76 | 20 2d 20 41 65 20 2d 20 | MM + Ev| - Ae - |
|00004da0| 41 33 3b 0a 58 0a 58 20 | 20 20 2f 2a 20 43 6f 72 |A3;.X.X | /* Cor|
|00004db0| 72 65 63 74 69 6f 6e 20 | 66 6f 72 20 74 68 65 20 |rection |for the |
|00004dc0| 65 71 75 61 74 69 6f 6e | 20 6f 66 20 74 68 65 20 |equation| of the |
|00004dd0| 63 65 6e 74 72 65 20 2a | 2f 0a 58 20 20 20 6d 45 |centre *|/.X mE|
|00004de0| 63 20 3d 20 36 2e 32 38 | 38 36 20 2a 20 73 69 6e |c = 6.28|86 * sin|
|00004df0| 28 74 6f 72 61 64 28 4d | 6d 50 29 29 3b 0a 58 0a |(torad(M|mP));.X.|
|00004e00| 58 20 20 20 2f 2a 20 41 | 6e 6f 74 68 65 72 20 63 |X /* A|nother c|
|00004e10| 6f 72 72 65 63 74 69 6f | 6e 20 74 65 72 6d 20 2a |orrectio|n term *|
|00004e20| 2f 0a 58 20 20 20 41 34 | 20 3d 20 30 2e 32 31 34 |/.X A4| = 0.214|
|00004e30| 20 2a 20 73 69 6e 28 74 | 6f 72 61 64 28 32 20 2a | * sin(t|orad(2 *|
|00004e40| 20 4d 6d 50 29 29 3b 0a | 58 0a 58 20 20 20 2f 2a | MmP));.|X.X /*|
|00004e50| 20 43 6f 72 72 65 63 74 | 65 64 20 6c 6f 6e 67 69 | Correct|ed longi|
|00004e60| 74 75 64 65 20 2a 2f 0a | 58 20 20 20 6c 50 20 3d |tude */.|X lP =|
|00004e70| 20 6d 6c 20 2b 20 45 76 | 20 2b 20 6d 45 63 20 2d | ml + Ev| + mEc -|
|00004e80| 20 41 65 20 2b 20 41 34 | 3b 0a 58 0a 58 20 20 20 | Ae + A4|;.X.X |
|00004e90| 2f 2a 20 56 61 72 69 61 | 74 69 6f 6e 20 2a 2f 0a |/* Varia|tion */.|
|00004ea0| 58 20 20 20 56 20 3d 20 | 30 2e 36 35 38 33 20 2a |X V = |0.6583 *|
|00004eb0| 20 73 69 6e 28 74 6f 72 | 61 64 28 32 20 2a 20 28 | sin(tor|ad(2 * (|
|00004ec0| 6c 50 20 2d 20 4c 61 6d | 62 64 61 73 75 6e 29 29 |lP - Lam|bdasun))|
|00004ed0| 29 3b 0a 58 0a 58 20 20 | 20 2f 2a 20 31 20 6c 6f |);.X.X | /* 1 lo|
|00004ee0| 6e 67 69 74 75 64 65 20 | 2a 2f 0a 58 20 20 20 6c |ngitude |*/.X l|
|00004ef0| 50 50 20 3d 20 6c 50 20 | 2b 20 56 3b 0a 58 0a 58 |PP = lP |+ V;.X.X|
|00004f00| 20 20 20 2f 2a 20 43 6f | 72 72 65 63 74 65 64 20 | /* Co|rrected |
|00004f10| 6c 6f 6e 67 69 74 75 64 | 65 20 6f 66 20 74 68 65 |longitud|e of the|
|00004f20| 20 6e 6f 64 65 20 2a 2f | 0a 58 20 20 20 4e 50 20 | node */|.X NP |
|00004f30| 3d 20 4d 4e 20 2d 20 30 | 2e 31 36 20 2a 20 73 69 |= MN - 0|.16 * si|
|00004f40| 6e 28 74 6f 72 61 64 28 | 4d 29 29 3b 0a 58 0a 58 |n(torad(|M));.X.X|
|00004f50| 20 20 20 2f 2a 20 59 20 | 69 6e 63 6c 69 6e 61 74 | /* Y |inclinat|
|00004f60| 69 6f 6e 20 63 6f 6f 72 | 64 69 6e 61 74 65 20 2a |ion coor|dinate *|
|00004f70| 2f 0a 58 20 20 20 79 20 | 3d 20 73 69 6e 28 74 6f |/.X y |= sin(to|
|00004f80| 72 61 64 28 6c 50 50 20 | 2d 20 4e 50 29 29 20 2a |rad(lPP |- NP)) *|
|00004f90| 20 63 6f 73 28 74 6f 72 | 61 64 28 6d 69 6e 63 29 | cos(tor|ad(minc)|
|00004fa0| 29 3b 0a 58 0a 58 20 20 | 20 2f 2a 20 58 20 69 6e |);.X.X | /* X in|
|00004fb0| 63 6c 69 6e 61 74 69 6f | 6e 20 63 6f 6f 72 64 69 |clinatio|n coordi|
|00004fc0| 6e 61 74 65 20 2a 2f 0a | 58 20 20 20 78 20 3d 20 |nate */.|X x = |
|00004fd0| 63 6f 73 28 74 6f 72 61 | 64 28 6c 50 50 20 2d 20 |cos(tora|d(lPP - |
|00004fe0| 4e 50 29 29 3b 0a 58 0a | 58 20 20 20 2f 2a 20 45 |NP));.X.|X /* E|
|00004ff0| 63 6c 69 70 74 69 63 20 | 6c 6f 6e 67 69 74 75 64 |cliptic |longitud|
|00005000| 65 20 2a 2f 0a 58 20 20 | 20 4c 61 6d 62 64 61 6d |e */.X | Lambdam|
|00005010| 6f 6f 6e 20 3d 20 74 6f | 64 65 67 28 61 74 61 6e |oon = to|deg(atan|
|00005020| 32 28 79 2c 20 78 29 29 | 3b 0a 58 20 20 20 4c 61 |2(y, x))|;.X La|
|00005030| 6d 62 64 61 6d 6f 6f 6e | 20 2b 3d 20 4e 50 3b 0a |mbdamoon| += NP;.|
|00005040| 58 0a 58 20 20 20 2f 2a | 20 43 61 6c 63 75 6c 61 |X.X /*| Calcula|
|00005050| 74 69 6f 6e 20 6f 66 20 | 74 68 65 20 70 68 61 73 |tion of |the phas|
|00005060| 65 20 6f 66 20 74 68 65 | 20 4d 6f 6f 6e 20 2a 2f |e of the| Moon */|
|00005070| 0a 58 0a 58 20 20 20 2f | 2a 20 41 67 65 20 6f 66 |.X.X /|* Age of|
|00005080| 20 74 68 65 20 4d 6f 6f | 6e 20 69 6e 20 64 65 67 | the Moo|n in deg|
|00005090| 72 65 65 73 20 2a 2f 0a | 58 20 20 20 4d 6f 6f 6e |rees */.|X Moon|
|000050a0| 41 67 65 20 3d 20 6c 50 | 50 20 2d 20 4c 61 6d 62 |Age = lP|P - Lamb|
|000050b0| 64 61 73 75 6e 3b 0a 58 | 0a 58 20 20 20 2f 2a 20 |dasun;.X|.X /* |
|000050c0| 50 68 61 73 65 20 6f 66 | 20 74 68 65 20 4d 6f 6f |Phase of| the Moo|
|000050d0| 6e 20 2a 2f 0a 58 20 20 | 20 4d 6f 6f 6e 50 68 61 |n */.X | MoonPha|
|000050e0| 73 65 20 3d 20 28 31 20 | 2d 20 63 6f 73 28 74 6f |se = (1 |- cos(to|
|000050f0| 72 61 64 28 4d 6f 6f 6e | 41 67 65 29 29 29 20 2f |rad(Moon|Age))) /|
|00005100| 20 32 3b 0a 58 0a 58 20 | 20 20 2f 2a 20 43 61 6c | 2;.X.X | /* Cal|
|00005110| 63 75 6c 61 74 65 20 64 | 69 73 74 61 6e 63 65 20 |culate d|istance |
|00005120| 6f 66 20 6d 6f 6f 6e 20 | 66 72 6f 6d 20 74 68 65 |of moon |from the|
|00005130| 20 63 65 6e 74 72 65 20 | 6f 66 20 74 68 65 20 45 | centre |of the E|
|00005140| 61 72 74 68 20 2a 2f 0a | 58 0a 58 20 20 20 4d 6f |arth */.|X.X Mo|
|00005150| 6f 6e 44 69 73 74 20 3d | 20 28 6d 73 6d 61 78 20 |onDist =| (msmax |
|00005160| 2a 20 28 31 20 2d 20 6d | 65 63 63 20 2a 20 6d 65 |* (1 - m|ecc * me|
|00005170| 63 63 29 29 20 2f 0a 58 | 20 20 20 20 20 20 28 31 |cc)) /.X| (1|
|00005180| 20 2b 20 6d 65 63 63 20 | 2a 20 63 6f 73 28 74 6f | + mecc |* cos(to|
|00005190| 72 61 64 28 4d 6d 50 20 | 2b 20 6d 45 63 29 29 29 |rad(MmP |+ mEc)))|
|000051a0| 3b 0a 58 0a 58 20 20 20 | 20 20 20 20 20 2f 2a 20 |;.X.X | /* |
|000051b0| 43 61 6c 63 75 6c 61 74 | 65 20 4d 6f 6f 6e 27 73 |Calculat|e Moon's|
|000051c0| 20 61 6e 67 75 6c 61 72 | 20 64 69 61 6d 65 74 65 | angular| diamete|
|000051d0| 72 20 2a 2f 0a 58 0a 58 | 20 20 20 4d 6f 6f 6e 44 |r */.X.X| MoonD|
|000051e0| 46 72 61 63 20 3d 20 4d | 6f 6f 6e 44 69 73 74 20 |Frac = M|oonDist |
|000051f0| 2f 20 6d 73 6d 61 78 3b | 0a 58 20 20 20 4d 6f 6f |/ msmax;|.X Moo|
|00005200| 6e 41 6e 67 20 3d 20 6d | 61 6e 67 73 69 7a 20 2f |nAng = m|angsiz /|
|00005210| 20 4d 6f 6f 6e 44 46 72 | 61 63 3b 0a 58 0a 58 20 | MoonDFr|ac;.X.X |
|00005220| 20 20 69 66 28 70 70 68 | 61 73 65 29 20 20 20 2a | if(pph|ase) *|
|00005230| 70 70 68 61 73 65 20 3d | 20 4d 6f 6f 6e 50 68 61 |pphase =| MoonPha|
|00005240| 73 65 3b 0a 58 20 20 20 | 69 66 28 6d 61 67 65 29 |se;.X |if(mage)|
|00005250| 20 20 20 20 20 2a 6d 61 | 67 65 20 3d 20 73 79 6e | *ma|ge = syn|
|00005260| 6d 6f 6e 74 68 20 2a 20 | 28 66 69 78 61 6e 67 6c |month * |(fixangl|
|00005270| 65 28 4d 6f 6f 6e 41 67 | 65 29 20 2f 20 33 36 30 |e(MoonAg|e) / 360|
|00005280| 2e 30 29 3b 0a 58 20 20 | 20 69 66 28 64 69 73 74 |.0);.X | if(dist|
|00005290| 29 20 20 20 20 20 2a 64 | 69 73 74 20 3d 20 4d 6f |) *d|ist = Mo|
|000052a0| 6f 6e 44 69 73 74 3b 0a | 58 20 20 20 69 66 28 61 |onDist;.|X if(a|
|000052b0| 6e 67 64 69 61 29 20 20 | 20 2a 61 6e 67 64 69 61 |ngdia) | *angdia|
|000052c0| 20 3d 20 4d 6f 6f 6e 41 | 6e 67 3b 0a 58 20 20 20 | = MoonA|ng;.X |
|000052d0| 69 66 28 73 75 64 69 73 | 74 29 20 20 20 2a 73 75 |if(sudis|t) *su|
|000052e0| 64 69 73 74 20 3d 20 53 | 75 6e 44 69 73 74 3b 0a |dist = S|unDist;.|
|000052f0| 58 20 20 20 69 66 28 73 | 75 61 6e 67 64 69 61 29 |X if(s|uangdia)|
|00005300| 20 2a 73 75 61 6e 67 64 | 69 61 20 3d 20 53 75 6e | *suangd|ia = Sun|
|00005310| 41 6e 67 3b 0a 58 20 20 | 20 72 65 74 75 72 6e 20 |Ang;.X | return |
|00005320| 66 69 78 61 6e 67 6c 65 | 28 4d 6f 6f 6e 41 67 65 |fixangle|(MoonAge|
|00005330| 29 20 2f 20 33 36 30 2e | 30 3b 0a 58 7d 0a 58 0a |) / 360.|0;.X}.X.|
|00005340| 58 2f 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |X/******|********|
|00005350| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00005360| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00005370| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00005380| 2a 2f 0a 58 2f 2a 20 20 | 20 20 20 20 20 20 20 20 |*/.X/* | |
|00005390| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000053a0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000053b0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000053c0| 20 20 20 2a 2f 0a 58 2f | 2a 20 20 4d 6f 6f 6e 50 | */.X/|* MoonP|
|000053d0| 68 61 73 65 20 20 20 20 | 20 20 20 20 20 20 20 20 |hase | |
|000053e0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000053f0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00005400| 20 20 20 20 20 20 2a 2f | 0a 58 2f 2a 20 20 20 20 | */|.X/* |
|00005410| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00005420| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00005430| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00005440| 20 20 20 20 20 20 20 20 | 20 2a 2f 0a 58 2f 2a 20 | | */.X/* |
|00005450| 20 49 6e 74 65 72 66 61 | 63 65 20 72 6f 75 74 69 | Interfa|ce routi|
|00005460| 6e 65 20 64 65 61 6c 69 | 6e 67 20 69 6e 20 52 65 |ne deali|ng in Re|
|00005470| 6d 69 6e 64 20 72 65 70 | 72 65 73 65 6e 74 61 74 |mind rep|resentat|
|00005480| 69 6f 6e 73 2e 20 20 20 | 20 20 20 20 2a 2f 0a 58 |ions. | */.X|
|00005490| 2f 2a 20 20 47 69 76 65 | 6e 20 61 20 6c 6f 63 61 |/* Give|n a loca|
|000054a0| 6c 20 64 61 74 65 20 61 | 6e 64 20 74 69 6d 65 2c |l date a|nd time,|
|000054b0| 20 72 65 74 75 72 6e 73 | 20 74 68 65 20 6d 6f 6f | returns| the moo|
|000054c0| 6e 20 70 68 61 73 65 20 | 61 74 20 20 20 20 20 2a |n phase |at *|
|000054d0| 2f 0a 58 2f 2a 20 20 74 | 68 61 74 20 64 61 74 65 |/.X/* t|hat date|
|000054e0| 20 61 6e 64 20 74 69 6d | 65 20 61 73 20 61 20 6e | and tim|e as a n|
|000054f0| 75 6d 62 65 72 20 66 72 | 6f 6d 20 30 20 74 6f 20 |umber fr|om 0 to |
|00005500| 33 36 30 2e 20 20 20 20 | 20 20 20 20 20 20 20 20 |360. | |
|00005510| 20 20 2a 2f 0a 58 2f 2a | 20 20 20 20 20 20 20 20 | */.X/*| |
|00005520| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00005530| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00005540| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00005550| 20 20 20 20 20 2a 2f 0a | 58 2f 2a 2a 2a 2a 2a 2a | */.|X/******|
|00005560| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00005570| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00005580| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00005590| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2f 0a 58 23 69 66 64 |********|*/.X#ifd|
|000055a0| 65 66 20 48 41 56 45 5f | 50 52 4f 54 4f 53 0a 58 |ef HAVE_|PROTOS.X|
|000055b0| 50 55 42 4c 49 43 20 69 | 6e 74 20 4d 6f 6f 6e 50 |PUBLIC i|nt MoonP|
|000055c0| 68 61 73 65 28 69 6e 74 | 20 64 61 74 65 2c 20 69 |hase(int| date, i|
|000055d0| 6e 74 20 74 69 6d 65 29 | 0a 58 23 65 6c 73 65 0a |nt time)|.X#else.|
|000055e0| 58 69 6e 74 20 4d 6f 6f | 6e 50 68 61 73 65 28 64 |Xint Moo|nPhase(d|
|000055f0| 61 74 65 2c 20 74 69 6d | 65 29 0a 58 69 6e 74 20 |ate, tim|e).Xint |
|00005600| 64 61 74 65 2c 20 74 69 | 6d 65 3b 0a 58 23 65 6e |date, ti|me;.X#en|
|00005610| 64 69 66 0a 58 7b 0a 58 | 20 20 20 69 6e 74 20 75 |dif.X{.X| int u|
|00005620| 74 63 64 2c 20 75 74 63 | 74 3b 0a 58 20 20 20 69 |tcd, utc|t;.X i|
|00005630| 6e 74 20 79 2c 20 6d 2c | 20 64 3b 0a 58 20 20 20 |nt y, m,| d;.X |
|00005640| 64 6f 75 62 6c 65 20 6a | 64 2c 20 6d 70 3b 0a 58 |double j|d, mp;.X|
|00005650| 0a 58 20 20 20 2f 2a 20 | 43 6f 6e 76 65 72 74 20 |.X /* |Convert |
|00005660| 66 72 6f 6d 20 6c 6f 63 | 61 6c 20 74 6f 20 55 54 |from loc|al to UT|
|00005670| 43 20 2a 2f 0a 58 20 20 | 20 4c 6f 63 61 6c 54 6f |C */.X | LocalTo|
|00005680| 55 54 43 28 64 61 74 65 | 2c 20 74 69 6d 65 2c 20 |UTC(date|, time, |
|00005690| 26 75 74 63 64 2c 20 26 | 75 74 63 74 29 3b 0a 58 |&utcd, &|utct);.X|
|000056a0| 0a 58 20 20 20 2f 2a 20 | 43 6f 6e 76 65 72 74 20 |.X /* |Convert |
|000056b0| 66 72 6f 6d 20 52 65 6d | 69 6e 64 20 72 65 70 72 |from Rem|ind repr|
|000056c0| 65 73 65 6e 74 61 74 69 | 6f 6e 20 74 6f 20 79 65 |esentati|on to ye|
|000056d0| 61 72 2f 6d 6f 6e 2f 64 | 61 79 20 2a 2f 0a 58 20 |ar/mon/d|ay */.X |
|000056e0| 20 20 46 72 6f 6d 4a 75 | 6c 69 61 6e 28 75 74 63 | FromJu|lian(utc|
|000056f0| 64 2c 20 26 79 2c 20 26 | 6d 2c 20 26 64 29 3b 0a |d, &y, &|m, &d);.|
|00005700| 58 0a 58 20 20 20 2f 2a | 20 43 6f 6e 76 65 72 74 |X.X /*| Convert|
|00005710| 20 74 6f 20 61 20 74 72 | 75 65 20 4a 75 6c 69 61 | to a tr|ue Julia|
|00005720| 6e 20 64 61 74 65 20 2d | 2d 20 73 6f 72 72 79 20 |n date -|- sorry |
|00005730| 66 6f 72 20 74 68 65 20 | 6e 61 6d 65 20 63 6c 61 |for the |name cla|
|00005740| 73 68 65 73 21 20 2a 2f | 0a 58 20 20 20 6a 64 20 |shes! */|.X jd |
|00005750| 3d 20 6a 74 69 6d 65 28 | 79 2c 20 6d 2c 20 64 2c |= jtime(|y, m, d,|
|00005760| 20 28 75 74 63 74 20 2f | 20 36 30 29 2c 20 28 75 | (utct /| 60), (u|
|00005770| 74 63 74 20 25 20 36 30 | 29 2c 20 30 29 3b 20 20 |tct % 60|), 0); |
|00005780| 20 0a 58 0a 58 20 20 20 | 2f 2a 20 43 61 6c 63 75 | .X.X |/* Calcu|
|00005790| 6c 61 74 65 20 6d 6f 6f | 6e 20 70 68 61 73 65 20 |late moo|n phase |
|000057a0| 2a 2f 0a 58 20 20 20 6d | 70 20 3d 20 33 36 30 2e |*/.X m|p = 360.|
|000057b0| 30 20 2a 20 70 68 61 73 | 65 28 6a 64 2c 20 4e 55 |0 * phas|e(jd, NU|
|000057c0| 4c 4c 2c 20 4e 55 4c 4c | 2c 20 4e 55 4c 4c 2c 20 |LL, NULL|, NULL, |
|000057d0| 4e 55 4c 4c 2c 20 4e 55 | 4c 4c 2c 20 4e 55 4c 4c |NULL, NU|LL, NULL|
|000057e0| 29 3b 0a 58 20 20 20 72 | 65 74 75 72 6e 20 28 69 |);.X r|eturn (i|
|000057f0| 6e 74 29 20 6d 70 3b 0a | 58 7d 0a 58 0a 58 2f 2a |nt) mp;.|X}.X.X/*|
|00005800| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00005810| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00005820| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00005830| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2f 0a |********|******/.|
|00005840| 58 2f 2a 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |X/* | |
|00005850| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00005860| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00005870| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00005880| 2a 2f 0a 58 2f 2a 20 20 | 48 75 6e 74 50 68 61 73 |*/.X/* |HuntPhas|
|00005890| 65 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |e | |
|000058a0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000058b0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000058c0| 20 20 20 2a 2f 0a 58 2f | 2a 20 20 20 20 20 20 20 | */.X/|* |
|000058d0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000058e0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000058f0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00005900| 20 20 20 20 20 20 2a 2f | 0a 58 2f 2a 20 20 47 69 | */|.X/* Gi|
|00005910| 76 65 6e 20 61 20 73 74 | 61 72 74 69 6e 67 20 64 |ven a st|arting d|
|00005920| 61 74 65 20 61 6e 64 20 | 74 69 6d 65 20 61 6e 64 |ate and |time and|
|00005930| 20 61 20 74 61 72 67 65 | 74 20 70 68 61 73 65 2c | a targe|t phase,|
|00005940| 20 66 69 6e 64 20 20 20 | 20 2a 2f 0a 58 2f 2a 20 | find | */.X/* |
|00005950| 20 74 68 65 20 66 69 72 | 73 74 20 64 61 74 65 20 | the fir|st date |
|00005960| 6f 6e 20 6f 72 20 61 66 | 74 65 72 20 74 68 65 20 |on or af|ter the |
|00005970| 73 74 61 72 74 69 6e 67 | 20 64 61 74 65 20 61 6e |starting| date an|
|00005980| 64 20 74 69 6d 65 20 77 | 68 65 6e 20 2a 2f 0a 58 |d time w|hen */.X|
|00005990| 2f 2a 20 20 74 68 65 20 | 6d 6f 6f 6e 20 68 69 74 |/* the |moon hit|
|000059a0| 73 20 74 68 65 20 73 70 | 65 63 69 66 69 65 64 20 |s the sp|ecified |
|000059b0| 70 68 61 73 65 2e 20 20 | 50 68 61 73 65 20 6d 75 |phase. |Phase mu|
|000059c0| 73 74 20 62 65 20 66 72 | 6f 6d 20 20 20 20 20 2a |st be fr|om *|
|000059d0| 2f 0a 58 2f 2a 20 20 30 | 20 74 6f 20 33 20 66 6f |/.X/* 0| to 3 fo|
|000059e0| 72 20 6e 65 77 2c 20 31 | 73 74 71 2c 20 66 75 6c |r new, 1|stq, ful|
|000059f0| 6c 2c 20 33 72 64 71 20 | 20 20 20 20 20 20 20 20 |l, 3rdq | |
|00005a00| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00005a10| 20 20 2a 2f 0a 58 2f 2a | 20 20 20 20 20 20 20 20 | */.X/*| |
|00005a20| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00005a30| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00005a40| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00005a50| 20 20 20 20 20 2a 2f 0a | 58 2f 2a 2a 2a 2a 2a 2a | */.|X/******|
|00005a60| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00005a70| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00005a80| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2a 2a 2a 2a 2a 2a 2a |********|********|
|00005a90| 2a 2a 2a 2a 2a 2a 2a 2a | 2a 2f 0a 58 23 69 66 64 |********|*/.X#ifd|
|00005aa0| 65 66 20 48 41 56 45 5f | 50 52 4f 54 4f 53 0a 58 |ef HAVE_|PROTOS.X|
|00005ab0| 50 55 42 4c 49 43 20 76 | 6f 69 64 20 48 75 6e 74 |PUBLIC v|oid Hunt|
|00005ac0| 50 68 61 73 65 28 69 6e | 74 20 73 74 61 72 74 64 |Phase(in|t startd|
|00005ad0| 61 74 65 2c 20 69 6e 74 | 20 73 74 61 72 74 74 69 |ate, int| startti|
|00005ae0| 6d 2c 20 69 6e 74 20 70 | 68 61 73 2c 20 69 6e 74 |m, int p|has, int|
|00005af0| 20 2a 64 61 74 65 2c 20 | 69 6e 74 20 2a 74 69 6d | *date, |int *tim|
|00005b00| 65 29 0a 58 23 65 6c 73 | 65 0a 58 76 6f 69 64 20 |e).X#els|e.Xvoid |
|00005b10| 48 75 6e 74 50 68 61 73 | 65 28 73 74 61 72 74 64 |HuntPhas|e(startd|
|00005b20| 61 74 65 2c 20 73 74 61 | 72 74 74 69 6d 2c 20 70 |ate, sta|rttim, p|
|00005b30| 68 61 73 2c 20 64 61 74 | 65 2c 20 74 69 6d 65 29 |has, dat|e, time)|
|00005b40| 0a 58 69 6e 74 20 73 74 | 61 72 74 64 61 74 65 2c |.Xint st|artdate,|
|00005b50| 20 73 74 61 72 74 74 69 | 6d 2c 20 70 68 61 73 2c | startti|m, phas,|
|00005b60| 20 2a 64 61 74 65 2c 20 | 2a 74 69 6d 65 3b 0a 58 | *date, |*time;.X|
|00005b70| 23 65 6e 64 69 66 0a 58 | 7b 0a 58 20 20 20 69 6e |#endif.X|{.X in|
|00005b80| 74 20 75 74 63 64 2c 20 | 75 74 63 74 3b 0a 58 20 |t utcd, |utct;.X |
|00005b90| 20 20 69 6e 74 20 79 2c | 20 6d 2c 20 64 3b 0a 58 | int y,| m, d;.X|
|00005ba0| 20 20 20 69 6e 74 20 68 | 2c 20 6d 69 6e 2c 20 73 | int h|, min, s|
|00005bb0| 3b 0a 58 20 20 20 69 6e | 74 20 64 31 2c 20 74 31 |;.X in|t d1, t1|
|00005bc0| 3b 0a 58 20 20 20 64 6f | 75 62 6c 65 20 6b 31 2c |;.X do|uble k1,|
|00005bd0| 20 6b 32 2c 20 6a 64 2c | 20 6a 64 6f 72 69 67 3b | k2, jd,| jdorig;|
|00005be0| 0a 58 20 20 20 64 6f 75 | 62 6c 65 20 6e 74 31 2c |.X dou|ble nt1,|
|00005bf0| 20 6e 74 32 3b 0a 58 20 | 20 20 2f 2a 20 43 6f 6e | nt2;.X | /* Con|
|00005c00| 76 65 72 74 20 66 72 6f | 6d 20 6c 6f 63 61 6c 20 |vert fro|m local |
|00005c10| 74 6f 20 55 54 43 20 2a | 2f 0a 58 20 20 20 4c 6f |to UTC *|/.X Lo|
|00005c20| 63 61 6c 54 6f 55 54 43 | 28 73 74 61 72 74 64 61 |calToUTC|(startda|
|00005c30| 74 65 2c 20 73 74 61 72 | 74 74 69 6d 2c 20 26 75 |te, star|ttim, &u|
|00005c40| 74 63 64 2c 20 26 75 74 | 63 74 29 3b 0a 58 0a 58 |tcd, &ut|ct);.X.X|
|00005c50| 20 20 20 2f 2a 20 43 6f | 6e 76 65 72 74 20 66 72 | /* Co|nvert fr|
|00005c60| 6f 6d 20 52 65 6d 69 6e | 64 20 72 65 70 72 65 73 |om Remin|d repres|
|00005c70| 65 6e 74 61 74 69 6f 6e | 20 74 6f 20 79 65 61 72 |entation| to year|
|00005c80| 2f 6d 6f 6e 2f 64 61 79 | 20 2a 2f 0a 58 20 20 20 |/mon/day| */.X |
|00005c90| 46 72 6f 6d 4a 75 6c 69 | 61 6e 28 75 74 63 64 2c |FromJuli|an(utcd,|
|00005ca0| 20 26 79 2c 20 26 6d 2c | 20 26 64 29 3b 0a 58 0a | &y, &m,| &d);.X.|
|00005cb0| 58 20 20 20 2f 2a 20 43 | 6f 6e 76 65 72 74 20 74 |X /* C|onvert t|
|00005cc0| 6f 20 61 20 74 72 75 65 | 20 4a 75 6c 69 61 6e 20 |o a true| Julian |
|00005cd0| 64 61 74 65 20 2d 2d 20 | 73 6f 72 72 79 20 66 6f |date -- |sorry fo|
|00005ce0| 72 20 74 68 65 20 6e 61 | 6d 65 20 63 6c 61 73 68 |r the na|me clash|
|00005cf0| 65 73 21 20 2a 2f 0a 58 | 20 20 20 6a 64 6f 72 69 |es! */.X| jdori|
|00005d00| 67 20 3d 20 6a 74 69 6d | 65 28 79 2c 20 6d 2c 20 |g = jtim|e(y, m, |
|00005d10| 64 2c 20 28 75 74 63 74 | 20 2f 20 36 30 29 2c 20 |d, (utct| / 60), |
|00005d20| 28 75 74 63 74 20 25 20 | 36 30 29 2c 20 30 29 3b |(utct % |60), 0);|
|00005d30| 20 20 20 0a 58 20 20 20 | 6a 64 20 3d 20 6a 64 6f | .X |jd = jdo|
|00005d40| 72 69 67 20 2d 20 34 35 | 3b 0a 58 20 20 20 6e 74 |rig - 45|;.X nt|
|00005d50| 31 20 3d 20 6d 65 61 6e | 70 68 61 73 65 28 6a 64 |1 = mean|phase(jd|
|00005d60| 2c 20 30 2e 30 2c 20 26 | 6b 31 29 3b 0a 58 20 20 |, 0.0, &|k1);.X |
|00005d70| 20 77 68 69 6c 65 28 31 | 29 20 7b 0a 58 20 20 20 | while(1|) {.X |
|00005d80| 20 20 20 6a 64 20 2b 3d | 20 73 79 6e 6d 6f 6e 74 | jd +=| synmont|
|00005d90| 68 3b 0a 58 20 20 20 20 | 20 20 6e 74 32 20 3d 20 |h;.X | nt2 = |
|00005da0| 6d 65 61 6e 70 68 61 73 | 65 28 6a 64 2c 20 30 2e |meanphas|e(jd, 0.|
|00005db0| 30 2c 20 26 6b 32 29 3b | 0a 58 20 20 20 20 20 20 |0, &k2);|.X |
|00005dc0| 69 66 20 28 6e 74 31 20 | 3c 3d 20 6a 64 6f 72 69 |if (nt1 |<= jdori|
|00005dd0| 67 20 26 26 20 6e 74 32 | 20 3e 20 6a 64 6f 72 69 |g && nt2| > jdori|
|00005de0| 67 29 20 62 72 65 61 6b | 3b 0a 58 20 20 20 20 20 |g) break|;.X |
|00005df0| 20 6e 74 31 20 3d 20 6e | 74 32 3b 0a 58 20 20 20 | nt1 = n|t2;.X |
|00005e00| 20 20 20 6b 31 20 3d 20 | 6b 32 3b 0a 58 20 20 20 | k1 = |k2;.X |
|00005e10| 7d 0a 58 20 20 20 6a 64 | 20 3d 20 74 72 75 65 70 |}.X jd| = truep|
|00005e20| 68 61 73 65 28 6b 31 2c | 20 70 68 61 73 2f 34 2e |hase(k1,| phas/4.|
|00005e30| 30 29 3b 0a 58 20 20 20 | 69 66 20 28 6a 64 20 3c |0);.X |if (jd <|
|00005e40| 20 6a 64 6f 72 69 67 29 | 20 6a 64 20 3d 20 74 72 | jdorig)| jd = tr|
|00005e50| 75 65 70 68 61 73 65 28 | 6b 32 2c 20 70 68 61 73 |uephase(|k2, phas|
|00005e60| 2f 34 2e 30 29 3b 0a 58 | 0a 58 20 20 20 2f 2a 20 |/4.0);.X|.X /* |
|00005e70| 43 6f 6e 76 65 72 74 20 | 62 61 63 6b 20 74 6f 20 |Convert |back to |
|00005e80| 52 65 6d 69 6e 64 20 66 | 6f 72 6d 61 74 20 2a 2f |Remind f|ormat */|
|00005e90| 0a 58 20 20 20 6a 79 65 | 61 72 28 6a 64 2c 20 26 |.X jye|ar(jd, &|
|00005ea0| 79 2c 20 26 6d 2c 20 26 | 64 29 3b 0a 58 20 20 20 |y, &m, &|d);.X |
|00005eb0| 6a 68 6d 73 28 6a 64 2c | 20 26 68 2c 20 26 6d 69 |jhms(jd,| &h, &mi|
|00005ec0| 6e 2c 20 26 73 29 3b 0a | 58 0a 58 20 20 20 64 31 |n, &s);.|X.X d1|
|00005ed0| 20 3d 20 4a 75 6c 69 61 | 6e 28 79 2c 20 6d 2c 20 | = Julia|n(y, m, |
|00005ee0| 64 29 3b 0a 58 20 20 20 | 74 31 20 3d 20 68 2a 36 |d);.X |t1 = h*6|
|00005ef0| 30 20 2b 20 6d 69 6e 3b | 0a 58 20 20 20 55 54 43 |0 + min;|.X UTC|
|00005f00| 54 6f 4c 6f 63 61 6c 28 | 64 31 2c 20 74 31 2c 20 |ToLocal(|d1, t1, |
|00005f10| 64 61 74 65 2c 20 74 69 | 6d 65 29 3b 0a 58 7d 0a |date, ti|me);.X}.|
|00005f20| 45 4e 44 5f 4f 46 5f 46 | 49 4c 45 0a 20 20 69 66 |END_OF_F|ILE. if|
|00005f30| 20 74 65 73 74 20 32 32 | 34 36 31 20 2d 6e 65 20 | test 22|461 -ne |
|00005f40| 60 77 63 20 2d 63 20 3c | 27 6d 6f 6f 6e 2e 63 27 |`wc -c <|'moon.c'|
|00005f50| 60 3b 20 74 68 65 6e 0a | 20 20 20 20 65 63 68 6f |`; then.| echo|
|00005f60| 20 73 68 61 72 3a 20 5c | 22 27 6d 6f 6f 6e 2e 63 | shar: \|"'moon.c|
|00005f70| 27 5c 22 20 75 6e 70 61 | 63 6b 65 64 20 77 69 74 |'\" unpa|cked wit|
|00005f80| 68 20 77 72 6f 6e 67 20 | 73 69 7a 65 21 0a 20 20 |h wrong |size!. |
|00005f90| 66 69 0a 20 20 23 20 65 | 6e 64 20 6f 66 20 27 6d |fi. # e|nd of 'm|
|00005fa0| 6f 6f 6e 2e 63 27 0a 66 | 69 0a 65 63 68 6f 20 73 |oon.c'.f|i.echo s|
|00005fb0| 68 61 72 3a 20 45 6e 64 | 20 6f 66 20 61 72 63 68 |har: End| of arch|
|00005fc0| 69 76 65 20 34 20 5c 28 | 6f 66 20 34 5c 29 2e 0a |ive 4 \(|of 4\)..|
|00005fd0| 63 70 20 2f 64 65 76 2f | 6e 75 6c 6c 20 61 72 6b |cp /dev/|null ark|
|00005fe0| 34 69 73 64 6f 6e 65 0a | 4d 49 53 53 49 4e 47 3d |4isdone.|MISSING=|
|00005ff0| 22 22 0a 66 6f 72 20 49 | 20 69 6e 20 31 20 32 20 |"".for I| in 1 2 |
|00006000| 33 20 34 20 3b 20 64 6f | 0a 20 20 20 20 69 66 20 |3 4 ; do|. if |
|00006010| 74 65 73 74 20 21 20 2d | 66 20 61 72 6b 24 7b 49 |test ! -|f ark${I|
|00006020| 7d 69 73 64 6f 6e 65 20 | 3b 20 74 68 65 6e 0a 09 |}isdone |; then..|
|00006030| 4d 49 53 53 49 4e 47 3d | 22 24 7b 4d 49 53 53 49 |MISSING=|"${MISSI|
|00006040| 4e 47 7d 20 24 7b 49 7d | 22 0a 20 20 20 20 66 69 |NG} ${I}|". fi|
|00006050| 0a 64 6f 6e 65 0a 69 66 | 20 74 65 73 74 20 22 24 |.done.if| test "$|
|00006060| 7b 4d 49 53 53 49 4e 47 | 7d 22 20 3d 20 22 22 20 |{MISSING|}" = "" |
|00006070| 3b 20 74 68 65 6e 0a 20 | 20 20 20 65 63 68 6f 20 |; then. | echo |
|00006080| 59 6f 75 20 68 61 76 65 | 20 75 6e 70 61 63 6b 65 |You have| unpacke|
|00006090| 64 20 61 6c 6c 20 34 20 | 61 72 63 68 69 76 65 73 |d all 4 |archives|
|000060a0| 2e 0a 20 20 20 20 72 6d | 20 2d 66 20 61 72 6b 5b |.. rm| -f ark[|
|000060b0| 31 2d 39 5d 69 73 64 6f | 6e 65 0a 65 6c 73 65 0a |1-9]isdo|ne.else.|
|000060c0| 20 20 20 20 65 63 68 6f | 20 59 6f 75 20 73 74 69 | echo| You sti|
|000060d0| 6c 6c 20 6d 75 73 74 20 | 75 6e 70 61 63 6b 20 74 |ll must |unpack t|
|000060e0| 68 65 20 66 6f 6c 6c 6f | 77 69 6e 67 20 61 72 63 |he follo|wing arc|
|000060f0| 68 69 76 65 73 3a 0a 20 | 20 20 20 65 63 68 6f 20 |hives:. | echo |
|00006100| 22 20 20 20 20 20 20 20 | 20 22 20 24 7b 4d 49 53 |" | " ${MIS|
|00006110| 53 49 4e 47 7d 0a 66 69 | 0a 65 78 69 74 20 30 0a |SING}.fi|.exit 0.|
|00006120| 65 78 69 74 20 30 20 23 | 20 4a 75 73 74 20 69 6e |exit 0 #| Just in|
|00006130| 20 63 61 73 65 2e 2e 2e | 0a | case...|. |
+--------+-------------------------+-------------------------+--------+--------+
