Незначительная ошибка в расчете внеземного излучения (метод "asce")?

Question

Пока смотрю на код для расчетов внеземного излучения get_extra_radiation() и перепроверяя это с литературой, я заметил, что для обоих методов "asce" и "spencer" solarposition._calculate_simple_day_angle(doy) функция используется. Эта функция в основном просто вычисляет (2. * np.pi / 365.) * (doy- 1), что правильно для method='spencer'.

Но я думаю, что для method='asce' это должно быть просто (2. * np.pi * doy / 365.) - без "- 1" - как описано на странице 9 в "JA Duffie and WA Beckman, " Солнечная инженерия тепловых процессов, 3-е издание " J. Wiley and Sons, New York (2006)"

Answer 1

У меня просто случайно появилась эта книга, поэтому я проверил ее на странице 9, и, как уже отмечалось, есть разница pvlib-python код. Одна формула использует n, другой использует n-1, Эта разница, вероятно, несущественна, но, возможно, ее следует исправить, чтобы она была правильной. Вопрос о GitHub будет лучшим способом поднять это.