1 |
jahn |
1.1 |
function radtrans_jd(i,j,k) |
2 |
|
|
c |
3 |
|
|
c |
4 |
|
|
c This function converts a calendar date to the corresponding Julian |
5 |
|
|
c day starting at noon on the calendar date. The algorithm used is |
6 |
|
|
c from Van Flandern and Pulkkinen, Ap. J. Supplement Series 41, |
7 |
|
|
c November 1979, p. 400. |
8 |
|
|
c |
9 |
|
|
c |
10 |
|
|
c Arguments |
11 |
|
|
c |
12 |
|
|
c Name Type I/O Description |
13 |
|
|
c ---- ---- --- ----------- |
14 |
|
|
c i I*4 I Year - e.g. 1970 |
15 |
|
|
c j I*4 I Month - (1-12) |
16 |
|
|
c k I*4 I Day - (1-31) |
17 |
|
|
c jd I*4 O Julian day |
18 |
|
|
c |
19 |
|
|
c external references |
20 |
|
|
c ------------------- |
21 |
|
|
c none |
22 |
|
|
c |
23 |
|
|
c |
24 |
|
|
c Written by Frederick S. Patt, GSC, November 4, 1992 |
25 |
|
|
c |
26 |
|
|
c |
27 |
|
|
INTEGER radtrans_jd,i,j,k |
28 |
|
|
|
29 |
|
|
radtrans_jd = 367*i - 7*(i+(j+9)/12)/4 + 275*j/9 + k + 1721014 |
30 |
|
|
|
31 |
|
|
c This additional calculation is needed only for dates outside of the |
32 |
|
|
c period March 1, 1900 to February 28, 2100 |
33 |
|
|
c radtrans_jd = radtrans_jd + 15 - 3*((i+(j-9)/7)/100+1)/4 |
34 |
|
|
|
35 |
|
|
return |
36 |
|
|
end |