日の出日没計算、日本の空港一覧でやってみよう

これ：
凝ったことはなにもしておりませぬ(numpy なしの生まれたまま python には変更してあります)：
  1 #! /usr/bin/python2.7
  2 # -*- coding: utf-8 -*-
  3 # http://www.hoshi-lab.info/env/solar-j.html の javascript を移植して「空港での…」に
  4 # 仕立て上げたもの、です。凝ったことは何もしてません。
  5 #
  6 import math
  7 
  8 # sin function using degree
  9 def _sind(d):
 10     return math.sin(math.radians(d))
 11 
 12 # cos function using degree
 13 def _cosd(d):
 14     return math.cos(math.radians(d))
 15 
 16 # tan function using degree
 17 def _tand(d):
 18     return math.tan(math.radians(d))
 19 
 20 # calculate Julius year (year from 2000/1/1, for variable "t")
 21 def jy(yy, mm, dd, h, m, s, i): # yy/mm/dd h:m:s, i: time difference
 22     yy -= 2000
 23     if (mm <= 2):
 24         mm += 12
 25         yy -= 1
 26 
 27     k = 365 * yy + 30 * mm + dd - 33.5 - i / 24.0 + math.floor(3 * (mm + 1) / 5.0) \
 28         + math.floor(yy / 4.0) - math.floor(yy / 100.0) + math.floor(yy / 400.0)
 29     k += ((s / 60.0 + m) / 60 + h) / 24.0 # plus time
 30     k += (65 + yy) / 86400.0 # plus delta T
 31     return k / 365.25
 32 
 33 # solar position1 (celestial longitude, degree)
 34 def spls(t): # t: Julius year
 35     L = 280.4603 + 360.00769 * t \
 36         + (1.9146 - 0.00005 * t) * _sind(357.538 + 359.991 * t) \
 37         + 0.0200 * _sind(355.05 +  719.981 * t) \
 38         + 0.0048 * _sind(234.95 +   19.341 * t) \
 39         + 0.0020 * _sind(247.1  +  329.640 * t) \
 40         + 0.0018 * _sind(297.8  + 4452.67  * t) \
 41         + 0.0018 * _sind(251.3  +    0.20  * t) \
 42         + 0.0015 * _sind(343.2  +  450.37  * t) \
 43         + 0.0013 * _sind( 81.4  +  225.18  * t) \
 44         + 0.0008 * _sind(132.5  +  659.29  * t) \
 45         + 0.0007 * _sind(153.3  +   90.38  * t) \
 46         + 0.0007 * _sind(206.8  +   30.35  * t) \
 47         + 0.0006 * _sind( 29.8  +  337.18  * t) \
 48         + 0.0005 * _sind(207.4  +    1.50  * t) \
 49         + 0.0005 * _sind(291.2  +   22.81  * t) \
 50         + 0.0004 * _sind(234.9  +  315.56  * t) \
 51         + 0.0004 * _sind(157.3  +  299.30  * t) \
 52         + 0.0004 * _sind( 21.1  +  720.02  * t) \
 53         + 0.0003 * _sind(352.5  + 1079.97  * t) \
 54         + 0.0003 * _sind(329.7  +   44.43  * t)
 55     while (L >= 360):
 56         L -= 360
 57     while (L < 0):
 58         L += 360
 59     return L
 60 
 61 # solar position2 (distance, AU)
 62 def spds(t): # t: Julius year
 63     r = (0.007256 - 0.0000002 * t) * _sind(267.54 + 359.991 * t) \
 64         + 0.000091 * _sind(265.1 +  719.98 * t) \
 65         + 0.000030 * _sind( 90.0) \
 66         + 0.000013 * _sind( 27.8 + 4452.67 * t) \
 67         + 0.000007 * _sind(254   +  450.4  * t) \
 68         + 0.000007 * _sind(156   +  329.6  * t)
 69     r = math.pow(10, r)
 70     return r
 71 
 72 # solar position3 (declination, degree)
 73 def spal(t): # t: Julius year
 74     ls = spls(t) ;
 75     ep = 23.439291 - 0.000130042 * t
 76     al = math.degrees(math.atan(_tand(ls) * _cosd(ep)))
 77     if ((ls >= 0) and (ls < 180)):
 78         while (al < 0):
 79             al += 180
 80         while (al >= 180):
 81             al -= 180
 82     else:
 83         while (al < 180):
 84             al += 180
 85         while (al >= 360):
 86             al -= 180
 87     return al ;
 88 
 89 # solar position4 (the right ascension, degree)
 90 def spdl(t): # t: Julius year
 91     ls = spls(t)
 92     ep = 23.439291 - 0.000130042 * t
 93     dl = math.degrees(math.asin(_sind(ls) * _sind(ep)))
 94     return dl
 95 
 96 # Calculate sidereal hour (degree)
 97 def sh(t, h, m, s, l, i):
 98     # t: julius year, h: hour, m: minute, s: second,
 99     # l: longitude, i: time difference
100     d = ((s / 60.0 + m) / 60.0 + h) / 24.0 # elapsed hour (from 0:00 a.m.)
101     th = 100.4606 + 360.007700536 * t + 0.00000003879 * t * t - 15 * i
102     th += l + 360 * d
103     while (th >= 360):
104         th -= 360
105     while (th < 0):
106         th += 360
107     return th
108 
109 # Calculating the seeming horizon altitude "sa"(degree)
110 def eandp(alt, ds): # subfunction for altitude and parallax
111     e = 0.035333333 * math.sqrt(alt)
112     p = 0.002442818 / ds
113     return p - e
114 
115 def sa(alt, ds): # alt: altitude (m), ds: solar distance (AU)
116     s = 0.266994444 / ds
117     r = 0.585555555
118     k = eandp(alt,ds) - s - r
119     return k
120 
121 # Calculating solar alititude (degree) {
122 def soal(la, th, al, dl):
123     # la: latitude, th: sidereal hour,
124     # al: solar declination, dl: right ascension
125     h = _sind(dl) * _sind(la) + _cosd(dl) * _cosd(la) * _cosd(th - al)
126     h = math.degrees(math.asin(h))
127     return h
128 
129 # Calculating solar direction (degree) {
130 def sodr(la, th, al, dl):
131     # la: latitude, th: sidereal hour,
132     # al: solar declination, dl: right ascension
133     t = th - al
134     dc = -_cosd(dl) * _sind(t)
135     dm = _sind(dl) * _sind(la) - _cosd(dl) * _cosd(la) * _cosd(t)
136     dr = math.degrees(math.atan2(dc, dm))
137 
138     if (dr < 0):
139         dr += 360
140     return dr
141 
142 def calc(yy, mm, dd, i, lon, lat, alt=0, hh=-1, m=-1):
143     """
144     return list:
145         item 0: astronomical twilight start time
146         item 1: voyage twilight start time
147         item 2: citizen twilight start time
148         item 3: sun rise time, and solar direction (in degrees)
149         item 4: meridian time, and solar altitude (in degrees)
150         item 5: sun set time, and solar direction (in degrees)
151         item 6: citizen twilight end time
152         item 7: voyage twilight end time
153         item 8: astronomical twilight end time
154     """
155     result = []
156 
157     t = jy(yy, mm, dd - 1, 23, 59, 0, i)
158     th = sh(t, 23, 59, 0, lon, i)
159     ds = spds(t)
160     ls = spls(t)
161     alp = spal(t)
162     dlt = spdl(t)
163     pht = soal(lat, th, alp, dlt)
164     pdr = sodr(lat, th, alp, dlt)
165 
166     for hh in range(0, 24):
167         for m in range(0, 60):
168             t = jy(yy, mm, dd, hh, m, 0, i)
169             th = sh(t, hh, m, 0, lon, i)
170             ds = spds(t)
171             ls = spls(t)
172             alp = spal(t)
173             dlt = spdl(t)
174             ht = soal(lat, th, alp, dlt)
175             dr = sodr(lat, th, alp, dlt)
176             tt = eandp(alt, ds)
177             t1 = tt - 18
178             t2 = tt - 12
179             t3 = tt - 6
180             t4 = sa(alt, ds)
181 
182             #key = "%d-%02d-%02d %02d:%02d:00+%02d" % (yy, mm, dd, hh, m, i)
183             key = "%02d:%02d+%02d" % (hh, m, i)
184 
185             if ((pht < t1) and (ht > t1)):
186                 result.append((key, ""))
187             if ((pht < t2) and (ht > t2)):
188                 result.append((key, ""))
189             if ((pht < t3) and (ht > t3)):
190                 result.append((key, ""))
191             if ((pht < t4) and (ht > t4)):
192                 result.append((key, math.floor(dr)))
193             if ((pdr < 180) and (dr > 180)):
194                 result.append((key, math.floor(ht)))
195             if ((pht > t4) and (ht < t4)):
196                 result.append((key, math.floor(dr)))
197             if ((pht > t3) and (ht < t3)):
198                 result.append((key, ""))
199             if ((pht > t2) and (ht < t2)):
200                 result.append((key, ""))
201             if ((pht > t1) and (ht < t1)):
202                 result.append((key, ""))
203             pht = ht
204             pdr = dr
205 
206     #f.result.value = ans
207     return result
208 
209 
210 ALL_AIRPORTS = [
211     ('RJCW', {'City': 'Wakkanai', 'Name': 'Wakkanai', 'Altitude': 30, 'Longitude': 141.800833, 'Latitude': 45.404167, 'Timezone': 9.0}),
212     ('RJER', {'City': 'Rishiri Island', 'Name': 'Rishiri', 'Altitude': 112, 'Longitude': 141.186431, 'Latitude': 45.242006, 'Timezone': 9.0}),
213     ('RJEB', {'City': 'Monbetsu', 'Name': 'Monbetsu', 'Altitude': 80, 'Longitude': 143.404028, 'Latitude': 44.303914, 'Timezone': 9.0}),
214     ('RJCM', {'City': 'Memanbetsu', 'Name': 'Memanbetsu', 'Altitude': 135, 'Longitude': 144.164053, 'Latitude': 43.880606, 'Timezone': 9.0}),
215     ('RJEC', {'City': 'Asahikawa', 'Name': 'Asahikawa', 'Altitude': 721, 'Longitude': 142.4475, 'Latitude': 43.670833, 'Timezone': 9.0}),
216     ('RJCN', {'City': 'Nakashibetsu', 'Name': 'Nakashibetsu', 'Altitude': 234, 'Longitude': 144.96, 'Latitude': 43.5775, 'Timezone': 9.0}),
217     ('RJCO', {'City': 'Sapporo', 'Name': 'Okadama Airport', 'Altitude': 25, 'Longitude': 141.38, 'Latitude': 43.1161, 'Timezone': 9.0}),
218     ('RJCK', {'City': 'Kushiro', 'Name': 'Kushiro Airport', 'Altitude': 327, 'Longitude': 144.193, 'Latitude': 43.041, 'Timezone': 9.0}),
219     ('RJCT', {'City': 'Tokachi', 'Name': 'Tokachi', 'Altitude': 281, 'Longitude': 143.158475, 'Latitude': 42.890544, 'Timezone': 9.0}),
220     ('RJCJ', {'City': 'Chitose', 'Name': 'Chitose', 'Altitude': 87, 'Longitude': 141.666447, 'Latitude': 42.794475, 'Timezone': 9.0}),
221     ('RJCC', {'City': 'Sapporo', 'Name': 'New Chitose', 'Altitude': 82, 'Longitude': 141.692283, 'Latitude': 42.7752, 'Timezone': 9.0}),
222     ('RJCB', {'City': 'Obihiro', 'Name': 'Obihiro', 'Altitude': 505, 'Longitude': 143.217222, 'Latitude': 42.733333, 'Timezone': 9.0}),
223     ('RJEO', {'City': 'Okushiri', 'Name': 'Okushiri', 'Altitude': 161, 'Longitude': 139.432911, 'Latitude': 42.071667, 'Timezone': 9.0}),
224     ('RJCH', {'City': 'Hakodate', 'Name': 'Hakodate', 'Altitude': 151, 'Longitude': 140.821944, 'Latitude': 41.77, 'Timezone': 9.0}),
225     ('RJSA', {'City': 'Aomori', 'Name': 'Aomori', 'Altitude': 664, 'Longitude': 140.690833, 'Latitude': 40.734722, 'Timezone': 9.0}),
226     ('RJSM', {'City': 'Misawa', 'Name': 'Misawa Ab', 'Altitude': 119, 'Longitude': 141.368364, 'Latitude': 40.703222, 'Timezone': 9.0}),
227     ('RJSH', {'City': 'Hachinoe', 'Name': 'Hachinohe', 'Altitude': 152, 'Longitude': 141.466325, 'Latitude': 40.556447, 'Timezone': 9.0}),
228     ('RJSR', {'City': 'Odate Noshiro', 'Name': 'Odate Noshiro Airport', 'Altitude': 292, 'Longitude': 140.371, 'Latitude': 40.1919, 'Timezone': 9.0}),
229     ('RJSK', {'City': 'Akita', 'Name': 'Akita', 'Altitude': 313, 'Longitude': 140.218611, 'Latitude': 39.615556, 'Timezone': 9.0}),
230     ('RJSI', {'City': 'Hanamaki', 'Name': 'Hanamaki', 'Altitude': 297, 'Longitude': 141.135278, 'Latitude': 39.428611, 'Timezone': 9.0}),
231     ('RJSY', {'City': 'Shonai', 'Name': 'Shonai Airport', 'Altitude': 86, 'Longitude': 139.787222, 'Latitude': 38.812222, 'Timezone': 9.0}),
232     ('RJSC', {'City': 'Yamagata', 'Name': 'Yamagata', 'Altitude': 353, 'Longitude': 140.371331, 'Latitude': 38.411894, 'Timezone': 9.0}),
233     ('RJST', {'City': 'Matsushima', 'Name': 'Matsushima', 'Altitude': 7, 'Longitude': 141.219572, 'Latitude': 38.404919, 'Timezone': 9.0}),
234     ('RJSS', {'City': 'Sendai', 'Name': 'Sendai', 'Altitude': 15, 'Longitude': 140.916944, 'Latitude': 38.139722, 'Timezone': 9.0}),
235     ('RJSD', {'City': 'Sado', 'Name': 'Sado', 'Altitude': 88, 'Longitude': 138.413928, 'Latitude': 38.060181, 'Timezone': 9.0}),
236     ('RJSN', {'City': 'Niigata', 'Name': 'Niigata', 'Altitude': 1, 'Longitude': 139.0646, 'Latitude': 37.5711, 'Timezone': 9.0}),
237     ('RJNW', {'City': 'Wajima', 'Name': 'Noto', 'Altitude': 718, 'Longitude': 136.961853, 'Latitude': 37.293097, 'Timezone': 9.0}),
238     ('RJSF', {'City': 'Fukushima', 'Name': 'Fukushima Airport', 'Altitude': 1221, 'Longitude': 140.431, 'Latitude': 37.2274, 'Timezone': 9.0}),
239     ('RJNT', {'City': 'Toyama', 'Name': 'Toyama', 'Altitude': 95, 'Longitude': 137.1875, 'Latitude': 36.648333, 'Timezone': 9.0}),
240     ('RJNK', {'City': 'Kanazawa', 'Name': 'Komatsu', 'Altitude': 36, 'Longitude': 136.406544, 'Latitude': 36.394611, 'Timezone': 9.0}),
241     ('RJNO', {'City': 'Oki Island', 'Name': 'Oki', 'Altitude': 311, 'Longitude': 133.324844, 'Latitude': 36.181125, 'Timezone': 9.0}),
242     ('RJAH', {'City': 'Ibaraki', 'Name': 'Hyakuri', 'Altitude': 105, 'Longitude': 140.415444, 'Latitude': 36.181083, 'Timezone': 9.0}),
243     ('RJAF', {'City': 'Matsumoto', 'Name': 'Matsumoto', 'Altitude': 2182, 'Longitude': 137.922669, 'Latitude': 36.166758, 'Timezone': 9.0}),
244     ('RJNF', {'City': 'Fukui', 'Name': 'Fukui', 'Altitude': 19, 'Longitude': 136.223922, 'Latitude': 36.142847, 'Timezone': 9.0}),
245     ('RJTJ', {'City': 'Iruma', 'Name': 'Iruma', 'Altitude': 295, 'Longitude': 139.410556, 'Latitude': 35.841944, 'Timezone': 9.0}),
246     ('RJTL', {'City': 'Shimofusa', 'Name': 'Shimofusa', 'Altitude': 98, 'Longitude': 140.011111, 'Latitude': 35.798944, 'Timezone': 9.0}),
247     ('RJAA', {'City': 'Tokyo', 'Name': 'Narita Intl', 'Altitude': 141, 'Longitude': 140.386389, 'Latitude': 35.764722, 'Timezone': 9.0}),
248     ('RJTY', {'City': 'Yokota', 'Name': 'Yokota Ab', 'Altitude': 463, 'Longitude': 139.348483, 'Latitude': 35.748492, 'Timezone': 9.0}),
249     ('RJTT', {'City': 'Tokyo', 'Name': 'Tokyo Intl', 'Altitude': 35, 'Longitude': 139.779694, 'Latitude': 35.552258, 'Timezone': 9.0}),
250     ('RJOR', {'City': 'Tottori', 'Name': 'Tottori', 'Altitude': 65, 'Longitude': 134.166553, 'Latitude': 35.530069, 'Timezone': 9.0}),
251     ('RJTR', {'City': 'Zama', 'Name': 'Kastner Aaf', 'Altitude': 360, 'Longitude': 139.393675, 'Latitude': 35.513769, 'Timezone': 9.0}),
252     ('RJBT', {'City': 'Toyooka', 'Name': 'Tajima Airport', 'Altitude': 578, 'Longitude': 134.786944, 'Latitude': 35.512778, 'Timezone': 9.0}),
253     ('RJOH', {'City': 'Miho', 'Name': 'Yonago Kitaro', 'Altitude': 20, 'Longitude': 133.236389, 'Latitude': 35.492222, 'Timezone': 9.0}),
254     ('RJTA', {'City': 'Atsugi', 'Name': 'Atsugi Naf', 'Altitude': 205, 'Longitude': 139.450167, 'Latitude': 35.454611, 'Timezone': 9.0}),
255     ('RJOC', {'City': 'Izumo', 'Name': 'Izumo', 'Altitude': 15, 'Longitude': 132.89, 'Latitude': 35.413611, 'Timezone': 9.0}),
256     ('RJTK', {'City': 'Kisarazu', 'Name': 'Kisarazu', 'Altitude': 10, 'Longitude': 139.909936, 'Latitude': 35.398272, 'Timezone': 9.0}),
257     ('RJNG', {'City': 'Gifu', 'Name': 'Gifu', 'Altitude': 128, 'Longitude': 136.869667, 'Latitude': 35.394078, 'Timezone': 9.0}),
258     ('RJNA', {'City': 'Nagoya', 'Name': 'Nagoya Airport', 'Altitude': 52, 'Longitude': 136.924, 'Latitude': 35.255, 'Timezone': 9.0}),
259     ('RJTE', {'City': 'Tateyama', 'Name': 'Tateyama', 'Altitude': 10, 'Longitude': 139.829208, 'Latitude': 34.987053, 'Timezone': 9.0}),
260     ('RJGG', {'City': 'Nagoya', 'Name': 'Chubu Centrair Intl', 'Altitude': 15, 'Longitude': 136.805408, 'Latitude': 34.858414, 'Timezone': 9.0}),
261     ('RJNY', {'City': 'Yaizu', 'Name': 'Shizuhama', 'Altitude': 23, 'Longitude': 138.298056, 'Latitude': 34.812778, 'Timezone': 9.0}),
262     ('RJNS', {'City': 'Shizuoka', 'Name': 'Mt. Fuji Shizuoka Airport', 'Altitude': 433, 'Longitude': 138.189444, 'Latitude': 34.796111, 'Timezone': 9.0}),
263     ('RJOO', {'City': 'Osaka', 'Name': 'Osaka Intl', 'Altitude': 50, 'Longitude': 135.438222, 'Latitude': 34.785528, 'Timezone': 9.0}),
264     ('RJTO', {'City': 'Oshima', 'Name': 'Oshima', 'Altitude': 130, 'Longitude': 139.360306, 'Latitude': 34.782033, 'Timezone': 9.0}),
265     ('RJOB', {'City': 'Okayama', 'Name': 'Okayama', 'Altitude': 806, 'Longitude': 133.855278, 'Latitude': 34.756944, 'Timezone': 9.0}),
266     ('RJNH', {'City': 'Hamamatsu', 'Name': 'Hamamatsu', 'Altitude': 150, 'Longitude': 137.703083, 'Latitude': 34.750239, 'Timezone': 9.0}),
267     ('RJOW', {'City': 'Iwami', 'Name': 'Iwami Airport', 'Altitude': 184, 'Longitude': 131.79, 'Latitude': 34.6764, 'Timezone': 9.0}),
268     ('RJBE', {'City': 'Kobe', 'Name': 'Kobe', 'Altitude': 22, 'Longitude': 135.223889, 'Latitude': 34.632778, 'Timezone': 9.0}),
269     ('RJOY', {'City': 'Osaka', 'Name': 'Yao', 'Altitude': 39, 'Longitude': 135.602944, 'Latitude': 34.596311, 'Timezone': 9.0}),
270     ('RJBK', {'City': 'Kohnan', 'Name': 'Kohnan', 'Altitude': 3, 'Longitude': 133.933225, 'Latitude': 34.590836, 'Timezone': 9.0}),
271     ('RJOE', {'City': 'Akeno', 'Name': 'Akeno Aero', 'Altitude': 12, 'Longitude': 136.672222, 'Latitude': 34.533333, 'Timezone': 9.0}),
272     ('RJOA', {'City': 'Hiroshima', 'Name': 'Hiroshima', 'Altitude': 1088, 'Longitude': 132.919444, 'Latitude': 34.436111, 'Timezone': 9.0}),
273     ('RJBB', {'City': 'Osaka', 'Name': 'Kansai', 'Altitude': 49, 'Longitude': 135.244167, 'Latitude': 34.4347222, 'Timezone': 9.0}),
274     ('RJBH', {'City': 'Hiroshima', 'Name': 'Hiroshima-Nishi', 'Altitude': 9, 'Longitude': 132.408, 'Latitude': 34.367, 'Timezone': 9.0}),
275     ('RJAN', {'City': 'Niijima', 'Name': 'Niijima Airport', 'Altitude': 133, 'Longitude': 139.268611, 'Latitude': 34.366944, 'Timezone': 9.0}),
276     ('RJDT', {'City': 'Tsushima', 'Name': 'Tsushima', 'Altitude': 213, 'Longitude': 129.33055, 'Latitude': 34.284889, 'Timezone': 9.0}),
277     ('RJOT', {'City': 'Takamatsu', 'Name': 'Takamatsu', 'Altitude': 607, 'Longitude': 134.015556, 'Latitude': 34.214167, 'Timezone': 9.0}),
278     ('RJOI', {'City': 'Iwakuni', 'Name': 'Iwakuni Mcas', 'Altitude': 7, 'Longitude': 132.23575, 'Latitude': 34.14386, 'Timezone': 9.0}),
279     ('RJOS', {'City': 'Tokushima', 'Name': 'Tokushima', 'Altitude': 26, 'Longitude': 134.606639, 'Latitude': 34.132808, 'Timezone': 9.0}),
280     ('RJTQ', {'City': 'Miyakejima', 'Name': 'Miyakejima Airport', 'Altitude': 67, 'Longitude': 139.56, 'Latitude': 34.0736, 'Timezone': 9.0}),
281     ('RJOZ', {'City': 'Ozuki', 'Name': 'Ozuki', 'Altitude': 13, 'Longitude': 131.052144, 'Latitude': 34.045322, 'Timezone': 9.0}),
282     ('RJOF', {'City': 'Hofu', 'Name': 'Hofu', 'Altitude': 7, 'Longitude': 131.549194, 'Latitude': 34.034667, 'Timezone': 9.0}),
283     ('RJDC', {'City': 'Yamaguchi', 'Name': 'Yamaguchi Ube', 'Altitude': 23, 'Longitude': 131.278611, 'Latitude': 33.93, 'Timezone': 9.0}),
284     ('RJFA', {'City': 'Ashiya', 'Name': 'Ashiya', 'Altitude': 98, 'Longitude': 130.653, 'Latitude': 33.883083, 'Timezone': 9.0}),
285     ('RJFR', {'City': 'Kitakyushu', 'Name': 'New Kitakyushu', 'Altitude': 21, 'Longitude': 131.034689, 'Latitude': 33.845942, 'Timezone': 9.0}),
286     ('RJOM', {'City': 'Matsuyama', 'Name': 'Matsuyama', 'Altitude': 25, 'Longitude': 132.699722, 'Latitude': 33.827222, 'Timezone': 9.0}),
287     ('RJDB', {'City': 'Iki', 'Name': 'Iki', 'Altitude': 41, 'Longitude': 129.785417, 'Latitude': 33.749027, 'Timezone': 9.0}),
288     ('RJFZ', {'City': 'Tsuiki', 'Name': 'Tsuiki', 'Altitude': 55, 'Longitude': 131.040278, 'Latitude': 33.685, 'Timezone': 9.0}),
289     ('RJBD', {'City': 'Nanki-shirahama', 'Name': 'Nanki Shirahama', 'Altitude': 298, 'Longitude': 135.364444, 'Latitude': 33.662222, 'Timezone': 9.0}),
290     ('RJFF', {'City': 'Fukuoka', 'Name': 'Fukuoka', 'Altitude': 32, 'Longitude': 130.450686, 'Latitude': 33.585942, 'Timezone': 9.0}),
291     ('RJOK', {'City': 'Kochi', 'Name': 'Kochi', 'Altitude': 42, 'Longitude': 133.669444, 'Latitude': 33.546111, 'Timezone': 9.0}),
292     ('RJFO', {'City': 'Oita', 'Name': 'Oita', 'Altitude': 19, 'Longitude': 131.737222, 'Latitude': 33.479444, 'Timezone': 9.0}),
293     ('RJFS', {'City': 'Saga', 'Name': 'Saga Airport', 'Altitude': 6, 'Longitude': 130.302, 'Latitude': 33.1497, 'Timezone': 9.0}),
294     ('RJTH', {'City': 'Hachijojima', 'Name': 'Hachijojima', 'Altitude': 303, 'Longitude': 139.785833, 'Latitude': 33.115, 'Timezone': 9.0}),
295     ('RJFU', {'City': 'Nagasaki', 'Name': 'Nagasaki', 'Altitude': 15, 'Longitude': 129.913611, 'Latitude': 32.916944, 'Timezone': 9.0}),
296     ('RJFT', {'City': 'Kumamoto', 'Name': 'Kumamoto', 'Altitude': 642, 'Longitude': 130.85505, 'Latitude': 32.837319, 'Timezone': 9.0}),
297     ('RJFE', {'City': 'Fukue', 'Name': 'Fukue', 'Altitude': 273, 'Longitude': 128.832808, 'Latitude': 32.666269, 'Timezone': 9.0}),
298     ('RJDA', {'City': 'Amakusa', 'Name': 'Amakusa Airfield', 'Altitude': 340, 'Longitude': 130.158889, 'Latitude': 32.482222, 'Timezone': 9.0}),
299     ('RJFN', {'City': 'Nyutabaru', 'Name': 'Nyutabaru', 'Altitude': 259, 'Longitude': 131.451389, 'Latitude': 32.083611, 'Timezone': 9.0}),
300     ('RJFM', {'City': 'Miyazaki', 'Name': 'Miyazaki', 'Altitude': 20, 'Longitude': 131.448611, 'Latitude': 31.877222, 'Timezone': 9.0}),
301     ('RJFK', {'City': 'Kagoshima', 'Name': 'Kagoshima', 'Altitude': 906, 'Longitude': 130.719408, 'Latitude': 31.803397, 'Timezone': 9.0}),
302     ('RJFY', {'City': 'Kanoya', 'Name': 'Kanoya', 'Altitude': 214, 'Longitude': 130.845456, 'Latitude': 31.367608, 'Timezone': 9.0}),
303     ('RJFG', {'City': 'Tanegashima', 'Name': 'New Tanegashima', 'Altitude': 768, 'Longitude': 130.991231, 'Latitude': 30.605067, 'Timezone': 9.0}),
304     ('RJFC', {'City': 'Yakushima', 'Name': 'Yakushima', 'Altitude': 124, 'Longitude': 130.659017, 'Latitude': 30.385569, 'Timezone': 9.0}),
305     ('RJKA', {'City': 'Amami', 'Name': 'Amami', 'Altitude': 27, 'Longitude': 129.712542, 'Latitude': 28.430633, 'Timezone': 9.0}),
306     ('RJKI', {'City': 'Kikai', 'Name': 'Kikai Airport', 'Altitude': 15, 'Longitude': 129.928056, 'Latitude': 28.321389, 'Timezone': 9.0}),
307     ('RJKN', {'City': 'Tokunoshima', 'Name': 'Tokunoshima', 'Altitude': 17, 'Longitude': 128.881253, 'Latitude': 27.836381, 'Timezone': 9.0}),
308     ('RJKB', {'City': 'Okierabu', 'Name': 'Okierabu', 'Altitude': 101, 'Longitude': 128.700903, 'Latitude': 27.425522, 'Timezone': 9.0}),
309     ('RORY', {'City': 'Yoron', 'Name': 'Yoron', 'Altitude': 52, 'Longitude': 128.401517, 'Latitude': 27.043964, 'Timezone': 9.0}),
310     ('RODE', {'City': 'Iejima', 'Name': 'Ie Shima Aux Ab', 'Altitude': 184, 'Longitude': 127.761775, 'Latitude': 26.728775, 'Timezone': 9.0}),
311     ('RORA', {'City': 'Aguni', 'Name': 'Aguni Airport', 'Altitude': 38, 'Longitude': 127.240278, 'Latitude': 26.592778, 'Timezone': 9.0}),
312     ('ROKJ', {'City': 'Kumejima', 'Name': 'Kumejima', 'Altitude': 23, 'Longitude': 126.713806, 'Latitude': 26.363506, 'Timezone': 9.0}),
313     ('RODN', {'City': 'Kadena', 'Name': 'Kadena Ab', 'Altitude': 143, 'Longitude': 127.767633, 'Latitude': 26.355612, 'Timezone': 9.0}),
314     ('ROTM', {'City': 'Futema', 'Name': 'Futenma Mcas', 'Altitude': 247, 'Longitude': 127.756494, 'Latitude': 26.274275, 'Timezone': 9.0}),
315     ('ROAH', {'City': 'Okinawa', 'Name': 'Naha', 'Altitude': 12, 'Longitude': 127.645869, 'Latitude': 26.195814, 'Timezone': 9.0}),
316     ('ROKR', {'City': 'Kerama', 'Name': 'Kerama Airport', 'Altitude': 156, 'Longitude': 127.293333, 'Latitude': 26.168333, 'Timezone': 9.0}),
317     ('RORK', {'City': 'Kitadaito', 'Name': 'Kitadaito', 'Altitude': 80, 'Longitude': 131.326944, 'Latitude': 25.944722, 'Timezone': 9.0}),
318     ('ROMD', {'City': 'Minami Daito', 'Name': 'Minami Daito', 'Altitude': 167, 'Longitude': 131.263494, 'Latitude': 25.846533, 'Timezone': 9.0}),
319     ('RORS', {'City': 'Shimojishima', 'Name': 'Shimojishima', 'Altitude': 54, 'Longitude': 125.144722, 'Latitude': 24.826667, 'Timezone': 9.0}),
320     ('RJAW', {'City': 'Iwojima', 'Name': 'Iwo Jima', 'Altitude': 384, 'Longitude': 141.322722, 'Latitude': 24.784, 'Timezone': 10.0}),
321     ('ROMY', {'City': 'Miyako', 'Name': 'Miyako', 'Altitude': 150, 'Longitude': 125.295111, 'Latitude': 24.782833, 'Timezone': 9.0}),
322     ('RORT', {'City': 'Tarama', 'Name': 'Tarama', 'Altitude': 36, 'Longitude': 124.675278, 'Latitude': 24.653889, 'Timezone': 9.0}),
323     ('ROYN', {'City': 'Yonaguni Jima', 'Name': 'Yonaguni', 'Altitude': 70, 'Longitude': 122.977778, 'Latitude': 24.466944, 'Timezone': 9.0}),
324     ('ROIG', {'City': 'Ishigaki', 'Name': 'Ishigaki', 'Altitude': 93, 'Longitude': 124.186983, 'Latitude': 24.344525, 'Timezone': 9.0}),
325     ('RJAM', {'City': 'Minami Tori Shima', 'Name': 'Minami Torishima', 'Altitude': 22, 'Longitude': 153.979119, 'Latitude': 24.289697, 'Timezone': 10.0}),
326 ]
327 #ALL_AIRPORTS_DICT = dict(ALL_AIRPORTS)
328 
329 if __name__ == '__main__':
330     import sys, codecs
331     import time
332 
333     fo = codecs.getwriter('utf-8')(open("jp_airports_sunrise_sunset.html", "w"))
334     
335     tm = time.localtime()
336 
337     fo.write("""\
338 <html dir="ltr" lang="ja">
339 <head>
340 <meta charset="UTF-8">
341 <meta name="robots" content="noindex,nofollow,noodp,noydir">
342     <style type="text/css">
343     <!--
344     body { margin: 2em; }
345     h1 {
346         font-size: 2em;
347         padding-bottom: 0.2em;
348         border-bottom: 5px #770000 solid;
349     }
350     .icao4letter { font-weight: 900; background-color: #eef; }
351     .sunrise { background-color: yellow; }
352     .sunset { background-color: orange; }
353     -->
354     </style>
355 </head>
356 <body>
357 """)
358     fo.write(u"""\
359 <h1>日本の空港の日の出日没時刻(%d-%02d-%02d)</h1>
360 <p>
361 <a href="http://hhsprings.pinoko.jp/site-hhs/2015/02/%%E6%%97%%A5%%E5%%87%%BA%%E6%%97%%A5%%E6%%B2%%A1%%E8%%A8%%88%%E7%%AE%%97%%E3%%80%%81%%E3%%82%%84%%E3%%81%%A3%%E3%%81%%A6%%E3%%81%%BF%%E3%%82%%88%%E3%%81%%86/" target="_blank">日出日没計算、やってみよう</a>の実践、です。CGIでもなんでもないただのバッチで毎日静的ページを作っているだけでございます。半年後を知りたい、のようなことには対応できませんです。
362 
363 <p>
364 <a href="http://www.hoshi-lab.info/env/solar-j.html">
365 近畿大学　星　岳彦 さん作成の日の出と日の入りの計算
366 </a>
367 と同じことをしているはずです、多分。Pythonに移植したままのテストプログラムの結果
368 そのもので、あまり正確なテストはしておりませんし、また、もともと1分の精度は
369 ありませんし、本来地形の影響等々様々な原因で日の出日没時間は現実とはずれます。
370 
371 ですので、本データを使用したいかなる損害についても、当方関知しませんし、「そういうものだ」
372 と思ってくださいますよう。
373 <p>
374 """ % (tm.tm_year, tm.tm_mon, tm.tm_mday))
375     fo.write("""\
376 <table border='1'>
377 """)
378     for icao4letter, airpt_info in ALL_AIRPORTS:
379         airpt_info["ICAO4letter"] = icao4letter
380         r = calc(tm.tm_year, tm.tm_mon, tm.tm_mday,
381                  airpt_info["Timezone"],
382                  airpt_info["Longitude"],
383                  airpt_info["Latitude"],
384                  airpt_info["Altitude"] * 0.3048)
385         fo.write("""\
386 <tr>
387 <tr><td rowspan="5" class="icao4letter">{ICAO4letter}</td><td>Name</td><td>{Name}</td></tr>
388 <tr><td>City</td><td>{City}</td></tr>
389 <tr><td>Lon / Lat</td><td>{Longitude} / {Latitude}</td></tr>
390 <tr><td>Altitude</td><td>{Altitude} [feets]</td></tr>
391 <tr>
392 """.format(**airpt_info))
393         fo.write("""\
394 <td></td>
395 <td>
396 <table border="1">
397 <tr><td>astronomical twilight start</td><td>%s</td><td>  </td></tr>
398 <tr><td>nautical twilight start    </td><td>%s</td><td>  </td></tr>
399 <tr><td>citizen twilight start     </td><td>%s</td><td>  </td></tr>
400 <tr class="sunrise"><td>sun rise                   </td><td>%s</td><td>solar direction (in degrees): %d</td></tr>
401 <tr><td>meridian                   </td><td>%s</td><td>solar altitude (in degrees): %d</td></tr>
402 <tr class="sunset"><td>sun set                    </td><td>%s</td><td>solar direction (in degrees): %d</td></tr>
403 <tr><td>citizen twilight end       </td><td>%s</td><td>  </td></tr>
404 <tr><td>nautical twilight end      </td><td>%s</td><td>  </td></tr>
405 <tr><td>astronomical twilight end  </td><td>%s</td><td>  </td></tr>
406 </table>
407 </td>
408 </tr>
409 </tr>
410 
411 """ % (
412                 r[0][0], #astronomical twilight start
413                 r[1][0], #nautical twilight start
414                 r[2][0], #citizen twilight start
415                 r[3][0], r[3][1], #sun rise
416                 r[4][0], r[4][1], #meridian
417                 r[5][0], r[5][1], #sun set
418                 r[6][0], #citizen twilight end
419                 r[7][0], #nautical twilight end
420                 r[8][0], #astronomical twilight end
421                 ))
422     fo.write("""\
423 </table>
424 </body>
425 </html>""")
426     fo.close()
2015年2月
日	月	火	水	木	金	土
« 1月				3月 »
1	2	3	4	5	6	7
8	9	10	11	12	13	14
15	16	17	18	19	20	21
22	23	24	25	26	27	28
Related Posts
