Закат на языке Си - PullRequest
Новая
       32

Закат на языке Си

0 голосов
/ 04 октября 2018

Я хочу рассчитать закат для конкретной локализации (широта, долгота) в C.

Я смоделировал на: http://edwilliams.org/sunrise_sunset_algorithm.htm

Для данной модели я попытался вычислить восход -это значение было хорошо.Я все еще получаю неправильное значение для заката: -11: -22.

Вот код: 

#include <stdio.h>
#include <math.h>
#define PI 3.1415926
#define ZENITH -.83

float calculateSunset(int year,int month,int day,float lat, float lng,int localOffset) {

//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;
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) {
    printf("the sun never rises on this location (on the specified date)");
    return 0;
}

if (cosH < -1) {
    printf("the sun never sets on this location (on the specified date)");
    return 0;
}

//7b. finish calculating H and convert into hours
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;

}

void printSunset() {

    float localT = calculateSunset(2018,10,4,51.8446,19.2094,2);
    double hours;
    float minutes = modf(localT,&hours)*60;
    printf("Sunset: ");
    printf("%.0f:%.0f",hours,minutes);   
    printf("\n"); 
}


int main()
{
  printSunset();
  return 0;
}

Кто-нибудь может мне помочь?Что я делаю не так?

1 Ответ

0 голосов
/ 04 октября 2018

Я думаю, что ваша основная проблема в том, что вы ожидаете, что fmod вернет положительный результат - это не так, fmod сохраняет знак, т.е. если вы принимаете fmod отрицательного значения, вы получите (правильный) отрицательный результат, еслиВы используете его на положительном значении (также правильный) результат будет положительным.Если вам требуется положительный результат, вы можете, например, просто добавить делитель, если результат отрицательный, т.е. 

float a = fmod(T, 360.0);
if(a < 0) { a = a + 360.0; }
...