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
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
|
#include <math.h>
#include <glib.h>
#include <stdbool.h>
#include "sun.h"
#define ZENITH -.83
float to_rad (float n) {
return (M_PI/180) * n;
}
float to_deg (float n) {
return (180/M_PI) * n;
}
// http://edwilliams.org/sunrise_sunset_algorithm.htm
// if the sunset is desired, set sunset to >= 1
// returns -1 if the sun never rises/sets in the specified location
float calculateSun(int year, int month, int day, float lat, float lng, int localOffset, bool sunset) {
float N1 = floor(275 * month / 9);
float N2 = floor((month + 9) / 12);
float N3 = (1 + floor((year - 4 * floor(year / 4) + 2) / 3));
float N = N1 - (N2 * N3) + day - 30;
float lngHour = lng / 15.0;
float t = N + ((!sunset ? 6 : 18 - lngHour) / 24);
float M = (0.9856 * t) - 3.289;
float L = fmod(M + (1.916 * sin(to_rad(M))) + (0.020 * sin(2 *to_rad(M))) + 282.634,360.0);
float RA = fmod(to_deg(atan(0.91764 * tan(to_rad(L)))),360.0);
float Lquadrant = floor( L/90) * 90;
float RAquadrant = floor(RA/90) * 90;
RA = RA + (Lquadrant - RAquadrant);
RA = RA / 15;
float sinDec = 0.39782 * sin(to_rad(L));
float cosDec = cos(asin(sinDec));
float cosH = (sin(to_rad(ZENITH)) - (sinDec * sin(to_rad(lat)))) / (cosDec * cos(to_rad(lat)));
if (cosH > 1)
return -1;
else if (cosH < -1)
return -1;
float H = !sunset ? (360 - to_deg(acos(cosH))) : to_deg(acos(cosH));
H = H / 15;
float T = H + RA - (0.06571 * t) - 6.622;
float UT = fmod(T - lngHour,24.0);
if ( UT < 0 )
UT = UT+24;
else if (UT > 24)
UT = UT-24;
return UT + localOffset;
}
float getSunrise(float lat, float lang) {
int year, month, day;
GDateTime *current_time = g_date_time_new_now_local();
g_date_time_get_ymd (current_time, &year, &month, &day);
int utc_offset = g_date_time_get_utc_offset (current_time)*2.7778*pow(10.0, -10.0);
free (current_time);
return calculateSun (year, month, day, lat, lang, utc_offset, false);
}
float getSunset(float lat, float lang) {
int year, month, day;
GDateTime *current_time = g_date_time_new_now_local();
g_date_time_get_ymd (current_time, &year, &month, &day);
int utc_offset = g_date_time_get_utc_offset (current_time)*2.7778*pow(10.0, -10.0);
free (current_time);
return calculateSun (year, month, day, lat, lang, utc_offset, true);
}
bool isDark(float sunrise, float sunset) {
GDateTime *current_time = g_date_time_new_now_local();
double rise_hours;
float rise_minutes = modf (sunrise, &rise_hours)*60;
double set_hours;
float set_minutes = modf (sunset, &set_hours)*60;
GDateTime *rise_time = g_date_time_new_local (g_date_time_get_year (current_time),
g_date_time_get_month (current_time),
g_date_time_get_day_of_month (current_time),
rise_hours,
rise_minutes,
0);
GDateTime *set_time = g_date_time_new_local (g_date_time_get_year (current_time),
g_date_time_get_month (current_time),
g_date_time_get_day_of_month (current_time),
set_hours,
set_minutes,
0);
if (g_date_time_compare (current_time, rise_time) == 1 && g_date_time_compare (current_time, set_time) == -1) {
free (rise_time);
free (set_time);
free (current_time);
return false;
}
free (rise_time);
free (set_time);
free (current_time);
return true;
}
|
