Mercurial > emacs
annotate lisp/calendar/lunar.el @ 13896:d23e1ac7b432
Doc fixes.
(ff-get-file): Use get-file-buffer to avoid false match.
(ff-switch-to-buffer): Rename arg FILE to BUFFER-OR-NAME.
author | Richard M. Stallman <rms@gnu.org> |
---|---|
date | Sun, 31 Dec 1995 00:10:18 +0000 |
parents | 9155a9ab5de9 |
children | 83f275dcd93a |
rev | line source |
---|---|
957 | 1 ;;; lunar.el --- calendar functions for phases of the moon. |
2 | |
13044
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
3 ;; Copyright (C) 1992, 1993, 1995 Free Software Foundation, Inc. |
957 | 4 |
5 ;; Author: Edward M. Reingold <reingold@cs.uiuc.edu> | |
2247
2c7997f249eb
Add or correct keywords
Eric S. Raymond <esr@snark.thyrsus.com>
parents:
957
diff
changeset
|
6 ;; Keywords: calendar |
2c7997f249eb
Add or correct keywords
Eric S. Raymond <esr@snark.thyrsus.com>
parents:
957
diff
changeset
|
7 ;; Human-Keywords: moon, lunar phases, calendar, diary |
957 | 8 |
9 ;; This file is part of GNU Emacs. | |
10 | |
6736
3e1323443b1a
Fix copying conditions for current GPL version.
Richard M. Stallman <rms@gnu.org>
parents:
6334
diff
changeset
|
11 ;; GNU Emacs is free software; you can redistribute it and/or modify |
3e1323443b1a
Fix copying conditions for current GPL version.
Richard M. Stallman <rms@gnu.org>
parents:
6334
diff
changeset
|
12 ;; it under the terms of the GNU General Public License as published by |
3e1323443b1a
Fix copying conditions for current GPL version.
Richard M. Stallman <rms@gnu.org>
parents:
6334
diff
changeset
|
13 ;; the Free Software Foundation; either version 2, or (at your option) |
3e1323443b1a
Fix copying conditions for current GPL version.
Richard M. Stallman <rms@gnu.org>
parents:
6334
diff
changeset
|
14 ;; any later version. |
957 | 15 |
6736
3e1323443b1a
Fix copying conditions for current GPL version.
Richard M. Stallman <rms@gnu.org>
parents:
6334
diff
changeset
|
16 ;; GNU Emacs is distributed in the hope that it will be useful, |
3e1323443b1a
Fix copying conditions for current GPL version.
Richard M. Stallman <rms@gnu.org>
parents:
6334
diff
changeset
|
17 ;; but WITHOUT ANY WARRANTY; without even the implied warranty of |
3e1323443b1a
Fix copying conditions for current GPL version.
Richard M. Stallman <rms@gnu.org>
parents:
6334
diff
changeset
|
18 ;; MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
3e1323443b1a
Fix copying conditions for current GPL version.
Richard M. Stallman <rms@gnu.org>
parents:
6334
diff
changeset
|
19 ;; GNU General Public License for more details. |
3e1323443b1a
Fix copying conditions for current GPL version.
Richard M. Stallman <rms@gnu.org>
parents:
6334
diff
changeset
|
20 |
3e1323443b1a
Fix copying conditions for current GPL version.
Richard M. Stallman <rms@gnu.org>
parents:
6334
diff
changeset
|
21 ;; You should have received a copy of the GNU General Public License |
3e1323443b1a
Fix copying conditions for current GPL version.
Richard M. Stallman <rms@gnu.org>
parents:
6334
diff
changeset
|
22 ;; along with GNU Emacs; see the file COPYING. If not, write to |
3e1323443b1a
Fix copying conditions for current GPL version.
Richard M. Stallman <rms@gnu.org>
parents:
6334
diff
changeset
|
23 ;; the Free Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA. |
957 | 24 |
25 ;;; Commentary: | |
26 | |
27 ;; This collection of functions implements lunar phases for calendar.el and | |
28 ;; diary.el. | |
29 | |
30 ;; Based on ``Astronomical Formulae for Calculators,'' 3rd ed., by Jean Meeus, | |
13044
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
31 ;; Willmann-Bell, Inc., 1985 and ``Astronomical Algorithms'' by Jean Meeus, |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
32 ;; Willmann-Bell, Inc., 1991. |
957 | 33 ;; |
34 ;; WARNING: The calculations will be accurate only to within a few minutes. | |
35 | |
36 ;; The author would be delighted to have an astronomically more sophisticated | |
37 ;; person rewrite the code for the lunar calculations in this file! | |
38 | |
39 ;; Comments, corrections, and improvements should be sent to | |
40 ;; Edward M. Reingold Department of Computer Science | |
41 ;; (217) 333-6733 University of Illinois at Urbana-Champaign | |
42 ;; reingold@cs.uiuc.edu 1304 West Springfield Avenue | |
43 ;; Urbana, Illinois 61801 | |
44 | |
45 ;;; Code: | |
46 | |
47 (if (fboundp 'atan) | |
48 (require 'lisp-float-type) | |
49 (error "Lunar calculations impossible since floating point is unavailable.")) | |
50 | |
51 (require 'solar) | |
52 | |
53 (defun lunar-phase-list (month year) | |
54 "List of lunar phases for three months starting with Gregorian MONTH, YEAR." | |
55 (let ((end-month month) | |
56 (end-year year) | |
57 (start-month month) | |
58 (start-year year)) | |
59 (increment-calendar-month end-month end-year 3) | |
60 (increment-calendar-month start-month start-year -1) | |
61 (let* ((end-date (list (list end-month 1 end-year))) | |
62 (start-date (list (list start-month | |
63 (calendar-last-day-of-month | |
64 start-month start-year) | |
65 start-year))) | |
66 (index (* 4 | |
67 (truncate | |
68 (* 12.3685 | |
69 (+ year | |
70 ( / (calendar-day-number (list month 1 year)) | |
71 366.0) | |
72 -1900))))) | |
73 (new-moon (lunar-phase index)) | |
74 (list)) | |
75 (while (calendar-date-compare new-moon end-date) | |
76 (if (calendar-date-compare start-date new-moon) | |
77 (setq list (append list (list new-moon)))) | |
78 (setq index (1+ index)) | |
79 (setq new-moon (lunar-phase index))) | |
80 list))) | |
81 | |
82 (defun lunar-phase (index) | |
83 "Local date and time of lunar phase INDEX. | |
84 Integer below INDEX/4 gives the lunation number, counting from Jan 1, 1900; | |
85 remainder mod 4 gives the phase: 0 new moon, 1 first quarter, 2 full moon, | |
86 3 last quarter." | |
7760
4edcac57a8f2
(calendar-mod): Remove; it was equivalent to `mod'. All callers changed.
Paul Eggert <eggert@twinsun.com>
parents:
6736
diff
changeset
|
87 (let* ((phase (mod index 4)) |
957 | 88 (index (/ index 4.0)) |
89 (time (/ index 1236.85)) | |
90 (date (+ (calendar-absolute-from-gregorian '(1 0.5 1900)) | |
91 0.75933 | |
92 (* 29.53058868 index) | |
93 (* 0.0001178 time time) | |
94 (* -0.000000155 time time time) | |
95 (* 0.00033 | |
96 (solar-sin-degrees (+ 166.56 | |
97 (* 132.87 time) | |
98 (* -0.009173 time time)))))) | |
4520 | 99 (sun-anomaly (mod |
957 | 100 (+ 359.2242 |
101 (* 29.105356 index) | |
102 (* -0.0000333 time time) | |
103 (* -0.00000347 time time time)) | |
104 360.0)) | |
4520 | 105 (moon-anomaly (mod |
957 | 106 (+ 306.0253 |
107 (* 385.81691806 index) | |
108 (* 0.0107306 time time) | |
109 (* 0.00001236 time time time)) | |
110 360.0)) | |
4520 | 111 (moon-lat (mod |
957 | 112 (+ 21.2964 |
113 (* 390.67050646 index) | |
114 (* -0.0016528 time time) | |
115 (* -0.00000239 time time time)) | |
116 360.0)) | |
117 (adjustment | |
118 (if (memq phase '(0 2)) | |
119 (+ (* (- 0.1734 (* 0.000393 time)) | |
120 (solar-sin-degrees sun-anomaly)) | |
121 (* 0.0021 (solar-sin-degrees (* 2 sun-anomaly))) | |
122 (* -0.4068 (solar-sin-degrees moon-anomaly)) | |
123 (* 0.0161 (solar-sin-degrees (* 2 moon-anomaly))) | |
124 (* -0.0004 (solar-sin-degrees (* 3 moon-anomaly))) | |
125 (* 0.0104 (solar-sin-degrees (* 2 moon-lat))) | |
126 (* -0.0051 (solar-sin-degrees (+ sun-anomaly moon-anomaly))) | |
127 (* -0.0074 (solar-sin-degrees (- sun-anomaly moon-anomaly))) | |
128 (* 0.0004 (solar-sin-degrees (+ (* 2 moon-lat) sun-anomaly))) | |
129 (* -0.0004 (solar-sin-degrees (- (* 2 moon-lat) sun-anomaly))) | |
130 (* -0.0006 (solar-sin-degrees | |
131 (+ (* 2 moon-lat) moon-anomaly))) | |
132 (* 0.0010 (solar-sin-degrees (- (* 2 moon-lat) moon-anomaly))) | |
133 (* 0.0005 (solar-sin-degrees | |
134 (+ (* 2 moon-anomaly) sun-anomaly)))) | |
135 (+ (* (- 0.1721 (* 0.0004 time)) | |
136 (solar-sin-degrees sun-anomaly)) | |
137 (* 0.0021 (solar-sin-degrees (* 2 sun-anomaly))) | |
138 (* -0.6280 (solar-sin-degrees moon-anomaly)) | |
139 (* 0.0089 (solar-sin-degrees (* 2 moon-anomaly))) | |
140 (* -0.0004 (solar-sin-degrees (* 3 moon-anomaly))) | |
141 (* 0.0079 (solar-sin-degrees (* 2 moon-lat))) | |
142 (* -0.0119 (solar-sin-degrees (+ sun-anomaly moon-anomaly))) | |
143 (* -0.0047 (solar-sin-degrees (- sun-anomaly moon-anomaly))) | |
144 (* 0.0003 (solar-sin-degrees (+ (* 2 moon-lat) sun-anomaly))) | |
145 (* -0.0004 (solar-sin-degrees (- (* 2 moon-lat) sun-anomaly))) | |
146 (* -0.0006 (solar-sin-degrees (+ (* 2 moon-lat) moon-anomaly))) | |
147 (* 0.0021 (solar-sin-degrees (- (* 2 moon-lat) moon-anomaly))) | |
148 (* 0.0003 (solar-sin-degrees | |
149 (+ (* 2 moon-anomaly) sun-anomaly))) | |
150 (* 0.0004 (solar-sin-degrees | |
151 (- sun-anomaly (* 2 moon-anomaly)))) | |
152 (* -0.0003 (solar-sin-degrees | |
153 (+ (* 2 sun-anomaly) moon-anomaly)))))) | |
154 (adj (+ 0.0028 | |
155 (* -0.0004 (solar-cosine-degrees | |
156 sun-anomaly)) | |
157 (* 0.0003 (solar-cosine-degrees | |
158 moon-anomaly)))) | |
159 (adjustment (cond ((= phase 1) (+ adjustment adj)) | |
160 ((= phase 2) (- adjustment adj)) | |
161 (t adjustment))) | |
162 (date (+ date adjustment)) | |
5213
e080f780f381
(lunar-phase): Add calendar-time-zone to solar ephemeris correction.
Paul Eggert <eggert@twinsun.com>
parents:
4520
diff
changeset
|
163 (date (+ date (/ (- calendar-time-zone |
e080f780f381
(lunar-phase): Add calendar-time-zone to solar ephemeris correction.
Paul Eggert <eggert@twinsun.com>
parents:
4520
diff
changeset
|
164 (solar-ephemeris-correction |
3871
a9f9a058567f
* lunar.el (lunar-phase): Use time conversion from solar.el
Jim Blandy <jimb@redhat.com>
parents:
2247
diff
changeset
|
165 (extract-calendar-year |
a9f9a058567f
* lunar.el (lunar-phase): Use time conversion from solar.el
Jim Blandy <jimb@redhat.com>
parents:
2247
diff
changeset
|
166 (calendar-gregorian-from-absolute |
5213
e080f780f381
(lunar-phase): Add calendar-time-zone to solar ephemeris correction.
Paul Eggert <eggert@twinsun.com>
parents:
4520
diff
changeset
|
167 (truncate date))))) |
e080f780f381
(lunar-phase): Add calendar-time-zone to solar ephemeris correction.
Paul Eggert <eggert@twinsun.com>
parents:
4520
diff
changeset
|
168 60.0 24.0))) |
957 | 169 (time (* 24 (- date (truncate date)))) |
7777
c48a233494e1
(lunar-phase): Revised to use the rewritten and new fcns.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
7760
diff
changeset
|
170 (date (calendar-gregorian-from-absolute (truncate date))) |
13044
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
171 (adj (dst-adjust-time date time))) |
7777
c48a233494e1
(lunar-phase): Revised to use the rewritten and new fcns.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
7760
diff
changeset
|
172 (list (car adj) (apply 'solar-time-string (cdr adj)) phase))) |
957 | 173 |
174 (defun lunar-phase-name (phase) | |
175 "Name of lunar PHASE. | |
176 0 = new moon, 1 = first quarter, 2 = full moon, 3 = last quarter." | |
177 (cond ((= 0 phase) "New Moon") | |
178 ((= 1 phase) "First Quarter Moon") | |
179 ((= 2 phase) "Full Moon") | |
180 ((= 3 phase) "Last Quarter Moon"))) | |
181 | |
182 (defun calendar-phases-of-moon () | |
183 "Create a buffer with the lunar phases for the current calendar window." | |
184 (interactive) | |
185 (message "Computing phases of the moon...") | |
186 (let ((m1 displayed-month) | |
187 (y1 displayed-year) | |
188 (m2 displayed-month) | |
9713
d62e29b1d7a2
Give lunar phases buffer a symbolic name; use new form of calendar-read-date.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
7777
diff
changeset
|
189 (y2 displayed-year)) |
957 | 190 (increment-calendar-month m1 y1 -1) |
191 (increment-calendar-month m2 y2 1) | |
192 (set-buffer (get-buffer-create lunar-phases-buffer)) | |
193 (setq buffer-read-only nil) | |
194 (calendar-set-mode-line | |
5697
47563495ca14
(calendar-phases-of-moon): Fix mode line.
Richard M. Stallman <rms@gnu.org>
parents:
5213
diff
changeset
|
195 (if (= y1 y2) |
47563495ca14
(calendar-phases-of-moon): Fix mode line.
Richard M. Stallman <rms@gnu.org>
parents:
5213
diff
changeset
|
196 (format "Phases of the Moon from %s to %s, %d%%-" |
47563495ca14
(calendar-phases-of-moon): Fix mode line.
Richard M. Stallman <rms@gnu.org>
parents:
5213
diff
changeset
|
197 (calendar-month-name m1) (calendar-month-name m2) y2) |
47563495ca14
(calendar-phases-of-moon): Fix mode line.
Richard M. Stallman <rms@gnu.org>
parents:
5213
diff
changeset
|
198 (format "Phases of the Moon from %s, %d to %s, %d%%-" |
47563495ca14
(calendar-phases-of-moon): Fix mode line.
Richard M. Stallman <rms@gnu.org>
parents:
5213
diff
changeset
|
199 (calendar-month-name m1) y1 (calendar-month-name m2) y2))) |
957 | 200 (erase-buffer) |
201 (insert | |
202 (mapconcat | |
203 '(lambda (x) | |
204 (let ((date (car x)) | |
205 (time (car (cdr x))) | |
206 (phase (car (cdr (cdr x))))) | |
207 (concat (calendar-date-string date) | |
208 ": " | |
209 (lunar-phase-name phase) | |
210 " " | |
211 time))) | |
212 (lunar-phase-list m1 y1) "\n")) | |
213 (goto-char (point-min)) | |
214 (set-buffer-modified-p nil) | |
215 (setq buffer-read-only t) | |
216 (display-buffer lunar-phases-buffer) | |
217 (message "Computing phases of the moon...done"))) | |
218 | |
219 ;;;###autoload | |
220 (defun phases-of-moon (&optional arg) | |
221 "Display the quarters of the moon for last month, this month, and next month. | |
222 If called with an optional prefix argument, prompts for month and year. | |
223 | |
224 This function is suitable for execution in a .emacs file." | |
225 (interactive "P") | |
226 (save-excursion | |
9748
45559582aa9d
Fix use of noday option for calendar-read-date.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9713
diff
changeset
|
227 (let* ((date (if arg |
9713
d62e29b1d7a2
Give lunar phases buffer a symbolic name; use new form of calendar-read-date.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
7777
diff
changeset
|
228 (calendar-read-date t) |
d62e29b1d7a2
Give lunar phases buffer a symbolic name; use new form of calendar-read-date.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
7777
diff
changeset
|
229 (calendar-current-date))) |
d62e29b1d7a2
Give lunar phases buffer a symbolic name; use new form of calendar-read-date.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
7777
diff
changeset
|
230 (displayed-month (extract-calendar-month date)) |
d62e29b1d7a2
Give lunar phases buffer a symbolic name; use new form of calendar-read-date.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
7777
diff
changeset
|
231 (displayed-year (extract-calendar-year date))) |
957 | 232 (calendar-phases-of-moon)))) |
233 | |
234 (defun diary-phases-of-moon () | |
235 "Moon phases diary entry." | |
236 (let* ((index (* 4 | |
237 (truncate | |
238 (* 12.3685 | |
239 (+ (extract-calendar-year date) | |
240 ( / (calendar-day-number date) | |
241 366.0) | |
242 -1900))))) | |
243 (phase (lunar-phase index))) | |
244 (while (calendar-date-compare phase (list date)) | |
245 (setq index (1+ index)) | |
246 (setq phase (lunar-phase index))) | |
247 (if (calendar-date-equal (car phase) date) | |
248 (concat (lunar-phase-name (car (cdr (cdr phase)))) " " | |
249 (car (cdr phase)))))) | |
250 | |
13044
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
251 |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
252 ;; For the Chinese calendar the calculations for the new moon need to be more |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
253 ;; accurate than those above, so we use more terms in the approximation. |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
254 |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
255 (defun lunar-new-moon-time (k) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
256 "Astronomical (Julian) day number of K th new moon." |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
257 (let* ((T (/ k 1236.85)) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
258 (T2 (* T T)) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
259 (T3 (* T T T)) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
260 (T4 (* T2 T2)) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
261 (JDE (+ 2451550.09765 |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
262 (* 29.530588853 k) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
263 (* 0.0001337 T2) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
264 (* -0.000000150 T3) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
265 (* 0.00000000073 T4))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
266 (E (- 1 (* 0.002516 T) (* 0.0000074 T2))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
267 (sun-anomaly (+ 2.5534 |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
268 (* 29.10535669 k) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
269 (* -0.0000218 T2) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
270 (* -0.00000011 T3))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
271 (moon-anomaly (+ 201.5643 |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
272 (* 385.81693528 k) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
273 (* 0.0107438 T2) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
274 (* 0.00001239 T3) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
275 (* -0.000000058 T4))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
276 (moon-argument (+ 160.7108 |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
277 (* 390.67050274 k) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
278 (* -0.0016341 T2) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
279 (* -0.00000227 T3) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
280 (* 0.000000011 T4))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
281 (omega (+ 124.7746 |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
282 (* -1.56375580 k) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
283 (* 0.0020691 T2) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
284 (* 0.00000215 T3))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
285 (A1 (+ 299.77 (* 0.107408 k) (* -0.009173 T2))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
286 (A2 (+ 251.88 (* 0.016321 k))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
287 (A3 (+ 251.83 (* 26.641886 k))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
288 (A4 (+ 349.42 (* 36.412478 k))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
289 (A5 (+ 84.66 (* 18.206239 k))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
290 (A6 (+ 141.74 (* 53.303771 k))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
291 (A7 (+ 207.14 (* 2.453732 k))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
292 (A8 (+ 154.84 (* 7.306860 k))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
293 (A9 (+ 34.52 (* 27.261239 k))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
294 (A10 (+ 207.19 (* 0.121824 k))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
295 (A11 (+ 291.34 (* 1.844379 k))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
296 (A12 (+ 161.72 (* 24.198154 k))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
297 (A13 (+ 239.56 (* 25.513099 k))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
298 (A14 (+ 331.55 (* 3.592518 k))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
299 (correction |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
300 (+ (* -0.40720 (solar-sin-degrees moon-anomaly)) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
301 (* 0.17241 E (solar-sin-degrees sun-anomaly)) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
302 (* 0.01608 (solar-sin-degrees (* 2 moon-anomaly))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
303 (* 0.01039 (solar-sin-degrees (* 2 moon-argument))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
304 (* 0.00739 E (solar-sin-degrees (- moon-anomaly sun-anomaly))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
305 (* -0.00514 E (solar-sin-degrees (+ moon-anomaly sun-anomaly))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
306 (* 0.00208 E E (solar-sin-degrees (* 2 sun-anomaly))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
307 (* -0.00111 (solar-sin-degrees |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
308 (- moon-anomaly (* 2 moon-argument)))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
309 (* -0.00057 (solar-sin-degrees |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
310 (+ moon-anomaly (* 2 moon-argument)))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
311 (* 0.00056 E (solar-sin-degrees |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
312 (+ (* 2 moon-anomaly) sun-anomaly))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
313 (* -0.00042 (solar-sin-degrees (* 3 moon-anomaly))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
314 (* 0.00042 E (solar-sin-degrees |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
315 (+ sun-anomaly (* 2 moon-argument)))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
316 (* 0.00038 E (solar-sin-degrees |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
317 (- sun-anomaly (* 2 moon-argument)))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
318 (* -0.00024 E (solar-sin-degrees |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
319 (- (* 2 moon-anomaly) sun-anomaly))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
320 (* -0.00017 (solar-sin-degrees omega)) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
321 (* -0.00007 (solar-sin-degrees |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
322 (+ moon-anomaly (* 2 sun-anomaly)))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
323 (* 0.00004 (solar-sin-degrees |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
324 (- (* 2 moon-anomaly) (* 2 moon-argument)))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
325 (* 0.00004 (solar-sin-degrees (* 3 sun-anomaly))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
326 (* 0.00003 (solar-sin-degrees (+ moon-anomaly sun-anomaly |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
327 (* -2 moon-argument)))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
328 (* 0.00003 (solar-sin-degrees |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
329 (+ (* 2 moon-anomaly) (* 2 moon-argument)))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
330 (* -0.00003 (solar-sin-degrees (+ moon-anomaly sun-anomaly |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
331 (* 2 moon-argument)))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
332 (* 0.00003 (solar-sin-degrees (- moon-anomaly sun-anomaly |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
333 (* -2 moon-argument)))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
334 (* -0.00002 (solar-sin-degrees (- moon-anomaly sun-anomaly |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
335 (* 2 moon-argument)))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
336 (* -0.00002 (solar-sin-degrees |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
337 (+ (* 3 moon-anomaly) sun-anomaly))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
338 (* 0.00002 (solar-sin-degrees (* 4 moon-anomaly))))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
339 (additional |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
340 (+ (* 0.000325 (solar-sin-degrees A1)) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
341 (* 0.000165 (solar-sin-degrees A2)) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
342 (* 0.000164 (solar-sin-degrees A3)) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
343 (* 0.000126 (solar-sin-degrees A4)) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
344 (* 0.000110 (solar-sin-degrees A5)) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
345 (* 0.000062 (solar-sin-degrees A6)) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
346 (* 0.000060 (solar-sin-degrees A7)) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
347 (* 0.000056 (solar-sin-degrees A8)) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
348 (* 0.000047 (solar-sin-degrees A9)) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
349 (* 0.000042 (solar-sin-degrees A10)) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
350 (* 0.000040 (solar-sin-degrees A11)) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
351 (* 0.000037 (solar-sin-degrees A12)) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
352 (* 0.000035 (solar-sin-degrees A13)) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
353 (* 0.000023 (solar-sin-degrees A14)))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
354 (newJDE (+ JDE correction additional))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
355 (+ newJDE |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
356 (- (solar-ephemeris-correction |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
357 (extract-calendar-year |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
358 (calendar-gregorian-from-absolute |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
359 (floor (calendar-absolute-from-astro newJDE)))))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
360 (/ calendar-time-zone 60.0 24.0)))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
361 |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
362 (defun lunar-new-moon-on-or-after (d) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
363 "Astronomical (Julian) day number of first new moon on or after astronomical |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
364 (Julian) day number d. The fractional part is the time of day. |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
365 |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
366 The date and time are local time, including any daylight savings rules, |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
367 as governed by the values of calendar-daylight-savings-starts, |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
368 calendar-daylight-savings-starts-time, calendar-daylight-savings-ends, |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
369 calendar-daylight-savings-ends-time, calendar-daylight-time-offset, and |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
370 calendar-time-zone." |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
371 (let* ((date (calendar-gregorian-from-absolute |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
372 (floor (calendar-absolute-from-astro d)))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
373 (year (+ (extract-calendar-year date) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
374 (/ (calendar-day-number date) 365.25))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
375 (k (floor (* (- year 2000.0) 12.3685))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
376 (date (lunar-new-moon-time k))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
377 (while (< date d) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
378 (setq k (1+ k)) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
379 (setq date (lunar-new-moon-time k))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
380 (let* ((a-date (calendar-absolute-from-astro date)) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
381 (time (* 24 (- a-date (truncate a-date)))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
382 (date (calendar-gregorian-from-absolute (truncate a-date))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
383 (adj (dst-adjust-time date time))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
384 (calendar-astro-from-absolute |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
385 (+ (calendar-absolute-from-gregorian (car adj)) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
386 (/ (car (cdr adj)) 24.0)))))) |
9155a9ab5de9
Added code to support Chinese calendar.
Edward M. Reingold <reingold@emr.cs.iit.edu>
parents:
9748
diff
changeset
|
387 |
957 | 388 (provide 'lunar) |
389 | |
390 ;;; lunar.el ends here |