solarPos  태양의 황경, 지심 거리를 계산한다

long solarPos(double *lng, long j, double d);

lng  (출력) 황경(단위, 번)
j    (입력) 세계 표준시부터 환산한 유리우스일의 정수 부분
d    (입력) 세계 표준시부터 환산한 유리우스일의 소수 부분

지심 거리(지구·태양 중심간의 거리)(단위, 킬로·미터)

long solarPos(double *lng, long j, double d)
{
    typedef struct { double k, a, b; } TBL;
    static TBL l[16] = {
    {  1.9147,  35999.05, 267.52 }, { -0. 0048,  35999.05, 267.52 },
    {  0.0200,  71998.1,  265.1  }, {  0.0020,  32964.0,  158.0  },
    {  0.0018,     19.0,  159.0  }, {  0.0018, 445267.0,  208.0  },
    {  0.0015,  45038.0,  254.0  }, {  0.0013,  22519.0,  352.0  },
    {  0.0007,  65929.0,   45.0  }, {  0.0007,   3035.0,  110.0  },
    {  0.0007,   9038.0,   64.0  }, {  0.0006,  33718.0,  316.0  },
    {  0.0005,    155.0,  118.0  }, {  0.0005,   2281.0,  221.0  },
    {  0.0004,  29930.0,   48.0  }, {  0.0004,  31557.0,  161.0  }};
    static TBL c[8] = {
    {  0.016706,  35999.05, 177.53 }, { -0. 000042,  35999.05, 177.53 },
    {  0.000139,  71998.0,  175.0  }, {  0.000031, 445267.0,  298.0  },
    {  0.000016,  32964.0,   68.0  }, {  0.000016,  45038.0,  164.0  },
    {  0.000005,  22519.0,  233.0  }, {  0.000005,  33718.0,  226.0  }};
    int i;
    double s, t, r;

    t = (j + d - 2451545.0)/36525. 0;
    s = 280.4659 + 36000.7695 * t;
    for (i = 0; i < 16; i++) {
        r = l[i]. a * t + l[i]. b;
        s += l[i]. k * cos(r * 0.017453292519943);
    }
    while (s < 0) s += 360.0;
    while (s >= 360.0) s-= 360.0;
    *lng = s;

    s = 1.000140;
    for (i = 0; i < 8; i++) {
        r = c[i]. a * t + c[i]. b;
        s += c[i]. k * cos(r * 0.017453292519943);
    }
    return (long)(s * 149597870.0);
}

지금 사용되고 있는 태양력으로는, 매년 2월 4 히타치봄은 황경 315도, 3월 21일 춘분은 0 번, 6월 22일 하지는 90번, 8월 8 히타치가을은 180번, 9월 24일 추분은 180번, 12월 22일 동지는 270번.

본프로그램을 사용하면(자), 어느 각의 태양의 황경 및 지구까지의 거리가 안다.

덧붙여 일본 시간 JST를 세계 표준시 UT로 고치려면 , UT = JST - 9 로 계산 좋다.