sets the initial brightness based on the time

2024-11-01 21:42:12 +13:00 · 2024-11-01 21:42:12 +13:00 · 041c51042d
commit 041c51042d
parent 50b724dabc
2 changed files with 190 additions and 25 deletions
--- a/sunriseSunset.h
+++ b/sunriseSunset.h
@ -0,0 +1,142 @@
+ 
+#include <math.h>
+//#define PI 3.1415926
+#define ZENITH -.83
+
+/* zenith calc
+offical      = 90 degrees 50'
+civil        = 96 degrees
+nautical     = 102 degrees
+astronomical = 108 degrees
+http://edwilliams.org/sunrise_sunset_algorithm.htm
+*/
+float calculateSunrise(int year,int month,int day,float lat, float lng,int localOffset, int daylightSavings) {
+    /*
+    localOffset will be <0 for western hemisphere and >0 for eastern hemisphere
+    daylightSavings should be 1 if it is in effect during the summer otherwise it should be 0
+    */
+    //1. first calculate the day of the year
+    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;
+
+    //2. convert the longitude to hour value and calculate an approximate time
+    float lngHour = lng / 15.0;      
+    float t = N + ((6 - lngHour) / 24);   //if rising time is desired:
+    //float t = N + ((18 - lngHour) / 24)   //if setting time is desired:
+
+    //3. calculate the Sun's mean anomaly   
+    float M = (0.9856 * t) - 3.289;
+
+    //4. calculate the Sun's true longitude
+    float L = fmod(M + (1.916 * sin((PI/180)*M)) + (0.020 * sin(2 *(PI/180) * M)) + 282.634,360.0);
+
+    //5a. calculate the Sun's right ascension      
+    float RA = fmod(180/PI*atan(0.91764 * tan((PI/180)*L)),360.0);
+
+    //5b. right ascension value needs to be in the same quadrant as L   
+    float Lquadrant  = floor( L/90) * 90;
+    float RAquadrant = floor(RA/90) * 90;
+    RA = RA + (Lquadrant - RAquadrant);
+
+    //5c. right ascension value needs to be converted into hours   
+    RA = RA / 15;
+
+    //6. calculate the Sun's declination
+    float sinDec = 0.39782 * sin((PI/180)*L);
+    float cosDec = cos(asin(sinDec));
+
+    //7a. calculate the Sun's local hour angle
+    float cosH = (sin((PI/180)*ZENITH) - (sinDec * sin((PI/180)*lat))) / (cosDec * cos((PI/180)*lat));
+    /*   
+    if (cosH >  1) 
+    the sun never rises on this location (on the specified date)
+    if (cosH < -1)
+    the sun never sets on this location (on the specified date)
+    */
+
+    //7b. finish calculating H and convert into hours
+    float H = 360 - (180/PI)*acos(cosH);   //   if if rising time is desired:
+    //float H = acos(cosH) //   if setting time is desired:      
+    H = H / 15;
+
+    //8. calculate local mean time of rising/setting      
+    float T = H + RA - (0.06571 * t) - 6.622;
+
+    //9. adjust back to UTC
+    float UT = fmod(T - lngHour,24.0);
+
+    //10. convert UT value to local time zone of latitude/longitude
+    return UT + localOffset + daylightSavings;
+
+    }
+    float calculateSunset(int year,int month,int day,float lat, float lng,int localOffset, int daylightSavings) {
+    /*
+    localOffset will be <0 for western hemisphere and >0 for eastern hemisphere
+    daylightSavings should be 1 if it is in effect during the summer otherwise it should be 0
+    */
+    //1. first calculate the day of the year
+    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;
+
+    //2. convert the longitude to hour value and calculate an approximate time
+    float lngHour = lng / 15.0;      
+    //float t = N + ((6 - lngHour) / 24);   //if rising time is desired:
+    float t = N + ((18 - lngHour) / 24);  //if setting time is desired:
+
+    //3. calculate the Sun's mean anomaly   
+    float M = (0.9856 * t) - 3.289;
+
+    //4. calculate the Sun's true longitude
+    float L = fmod(M + (1.916 * sin((PI/180)*M)) + (0.020 * sin(2 *(PI/180) * M)) + 282.634,360.0);
+
+    //5a. calculate the Sun's right ascension      
+    float RA = fmod(180/PI*atan(0.91764 * tan((PI/180)*L)),360.0);
+
+    //5b. right ascension value needs to be in the same quadrant as L   
+    float Lquadrant  = floor( L/90) * 90;
+    float RAquadrant = floor(RA/90) * 90;
+    RA = RA + (Lquadrant - RAquadrant);
+
+    //5c. right ascension value needs to be converted into hours   
+    RA = RA / 15;
+
+    //6. calculate the Sun's declination
+    float sinDec = 0.39782 * sin((PI/180)*L);
+    float cosDec = cos(asin(sinDec));
+
+    //7a. calculate the Sun's local hour angle
+    float cosH = (sin((PI/180)*ZENITH) - (sinDec * sin((PI/180)*lat))) / (cosDec * cos((PI/180)*lat));
+    /*   
+    if (cosH >  1) 
+    the sun never rises on this location (on the specified date)
+    if (cosH < -1)
+    the sun never sets on this location (on the specified date)
+    */
+
+    //7b. finish calculating H and convert into hours
+    //float H = 360 - (180/PI)*acos(cosH);   //   if if rising time is desired:
+    float H = (180/PI)*acos(cosH);//   if setting time is desired:      
+    H = H / 15;
+
+    //8. calculate local mean time of rising/setting      
+    float T = H + RA - (0.06571 * t) - 6.622;
+
+    //9. adjust back to UTC
+    float UT = fmod(T - lngHour,24.0);
+
+    //10. convert UT value to local time zone of latitude/longitude
+    return UT + localOffset + daylightSavings;
+
+    }
+/*
+    void printSunrise() {
+    float localT = calculateSunrise(args);
+    double hours;
+    float minutes = modf(localT,&hours)*60;
+    printf("%.0f:%.0f",hours,minutes);
+    }
+    */