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 |