calc_sun 0.1.1 → 1.2.6

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Files changed (36) hide show
  1. checksums.yaml +4 -4
  2. checksums.yaml.gz.sig +0 -0
  3. data.tar.gz.sig +0 -0
  4. data/.autotest +15 -0
  5. data/.gitignore +12 -0
  6. data/.source_index +0 -0
  7. data/.travis.yml +5 -0
  8. data/CODE_OF_CONDUCT.md +74 -0
  9. data/Gemfile +19 -28
  10. data/Guardfile +41 -0
  11. data/LICENSE.txt +21 -21
  12. data/Manifest.txt +20 -4
  13. data/README.rdoc +45 -37
  14. data/Rakefile +76 -11
  15. data/calc_sun.gemspec +10 -15
  16. data/dev.md +49 -0
  17. data/example/Almanac.html +1456 -0
  18. data/example/Makefile +10 -0
  19. data/example/adt_helper.rb +91 -0
  20. data/example/askgeo_query.rb +32 -0
  21. data/example/eq_center.md +15 -0
  22. data/example/sunriset.c +670 -0
  23. data/example/sunriset.exe +0 -0
  24. data/example/sunriset.rb +16 -28
  25. data/example/test_my_times.rb +134 -0
  26. data/ext/calc_sun/calc_sun.c +491 -224
  27. data/ext/side_time/extconf.rb +5 -0
  28. data/ext/side_time/side_time.c +117 -0
  29. data/lib/calc_sun.rb +1 -0
  30. data/lib/calc_sun/version.rb +3 -2
  31. data/lib/side_time/side_time.so +0 -0
  32. data/lib/side_time/version.rb +7 -0
  33. data/lib/sidereal_time.rb +4 -0
  34. metadata +37 -62
  35. metadata.gz.sig +0 -0
  36. data/test/calc_sun_test.rb +0 -230
data/calc_sun.gemspec CHANGED
@@ -8,26 +8,21 @@ Gem::Specification.new do |spec|
8
8
  spec.version = CalcSun::VERSION
9
9
  spec.author = 'Douglas Allen'
10
10
  spec.email = ['kb9agt@gmail.com']
11
-
12
11
  spec.summary = 'Calculates Sun Times ea. rise, noon, set'
13
- spec.description = 'supply the jd date, latitude, and longitude in the calls
14
- to each method t_rise(jd, lon, lat),
15
- t_mid_day(jd, lon, lat), and t_set(jd, lon, lat)'
12
+ spec.description = 'supply the date yyyy-mm-dd, latitude decimal, and
13
+ longitude decimal in the calls to each method
14
+ rise(date, lat, lon),
15
+ noon(date, lat, lon), and
16
+ set(date, lat, lon)'
16
17
  spec.homepage = 'https://github.com/DouglasAllen/calc_sun'
17
18
  spec.license = 'MIT'
19
+ require 'English'
20
+ spec.files = `git ls-files -z`.split("\x0").reject do |f|
21
+ f.match(%r{^(test|spec|features)/})
22
+ end
18
23
  require 'rake'
19
- spec.files = FileList['lib .rb',
20
- '[A-Z]*',
21
- 'test/ *'].to_a
22
- spec.extensions << 'ext/eot/extconf.rb'
24
+ spec.extensions << 'ext/calc_sun/extconf.rb'
23
25
  spec.test_files = spec.files.grep(%r{^(test|spec|features)/})
24
26
  spec.require_paths = ['lib']
25
27
  spec.platform = Gem::Platform::RUBY
26
- spec.add_development_dependency 'bundler', '~> 1.13'
27
- spec.add_development_dependency 'rake', '~> 10.0'
28
- spec.add_development_dependency 'rake-compiler', '~> 0.9.3'
29
- spec.add_development_dependency 'hoe', '~> 3.12.0'
30
- spec.add_development_dependency 'guard', '~> 2.6.1'
31
- spec.add_development_dependency 'guard-minitest', '~> 2.3.2'
32
- spec.add_development_dependency 'minitest', '~> 5.4.1'
33
28
  end
data/dev.md ADDED
@@ -0,0 +1,49 @@
1
+
2
+ # calc_sun
3
+
4
+ ## develop
5
+
6
+ ```code
7
+ $> autotest
8
+ ```
9
+
10
+ or
11
+
12
+ ```code
13
+ $> rake
14
+ ```
15
+
16
+ ## example
17
+
18
+ ```code
19
+ $> ruby example/sunriset.rb
20
+ ```
21
+
22
+ Change the lat, lon, and date in sunriset.rb for your location or preferences.
23
+
24
+ pull requests welcome
25
+
26
+ == LICENSE:
27
+
28
+ (The MIT License)
29
+
30
+ Copyright (c) 2016
31
+
32
+ Permission is hereby granted, free of charge, to any person obtaining
33
+ a copy of this software and associated documentation files (the
34
+ 'Software'), to deal in the Software without restriction, including
35
+ without limitation the rights to use, copy, modify, merge, publish,
36
+ distribute, sublicense, and/or sell copies of the Software, and to
37
+ permit persons to whom the Software is furnished to do so, subject to
38
+ the following conditions:
39
+
40
+ The above copyright notice and this permission notice shall be
41
+ included in all copies or substantial portions of the Software.
42
+
43
+ THE SOFTWARE IS PROVIDED 'AS IS', WITHOUT WARRANTY OF ANY KIND,
44
+ EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
45
+ MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.
46
+ IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY
47
+ CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT,
48
+ TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
49
+ SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
data/example/Almanac.html ADDED
@@ -0,0 +1,1456 @@
1
+ <!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
2
+ <!-- saved from url=(0033)http://www.celnav.de/longterm.htm -->
3
+ <html><head><meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
4
+
5
+
6
+
7
+ <title>Almanac</title>
8
+
9
+ <script language="JavaScript">
10
+
11
+ var dtr = Math.PI/180;
12
+
13
+ // Main function
14
+ function Main()
15
+ {
16
+ window.status = "Calculating... Please wait.";
17
+ ReadData();
18
+ TimeMeasures();
19
+ Nutation();
20
+ Aries();
21
+ Sun();
22
+ Moon();
23
+ Polaris();
24
+ MoonPhase();
25
+ Weekday();
26
+ Output();
27
+ window.status = "Done.";
28
+ }
29
+
30
+
31
+ //Auxiliary functions
32
+
33
+ //Sine of angles in degrees
34
+ function sind(x)
35
+ {
36
+ return Math.sin(dtr*x);
37
+ }
38
+
39
+ //Cosine of angles in degrees
40
+ function cosd(x)
41
+ {
42
+ return Math.cos(dtr*x);
43
+ }
44
+
45
+ //Tangent of angles in degrees
46
+ function tand(x)
47
+ {
48
+ return Math.tan(dtr*x);
49
+ }
50
+
51
+ //Normalize large angles
52
+ function norm_360_deg(x)
53
+ {
54
+ while (x<0) {
55
+ x+=360;
56
+ }
57
+ return x%360;
58
+ }
59
+
60
+
61
+ //Output Hour Angle
62
+ function OutHA(x)
63
+ {
64
+ with(Math)
65
+ {
66
+ var GHAdeg = floor(x);
67
+ var GHAmin = floor(60*(x-GHAdeg));
68
+ var GHAsec = round(3600*(x-GHAdeg-GHAmin/60));
69
+ if (GHAsec==60) {GHAsec=0; GHAmin+=1};
70
+ if (GHAmin==60) {GHAmin=0; GHAdeg+=1};
71
+ if (GHAdeg==0) GHAdeg="000";
72
+ else if (GHAdeg<10) GHAdeg="00"+GHAdeg;
73
+ else if (GHAdeg<100) GHAdeg="0"+GHAdeg;
74
+ if (GHAmin==0) GHAmin="00";
75
+ else if (GHAmin<10) GHAmin="0"+GHAmin;
76
+ if (GHAsec<10) GHAsec="0"+GHAsec;
77
+ return " "+GHAdeg+"°"+" "+GHAmin+"'"+" "+GHAsec+"''";
78
+ }
79
+ }
80
+
81
+ //Output Mean Obliqity of Ecliptic
82
+ function OutECL(x)
83
+ {
84
+ with(Math)
85
+ {
86
+ var ECLdeg = floor(x);
87
+ var ECLmin = floor(60*(x-ECLdeg));
88
+ var ECLsec = round(3600000*(x-ECLdeg-ECLmin/60))/1000;
89
+ if (ECLsec==60) {ECLsec=0; ECLmin+=1};
90
+ if (ECLmin==60) {ECLmin=0; ECLdeg+=1};
91
+ if (ECLmin==0) ECLmin="00";
92
+ else if (ECLmin<10) ECLmin="0"+ECLmin;
93
+ if (ECLsec<10) ECLsec="0"+ECLsec;
94
+ return " "+ECLdeg+"°"+" "+ECLmin+"'"+" "+ECLsec+"''";
95
+ }
96
+ }
97
+
98
+ //Output Sidereal Time
99
+ function OutSidTime(x)
100
+ {
101
+ with(Math)
102
+ {
103
+ var GMSTdecimal = x/15;
104
+ var GMSTh = floor(GMSTdecimal);
105
+ var GMSTmdecimal = 60*(GMSTdecimal-GMSTh);
106
+ var GMSTm = floor(GMSTmdecimal);
107
+ var GMSTsdecimal = 60*(GMSTmdecimal-GMSTm);
108
+ var GMSTs = round(1000*GMSTsdecimal)/1000;
109
+ if (GMSTs-floor(GMSTs)==0) GMSTs+=".000";
110
+ else if (10*GMSTs-floor(10*GMSTs)==0) GMSTs+="00";
111
+ else if (100*GMSTs-floor(100*GMSTs)==0) GMSTs+="0";
112
+ return " "+GMSTh+"h"+" "+GMSTm+"m"+" "+GMSTs+"s";
113
+ }
114
+ }
115
+
116
+
117
+ //Output Declination
118
+ function OutDec(x)
119
+ {
120
+ with(Math)
121
+ {
122
+ if (x<0) {var signDEC=-1; var name="S";}
123
+ else {var signDEC=1; var name="N";}
124
+ var DEC = abs(x);
125
+ var DECdeg = floor(DEC);
126
+ var DECmin = floor(60*(DEC-DECdeg));
127
+ var DECsec = round(3600*(DEC-DECdeg-DECmin/60));
128
+ if (DECsec==60) {DECsec=0; DECmin+=1};
129
+ if (DECmin==60) {DECmin=0; DECdeg+=1};
130
+ if (DECdeg==0) DECdeg="00";
131
+ else if (DECdeg<10) DECdeg="0"+DECdeg;
132
+ if (DECmin==0) DECmin="00";
133
+ else if (DECmin<10) DECmin="0"+DECmin;
134
+ if (DECsec<10) DECsec="0"+DECsec;
135
+ return " "+name+" "+DECdeg+"°"+" "+DECmin+"'"+" "+DECsec+"''";
136
+ }
137
+ }
138
+
139
+
140
+ //Output SD and HP
141
+ function OutSDHP(x)
142
+ {
143
+ with(Math)
144
+ {
145
+ x = round(10*x)/10;
146
+ if (x-floor(x)== 0) x+=".0";
147
+ return " "+x+"''";
148
+ }
149
+ }
150
+
151
+
152
+ // Input data conversion
153
+ function ReadData()
154
+ {
155
+ with(document)
156
+ {
157
+ if(InForm.year.value == "") alert("Missing year! Restart calculation.");
158
+ else year = eval(InForm.year.value);
159
+ if(InForm.month.value == "") alert("Missing month! Restart calculation.");
160
+ else month = eval(InForm.month.value);
161
+ if(month < 1 || month > 12) alert("Month out of range! Restart calculation.");
162
+ if(InForm.day.value == "") alert("Missing day! Restart calculation.");
163
+ else day = eval(InForm.day.value);
164
+ if(day < 1 || day > 31) alert("Day out of range! Restart calculation.");
165
+ var schj=0;
166
+ if(year/4-Math.floor(year/4) == 0) schj=1;
167
+ if(year/100-Math.floor(year/100) == 0) schj=0;
168
+ if(year/400-Math.floor(year/400) == 0) schj=1;
169
+ if(month == 2 && day > 28 && schj == 0) alert("February has only 28 days! Restart calculation.");
170
+ if(month == 2 && day > 29 && schj == 1) alert("February has only 29 days in a leap year! Restart calculation.");
171
+ if(month == 4 && day > 30) alert("April has only 30 days! Restart calculation.");
172
+ if(month == 6 && day > 30) alert("June has only 30 days! Restart calculation.");
173
+ if(month == 9 && day > 30) alert("September has only 30 days! Restart calculation.");
174
+ if(month == 11 && day > 30) alert("November has only 30 days! Restart calculation.");
175
+ if(InForm.hour.value == "") var hour = 0;
176
+ else var hour = eval(InForm.hour.value);
177
+ if(InForm.minute.value == "") var minute = 0;
178
+ else var minute = eval(InForm.minute.value);
179
+ if(InForm.second.value == "") var second = 0;
180
+ else var second = eval(InForm.second.value);
181
+ dayfraction = (hour + minute/60 + second/3600)/24;
182
+ if(dayfraction < 0 || dayfraction > 1) alert("Time out of range! Restart calculation.");
183
+ if(InForm.delta.value == "") deltaT = 0;
184
+ else deltaT = eval(InForm.delta.value);
185
+ }
186
+ }
187
+
188
+
189
+ // Calculating Julian date, century, and millennium
190
+ function TimeMeasures()
191
+ {
192
+ with(Math)
193
+ {
194
+ //Julian day (UT1)
195
+ if(month <= 2) {year -=1; month += 12;}
196
+ var A = floor(year/100);
197
+ var B = 2-A+floor(A/4);
198
+ JD0h = floor(365.25*(year+4716))+floor(30.6001*(month+1))+day+B-1524.5;
199
+ JD = JD0h+dayfraction;
200
+
201
+ //Julian centuries (GMT) since 2000 January 0.5
202
+ T = (JD-2451545)/36525;
203
+ T2 = T*T;
204
+ T3 = T*T2;
205
+ T4 = T*T3;
206
+ T5 = T*T4;
207
+
208
+ //Julian ephemeris day (TDT)
209
+ JDE = JD+deltaT/86400;
210
+
211
+ //Julian centuries (TDT) from 2000 January 0.5
212
+ TE = (JDE-2451545)/36525;
213
+ TE2 = TE*TE;
214
+ TE3 = TE*TE2;
215
+ TE4 = TE*TE3;
216
+ TE5 = TE*TE4;
217
+
218
+ //Julian millenniums (TDT) from 2000 January 0.5
219
+ Tau = 0.1*TE;
220
+ Tau2 = Tau*Tau;
221
+ Tau3 = Tau*Tau2;
222
+ Tau4 = Tau*Tau3;
223
+ Tau5 = Tau*Tau4;
224
+ }
225
+ }
226
+
227
+
228
+ //Nutation, obliquity of the ecliptic
229
+ function Nutation()
230
+ {
231
+ //IAU 1980 nutation theory:
232
+
233
+ //Mean anomaly of the moon
234
+ var Mm = 134.962981389+198.867398056*TE+norm_360_deg(477000*TE)+0.008697222222*TE2+TE3/56250;
235
+
236
+ //Mean anomaly of the sun
237
+ var M = 357.527723333+359.05034*TE+norm_360_deg(35640*TE)-0.0001602777778*TE2-TE3/300000;
238
+ SunMA = norm_360_deg(M);
239
+
240
+ //Mean distance of the moon from the ascending node
241
+ var F = 93.271910277+82.017538055*TE+norm_360_deg(483120*TE)-0.0036825*TE2+TE3/327272.7273;
242
+
243
+ //Mean elongation of the moon
244
+ var D = 297.850363055+307.11148*TE+norm_360_deg(444960*TE)-0.001914166667*TE2+TE3/189473.6842;
245
+
246
+ //Longitude of the ascending node of the moon
247
+ var omega = 125.044522222-134.136260833*TE-norm_360_deg(1800*TE)+0.002070833333*TE2+TE3/450000;
248
+
249
+ //Periodic terms for nutation
250
+ var nut = new Array(106);
251
+ nut[0] = " 0 0 0 0 1-171996-174.2 92025 8.9 ";
252
+ nut[1] = " 0 0 2-2 2 -13187 -1.6 5736-3.1 ";
253
+ nut[2] = " 0 0 2 0 2 -2274 -0.2 977-0.5 ";
254
+ nut[3] = " 0 0 0 0 2 2062 0.2 -895 0.5 ";
255
+ nut[4] = " 0-1 0 0 0 -1426 3.4 54-0.1 ";
256
+ nut[5] = " 1 0 0 0 0 712 0.1 -7 0.0 ";
257
+ nut[6] = " 0 1 2-2 2 -517 1.2 224-0.6 ";
258
+ nut[7] = " 0 0 2 0 1 -386 -0.4 200 0.0 ";
259
+ nut[8] = " 1 0 2 0 2 -301 0.0 129-0.1 ";
260
+ nut[9] = " 0-1 2-2 2 217 -0.5 -95 0.3 ";
261
+ nut[10] = "-1 0 0 2 0 158 0.0 -1 0.0 ";
262
+ nut[11] = " 0 0 2-2 1 129 0.1 -70 0.0 ";
263
+ nut[12] = "-1 0 2 0 2 123 0.0 -53 0.0 ";
264
+ nut[13] = " 1 0 0 0 1 63 0.1 -33 0.0 ";
265
+ nut[14] = " 0 0 0 2 0 63 0.0 -2 0.0 ";
266
+ nut[15] = "-1 0 2 2 2 -59 0.0 26 0.0 ";
267
+ nut[16] = "-1 0 0 0 1 -58 -0.1 32 0.0 ";
268
+ nut[17] = " 1 0 2 0 1 -51 0.0 27 0.0 ";
269
+ nut[18] = "-2 0 0 2 0 -48 0.0 1 0.0 ";
270
+ nut[19] = "-2 0 2 0 1 46 0.0 -24 0.0 ";
271
+ nut[20] = " 0 0 2 2 2 -38 0.0 16 0.0 ";
272
+ nut[21] = " 2 0 2 0 2 -31 0.0 13 0.0 ";
273
+ nut[22] = " 2 0 0 0 0 29 0.0 -1 0.0 ";
274
+ nut[23] = " 1 0 2-2 2 29 0.0 -12 0.0 ";
275
+ nut[24] = " 0 0 2 0 0 26 0.0 -1 0.0 ";
276
+ nut[25] = " 0 0 2-2 0 -22 0.0 0 0.0 ";
277
+ nut[26] = "-1 0 2 0 1 21 0.0 -10 0.0 ";
278
+ nut[27] = " 0 2 0 0 0 17 -0.1 0 0.0 ";
279
+ nut[28] = " 0 2 2-2 2 -16 0.1 7 0.0 ";
280
+ nut[29] = "-1 0 0 2 1 16 0.0 -8 0.0 ";
281
+ nut[30] = " 0 1 0 0 1 -15 0.0 9 0.0 ";
282
+ nut[31] = " 1 0 0-2 1 -13 0.0 7 0.0 ";
283
+ nut[32] = " 0-1 0 0 1 -12 0.0 6 0.0 ";
284
+ nut[33] = " 2 0-2 0 0 11 0.0 0 0.0 ";
285
+ nut[34] = "-1 0 2 2 1 -10 0.0 5 0.0 ";
286
+ nut[35] = " 1 0 2 2 2 -8 0.0 3 0.0 ";
287
+ nut[36] = " 0-1 2 0 2 -7 0.0 3 0.0 ";
288
+ nut[37] = " 0 0 2 2 1 -7 0.0 3 0.0 ";
289
+ nut[38] = " 1 1 0-2 0 -7 0.0 0 0.0 ";
290
+ nut[39] = " 0 1 2 0 2 7 0.0 -3 0.0 ";
291
+ nut[40] = "-2 0 0 2 1 -6 0.0 3 0.0 ";
292
+ nut[41] = " 0 0 0 2 1 -6 0.0 3 0.0 ";
293
+ nut[42] = " 2 0 2-2 2 6 0.0 -3 0.0 ";
294
+ nut[43] = " 1 0 0 2 0 6 0.0 0 0.0 ";
295
+ nut[44] = " 1 0 2-2 1 6 0.0 -3 0.0 ";
296
+ nut[45] = " 0 0 0-2 1 -5 0.0 3 0.0 ";
297
+ nut[46] = " 0-1 2-2 1 -5 0.0 3 0.0 ";
298
+ nut[47] = " 2 0 2 0 1 -5 0.0 3 0.0 ";
299
+ nut[48] = " 1-1 0 0 0 5 0.0 0 0.0 ";
300
+ nut[49] = " 1 0 0-1 0 -4 0.0 0 0.0 ";
301
+ nut[50] = " 0 0 0 1 0 -4 0.0 0 0.0 ";
302
+ nut[51] = " 0 1 0-2 0 -4 0.0 0 0.0 ";
303
+ nut[52] = " 1 0-2 0 0 4 0.0 0 0.0 ";
304
+ nut[53] = " 2 0 0-2 1 4 0.0 -2 0.0 ";
305
+ nut[54] = " 0 1 2-2 1 4 0.0 -2 0.0 ";
306
+ nut[55] = " 1 1 0 0 0 -3 0.0 0 0.0 ";
307
+ nut[56] = " 1-1 0-1 0 -3 0.0 0 0.0 ";
308
+ nut[57] = "-1-1 2 2 2 -3 0.0 1 0.0 ";
309
+ nut[58] = " 0-1 2 2 2 -3 0.0 1 0.0 ";
310
+ nut[59] = " 1-1 2 0 2 -3 0.0 1 0.0 ";
311
+ nut[60] = " 3 0 2 0 2 -3 0.0 1 0.0 ";
312
+ nut[61] = "-2 0 2 0 2 -3 0.0 1 0.0 ";
313
+ nut[62] = " 1 0 2 0 0 3 0.0 0 0.0 ";
314
+ nut[63] = "-1 0 2 4 2 -2 0.0 1 0.0 ";
315
+ nut[64] = " 1 0 0 0 2 -2 0.0 1 0.0 ";
316
+ nut[65] = "-1 0 2-2 1 -2 0.0 1 0.0 ";
317
+ nut[66] = " 0-2 2-2 1 -2 0.0 1 0.0 ";
318
+ nut[67] = "-2 0 0 0 1 -2 0.0 1 0.0 ";
319
+ nut[68] = " 2 0 0 0 1 2 0.0 -1 0.0 ";
320
+ nut[69] = " 3 0 0 0 0 2 0.0 0 0.0 ";
321
+ nut[70] = " 1 1 2 0 2 2 0.0 -1 0.0 ";
322
+ nut[71] = " 0 0 2 1 2 2 0.0 -1 0.0 ";
323
+ nut[72] = " 1 0 0 2 1 -1 0.0 0 0.0 ";
324
+ nut[73] = " 1 0 2 2 1 -1 0.0 1 0.0 ";
325
+ nut[74] = " 1 1 0-2 1 -1 0.0 0 0.0 ";
326
+ nut[75] = " 0 1 0 2 0 -1 0.0 0 0.0 ";
327
+ nut[76] = " 0 1 2-2 0 -1 0.0 0 0.0 ";
328
+ nut[77] = " 0 1-2 2 0 -1 0.0 0 0.0 ";
329
+ nut[78] = " 1 0-2 2 0 -1 0.0 0 0.0 ";
330
+ nut[79] = " 1 0-2-2 0 -1 0.0 0 0.0 ";
331
+ nut[80] = " 1 0 2-2 0 -1 0.0 0 0.0 ";
332
+ nut[81] = " 1 0 0-4 0 -1 0.0 0 0.0 ";
333
+ nut[82] = " 2 0 0-4 0 -1 0.0 0 0.0 ";
334
+ nut[83] = " 0 0 2 4 2 -1 0.0 0 0.0 ";
335
+ nut[84] = " 0 0 2-1 2 -1 0.0 0 0.0 ";
336
+ nut[85] = "-2 0 2 4 2 -1 0.0 1 0.0 ";
337
+ nut[86] = " 2 0 2 2 2 -1 0.0 0 0.0 ";
338
+ nut[87] = " 0-1 2 0 1 -1 0.0 0 0.0 ";
339
+ nut[88] = " 0 0-2 0 1 -1 0.0 0 0.0 ";
340
+ nut[89] = " 0 0 4-2 2 1 0.0 0 0.0 ";
341
+ nut[90] = " 0 1 0 0 2 1 0.0 0 0.0 ";
342
+ nut[91] = " 1 1 2-2 2 1 0.0 -1 0.0 ";
343
+ nut[92] = " 3 0 2-2 2 1 0.0 0 0.0 ";
344
+ nut[93] = "-2 0 2 2 2 1 0.0 -1 0.0 ";
345
+ nut[94] = "-1 0 0 0 2 1 0.0 -1 0.0 ";
346
+ nut[95] = " 0 0-2 2 1 1 0.0 0 0.0 ";
347
+ nut[96] = " 0 1 2 0 1 1 0.0 0 0.0 ";
348
+ nut[97] = "-1 0 4 0 2 1 0.0 0 0.0 ";
349
+ nut[98] = " 2 1 0-2 0 1 0.0 0 0.0 ";
350
+ nut[99] = " 2 0 0 2 0 1 0.0 0 0.0 ";
351
+ nut[100]= " 2 0 2-2 1 1 0.0 -1 0.0 ";
352
+ nut[101]= " 2 0-2 0 1 1 0.0 0 0.0 ";
353
+ nut[102]= " 1-1 0-2 0 1 0.0 0 0.0 ";
354
+ nut[103]= "-1 0 0 1 1 1 0.0 0 0.0 ";
355
+ nut[104]= "-1-1 0 2 1 1 0.0 0 0.0 ";
356
+ nut[105]= " 0 1 0 1 0 1 0.0 0 0.0 ";
357
+
358
+ //Reading periodic terms
359
+ var fMm, fM, fF, fD, f_omega, dp=0, de=0;
360
+
361
+ for (x=0; x<105; x++)
362
+ {
363
+ fMm = eval(nut[x].substring(0,2));
364
+ fM = eval(nut[x].substring(2,4));
365
+ fF = eval(nut[x].substring(4,6));
366
+ fD = eval(nut[x].substring(6,8));
367
+ f_omega = eval(nut[x].substring(8,10));
368
+ dp += (eval(nut[x].substring(10,17))+TE*eval(nut[x].substring(17,23)))*sind(fD*D+fM*M+fMm*Mm+fF*F+f_omega*omega);
369
+ de += (eval(nut[x].substring(23,29))+TE*eval(nut[x].substring(29,33)))*cosd(fD*D+fM*M+fMm*Mm+fF*F+f_omega*omega);
370
+ }
371
+
372
+ //Corrections (Herring, 1987)
373
+ /*
374
+ var corr = new Array(4);
375
+ corr[0] = " 0 0 0 0 1-725 417 213 224 ";
376
+ corr[1] = " 0 1 0 0 0 523 61 208 -24 ";
377
+ corr[2] = " 0 0 2-2 2 102-118 -41 -47 ";
378
+ corr[3] = " 0 0 2 0 2 -81 0 32 0 ";
379
+
380
+ for (x=0; x<4; x++)
381
+ {
382
+ fMm = eval(corr[x].substring(0,2));
383
+ fM = eval(corr[x].substring(2,4));
384
+ fF = eval(corr[x].substring(4,6));
385
+ fD = eval(corr[x].substring(6,8));
386
+ f_omega = eval(corr[x].substring(8,10));
387
+ dp += 0.1*(eval(corr[x].substring(10,14))*sind(fD*D+fM*M+fMm*Mm+fF*F+f_omega*omega)+eval(corr[x].substring(14,18))*cosd(fD*D+fM*M+fMm*Mm+fF*F+f_omega*omega));
388
+ de += 0.1*(eval(corr[x].substring(18,22))*cosd(fD*D+fM*M+fMm*Mm+fF*F+f_omega*omega)+eval(corr[x].substring(22,26))*sind(fD*D+fM*M+fMm*Mm+fF*F+f_omega*omega));
389
+ }
390
+ */
391
+
392
+ //Nutation in longitude
393
+ delta_psi = dp/36000000;
394
+
395
+ //Nutation in obliquity
396
+ delta_eps = de/36000000;
397
+
398
+ //Mean obliquity of the ecliptic
399
+ eps0 = (84381.448-46.815*TE-0.00059*TE2+0.001813*TE3)/3600;
400
+
401
+ //True obliquity of the ecliptic
402
+ eps = eps0+delta_eps;
403
+ }
404
+
405
+
406
+ //GHA Aries, GAST, GMST, equation of the equinoxes
407
+ function Aries()
408
+ {
409
+ //Mean GHA Aries
410
+ var GHAAmean = norm_360_deg(280.46061837+ 360.98564736629*(JD-2451545)+0.000387933*T2-T3/38710000);
411
+ ARIESmean = GHAAmean;
412
+ //GMST
413
+ SidTm = OutSidTime(GHAAmean);
414
+
415
+ //True GHA Aries
416
+ GHAAtrue = norm_360_deg(GHAAmean+delta_psi*cosd(eps));
417
+
418
+ //GAST
419
+ SidTa = OutSidTime(GHAAtrue);
420
+
421
+ //Equation of the equinoxes
422
+ EoE = 240*delta_psi*cosd(eps);
423
+ EoE = delta_psi*cosd(eps);
424
+ EoEout = Math.round(1000*EoE)/1000;
425
+ EoEout = " "+EoEout+"s";
426
+ }
427
+
428
+
429
+ //Calculations for the sun
430
+ function Sun()
431
+ {
432
+ with(Math)
433
+ {
434
+ //Periodic terms for the sun
435
+
436
+ //Longitude
437
+ var L0=175347046;
438
+ L0+=3341656*cos(4.6692568+6283.0758500*Tau);
439
+ L0+=34894*cos(4.62610+12566.15170*Tau);
440
+ L0+=3497*cos(2.7441+5753.3849*Tau);
441
+ L0+=3418*cos(2.8289+3.5231*Tau);
442
+ L0+=3136*cos(3.6277+77713.7715*Tau);
443
+ L0+=2676*cos(4.4181+7860.4194*Tau);
444
+ L0+=2343*cos(6.1352+3930.2097*Tau);
445
+ L0+=1324*cos(0.7425+11506.7698*Tau);
446
+ L0+=1273*cos(2.0371+529.6910*Tau);
447
+
448
+ L0+=1199*cos(1.1096+1577.3435*Tau);
449
+ L0+=990*cos(5.233+5884.927*Tau);
450
+ L0+=902*cos(2.045+26.298*Tau);
451
+ L0+=857*cos(3.508+398.149*Tau);
452
+ L0+=780*cos(1.179+5223.694*Tau);
453
+ L0+=753*cos(2.533+5507.553*Tau);
454
+ L0+=505*cos(4.583+18849.228*Tau);
455
+ L0+=492*cos(4.205+775.523*Tau);
456
+ L0+=357*cos(2.920+0.067*Tau);
457
+ L0+=317*cos(5.849+11790.629*Tau);
458
+
459
+ L0+=284*cos(1.899+796.298*Tau);
460
+ L0+=271*cos(0.315+10977.079*Tau);
461
+ L0+=243*cos(0.345+5486.778*Tau);
462
+ L0+=206*cos(4.806+2544.314*Tau);
463
+ L0+=205*cos(1.869+5573.143*Tau);
464
+ L0+=202*cos(2.458+6069.777*Tau);
465
+ L0+=156*cos(0.833+213.299*Tau);
466
+ L0+=132*cos(3.411+2942.463*Tau);
467
+ L0+=126*cos(1.083+20.775*Tau);
468
+ L0+=115*cos(0.645+0.980*Tau);
469
+
470
+ L0+=103*cos(0.636+4694.003*Tau);
471
+ L0+=102*cos(0.976+15720.839*Tau);
472
+ L0+=102*cos(4.267+7.114*Tau);
473
+ L0+=99*cos(6.21+2146.17*Tau);
474
+ L0+=98*cos(0.68+155.42*Tau);
475
+ L0+=86*cos(5.98+161000.69*Tau);
476
+ L0+=85*cos(1.30+6275.96*Tau);
477
+ L0+=85*cos(3.67+71430.70*Tau);
478
+ L0+=80*cos(1.81+17260.15*Tau);
479
+ L0+=79*cos(3.04+12036.46*Tau);
480
+
481
+ L0+=75*cos(1.76+5088.63*Tau);
482
+ L0+=74*cos(3.50+3154.69*Tau);
483
+ L0+=74*cos(4.68+801.82*Tau);
484
+ L0+=70*cos(0.83+9437.76*Tau);
485
+ L0+=62*cos(3.98+8827.39*Tau);
486
+ L0+=61*cos(1.82+7084.90*Tau);
487
+ L0+=57*cos(2.78+6286.60*Tau);
488
+ L0+=56*cos(4.39+14143.50*Tau);
489
+ L0+=56*cos(3.47+6279.55*Tau);
490
+ L0+=52*cos(0.19+12139.55*Tau);
491
+
492
+ L0+=52*cos(1.33+1748.02*Tau);
493
+ L0+=51*cos(0.28+5856.48*Tau);
494
+ L0+=49*cos(0.49+1194.45*Tau);
495
+ L0+=41*cos(5.37+8429.24*Tau);
496
+ L0+=41*cos(2.40+19651.05*Tau);
497
+ L0+=39*cos(6.17+10447.39*Tau);
498
+ L0+=37*cos(6.04+10213.29*Tau);
499
+ L0+=37*cos(2.57+1059.38*Tau);
500
+ L0+=36*cos(1.71+2352.87*Tau);
501
+ L0+=36*cos(1.78+6812.77*Tau);
502
+
503
+ L0+=33*cos(0.59+17789.85*Tau);
504
+ L0+=30*cos(0.44+83996.85*Tau);
505
+ L0+=30*cos(2.74+1349.87*Tau);
506
+ L0+=25*cos(3.16+4690.48*Tau);
507
+
508
+
509
+ var L1=628331966747;
510
+ L1+=206059*cos(2.678235+6283.075850*Tau);
511
+ L1+=4303*cos(2.6351+12566.1517*Tau);
512
+ L1+=425*cos(1.590+3.523*Tau);
513
+ L1+=119*cos(5.796+26.298*Tau);
514
+ L1+=109*cos(2.966+1577.344*Tau);
515
+ L1+=93*cos(2.59+18849.23*Tau);
516
+ L1+=72*cos(1.14+529.69*Tau);
517
+ L1+=68*cos(1.87+398.15*Tau);
518
+ L1+=67*cos(4.41+5507.55*Tau);
519
+
520
+ L1+=59*cos(2.89+5223.69*Tau);
521
+ L1+=56*cos(2.17+155.42*Tau);
522
+ L1+=45*cos(0.40+796.30*Tau);
523
+ L1+=36*cos(0.47+775.52*Tau);
524
+ L1+=29*cos(2.65+7.11*Tau);
525
+ L1+=21*cos(5.34+0.98*Tau);
526
+ L1+=19*cos(1.85+5486.78*Tau);
527
+ L1+=19*cos(4.97+213.30*Tau);
528
+ L1+=17*cos(2.99+6275.96*Tau);
529
+ L1+=16*cos(0.03+2544.31*Tau);
530
+
531
+ L1+=16*cos(1.43+2146.17*Tau);
532
+ L1+=15*cos(1.21+10977.08*Tau);
533
+ L1+=12*cos(2.83+1748.02*Tau);
534
+ L1+=12*cos(3.26+5088.63*Tau);
535
+ L1+=12*cos(5.27+1194.45*Tau);
536
+ L1+=12*cos(2.08+4694.00*Tau);
537
+ L1+=11*cos(0.77+553.57*Tau);
538
+ L1+=10*cos(1.30+6286.60*Tau);
539
+ L1+=10*cos(4.24+1349.87*Tau);
540
+ L1+=9*cos(2.70+242.73*Tau);
541
+
542
+ L1+=9*cos(5.64+951.72*Tau);
543
+ L1+=8*cos(5.30+2352.87*Tau);
544
+ L1+=6*cos(2.65+9437.76*Tau);
545
+ L1+=6*cos(4.67+4690.48*Tau);
546
+
547
+
548
+ var L2=52919;
549
+ L2+=8720*cos(1.0721+6283.0758*Tau);
550
+ L2+=309*cos(0.867+12566.152*Tau);
551
+ L2+=27*cos(0.05+3.52*Tau);
552
+ L2+=16*cos(5.19+26.30*Tau);
553
+ L2+=16*cos(3.68+155.42*Tau);
554
+ L2+=10*cos(0.76+18849.23*Tau);
555
+ L2+=9*cos(2.06+77713.77*Tau);
556
+ L2+=7*cos(0.83+775.52*Tau);
557
+ L2+=5*cos(4.66+1577.34*Tau);
558
+
559
+ L2+=4*cos(1.03+7.11*Tau);
560
+ L2+=4*cos(3.44+5573.14*Tau);
561
+ L2+=3*cos(5.14+796.30*Tau);
562
+ L2+=3*cos(6.05+5507.55*Tau);
563
+ L2+=3*cos(1.19+242.73*Tau);
564
+ L2+=3*cos(6.12+529.69*Tau);
565
+ L2+=3*cos(0.31+398.15*Tau);
566
+ L2+=3*cos(2.28+553.57*Tau);
567
+ L2+=2*cos(4.38+5223.69*Tau);
568
+ L2+=2*cos(3.75+0.98*Tau);
569
+
570
+
571
+ var L3=289*cos(5.844+6283.076*Tau);
572
+ L3+=35;
573
+ L3+=17*cos(5.49+12566.15*Tau);
574
+ L3+=3*cos(5.20+155.42*Tau);
575
+ L3+=1*cos(4.72+3.52*Tau);
576
+ L3+=1*cos(5.30+18849.23*Tau);
577
+ L3+=1*cos(5.97+242.73*Tau);
578
+
579
+
580
+ var L4=114*cos(3.142);
581
+ L4+=8*cos(4.13+6283.08*Tau);
582
+ L4+=1*cos(3.84+12566.15*Tau);
583
+
584
+
585
+ var L5 = 1*cos(3.14);
586
+
587
+ //Mean longitude of the sun
588
+ var Lsun_mean = norm_360_deg(280.4664567+360007.6982779*Tau+0.03032028*Tau2+Tau3/49931-Tau4/15299-Tau5/1988000);
589
+ SunML = Lsun_mean
590
+
591
+ //Heliocentric longitude
592
+ var Lhelioc = norm_360_deg((L0+L1*Tau+L2*Tau2+L3*Tau3+L4*Tau4+L5*Tau5)/1e8/dtr);
593
+ Hlon = Lhelioc
594
+
595
+ //Geocentric longitude
596
+ Lsun_true = norm_360_deg(Lhelioc+180-0.000025);
597
+
598
+ //Latitude
599
+ var B0=280*cos(3.199+84334.662*Tau);
600
+ B0+=102*cos(5.422+5507.553*Tau);
601
+ B0+=80*cos(3.88+5223.69*Tau);
602
+ B0+=44*cos(3.70+2352.87*Tau);
603
+ B0+=32*cos(4.00+1577.34*Tau);
604
+
605
+ var B1=9*cos(3.90+5507.55*Tau);
606
+ B1+=6*cos(1.73+5223.69*Tau);
607
+
608
+ //Heliocentric latitude
609
+ var B = (B0+B1*Tau)/1e8/dtr;
610
+ Hlat = B;
611
+
612
+ //Geocentric latitude
613
+ var beta = norm_360_deg(-B);
614
+ Glat = beta
615
+
616
+ //Corrections
617
+ Lsun_prime = norm_360_deg(Lhelioc+180-1.397*TE-0.00031*TE2);
618
+
619
+ beta = beta+0.000011*(cosd(Lsun_prime)-sind(Lsun_prime));
620
+
621
+
622
+ //Distance earth-sun
623
+ var R0=100013989;
624
+ R0+=1670700*cos(3.0984635+6283.0758500*Tau);
625
+ R0+=13956*cos(3.05525+12566.15170*Tau);
626
+ R0+=3084*cos(5.1985+77713.7715*Tau);
627
+ R0+=1628*cos(1.1739+5753.3849*Tau);
628
+ R0+=1576*cos(2.8469+7860.4194*Tau);
629
+ R0+=925*cos(5.453+11506.770*Tau);
630
+ R0+=542*cos(4.564+3930.210*Tau);
631
+ R0+=472*cos(3.661+5884.927*Tau);
632
+ R0+=346*cos(0.964+5507.553*Tau);
633
+
634
+ R0+=329*cos(5.900+5223.694*Tau);
635
+ R0+=307*cos(0.299+5573.143*Tau);
636
+ R0+=243*cos(4.273+11790.629*Tau);
637
+ R0+=212*cos(5.847+1577.344*Tau);
638
+ R0+=186*cos(5.022+10977.079*Tau);
639
+ R0+=175*cos(3.012+18849.228*Tau);
640
+ R0+=110*cos(5.055+5486.778*Tau);
641
+ R0+=98*cos(0.89+6069.78*Tau);
642
+ R0+=86*cos(5.69+15720.84*Tau);
643
+ R0+=86*cos(1.27+161000.69*Tau);
644
+
645
+ R0+=65*cos(0.27+17260.15*Tau);
646
+ R0+=63*cos(0.92+529.69*Tau);
647
+ R0+=57*cos(2.01+83996.85*Tau);
648
+ R0+=56*cos(5.24+71430.70*Tau);
649
+ R0+=49*cos(3.25+2544.31*Tau);
650
+ R0+=47*cos(2.58+775.52*Tau);
651
+ R0+=45*cos(5.54+9437.76*Tau);
652
+ R0+=43*cos(6.01+6275.96*Tau);
653
+ R0+=39*cos(5.36+4694.00*Tau);
654
+ R0+=38*cos(2.39+8827.39*Tau);
655
+
656
+ R0+=37*cos(0.83+19651.05*Tau);
657
+ R0+=37*cos(4.90+12139.55*Tau);
658
+ R0+=36*cos(1.67+12036.46*Tau);
659
+ R0+=35*cos(1.84+2942.46*Tau);
660
+ R0+=33*cos(0.24+7084.90*Tau);
661
+ R0+=32*cos(0.18+5088.63*Tau);
662
+ R0+=32*cos(1.78+398.15*Tau);
663
+ R0+=28*cos(1.21+6286.60*Tau);
664
+ R0+=28*cos(1.90+6279.55*Tau);
665
+ R0+=26*cos(4.59+10447.39*Tau);
666
+
667
+
668
+ var R1=103019*cos(1.107490+6283.075850*Tau);
669
+ R1+=1721*cos(1.0644+12566.1517*Tau);
670
+ R1+=702*cos(3.142);
671
+ R1+=32*cos(1.02+18849.23*Tau);
672
+ R1+=31*cos(2.84+5507.55*Tau);
673
+ R1+=25*cos(1.32+5223.69*Tau);
674
+ R1+=18*cos(1.42+1577.34*Tau);
675
+ R1+=10*cos(5.91+10977.08*Tau);
676
+ R1+=9*cos(1.42+6275.96*Tau);
677
+ R1+=9*cos(0.27+5486.78*Tau);
678
+
679
+
680
+ var R2=4359*cos(5.7846+6283.0758*Tau);
681
+ R2+=124*cos(5.579+12566.152*Tau);
682
+ R2+=12*cos(3.14);
683
+ R2+=9*cos(3.63+77713.77*Tau);
684
+ R2+=6*cos(1.87+5573.14*Tau);
685
+ R2+=3*cos(5.47+18849.23*Tau);
686
+
687
+
688
+ var R3=145*cos(4.273+6283.076*Tau);
689
+ R3+=7*cos(3.92+12566.15*Tau);
690
+
691
+
692
+ var R4 = 4*cos(2.56+6283.08*Tau);
693
+
694
+ R = (R0+R1*Tau+R2*Tau2+R3*Tau3+R4*Tau4)/1e8;
695
+
696
+
697
+ //Apparent longitude of the sun
698
+ lambda = norm_360_deg(Lsun_true+delta_psi-0.005691611/R);
699
+
700
+ //Right ascension of the sun, apparent
701
+ RAsun = norm_360_deg(atan2((sind(lambda)*cosd(eps)-tand(beta)*sind(eps)),cosd(lambda))/dtr);
702
+
703
+ //Sidereal hour angle of the sun, apparent
704
+ SHAsun = 360-RAsun;
705
+
706
+ //Declination of the sun, apparent
707
+ DECsun = asin(sind(beta)*cosd(eps)+cosd(beta)*sind(eps)*sind(lambda))/dtr;
708
+ Dsun = DECsun;
709
+
710
+ //GHA of the sun
711
+ GHAsun = norm_360_deg(GHAAtrue-RAsun);
712
+
713
+ //Semidiameter of the sun
714
+ SDsun = 959.63/R;
715
+
716
+ //Horizontal parallax of the sun
717
+ HPsun = 8.794/R;
718
+
719
+ //Equation of time
720
+ //EOT = 4*(Lsun_mean-0.0057183-0.0008-RAsun+delta_psi*cosd(eps));
721
+ Eot = 4*GHAsun+720-1440*dayfraction;
722
+ if (Eot>20) Eot-=1440;
723
+ if (Eot<-20) Eot+=1440;
724
+ }
725
+ }
726
+
727
+
728
+ //Calculation of ephemerides for the moon
729
+ function Moon()
730
+ {
731
+ with (Math)
732
+ {
733
+ // Mean longitude of the moon
734
+ Lmoon_mean = norm_360_deg(218.3164591+481267.88134236*TE-0.0013268*TE2+TE3/538841-TE4/65194000);
735
+
736
+ //Mean elongation of the moon
737
+ D = norm_360_deg(297.8502042+445267.1115168*TE-0.00163*TE2+TE3/545868-TE4/113065000);
738
+
739
+ // Mean anomaly of the sun
740
+ Msun_mean = norm_360_deg(357.5291092+35999.0502909*TE-0.0001536*TE2+TE3/24490000);
741
+
742
+ //Mean anomaly of the moon
743
+ Mmoon_mean = norm_360_deg(134.9634114+477198.8676313*TE+0.008997*TE2+TE3/69699-TE4/14712000);
744
+
745
+ //Mean distance of the moon from her ascending node
746
+ F = norm_360_deg(93.2720993+483202.0175273*TE-0.0034029*TE2-TE3/3526000+TE4/863310000);
747
+
748
+ //Corrections
749
+ A1 =119.75+131.849*TE;
750
+ A1 = 360*(A1/360 - floor(A1/360));
751
+ A2 = 53.09+479264.29*TE;
752
+ A2 = 360*(A2/360 - floor(A2/360));
753
+ A3 = 313.45+481266.484*TE;
754
+ A3 = 360*(A3/360 - floor(A3/360));
755
+
756
+ fE = 1-0.002516*TE-0.0000074*TE2;
757
+ fE2 = fE*fE;
758
+
759
+
760
+ //Periodic terms for the moon
761
+
762
+ //Longitude and distance
763
+ fD = new Array(60);
764
+ fMms = new Array(60);
765
+ fMmm = new Array(60);
766
+ fF = new Array(60);
767
+ coeffs = new Array(60);
768
+ coeffc = new Array(60);
769
+
770
+ fD[0]=0; fMms[0]=0; fMmm[0]=1; fF[0]=0; coeffs[0]=6288774; coeffc[0]=-20905355;
771
+ fD[1]=2; fMms[1]=0; fMmm[1]=-1; fF[1]=0; coeffs[1]=1274027; coeffc[1]=-3699111;
772
+ fD[2]=2; fMms[2]=0; fMmm[2]=0; fF[2]=0; coeffs[2]=658314; coeffc[2]=-2955968;
773
+ fD[3]=0; fMms[3]=0; fMmm[3]=2; fF[3]=0; coeffs[3]=213618; coeffc[3]=-569925;
774
+ fD[4]=0; fMms[4]=1; fMmm[4]=0; fF[4]=0; coeffs[4]=-185116; coeffc[4]=48888;
775
+ fD[5]=0; fMms[5]=0; fMmm[5]=0; fF[5]=2; coeffs[5]=-114332; coeffc[5]=-3149;
776
+ fD[6]=2; fMms[6]=0; fMmm[6]=-2; fF[6]=0; coeffs[6]=58793; coeffc[6]=246158;
777
+ fD[7]=2; fMms[7]=-1; fMmm[7]=-1; fF[7]=0; coeffs[7]=57066; coeffc[7]=-152138;
778
+ fD[8]=2; fMms[8]=0; fMmm[8]=1; fF[8]=0; coeffs[8]=53322; coeffc[8]=-170733;
779
+ fD[9]=2; fMms[9]=-1; fMmm[9]=0; fF[9]=0; coeffs[9]=45758; coeffc[9]=-204586;
780
+
781
+ fD[10]=0; fMms[10]=1; fMmm[10]=-1; fF[10]=0; coeffs[10]=-40923; coeffc[10]=-129620;
782
+ fD[11]=1; fMms[11]=0; fMmm[11]=0; fF[11]=0; coeffs[11]=-34720; coeffc[11]=108743;
783
+ fD[12]=0; fMms[12]=1; fMmm[12]=1; fF[12]=0; coeffs[12]=-30383; coeffc[12]=104755;
784
+ fD[13]=2; fMms[13]=0; fMmm[13]=0; fF[13]=-2; coeffs[13]=15327; coeffc[13]=10321;
785
+ fD[14]=0; fMms[14]=0; fMmm[14]=1; fF[14]=2; coeffs[14]=-12528; coeffc[14]=0;
786
+ fD[15]=0; fMms[15]=0; fMmm[15]=1; fF[15]=-2; coeffs[15]=10980; coeffc[15]=79661;
787
+ fD[16]=4; fMms[16]=0; fMmm[16]=-1; fF[16]=0; coeffs[16]=10675; coeffc[16]=-34782;
788
+ fD[17]=0; fMms[17]=0; fMmm[17]=3; fF[17]=0; coeffs[17]=10034; coeffc[17]=-23210;
789
+ fD[18]=4; fMms[18]=0; fMmm[18]=-2; fF[18]=0; coeffs[18]=8548; coeffc[18]=-21636;
790
+ fD[19]=2; fMms[19]=1; fMmm[19]=-1; fF[19]=0; coeffs[19]=-7888; coeffc[19]=24208;
791
+
792
+ fD[20]=2; fMms[20]=1; fMmm[20]=0; fF[20]=0; coeffs[20]=-6766; coeffc[20]=30824;
793
+ fD[21]=1; fMms[21]=0; fMmm[21]=-1; fF[21]=0; coeffs[21]=-5163; coeffc[21]=-8379;
794
+ fD[22]=1; fMms[22]=1; fMmm[22]=0; fF[22]=0; coeffs[22]=4987; coeffc[22]=-16675;
795
+ fD[23]=2; fMms[23]=-1; fMmm[23]=1; fF[23]=0; coeffs[23]=4036; coeffc[23]=-12831;
796
+ fD[24]=2; fMms[24]=0; fMmm[24]=2; fF[24]=0; coeffs[24]=3994; coeffc[24]=-10445;
797
+ fD[25]=4; fMms[25]=0; fMmm[25]=0; fF[25]=0; coeffs[25]=3861; coeffc[25]=-11650;
798
+ fD[26]=2; fMms[26]=0; fMmm[26]=-3; fF[26]=0; coeffs[26]=3665; coeffc[26]=14403;
799
+ fD[27]=0; fMms[27]=1; fMmm[27]=-2; fF[27]=0; coeffs[27]=-2689; coeffc[27]=-7003;
800
+ fD[28]=2; fMms[28]=0; fMmm[28]=-1; fF[28]=2; coeffs[28]=-2602; coeffc[28]=0;
801
+ fD[29]=2; fMms[29]=-1; fMmm[29]=-2; fF[29]=0; coeffs[29]=2390; coeffc[29]=10056;
802
+
803
+ fD[30]=1; fMms[30]=0; fMmm[30]=1; fF[30]=0; coeffs[30]=-2348; coeffc[30]=6322;
804
+ fD[31]=2; fMms[31]=-2; fMmm[31]=0; fF[31]=0; coeffs[31]=2236; coeffc[31]=-9884;
805
+ fD[32]=0; fMms[32]=1; fMmm[32]=2; fF[32]=0; coeffs[32]=-2120; coeffc[32]=5751;
806
+ fD[33]=0; fMms[33]=2; fMmm[33]=0; fF[33]=0; coeffs[33]=-2069; coeffc[33]=0;
807
+ fD[34]=2; fMms[34]=-2; fMmm[34]=-1; fF[34]=0; coeffs[34]=2048; coeffc[34]=-4950;
808
+ fD[35]=2; fMms[35]=0; fMmm[35]=1; fF[35]=-2; coeffs[35]=-1773; coeffc[35]=4130;
809
+ fD[36]=2; fMms[36]=0; fMmm[36]=0; fF[36]=2; coeffs[36]=-1595; coeffc[36]=0;
810
+ fD[37]=4; fMms[37]=-1; fMmm[37]=-1; fF[37]=0; coeffs[37]=1215; coeffc[37]=-3958;
811
+ fD[38]=0; fMms[38]=0; fMmm[38]=2; fF[38]=2; coeffs[38]=-1110; coeffc[38]=0;
812
+ fD[39]=3; fMms[39]=0; fMmm[39]=-1; fF[39]=0; coeffs[39]=-892; coeffc[39]=3258;
813
+
814
+ fD[40]=2; fMms[40]=1; fMmm[40]=1; fF[40]=0; coeffs[40]=-810; coeffc[40]=2616;
815
+ fD[41]=4; fMms[41]=-1; fMmm[41]=-2; fF[41]=0; coeffs[41]=759; coeffc[41]=-1897;
816
+ fD[42]=0; fMms[42]=2; fMmm[42]=-1; fF[42]=0; coeffs[42]=-713; coeffc[42]=-2117;
817
+ fD[43]=2; fMms[43]=2; fMmm[43]=-1; fF[43]=0; coeffs[43]=-700; coeffc[43]=2354;
818
+ fD[44]=2; fMms[44]=1; fMmm[44]=-2; fF[44]=0; coeffs[44]=691; coeffc[44]=0;
819
+ fD[45]=2; fMms[45]=-1; fMmm[45]=0; fF[45]=-2; coeffs[45]=596; coeffc[45]=0;
820
+ fD[46]=4; fMms[46]=0; fMmm[46]=1; fF[46]=0; coeffs[46]=549; coeffc[46]=-1423;
821
+ fD[47]=0; fMms[47]=0; fMmm[47]=4; fF[47]=0; coeffs[47]=537; coeffc[47]=-1117;
822
+ fD[48]=4; fMms[48]=-1; fMmm[48]=0; fF[48]=0; coeffs[48]=520; coeffc[48]=-1571;
823
+ fD[49]=1; fMms[49]=0; fMmm[49]=-2; fF[49]=0; coeffs[49]=-487; coeffc[49]=-1739;
824
+
825
+ fD[50]=2; fMms[50]=1; fMmm[50]=0; fF[50]=-2; coeffs[50]=-399; coeffc[50]=0;
826
+ fD[51]=0; fMms[51]=0; fMmm[51]=2; fF[51]=-2; coeffs[51]=-381; coeffc[51]=-4421;
827
+ fD[52]=1; fMms[52]=1; fMmm[52]=1; fF[52]=0; coeffs[52]=351; coeffc[52]=0;
828
+ fD[53]=3; fMms[53]=0; fMmm[53]=-2; fF[53]=0; coeffs[53]=-340; coeffc[53]=0;
829
+ fD[54]=4; fMms[54]=0; fMmm[54]=-3; fF[54]=0; coeffs[54]=330; coeffc[54]=0;
830
+ fD[55]=2; fMms[55]=-1; fMmm[55]=2; fF[55]=0; coeffs[55]=327; coeffc[55]=0;
831
+ fD[56]=0; fMms[56]=2; fMmm[56]=1; fF[56]=0; coeffs[56]=-323; coeffc[56]=1165;
832
+ fD[57]=1; fMms[57]=1; fMmm[57]=-1; fF[57]=0; coeffs[57]=299; coeffc[57]=0;
833
+ fD[58]=2; fMms[58]=0; fMmm[58]=3; fF[58]=0; coeffs[58]=294; coeffc[58]=0;
834
+ fD[59]=2; fMms[59]=0; fMmm[59]=-1; fF[59]=-2; coeffs[59]=0; coeffc[59]=8752;
835
+
836
+ //Latitude
837
+ fD2 = new Array(60);
838
+ fMms2= new Array(60);
839
+ fMmm2 = new Array(60);
840
+ fF2 = new Array(60);
841
+ coeffs2 = new Array(60);
842
+
843
+ fD2[0]=0; fMms2[0]=0; fMmm2[0]=0; fF2[0]=1; coeffs2[0]=5128122;
844
+ fD2[1]=0; fMms2[1]=0; fMmm2[1]=1; fF2[1]=1; coeffs2[1]=280602;
845
+ fD2[2]=0; fMms2[2]=0; fMmm2[2]=1; fF2[2]=-1; coeffs2[2]=277693;
846
+ fD2[3]=2; fMms2[3]=0; fMmm2[3]=0; fF2[3]=-1; coeffs2[3]=173237;
847
+ fD2[4]=2; fMms2[4]=0; fMmm2[4]=-1; fF2[4]=1; coeffs2[4]=55413;
848
+ fD2[5]=2; fMms2[5]=0; fMmm2[5]=-1; fF2[5]=-1; coeffs2[5]=46271;
849
+ fD2[6]=2; fMms2[6]=0; fMmm2[6]=0; fF2[6]=1; coeffs2[6]=32573;
850
+ fD2[7]=0; fMms2[7]=0; fMmm2[7]=2; fF2[7]=1; coeffs2[7]=17198;
851
+ fD2[8]=2; fMms2[8]=0; fMmm2[8]=1; fF2[8]=-1; coeffs2[8]=9266;
852
+ fD2[9]=0; fMms2[9]=0; fMmm2[9]=2; fF2[9]=-1; coeffs2[9]=8822;
853
+
854
+ fD2[10]=2; fMms2[10]=-1; fMmm2[10]=0; fF2[10]=-1; coeffs2[10]=8216;
855
+ fD2[11]=2; fMms2[11]=0; fMmm2[11]=-2; fF2[11]=-1; coeffs2[11]=4324;
856
+ fD2[12]=2; fMms2[12]=0; fMmm2[12]=1; fF2[12]=1; coeffs2[12]=4200;
857
+ fD2[13]=2; fMms2[13]=1; fMmm2[13]=0; fF2[13]=-1; coeffs2[13]=-3359;
858
+ fD2[14]=2; fMms2[14]=-1; fMmm2[14]=-1; fF2[14]=1; coeffs2[14]=2463;
859
+ fD2[15]=2; fMms2[15]=-1; fMmm2[15]=0; fF2[15]=1; coeffs2[15]=2211;
860
+ fD2[16]=2; fMms2[16]=-1; fMmm2[16]=-1; fF2[16]=-1; coeffs2[16]=2065;
861
+ fD2[17]=0; fMms2[17]=1; fMmm2[17]=-1; fF2[17]=-1; coeffs2[17]=-1870;
862
+ fD2[18]=4; fMms2[18]=0; fMmm2[18]=-1; fF2[18]=-1; coeffs2[18]=1828;
863
+ fD2[19]=0; fMms2[19]=1; fMmm2[19]=0; fF2[19]=1; coeffs2[19]=-1794;
864
+
865
+ fD2[20]=0; fMms2[20]=0; fMmm2[20]=0; fF2[20]=3; coeffs2[20]=-1749;
866
+ fD2[21]=0; fMms2[21]=1; fMmm2[21]=-1; fF2[21]=1; coeffs2[21]=-1565;
867
+ fD2[22]=1; fMms2[22]=0; fMmm2[22]=0; fF2[22]=1; coeffs2[22]=-1491;
868
+ fD2[23]=0; fMms2[23]=1; fMmm2[23]=1; fF2[23]=1; coeffs2[23]=-1475;
869
+ fD2[24]=0; fMms2[24]=1; fMmm2[24]=1; fF2[24]=-1; coeffs2[24]=-1410;
870
+ fD2[25]=0; fMms2[25]=1; fMmm2[25]=0; fF2[25]=-1; coeffs2[25]=-1344;
871
+ fD2[26]=1; fMms2[26]=0; fMmm2[26]=0; fF2[26]=-1; coeffs2[26]=-1335;
872
+ fD2[27]=0; fMms2[27]=0; fMmm2[27]=3; fF2[27]=1; coeffs2[27]=1107;
873
+ fD2[28]=4; fMms2[28]=0; fMmm2[28]=0; fF2[28]=-1; coeffs2[28]=1021;
874
+ fD2[29]=4; fMms2[29]=0; fMmm2[29]=-1; fF2[29]=1; coeffs2[29]=833;
875
+
876
+ fD2[30]=0; fMms2[30]=0; fMmm2[30]=1; fF2[30]=-3; coeffs2[30]=777;
877
+ fD2[31]=4; fMms2[31]=0; fMmm2[31]=-2; fF2[31]=1; coeffs2[31]=671;
878
+ fD2[32]=2; fMms2[32]=0; fMmm2[32]=0; fF2[32]=-3; coeffs2[32]=607;
879
+ fD2[33]=2; fMms2[33]=0; fMmm2[33]=2; fF2[33]=-1; coeffs2[33]=596;
880
+ fD2[34]=2; fMms2[34]=-1; fMmm2[34]=1; fF2[34]=-1; coeffs2[34]=491;
881
+ fD2[35]=2; fMms2[35]=0; fMmm2[35]=-2; fF2[35]=1; coeffs2[35]=-451;
882
+ fD2[36]=0; fMms2[36]=0; fMmm2[36]=3; fF2[36]=-1; coeffs2[36]=439;
883
+ fD2[37]=2; fMms2[37]=0; fMmm2[37]=2; fF2[37]=1; coeffs2[37]=422;
884
+ fD2[38]=2; fMms2[38]=0; fMmm2[38]=-3; fF2[38]=-1; coeffs2[38]=421;
885
+ fD2[39]=2; fMms2[39]=1; fMmm2[39]=-1; fF2[39]=1; coeffs2[39]=-366;
886
+
887
+ fD2[40]=2; fMms2[40]=1; fMmm2[40]=0; fF2[40]=1; coeffs2[40]=-351;
888
+ fD2[41]=4; fMms2[41]=0; fMmm2[41]=0; fF2[41]=1; coeffs2[41]=331;
889
+ fD2[42]=2; fMms2[42]=-1; fMmm2[42]=1; fF2[42]=1; coeffs2[42]=315;
890
+ fD2[43]=2; fMms2[43]=-2; fMmm2[43]=0; fF2[43]=-1; coeffs2[43]=302;
891
+ fD2[44]=0; fMms2[44]=0; fMmm2[44]=1; fF2[44]=3; coeffs2[44]=-283;
892
+ fD2[45]=2; fMms2[45]=1; fMmm2[45]=1; fF2[45]=-1; coeffs2[45]=-229;
893
+ fD2[46]=1; fMms2[46]=1; fMmm2[46]=0; fF2[46]=-1; coeffs2[46]=223;
894
+ fD2[47]=1; fMms2[47]=1; fMmm2[47]=0; fF2[47]=1; coeffs2[47]=223;
895
+ fD2[48]=0; fMms2[48]=1; fMmm2[48]=-2; fF2[48]=-1; coeffs2[48]=-220;
896
+ fD2[49]=2; fMms2[49]=1; fMmm2[49]=-1; fF2[49]=-1; coeffs2[49]=-220;
897
+
898
+ fD2[50]=1; fMms2[50]=0; fMmm2[50]=1; fF2[50]=1; coeffs2[50]=-185;
899
+ fD2[51]=2; fMms2[51]=-1; fMmm2[51]=-2; fF2[51]=-1; coeffs2[51]=181;
900
+ fD2[52]=0; fMms2[52]=1; fMmm2[52]=2; fF2[52]=1; coeffs2[52]=-177;
901
+ fD2[53]=4; fMms2[53]=0; fMmm2[53]=-2; fF2[53]=-1; coeffs2[53]=176;
902
+ fD2[54]=4; fMms2[54]=-1; fMmm2[54]=-1; fF2[54]=-1; coeffs2[54]=166;
903
+ fD2[55]=1; fMms2[55]=0; fMmm2[55]=1; fF2[55]=-1; coeffs2[55]=-164;
904
+ fD2[56]=4; fMms2[56]=0; fMmm2[56]=1; fF2[56]=-1; coeffs2[56]=132;
905
+ fD2[57]=1; fMms2[57]=0; fMmm2[57]=-1; fF2[57]=-1; coeffs2[57]=-119;
906
+ fD2[58]=4; fMms2[58]=-1; fMmm2[58]=0; fF2[58]=-1; coeffs2[58]=115;
907
+ fD2[59]=2; fMms2[59]=-2; fMmm2[59]=0; fF2[59]=1; coeffs2[59]=107;
908
+
909
+
910
+ var sumL = 0; var sumr = 0; var sumB = 0;
911
+
912
+ for (x=0; x<60; x++)
913
+ {
914
+ var f = 1;
915
+ if(abs(fMms[x])==1) f=fE;
916
+ if(abs(fMms[x])==2) f=fE2;
917
+ sumL += f*(coeffs[x]*sind(fD[x]*D+fMms[x]*Msun_mean+fMmm[x]*Mmoon_mean+fF[x]*F));
918
+ sumr += f*(coeffc[x]*cosd(fD[x]*D+fMms[x]*Msun_mean+fMmm[x]*Mmoon_mean+fF[x]*F));
919
+ f = 1;
920
+ if(abs(fMms2[x])==1) f=fE;
921
+ if(abs(fMms2[x])==2) f=fE2;
922
+ sumB += f*(coeffs2[x]*sind(fD2[x]*D+fMms2[x]*Msun_mean+fMmm2[x]*Mmoon_mean+fF2[x]*F));
923
+ }
924
+
925
+ //Corrections
926
+ sumL = sumL+3958*sind(A1)+1962*sind(Lmoon_mean-F)+318*sind(A2);
927
+ sumB = sumB-2235*sind(Lmoon_mean)+382*sind(A3)+175*sind(A1-F)+175*sind(A1+F)+127*sind(Lmoon_mean-Mmoon_mean)-115*sind(Lmoon_mean+Mmoon_mean);
928
+
929
+ //Longitude of the moon
930
+ lambdaMm = norm_360_deg(Lmoon_mean+sumL/1000000);
931
+
932
+ //Latitude of the moon
933
+ betaM = sumB/1000000;
934
+
935
+ //Distance earth-moon
936
+ dEM = 385000.56+sumr/1000;
937
+
938
+ //Apparent longitude of the moon
939
+ lambdaMapp = lambdaMm+delta_psi;
940
+
941
+ //Right ascension of the moon, apparent
942
+ RAmoon = norm_360_deg(atan2((sind(lambdaMapp)*cosd(eps)-tand(betaM)*sind(eps)),cosd(lambdaMapp))/dtr);
943
+
944
+ //Sidereal hour angle of the moon, apparent
945
+ SHAmoon = 360-RAmoon;
946
+
947
+ //Declination of the moon
948
+ DECmoon = asin(sind(betaM)*cosd(eps)+cosd(betaM)*sind(eps)*sind(lambdaMapp))/dtr;
949
+ Dmoon = DECmoon;
950
+
951
+ //GHA of the moon
952
+ GHAmoon = norm_360_deg(GHAAtrue-RAmoon);
953
+
954
+ //Horizontal parallax of the moon
955
+ HPmoon = 3600*asin(6378.14/dEM)/dtr;
956
+
957
+ //Semidiameter of the moon
958
+ SDmoon = 3600*asin(1738/dEM)/dtr;
959
+
960
+ //Geocentric angular distance between moon and sun
961
+ LDist = acos(sind(Dmoon)*sind(Dsun)+cosd(Dmoon)*cosd(Dsun)*cosd(RAmoon-RAsun))/dtr;
962
+
963
+ //Phase of the moon
964
+ var i = lambdaMapp-lambda;
965
+ k = 100*(1-cosd(i))/2;
966
+ k = round(10*k)/10;
967
+ }
968
+ }
969
+
970
+
971
+ //Ephemerides of Polaris
972
+ function Polaris()
973
+ {
974
+ //Equatorial coordinates of Polaris at 2000.0 (mean equinox and equator 2000.0)
975
+ var RApol0 = 37.95293333;
976
+ var DECpol0 = 89.26408889;
977
+
978
+ //Proper motion per year
979
+ var dRApol = 2.98155/3600;
980
+ var dDECpol = -0.0152/3600;
981
+
982
+ //Equatorial coordinates at Julian Date T (mean equinox and equator 2000.0)
983
+ var RApol1 = RApol0+100*TE*dRApol;
984
+ var DECpol1 = DECpol0+100*TE*dDECpol;
985
+
986
+ //Mean obliquity of ecliptic at 2000.0 in degrees
987
+ var eps0_2000 = 23.439291111;
988
+
989
+ //Transformation to ecliptic coordinates in radians (mean equinox and equator 2000.0)
990
+ with(Math)
991
+ {
992
+ var lambdapol1 = atan2((sind(RApol1)*cosd(eps0_2000)+tand(DECpol1)*sind(eps0_2000)),cosd(RApol1));
993
+ var betapol1 = asin(sind(DECpol1)*cosd(eps0_2000)-cosd(DECpol1)*sind(eps0_2000)*sind(RApol1));
994
+ }
995
+
996
+ //Precession
997
+ var eta = (47.0029*TE-0.03302*TE2+0.00006*TE3)*dtr/3600;
998
+ var PI0 = (174.876384-(869.8089*TE+0.03536*TE2)/3600)*dtr;
999
+ var p0 = (5029.0966*TE+1.11113*TE2-0.0000006*TE3)*dtr/3600;
1000
+ with(Math)
1001
+ {
1002
+ var A1 = cos(eta)*cos(betapol1)*sin(PI0-lambdapol1)-sin(eta)*sin(betapol1);
1003
+ var B1 = cos(betapol1)*cos(PI0-lambdapol1);
1004
+ var C1 = cos(eta)*sin(betapol1)+sin(eta)*cos(betapol1)*sin(PI0-lambdapol1);
1005
+ var lambdapol2 = p0+PI0-atan2(A1,B1);
1006
+ var betapol2 = asin(C1);
1007
+ }
1008
+
1009
+ //Nutation in longitude
1010
+ lambdapol2 += dtr*delta_psi;
1011
+
1012
+ //Aberration
1013
+ var kappa = dtr*20.49552/3600;
1014
+ pi0 = dtr*(102.93735+1.71953*TE+0.00046*TE2);
1015
+ var e = 0.016708617-0.000042037*TE-0.0000001236*TE2;
1016
+ with(Math)
1017
+ {
1018
+ dlambdapol = (e*kappa*cos(pi0-lambdapol2)-kappa*cos(dtr*Lsun_true-lambdapol2))/cos(betapol2);
1019
+ dbetapol = -kappa*sin(betapol2)*(sin(dtr*Lsun_true-lambdapol2)-e*sin(pi0-lambdapol2));
1020
+ }
1021
+ lambdapol2 += dlambdapol;
1022
+ betapol2 += dbetapol;
1023
+
1024
+ //Transformation back to equatorial coordinates in radians
1025
+ with(Math)
1026
+ {
1027
+ var RApol2 = atan2((sin(lambdapol2)*cosd(eps)-tan(betapol2)*sind(eps)),cos(lambdapol2));
1028
+ var DECpol2 = asin(sin(betapol2)*cosd(eps)+cos(betapol2)*sind(eps)*sin(lambdapol2));
1029
+ }
1030
+
1031
+ //Finals
1032
+ GHApol = GHAAtrue-RApol2/dtr;
1033
+ GHApol = norm_360_deg(GHApol);
1034
+ SHApol = 360-RApol2/dtr;
1035
+ SHApol = norm_360_deg(SHApol);
1036
+ DECpol = DECpol2/dtr;
1037
+ }
1038
+
1039
+
1040
+ //Calculation of the phase of the moon
1041
+ function MoonPhase()
1042
+ {
1043
+ var x = lambdaMapp-lambda;
1044
+ x = norm_360_deg(x);
1045
+ if(x>0 && x<180) quarter = " (+)";
1046
+ if(x>180 && x<360) quarter = " (-)";
1047
+ if (k==0) quarter = ", new";
1048
+ if (k==100) quarter = ", full";
1049
+ }
1050
+
1051
+ //Day of the week
1052
+ function Weekday()
1053
+ {
1054
+ JD0h += 1.5;
1055
+ var res = JD0h-7*Math.floor(JD0h/7);
1056
+ if (res == 0) DoW = " SUN";
1057
+ if (res == 1) DoW = " MON";
1058
+ if (res == 2) DoW = " TUE";
1059
+ if (res == 3) DoW = " WED";
1060
+ if (res == 4) DoW = " THU";
1061
+ if (res == 5) DoW = " FRI";
1062
+ if (res == 6) DoW = " SAT";
1063
+ }
1064
+
1065
+ // Data output
1066
+ function Output()
1067
+ {
1068
+ //Sun
1069
+ GHASun = OutHA(GHAsun);
1070
+ SHAsun = OutHA(SHAsun);
1071
+ DECSun = OutDec(DECsun);
1072
+ SDsun = OutSDHP(SDsun);
1073
+ HPsun = OutSDHP(HPsun);
1074
+
1075
+ //Moon
1076
+ GHAmoon = OutHA(GHAmoon);
1077
+ SHAmoon = OutHA(SHAmoon);
1078
+ DECmoon = OutDec(DECmoon);
1079
+ SDmoon = OutSDHP(SDmoon);
1080
+ HPmoon = OutSDHP(HPmoon);
1081
+
1082
+ //Aries
1083
+ AR = OutHA(GHAAtrue);
1084
+
1085
+ //Polaris
1086
+ GHApolaris = OutHA(GHApol);
1087
+ SHApolaris = OutHA(SHApol);
1088
+ DECpolaris = OutDec(DECpol);
1089
+
1090
+ //Obliquity of Ecliptic
1091
+ OoE = OutECL(eps0);
1092
+ tOoE = OutECL(eps);
1093
+
1094
+ with(Math)
1095
+ {
1096
+ //Equation of time
1097
+ if (Eot<0) var sign="-";
1098
+ else var sign="+";
1099
+ Eot = abs(Eot);
1100
+ var EOTmin = floor(Eot);
1101
+ var EOTsec = round(600*(Eot-EOTmin))/10;
1102
+ if (EOTsec-floor(EOTsec)==0) EOTsec+=".0";
1103
+ if (EOTmin==0) EOT = " "+sign+" "+EOTsec+"s";
1104
+ else EOT = " "+sign+" "+EOTmin+"m "+EOTsec+"s";
1105
+ Eot = sign+Eot;
1106
+
1107
+ //Illumination
1108
+ var illum=" "+k+"%"+quarter;
1109
+
1110
+ //Lunar Distance of Sun
1111
+ LDist = OutHA(LDist);
1112
+ }
1113
+
1114
+ with(document)
1115
+ {
1116
+ OutForm.SunMA.value = SunMA;
1117
+ OutForm.SunML.value = SunML;
1118
+ OutForm.SunRA.value = RAsun;
1119
+ OutForm.SunDEC.value = DECsun;
1120
+ OutForm.SunGHA.value = GHAsun;
1121
+ OutForm.HLon.value = Hlon;
1122
+ OutForm.GLon.value = Lsun_true;
1123
+ OutForm.HLat.value = Hlat;
1124
+ OutForm.GLat.value = Glat;
1125
+ OutForm.EOT.value = Eot;
1126
+ OutForm.SidTapp.value = GHAAtrue;
1127
+ OutForm.SidTmean.value = ARIESmean;
1128
+ OutForm.EoEquin.value = EoE;
1129
+ OutForm.dpsi.value = delta_psi;//Math.round(3600000*delta_psi)/1000+"''";
1130
+ OutForm.deps.value = delta_eps;//Math.round(3600000*delta_eps)/1000+"''";
1131
+ OutForm.obliq.value = eps0;
1132
+ OutForm.trueobliq.value = eps;
1133
+ OutForm.JulianDay.value = JD//" "+Math.round(1000000*JD)/1000000;
1134
+ OutForm.JulianEphemDay.value = JDE//" "+Math.round(1000000*JDE)/1000000;
1135
+ OutForm.LunarDist.value = " "+LDist;
1136
+ OutForm.DayOfWeek.value = " "+DoW;
1137
+ }
1138
+ }
1139
+
1140
+ </script>
1141
+ </head>
1142
+
1143
+ <body bgcolor="#eeeeee">
1144
+ <p><font color="#500000" size="4">&nbsp;LONG-TERM ALMANAC FOR SUN, MOON, AND POLARIS&nbsp;&nbsp;V1.11</font>
1145
+ <br>
1146
+ <font color="#500000" size="3">&nbsp;<i>Copyright © 2001−2016 Henning Umland</i></font></p>
1147
+
1148
+ <form name="InForm">
1149
+ <thead>
1150
+ <tr>
1151
+ <th><font color="#000080" size="4"><b>&nbsp;Date:&nbsp;&nbsp;</b></font></th>
1152
+ <td>&nbsp; &nbsp; &nbsp;&nbsp;&nbsp;</td>
1153
+ <td>&nbsp; &nbsp; &nbsp;&nbsp;&nbsp;</td>
1154
+ <td>&nbsp; &nbsp; &nbsp;&nbsp;&nbsp;</td>
1155
+ <td>&nbsp; &nbsp; &nbsp;&nbsp;&nbsp;</td>
1156
+ <td>&nbsp; &nbsp; &nbsp;&nbsp;&nbsp;</td>
1157
+ <td>&nbsp; &nbsp; &nbsp;&nbsp;&nbsp;</td>
1158
+ <td>&nbsp; &nbsp; &nbsp;&nbsp;&nbsp;</td>
1159
+ <th><font color="#000080" size="4"><b>Time (UT1):&nbsp;</b></font></th>
1160
+ </tr>
1161
+ </thead>
1162
+ <table border="0">
1163
+ <tbody>
1164
+ <tr>
1165
+ <th align="center"><font color="#000080" size="4">y</b></font></th>
1166
+ <td><input size="4" name="year"></td>
1167
+ <th align="center">m</th>
1168
+ <td><input size="2" name="month"></td>
1169
+ <th align="center">d</th>
1170
+ <td><input size="2" name="day"></td>
1171
+ <td>&nbsp; &nbsp; &nbsp;</td>
1172
+ <th align="center">h</th>
1173
+ <td><input size="2" name="hour"></td>
1174
+ <th align="center">m</th>
1175
+ <td><input size="2" name="minute"></td>
1176
+ <th align="center"><font color="#000080" size="4">s</b></font></th>
1177
+ <td><input size="8" name="second"></td>
1178
+ <td>&nbsp; &nbsp; &nbsp;</td>
1179
+ <td><font color="#000080" size="4"><b>&nbsp;ΔT:&nbsp;</b></font></td>
1180
+ <td><input size="3" name="delta"></td>
1181
+ <td>&nbsp; &nbsp; &nbsp;&nbsp;&nbsp;</td>
1182
+ <td><input onclick="Main()" value="Calculate" type="button">&nbsp;</td>
1183
+ <td><input value="Clear" type="reset"></td>
1184
+ </tr>
1185
+
1186
+ <tr>
1187
+ <td colspan="15">---------------------------------------------------------------------------------------------------------------------</td>
1188
+ </tr>
1189
+ </tbody>
1190
+ </table>
1191
+ </form>
1192
+
1193
+ <form name="OutForm">
1194
+ <table border="0" cellspacing="4">
1195
+ <tbody>
1196
+ <tr>
1197
+ <td>&nbsp;</td>
1198
+ <td align="center">mean anomaly</td>
1199
+ <td></td>
1200
+ <td align="center">mean longitude</td>
1201
+ <td></td>
1202
+ <td align="center">RA</td>
1203
+ <td></td>
1204
+ <td align="center">DEC</td>
1205
+ <td></td>
1206
+ <td align="center">GHA</td>
1207
+ <td></td>
1208
+ </tr>
1209
+ <tr>
1210
+ <td><font color="#500000" size="4"><b>Sun </b></font></td>
1211
+ <td><input size="20" name="SunMA"></td>
1212
+ <td></td>
1213
+ <td><input size="20" name="SunML"></td>
1214
+ <td></td>
1215
+ <td><input size="20" name="SunRA"></td>
1216
+ <td></td>
1217
+ <td><input size="20" name="SunDEC"></td>
1218
+ <td></td>
1219
+ <td><input size="20" name="SunGHA"></td>
1220
+ <td></td>
1221
+ </tr>
1222
+ <tr>
1223
+ <td colspan="12"></td>
1224
+ </tr>
1225
+ <tr>
1226
+ <td>&nbsp;</td>
1227
+ <td align="center">HLon</td>
1228
+ <td></td>
1229
+ <td align="center">GLon</td>
1230
+ <td></td>
1231
+ <td align="center">HLat</td>
1232
+ <td></td>
1233
+ <td align="center">GLat</td>
1234
+ <td></td>
1235
+ <td align="center">Equation of Time</td>
1236
+ </tr>
1237
+ <tr>
1238
+ <td></td>
1239
+ <td><input size="20" name="HLon"></td>
1240
+ <td></td>
1241
+ <td><input size="20" name="GLon"></td>
1242
+ <td></td>
1243
+ <td><input size="20" name="HLat"></td>
1244
+ <td></td>
1245
+ <td><input size="20" name="GLat"></td>
1246
+ <td></td>
1247
+ <td><input size="20" name="EOT"></td>
1248
+ </tr>
1249
+ <tr>
1250
+ <td colspan="12"></td>
1251
+ </tr>
1252
+
1253
+ <tr>
1254
+ <td colspan="12"></td>
1255
+ </tr>
1256
+ <tr>
1257
+ <td colspan="12"></td>
1258
+ </tr>
1259
+ <tr>
1260
+ <td colspan="12"></td>
1261
+ </tr>
1262
+ <tr><td></td></tr>
1263
+ <tr>
1264
+ <td>&nbsp;</td>
1265
+ <td align="center">GAST</td>
1266
+ <td></td>
1267
+ <td align="center">GMST</td>
1268
+ <td></td>
1269
+ <td align="center">Eq. of Equinoxes</td>
1270
+ <td></td>
1271
+ <td align="center">Δψ</td>
1272
+ <td></td>
1273
+ <td align="center">Δε</td>
1274
+ <td></td>
1275
+ <td align="center">Mean Obl. of Ecl.</td>
1276
+ </tr>
1277
+ <tr>
1278
+ <td><font color="#500000" size="4"><b>Misc.</b></font></td>
1279
+ <td><input size="20" name="SidTapp"></td>
1280
+ <td></td>
1281
+ <td><input size="20" name="SidTmean"></td>
1282
+ <td></td>
1283
+ <td><input size="20" name="EoEquin"></td>
1284
+ <td></td>
1285
+ <td><input size="20" name="dpsi"></td>
1286
+ <td></td>
1287
+ <td><input size="20" name="deps"></td>
1288
+ <td></td>
1289
+ <td><input size="20" name="obliq"></td>
1290
+ </tr>
1291
+ <tr>
1292
+ <td colspan="12"></td>
1293
+ </tr>
1294
+ <tr><td></td></tr>
1295
+ <tr>
1296
+ <td>&nbsp;</td>
1297
+ <td align="center">JD</td>
1298
+ <td></td>
1299
+ <td align="center">JDE</td>
1300
+ <td></td>
1301
+ <td align="center">Lunar Dist. of Sun</td>
1302
+ <td></td>
1303
+ <td align="center">Weekday</td>
1304
+ <td></td>
1305
+ <td align="center"></td>
1306
+ <td></td>
1307
+ <td align="center">True Obl. of Ecl.</td>
1308
+ </tr>
1309
+ <tr>
1310
+ <td>&nbsp;</td>
1311
+ <td><input size="20" name="JulianDay"></td>
1312
+ <td></td>
1313
+ <td><input size="20" name="JulianEphemDay"></td>
1314
+ <td></td>
1315
+ <td><input size="20" name="LunarDist"></td>
1316
+ <td></td>
1317
+ <td><input size="20" name="DayOfWeek"></td>
1318
+ <td></td>
1319
+ <td></td>
1320
+ <td></td>
1321
+ <td><input size="20" name="trueobliq"></td>
1322
+ </tr>
1323
+ </tbody>
1324
+ </table>
1325
+ </form>
1326
+
1327
+ <br><br>
1328
+
1329
+ <p align="justify"><b>Description:</b> <br>This computer almanac is based
1330
+ upon formulas published in <i>Astronomical Algorithms</i> by Jean Meeus.
1331
+ The program calculates Greenwich hour angle (GHA), sidereal hour angle
1332
+ (SHA), and declination (Dec) for sun, moon, and polaris.<br><br>
1333
+ Further, the following quantities are provided:<br><br>
1334
+ Greenwich hour angle of the vernal point (Aries)<br>
1335
+ Geocentric semidiameter (SD) and equatorial horizontal parallax (HP) for sun
1336
+ and moon<br>
1337
+ Equation of time<br>
1338
+ Phase of the moon (illuminated fraction of the moon's disk)<br>
1339
+ Greenwich apparent sidereal time (GAST)<br>
1340
+ Greenwich mean sidereal time (GMST)<br>
1341
+ Equation of the equinoxes (= GAST−GMST)<br>
1342
+ Nutation in longitude (Δψ)<br>
1343
+ Nutation in obliquity (Δε)<br>
1344
+ Mean obliquity of the ecliptic<br>
1345
+ True obliquity of the ecliptic (= mean obliquity + Δε)<br>
1346
+ Julian date (JD)<br>
1347
+ Julian ephemeris date (JDE)<br>
1348
+ Geocentric lunar distance of the sun (center-center)<br>
1349
+ Day of the week<br><br>
1350
+ The almanac can be used for many decades, provided the ΔT value (= TT−UT1)
1351
+ for the given date is known. An accuracy of approx. ±1s is sufficient for most
1352
+ applications. Errors in ΔT have a much greater influence on the coordinates
1353
+ of the moon than on the other results. ΔT is obtained through the following
1354
+ formula:<br><br>
1355
+ ΔT = 32.184s + (TAI−UTC) − DUT1 <br><br>
1356
+ Current values for TAI−UTC and DUT1 (= UT1−UTC) are published on the web site
1357
+ of the <a href="http://maia.usno.navy.mil/"><i>IERS Rapid Service / Prediction
1358
+ Center</i></a> (IERS Bulletin A). <br><br>Reliable long-term predictions for
1359
+ ΔT are not possible. Here are some ΔT values of the past:
1360
+ <br><br>
1361
+ 1970.0:&nbsp;&nbsp;+40.2s<br>
1362
+ 1975.0:&nbsp;&nbsp;+45.5s<br>
1363
+ 1980.0:&nbsp;&nbsp;+50.5s<br>
1364
+ 1985.0:&nbsp;&nbsp;+54.3s<br>
1365
+ 1990.0:&nbsp;&nbsp;+56.9s<br>
1366
+ 1995.0:&nbsp;&nbsp;+60.8s<br>
1367
+ 2000.0:&nbsp;&nbsp;+63.8s<br>
1368
+ 2005.0:&nbsp;&nbsp;+64.7s<br>
1369
+ 2010.0:&nbsp;&nbsp;+66.1s<br>
1370
+ 2015.0:&nbsp;&nbsp;+67.6s<br>
1371
+ 2016.0:&nbsp;&nbsp;+68.1s<br><br>
1372
+
1373
+ The program regards any blank dialog box in the time input field as zero. A
1374
+ missing year, month, or day will result in an error message (program must be
1375
+ restarted). The number of the year must be entered in 4-digit format.<br><br>
1376
+ A "(+)" behind the value for the illuminated fraction of the moon's disk
1377
+ indicates a waxing moon, a "(−)" a waning moon.
1378
+ <br><br>
1379
+ SHA and GHA<sub>Aries</sub> refer to the true equinox of date.<br><br>
1380
+
1381
+ </p><b>Accuracy (approximate values):</b><br>
1382
+
1383
+ <table border="0">
1384
+ <tbody>
1385
+ <tr>
1386
+ <td>GHA, SHA, and Dec of the sun: </td>
1387
+ <td>±1"</td>
1388
+ </tr>
1389
+ <tr>
1390
+ <td>GHA and SHA of the moon: &nbsp;&nbsp;</td>
1391
+ <td>±10"</td>
1392
+ </tr>
1393
+ <tr>
1394
+ <td>Dec of the moon: &nbsp;&nbsp;</td>
1395
+ <td>±5"</td>
1396
+ </tr>
1397
+ <tr>
1398
+ <td>GHA and SHA of Polaris: </td>
1399
+ <td>±2"</td>
1400
+ </tr>
1401
+ <tr>
1402
+ <td>Dec of Polaris: </td>
1403
+ <td>±1"</td>
1404
+ </tr>
1405
+ <tr>
1406
+ <td>GHA Aries: &nbsp;&nbsp;</td>
1407
+ <td>±1"</td>
1408
+ </tr>
1409
+ <tr>
1410
+ <td>HP and SD: </td>
1411
+ <td>±0.1"</td>
1412
+ </tr>
1413
+ <tr>
1414
+ <td>Equation of Time: </td>
1415
+ <td>±0.1s</td>
1416
+ </tr>
1417
+ <tr>
1418
+ <td>GAST, GMST, and Equation of Equinoxes: &nbsp;&nbsp;</td>
1419
+ <td>±0.001s</td>
1420
+ </tr>
1421
+ <tr>
1422
+ <td>Nutation (Δψ and Δε): </td>
1423
+ <td>±0.001"</td>
1424
+ </tr>
1425
+ <tr>
1426
+ <td>Mean Obliquity of the Ecliptic: </td>
1427
+ <td>±0.001"</td>
1428
+ </tr>
1429
+ <tr>
1430
+ <td>Lunar distance of sun: </td>
1431
+ <td>±10"</td>
1432
+ </tr>
1433
+ </tbody>
1434
+ </table>
1435
+
1436
+ <br><br>
1437
+
1438
+ <hr>
1439
+ <p align="justify">LICENSE<br><br>
1440
+ This program is free software: you can redistribute it and/or modify
1441
+ it under the terms of the GNU General Public License as published by the Free Software
1442
+ Foundation, either version 3 of the license or any later version.<br><br>
1443
+ This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
1444
+ without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
1445
+ See the GNU General Public License
1446
+ (<a href="http://www.celnav.de/index.htm">http://www.gnu.org/licenses/</a>) for more details.
1447
+ <br><br>
1448
+ Henning Umland &nbsp;&nbsp;&nbsp;N 53° 20' 34'' &nbsp; E 9° 52' 00''
1449
+ <br><br>
1450
+ Check this web site for updated versions:
1451
+ <a href="http://www.celnav.de/index.htm">http://www.celnav.de/index.htm</a></p>
1452
+ <hr>
1453
+
1454
+
1455
+
1456
+ </body></html>