| 1 |
jahn |
1.1 |
C $Header$ |
| 2 |
|
|
C $Name$ |
| 3 |
|
|
|
| 4 |
|
|
#include "RADTRANS_OPTIONS.h" |
| 5 |
|
|
|
| 6 |
|
|
CBOP |
| 7 |
|
|
C !ROUTINE: RADTRANS_EPHPARMS |
| 8 |
|
|
|
| 9 |
|
|
C !INTERFACE: ========================================================== |
| 10 |
|
|
subroutine radtrans_ephparms (t, |
| 11 |
|
|
O xls, gs, xlm, omega) |
| 12 |
|
|
|
| 13 |
|
|
C !DESCRIPTION: |
| 14 |
|
|
c This subroutine computes ephemeris parameters used by other Mission |
| 15 |
|
|
c Operations routines: the solar mean longitude and mean anomaly, and |
| 16 |
|
|
c the lunar mean longitude and mean ascending node. It uses the model |
| 17 |
|
|
c referenced in The Astronomical Almanac for 1984, Section S |
| 18 |
|
|
c (Supplement) and documented and documented in Exact closed-form |
| 19 |
|
|
c geolocation algorithm for Earth survey sensors, by F.S. Patt and |
| 20 |
|
|
c W.W. Gregg, Int. Journal of Remote Sensing, 1993. These parameters |
| 21 |
|
|
c are used to compute the solar longitude and the nutation in |
| 22 |
|
|
c longitude and obliquity. |
| 23 |
|
|
c |
| 24 |
|
|
c Program written by: Frederick S. Patt |
| 25 |
|
|
c General Sciences Corporation |
| 26 |
|
|
c November 2, 1992 |
| 27 |
|
|
|
| 28 |
|
|
c Calling Arguments |
| 29 |
|
|
C !USES: =============================================================== |
| 30 |
|
|
IMPLICIT NONE |
| 31 |
|
|
|
| 32 |
|
|
C !INPUT PARAMETERS: =================================================== |
| 33 |
|
|
c t :: Time in days since January 1, 2000 at 12 hours UT |
| 34 |
|
|
_RL t |
| 35 |
|
|
c INTEGER myThid |
| 36 |
|
|
|
| 37 |
|
|
C !OUTPUT PARAMETERS: ================================================== |
| 38 |
|
|
c xls :: Mean solar longitude (degrees) |
| 39 |
|
|
c gs :: Mean solar anomaly (degrees) |
| 40 |
|
|
c xlm :: Mean lunar longitude (degrees) |
| 41 |
|
|
c omega :: Ascending node of mean lunar orbit (degrees) |
| 42 |
|
|
_RL xls, gs, xlm, omega |
| 43 |
|
|
CEOP |
| 44 |
|
|
|
| 45 |
|
|
C !LOCAL VARIABLES: ==================================================== |
| 46 |
|
|
c |
| 47 |
|
|
|
| 48 |
|
|
c Sun Mean Longitude |
| 49 |
|
|
xls = 280.46592D0 + 0.9856473516D0*t |
| 50 |
|
|
xls = mod(xls,360.0) |
| 51 |
|
|
|
| 52 |
|
|
c Sun Mean Anomaly |
| 53 |
|
|
gs = 357.52772D0 + 0.9856002831D0*t |
| 54 |
|
|
gs = mod(gs,360.0) |
| 55 |
|
|
|
| 56 |
|
|
c Moon Mean Longitude |
| 57 |
|
|
xlm = 218.31643D0 + 13.17639648D0*t |
| 58 |
|
|
xlm = mod(xlm,360.0) |
| 59 |
|
|
|
| 60 |
|
|
c Ascending Node of Moons Mean Orbit |
| 61 |
|
|
omega = 125.04452D0 - 0.0529537648D0*t |
| 62 |
|
|
omega = mod(omega,360.0) |
| 63 |
|
|
|
| 64 |
|
|
return |
| 65 |
|
|
end |
| 66 |
|
|
|