| 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 |
Parent Directory
| 
Revision Log
| 
Revision Graph