суть в следующем: для расчета интенсивности потока прямой солнечной радиации необходимо знать высоту солнца h.
sin(h) = sin(j)*cos(d) + cos(j)*cos(d)*cos(t)
где d - склонение
t - часовой угол
j - географическая широта
с j и t вопросов нет.
а склонение? его по таблицам, что ли брать? а если мне необходимо проследить изменение его при разной величине часового угла в течение суток?
есть конечно выражения связывающие экваториальные и горизонтальные координаты:
sin(d) = sin(j)*cos(z) - cos(j)*sin(z)*cos(A)
sin(A)=cos(d)sin(t)/cos(h)
z+h=90 град
z - зенитное расстояние
A  - азимут
но я сомневаюсь, что склонение считается по ним. так как в них не наблюдается изменение склонения в зависимости от времени года. а, насколько я понимаю,  она - зависимость - должна быть.
в общем как-то неясен этот вопрос...
thank

http://www.astro.uio.no/~bgranslo/aares/calculate.html - вполне оригинальные формулы

Это конечно все здорово: программы и таблицы. Особенно благодарен за ссылки на англоязычный сайт (сложновато для меня, а ошибиться не хочется) . НО
смысл в том чтобы используя формулы получить значения склонения, которые в свою очередь используются для расчета интенсивности солнечной радиации. И это все в динамике на протяжении суток и всего лета.
Естественно это не конечная цель. полученные данные применяются  в дальнейших расчетах.

поэтому: как найти эти самые исходные формулы?

...смысл в том чтобы используя формулы получить значения склонения, которые в свою очередь используются для ...
поэтому: как найти эти самые исходные формулы?

Смотря какая точность требуется. Если пренебречь эллиптичностью земной орбиты (а она невелика), то формула зависимости склонения от порядкового дня в году будет элементарной: delta = eps*sin((D-D0)/365.25), где  eps - угол наклона эклиптики к экватору (~ 23њ26'), D - порядковый номер текущего (т.е. для которого ищется склонение Солнца) дня в году,  D0 - порядковый номер дня весеннего равноденствия (примерно 81 - это можно уточнить). С порядковыми днями удобно работать через юлианские, но если особой точности не требуется, то есть разные пути радикального упрощения таких вычислений.

А для большей точности уже придется учитывать неравномерность орбитального движения Земли - в первом приближении движущееся во эклиптике Солнце заставить двигаться с переменной скорость (на среднюю скорость наложить соответствующую синусоиду), а затем эклиптические координаты перевести в экваториальные. И тогда, уже после перевода в горизонтальные, не забыть еще и рефракцию учесть...